High-throughput sequencing of black pepper root transcriptome

Assessing the biodiversity of rhizosphere and endophytic fungi in Knoxia valerianoides under continuous cropping conditions

BACKGROUND

Rhizosphere and endophytic fungi play important roles in plant health and crop productivity. However, their community dynamics during the continuous cropping of Knoxia valerianoides have rarely been reported. K. valerianoides is a perennial herb of the family Rubiaceae and has been used in herbal medicines for ages. Here, we used high-throughput sequencing technology Illumina MiSeq to study the structural and functional dynamics of the rhizosphere and endophytic fungi of K. valerianoides.

RESULTS

The findings indicate that continuous planting has led to an increase in the richness and diversity of rhizosphere fungi, while concomitantly resulting in a decrease in the richness and diversity of root fungi. The diversity of endophytic fungal communities in roots was lower than that of the rhizosphere fungi. Ascomycota and Basidiomycota were the dominant phyla detected during the continuous cropping of K. valerianoides. In addition, we found that root rot directly affected the structure and diversity of fungal communities in the rhizosphere and the roots of K. valerianoides. Consequently, both the rhizosphere and endophyte fungal communities of root rot-infected plants showed higher richness than the healthy plants. The relative abundance of Fusarium in two and three years old root rot-infected plants was significantly higher than the control, indicating that continuous planting negatively affected the health of K. valerianoides plants. Decision Curve Analysis showed that soil pH, organic matter (OM), available K, total K, soil sucrase (S_SC), soil catalase (S_CAT), and soil cellulase (S_CL) were significantly related (p < 0.05) to the fungal community dynamics.

CONCLUSIONS

The diversity of fungal species in the rhizosphere and root of K. valerianoides was reported for the first time. The fungal diversity of rhizosphere soil was higher than that of root endophytic fungi. The fungal diversity of root rot plants was higher than that of healthy plants. Soil pH, OM, available K, total K, S_CAT, S_SC, and S_CL were significantly related to the fungal diversity. The occurrence of root rot had an effect on the community structure and diversity of rhizosphere and root endophytic fungi.

Applications of some advanced sequencing, analytical, and computational approaches in medicinal plant research: a review

The potential active chemicals found in medicinal plants, which have long been employed as natural medicines, are abundant. Exploring the genes responsible for producing these compounds has given new insights into medicinal plant research. Previously, the authentication of medicinal plants was done via DNA marker sequencing. With the advancement of sequencing technology, several new techniques like next-generation sequencing, single molecule sequencing, and fourth-generation sequencing have emerged. These techniques enshrined the role of molecular approaches for medicinal plants because all the genes involved in the biosynthesis of medicinal compound(s) could be identified through RNA-seq analysis. In several research insights, transcriptome data have also been used for the identification of biosynthesis pathways. miRNAs in several medicinal plants and their role in the biosynthesis pathway as well as regulation of the disease-causing genes were also identified. In several research articles, an in silico study was also found to be effective in identifying the inhibitory effect of medicinal plant-based compounds against virus' gene(s). The use of advanced analytical methods like spectroscopy and chromatography in metabolite proofing of secondary metabolites has also been reported in several recent research findings. Furthermore, advancement in molecular and analytic methods will give new insight into studying the traditionally important medicinal plants that are still unexplored.

Identification and screening of fungicides against Piper nigrum basal Fusarium wilt disease in Hainan, China

Fusarium wilt has occurred in the main Piper nigrum cultivation regions, which seriously affects the yield and quality of P. nigrum. To identify the pathogen of this disease, the diseased roots were collected from a demonstration base in Hainan Province. The pathogen was obtained by tissue isolation method and confirmed by pathogenicity test. Based on the morphological observation, sequence analyses of TEF1-α nuclear gene, Fusarium solani was identified as the pathogen causing P. nigrum Fusarium wilt and induced symptoms on inoculated plants, including chlorosis, necrotic spots, wilt, drying, and root rot. The experiments for the antifungal activity showed that all the 11 fungicides selected in this study showed certain inhibitory effects on the colony growth of F. solani, where 2% kasugamycin AS, 45% prochloraz EW, 25 g·L-1 fludioxonil SC and 430 g·L-1 tebuconazole SC exhibited relative higher inhibitory effects with EC50 as 0.065, 0.205, 0.395, and 0.483 mg·L-1 , respectively, and were selected to perform SEM analysis and test in seeds in vitro. The SEM analysis showed that kasugamycin, prochloraz, fludioxonil, and tebuconazole might have exerted their antifungal effect by damaging F. solani mycelia or microconidia. These preparations were applied as a seed coating of P. nigrum Reyin-1. The kasugamycin treatment was most effective in reducing the harmful impact of F. solani on the seed germination. These results presented herein provide useful guidance for the effective control of P. nigrum Fusarium wilt.

A comprehensive review on genomic resources in medicinally and industrially important major spices for future breeding programs: Status, utility and challenges

In the global market, spices possess a high-value but low-volume commodities of commerce. The food industry depends largely on spices for taste, flavor, and therapeutic properties in replacement of cheap synthetic ones. The estimated growth rate for spices demand in the world is ∼3.19%. Since spices grow in limited geographical regions, India is one of the leading producer of spices, contributing 25-30 percent of total world trade. Hitherto, there has been no comprehensive review of the genomic resources of industrially important major medicinal spices to overcome major impediments in varietal improvement and management. This review focuses on currently available genomic resources of 24 commercially significant spices, namely, Ajwain, Allspice, Asafoetida, Black pepper, Cardamom large, Cardamom small, Celery, Chillies, Cinnamon, Clove, Coriander, Cumin, Curry leaf, Dill seed, Fennel, Fenugreek, Garlic, Ginger, Mint, Nutmeg, Saffron, Tamarind, Turmeric and Vanilla. The advent of low-cost sequencing machines has contributed immensely to the voluminous data generation of these spices, cracking the complex genomic architecture, marker discovery, and understanding comparative and functional genomics. This review of spice genomics resources concludes the perspective and way forward to provide footprints by uncovering genome assemblies, sequencing and re-sequencing projects, transcriptome-based studies, non-coding RNA-mediated regulation, organelles-based resources, developed molecular markers, web resources, databases and AI-directed resources in candidate spices for enhanced breeding potential in them. Further, their integration with molecular breeding could be of immense use in formulating a strategy to protect and expand the production of the spices due to increased global demand.

Elucidating the Mesocarp Drupe Transcriptome of Açai (Euterpe oleracea Mart.): An Amazonian Tree Palm Producer of Bioactive Compounds

Euterpe oleracea palm, endemic to the Amazon region, is well known for açai, a fruit violet beverage with nutritional and medicinal properties. During E. oleracea fruit ripening, anthocyanin accumulation is not related to sugar production, contrarily to grape and blueberry. Ripened fruits have a high content of anthocyanins, isoprenoids, fibers, and proteins, and are poor in sugars. E. oleracea is proposed as a new genetic model for metabolism partitioning in the fruit. Approximately 255 million single-end-oriented reads were generated on an Ion Proton NGS platform combining fruit cDNA libraries at four ripening stages. The de novo transcriptome assembly was tested using six assemblers and 46 different combinations of parameters, a pre-processing and a post-processing step. The multiple k-mer approach with TransABySS as an assembler and Evidential Gene as a post-processer have shown the best results, with an N50 of 959 bp, a read coverage mean of 70x, a BUSCO complete sequence recovery of 36% and an RBMT of 61%. The fruit transcriptome dataset included 22,486 transcripts representing 18 Mbp, of which a proportion of 87% had significant homology with other plant sequences. Approximately 904 new EST-SSRs were described, and were common and transferable to Phoenix dactylifera and Elaeis guineensis, two other palm trees. The global GO classification of transcripts showed similar categories to that in P. dactylifera and E. guineensis fruit transcriptomes. For an accurate annotation and functional description of metabolism genes, a bioinformatic pipeline was developed to precisely identify orthologs, such as one-to-one orthologs between species, and to infer multigenic family evolution. The phylogenetic inference confirmed an occurrence of duplication events in the Arecaceae lineage and the presence of orphan genes in E. oleracea. Anthocyanin and tocopherol pathways were annotated entirely. Interestingly, the anthocyanin pathway showed a high number of paralogs, similar to in grape, whereas the tocopherol pathway exhibited a low and conserved gene number and the prediction of several splicing forms. The release of this exhaustively annotated molecular dataset of E. oleracea constitutes a valuable tool for further studies in metabolism partitioning and opens new great perspectives to study fruit physiology with açai as a model.

Toxic Elemental Impurities in Herbal Weight Loss Supplements; A Study Using ICP-OES Microwave-Assisted Digestion

The tendency of using weight loss herbal preparations is continuously increasing, especially for the widespread consumption of junk food that is characterized by high calories. Weight loss herbal preparations are considered a type of food supplement product, and, as such, the regulations governing their quality control might be minimal. These products could be locally formulated in any country or internationally imported. Being non-controlled products, the herbal weight-loss products may contain high levels of elemental impurities that might exceed the permissible ranges. Moreover, these products contribute to the total daily intake (TDI) of such elements, which might represent concerns about their potential toxicological danger. In this research, the elemental contents in such products were investigated. The inductively coupled plasma with optical emission spectrometer (ICP-OES) was used to determine the levels of 15 elemental contents, namely, Na, K, Ca, Mg, Al, Cu, Fe, Li, Mn, As, Co, Cr, Cd, Ni and Pb. The results showed that seven micro-elements, namely Cd, Co, Ni, Cr, Pb, Li and Cu, were either not detectable or at a concentration much lower than their tolerable limits. However, all studied macro-elements (Na, K, Ca and Mg), together with Fe, were found at considerable, yet safe levels. On the other hand, Mn, Al and As contents showed perturbing levels in some of the studied products. Finally, a conclusion was highlighted for the necessity for stricter surveillance of such herbal products.

Comparative Transcriptome and Metabolome Analysis of Resistant and Susceptible Piper Species Upon Infection by the Oomycete Phytophthora Capsici

Phytophthora capsici is a destructive oomycete pathogen that causes devastating disease in black pepper, resulting in a significant decline in yield and economic losses. Piper nigrum (black pepper) is documented as susceptible to P. capsici, whereas its close relative Piper flaviflorum is known to be resistant. However, the molecular mechanism underlying the resistance of P. flaviflorum remains obscure. In this study, we conducted a comparative transcriptome and metabolome analysis between P. flaviflorum and P. nigrum upon P. capsici infection and found substantial differences in their gene expression profiles, with altered genes being significantly enriched in terms relating to plant-pathogen interaction, phytohormone signal transduction, and secondary metabolic pathways, including phenylpropanoid biosynthesis. Further metabolome analysis revealed the resistant P. flaviflorum to have a high background endogenous ABA reservoir and time-course-dependent accumulation of ABA and SA upon P. capsici inoculation, while the susceptible P. nigrum had a high background endogenous IAA reservoir and time-course-dependent accumulation of JA-Ile, the active form of JA. Investigation of the phenylpropanoid biosynthesis metabolome further indicated the resistant P. flaviflorum to have more accumulation of lignin precursors than the susceptible P. nigrum, resulting in a higher accumulation after inoculation. This study provides an overall characterization of biologically important pathways underlying the resistance of P. flaviflorum, which theoretically explains the advantage of using this species as rootstock for the management of oomycete pathogen in black pepper production.

Rapid Genome-Wide Location-Specific Polymorphic SSR Marker Discovery in Black Pepper by GBS Approach

Black pepper (Piper nigrum), the "King of Spices," is an economically important spice in India and is known for its medicinal and cultural values. SSRs, the tandem repeats of small DNA sequences, are often polymorphic in nature with diverse applications. For population structure, QTL/gene discovery, MAS, and diversity analysis, it is imperative to have their location specificity. The existing PinigSSRdb catalogs ~70K putative SSR markers but these are anonymous (unknown chromosomal location), based on 916 scaffolds rather than 26 chromosomes. Under this study, we generated ddRAD sequence data of 29 black pepper genotypes from all over India, being low-cost and most efficient technique for the identification of polymorphic markers. The major limitation of ddRAD with compromised/non-uniform coverage has been successfully overcome by taking advantage of chromosome-wise data availability. The latest black pepper genome assembly was used to extract genome-wide SSRs. A total of 276,230 genomic SSRs were mined distributed over 26 chromosomes, with relative density of 362.88 SSRs/Mb and average distance of 2.76 Kb between two SSRs. This assembly was also used to find the polymorphic SSRs in the generated GBS data of 29 black pepper genotypes utilizing rapid and cost-effective method giving 3,176 polymorphic SSRs, out of which 2015 were found to be hypervariable. The developed web-genomic resource, BlackP2MSATdb (http://webtom.cabgrid.res.in/blackp2msatdb/), is the largest and first reported web resource for genomic and polymorphic SSRs of black pepper, which is useful to develop varietal signature, coreset, physical map, QTL/gene identification, and MAS in endeavor of black pepper production.

Transcriptome and proteome of the corm, leaf and flower of Hypoxis hemerocallidea (African potato)

The corm of Hypoxis hemerocallidea, commonly known as the African potato, is used in traditional medicine to treat several medical conditions such as urinary infections, benign prostate hyperplasia, inflammatory conditions and testicular tumours. The metabolites contributing to the medicinal properties of H. hemerocallidea have been identified in several studies and, more recently, the active terpenoids of the plant were profiled. However, the biosynthetic pathways and the enzymes involved in the production of the terpene metabolites in H. hemerocallidea have not been characterised at a transcriptomic or proteomic level. In this study, total RNA extracted from the corm, leaf and flower tissues of H. hemerocallidea was sequenced on the Illumina HiSeq 2500 platform. A total of 143,549 transcripts were assembled de novo using Trinity and 107,131 transcripts were functionally annotated using the nr, GO, COG, KEGG and SWISS-PROT databases. Additionally, the proteome of the three tissues were sequenced using LC-MS/MS, revealing aspects of secondary metabolism and serving as data validation for the transcriptome. Functional annotation led to the identification of numerous terpene synthases such as nerolidol synthase, germacrene D synthase, and cycloartenol synthase amongst others. Annotations also revealed a transcript encoding the terpene synthase phytoalexin momilactone A synthase. Differential expression analysis using edgeR identified 946 transcripts differentially expressed between the three tissues and revealed that the leaf upregulates linalool synthase compared to the corm and the flower tissues. The transcriptome as well as the proteome of Hypoxis hemerocallidea presented here provide a foundation for future research.

Transcriptome Analysis of Chloris virgata, Which Shows the Fastest Germination and Growth in the Major Mongolian Grassland Plant

Plants in Mongolian grasslands are exposed to short, dry summers and long, cold winters. These plants should be prepared for fast germination and growth activity in response to the limited summer rainfall. The wild plant species adapted to the Mongolian grassland environment may allow us to explore useful genes, as a source of unique genetic codes for crop improvement. Here, we identified the Chloris virgata Dornogovi accession as the fastest germinating plant in major Mongolian grassland plants. It germinated just 5 h after treatment for germination initiation and showed rapid growth, especially in its early and young development stages. This indicates its high growth potential compared to grass crops such as rice and wheat. By assessing growth recovery after animal bite treatment (mimicked by cutting the leaves with scissors), we found that C. virgata could rapidly regenerate leaves after being damaged, suggesting high regeneration potential against grazing. To analyze the regulatory mechanism involved in the high growth potential of C. virgata, we performed RNA-seq-based transcriptome analysis and illustrated a comprehensive gene expression map of the species. Through de novo transcriptome assembly with the RNA-seq reads from whole organ samples of C. virgata at the germination stage (2 days after germination, DAG), early young development stage (8 DAG), young development stage (17 DAG), and adult development stage (28 DAG), we identified 21,589 unified transcripts (contigs) and found that 19,346 and 18,156 protein-coding transcripts were homologous to those in rice and Arabidopsis, respectively. The best-aligned sequences were annotated with gene ontology groups. When comparing the transcriptomes across developmental stages, we found an over-representation of genes involved in growth regulation in the early development stage in C. virgata. Plant development is tightly regulated by phytohormones such as brassinosteroids, gibberellic acid, abscisic acid, and strigolactones. Moreover, our transcriptome map demonstrated the expression profiles of orthologs involved in the biosynthesis of these phytohormones and their signaling networks. We discuss the possibility that C. virgata phytohormone signaling and biosynthesis genes regulate early germination and growth advantages. Comprehensive transcriptome information will provide a useful resource for gene discovery and facilitate a deeper understanding of the diversity of the regulatory systems that have evolved in C. virgata while adapting to severe environmental conditions.

Identification and characterization of piperine synthase from black pepper, Piper nigrum L

Black pepper (Piper nigrum L.) is the world's most popular spice and is also used as an ingredient in traditional medicine. Its pungent perception is due to the interaction of its major compound, piperine (1-piperoyl-piperidine) with the human TRPV-1 or vanilloid receptor. We now identify the hitherto concealed enzymatic formation of piperine from piperoyl coenzyme A and piperidine based on a differential RNA-Seq approach from developing black pepper fruits. This enzyme is described as piperine synthase (piperoyl-CoA:piperidine piperoyl transferase) and is a member of the BAHD-type of acyltransferases encoded by a gene that is preferentially expressed in immature fruits. A second BAHD-type enzyme, also highly expressed in immature black pepper fruits, has a rather promiscuous substrate specificity, combining diverse CoA-esters with aliphatic and aromatic amines with similar efficiencies, and was termed piperamide synthase. Recombinant piperine and piperamide synthases are members of a small gene family in black pepper. They can be used to facilitate the microbial production of a broad range of medicinally relevant aliphatic and aromatic piperamides based on a wide array of CoA-donors and amine-derived acceptors, offering widespread applications.

A piperic acid CoA ligase produces a putative precursor of piperine, the pungent principle from black pepper fruits

Black pepper (Piper nigrum L.) is known for its high content of piperine, a cinnamoyl amide derivative regarded as largely responsible for the pungent taste of this widely used spice. Despite its long history and worldwide use, the biosynthesis of piperine and related amides has been enigmatic up to now. In this report we describe a specific piperic acid CoA ligase from immature green fruits of P. nigrum. The corresponding enzyme was cloned and functionally expressed in E. coli. The recombinant enzyme displays a high specificity for piperic acid and does not accept the structurally related feruperic acid characterized by a similar C-2 extension of the general C6-C3 phenylpropanoid structure. The enzyme is also inactive with the standard set of hydroxycinnamic acids tested including caffeic acid, 4-coumaric acid, ferulic acid, and sinapic acid. Substrate specificity is corroborated by in silico modelling that suggests a perfect fit for the substrate piperic acid to the active site of the piperic acid CoA ligase. The CoA ligase gene shows its highest expression levels in immature green fruits, is also expressed in leaves and flowers, but not in roots. Virus-induced gene silencing provided some preliminary indications that the production of piperoyl-CoA is required for the biosynthesis of piperine in black pepper fruits.

Transcriptomic analysis of pepper plants provides insights into host responses to Fusarium solani infestation

Black pepper is an important commodity crop in Malaysia that generates millions of annual revenue for the country. However, black pepper yield is affected by slow decline disease caused by a soil-borne fungus Fusarium solani. RNA sequencing transcriptomics approach has been employed in this study to explore the differential gene expression in susceptible Piper nigrum L. and resistant Piper colubrinum Link. Gene expression comparative analysis of the two pepper species has yielded 2,361 differentially expressed genes (DEGs). Among them, higher expression of 1,426 DEGs was detected in resistant plant. These DEGs practically demonstrated the major branches of plant-pathogen interaction pathway (Path: ko04626). We selected five groups of defence-related DEGs for downstream qRT-PCR analysis. Cf-9, the gene responsible for recognizing fungal avirulence protein activity was found inexpressible in susceptible plant. However, this gene exhibited promising expression in resistant plant. Inactivation of Cf-9 could be the factor that causes susceptible plant fail in recognition of F. solani and subsequently delay activation of adaptive response to fungal invasion. This vital study advance the understanding of pepper plant defence in response to F. solani and aid in identifying potential solution to manage slow decline disease in black pepper cultivation.

J, A. S, Bhat KV, Gaikwad AB. Genome wide identification and characterization of microsatellite markers in black pepper (Piper nigrum): A valuable resource for boosting genomics applications

Black pepper is one of the most valued and widely used spices in the world and dominates multi-billion dollar global spices trade. India is amongst the major producers, consumers and exporters of black pepper. In spite of its commercial and cultural importance, black pepper has received meagre attention in terms of generation of genomic resources. Availability of markers distributed throughout the genome would facilitate and accelerate genetic studies, QTL identification, genetic enhancement and crop improvement in black pepper. In this perspective, the sequence information from the recently sequenced black pepper (Piper nigrum) genome has been used for identification and characterisation of Simple Sequence Repeats (SSRs). Total 69,126 SSRs were identified from assembled genomic sequence of P. nigrum. The SSR frequency was 158 per MB making it, one SSR for every 6.3 kb in the assembled genome. Among the different types of microsatellite repeat motifs, dinucleotides were the most abundant (48.6%), followed by trinucleotide (23.7%) and compound repeats (20.62%). A set of 85 SSRs were used for validation, of which 74 produced amplification products of expected size. Genetic diversity of 30 black pepper accessions using 50 SSRs revealed four distinct clusters. Further, the cross species transferability of the SSRs was checked in nine other Piper species. Out of 50 SSRs used, 19 and 31 SSRs were amplified in nine and seven species, respectively. Thus the identified SSRs may have application in other species of the genus Piper where genome sequence is not available yet. Present study reports the first NGS based genomic SSRs in black pepper and thus constitute a valuable resource for a whole fleet of applications in genetics and plant breeding studies such as genetic map construction, QTL identification, map-based gene cloning, marker-assisted selection and evolutionary studies in Piper nigrum and related species.

Identification of stable reference gene for transcript normalization in black pepper-Phytophthora capsici pathosystem

A systematic validation of reference genes is a pre-requisite for the proper normalization of gene transcripts. In the present study, the annotated sequences from black pepper (Piper nigrum L.) leaf transcriptome were used as reference genes namely actin (PnACT), glyceraldehyde phosphate dehydrogenase (PnGAPDH), β-tubulin (PnTUB), ubiquitin conjugating enzyme (PnUBCE), 18srRNA and elongation factor-1-α (PnElF) to identify the stable reference gene. We focused the selection of stable reference gene on important biotic stress (Phytophthora) with different algorithms (geNorm, NormFinder and BestKeeper) along with Reffinder which resulted in identification of PnGAPDH and PnUBCE as stable genes. Norm qPCR (R package) was also used to estimate the stability of the selected genes. We elucidated the expression patterns of a target gene PnBGLU which codes for 1,3 beta glucanase with most stable as well as least stable reference genes by which the importance of selecting the stable gene for gene expression studies in this system was emphasized. The mean expression levels of PnBGLU was significantly overestimated and misinterpreted when least stable reference gene was used as normalizer. The selected reference genes on further analysis of the expression dynamics of PnBGLU among resistant and susceptible genotypes showed PnGAPDH as the suitable reference gene for P. nigrum-P. capsici pathosystem.

Transcriptome analyses of Paris polyphylla var. chinensis, Ypsilandra thibetica, and Polygonatum kingianum characterize their steroidal saponin biosynthesis pathway

Steroidal saponins, one of the most diverse groups of plant-derived natural products, elicit biological and pharmacological activities; however, the genes involved in their biosynthesis and the corresponding biosynthetic pathway in monocotyledon plants remain unclear. This study aimed to identify genes involved in the biosynthesis of steroidal saponins by performing a comparative analysis among transcriptomes of Paris polyphylla var. chinensis (PPC), Ypsilandra thibetica (YT), and Polygonatum kingianum (PK). De novo transcriptome assemblies generated 57,537, 140,420, and 151,773 unigenes from PPC, YT, and PK, respectively, of which 56.54, 47.81, and 44.30% were successfully annotated, respectively. Among the transcriptomes for PPC, YT, and PK, we identified 194, 169, and 131; 17, 14, and 26; and, 80, 122, and 113 unigenes corresponding to terpenoid backbone biosynthesis; sesquiterpenoid and triterpenoid biosynthesis; and, steroid biosynthesis pathways, respectively. These genes are putatively involved in the biosynthesis of cholesterol that is the primary precursor of steroidal saponins. Phylogenetic analyses indicated that lanosterol synthase may be exclusive to dicotyledon plant species, and the cytochrome P450 unigenes were closely related to clusters CYP90B1 and CYP734A1, which are UDP-glycosyltransferases unigenes homologous with the UGT73 family. Thus, unigenes of β-glucosidase may be candidate genes for catalysis of later period modifications of the steroidal saponin skeleton. Our data provide evidence to support the hypothesis that monocotyledons biosynthesize steroidal saponins from cholesterol via the cycloartenol pathway.

Secondary Metabolic Profiles of Two Cultivars of Piper nigrum (Black Pepper) Resulting from Infection by Fusarium solani f. sp. piperis

Bragantina and Cingapura are the main black pepper (Piper nigrum L.) cultivars and the Pará state is the largest producer in Brazil with about 90% of national production, representing the third largest production in the world. The infection of Fusarium solani f. sp. piperis, the causal agent of Fusarium disease in black pepper, was monitored on the cultivars Bragantina (susceptible) and Cingapura (tolerant), during 45 days’ post infection (dpi). Gas Chromatography-Mass spectrometry (GC-MS) analysis of the volatile concentrates of both cultivars showed that the Bragantina responded with the production of higher contents of α-bisabolol at 21 dpi and a decrease of elemol, mostly at 30 dpi; while Cingapura displayed an decrease of δ-elemene production, except at 15 dpi. The phenolic content determined by the Folin Ciocalteu method showed an increase in the leaves of plants inoculated at 7 dpi (Bragantina) and 7–15 dpi (Cingapura); in the roots, the infection caused a phenolic content decrease in Bragantina cultivar at 45 dpi and an increase in the Cingapura cultivar at 15, 30 and 45 dpi. High Performance Liquid Chromatography-Mass spectrometry (HPLC-MS) analysis of the root extracts showed a qualitative variation of alkamides during infection. The results indicated that there is a possible relationship between secondary metabolites and tolerance against phytopathogens.

Localization of Low-Copy DNA Sequences on Mitotic Chromosomes by FISH

Fluorescence in situ hybridization (FISH) is a widely used method to localize DNA sequences on mitotic and meiotic chromosomes and interphase nuclei. It was developed in early 1980s and since then it has contributed to numerous studies and important discoveries. Over the decades, the protocol was modified for ease of use, allowing for localizing multiple probes simultaneously and increasing its sensitivity and specificity. Despite the continuous improvements, the ability to detect short single-copy sequences of only a few kilobases or less, such as genes, remains limited. Here, we provide a detailed protocol for detection of short, single- or low-copy sequences on plant mitotic metaphase chromosomes.

Proteomic profile of Piper tuberculatum (Piperaceae)

Abstract Piper tuberculatum (Piperaceae) is a species that accumulates especially amides as secondary metabolites and several biological activities was previously reported. In this article, we report a proteomic study of P. tuberculatum. Bidimensional electrophoresis (2D SDS-PAGE) and mass spectrometry (ESI-Q-TOF) were used in this study. Over a hundred spots and various peptides were identified in this species and the putative functions of these peptides related to defense mechanism as biotic and abiotic stress were assigned. The information presented extend the range of molecular information of P. tuberculatum.

Modeling the impact of climate change on wild Piper nigrum (Black Pepper) in Western Ghats, India using ecological niche models

The center of diversity of Piper nigrum L. (Black Pepper), one of the highly valued spice crops is reported to be from India. Black pepper is naturally distributed in India in the Western Ghats biodiversity hotspot and is the only known existing source of its wild germplasm in the world. We used ecological niche models to predict the potential distribution of wild P. nigrum in the present and two future climate change scenarios viz (A1B) and (A2A) for the year 2080. Three topographic and nine uncorrelated bioclim variables were used to develop the niche models. The environmental variables influencing the distribution of wild P. nigrum across different climate change scenarios were identified. We also assessed the direction and magnitude of the niche centroid shift and the change in niche breadth to estimate the impact of projected climate change on the distribution of P. nigrum. The study shows a niche centroid shift in the future climate scenarios. Both the projected future climate scenarios predicted a reduction in the habitat of P. nigrum in Southern Western Ghats, which harbors many wild accessions of P. nigrum. Our results highlight the impact of future climate change on P. nigrum and provide useful information for designing sound germplasm conservation strategies for P. nigrum.

De novo transcriptome sequencing of black pepper (Piper nigrum L.) and an analysis of genes involved in phenylpropanoid metabolism in response to Phytophthora capsici

BACKGROUND

Piper nigrum L., or "black pepper", is an economically important spice crop in tropical regions. Black pepper production is markedly affected by foot rot disease caused by Phytophthora capsici, and genetic improvement of black pepper is essential for combating foot rot diseases. However, little is known about the mechanism of anti- P. capsici in black pepper. The molecular mechanisms underlying foot rot susceptibility were studied by comparing transcriptome analysis between resistant (Piper flaviflorum) and susceptible (Piper nigrum cv. Reyin-1) black pepper species.

RESULTS

116,432 unigenes were acquired from six libraries (three replicates of resistant and susceptible black pepper samples), which were integrated by applying BLAST similarity searches and noted by adopting Kyoto Encyclopaedia of Genes and Gene Ontology (GO) genome orthology identifiers. The reference transcriptome was mapped using two sets of digital gene expression data. Using GO enrichment analysis for the differentially expressed genes, the majority of the genes associated with the phenylpropanoid biosynthesis pathway were identified in P. flaviflorum. In addition, the expression of genes revealed that after susceptible and resistant species were inoculated with P. capsici, the majority of genes incorporated in the phenylpropanoid metabolism pathway were up-regulated in both species. Among various treatments and organs, all the genes were up-regulated to a relatively high degree in resistant species. Phenylalanine ammonia lyase and peroxidase enzyme activity increased in susceptible and resistant species after inoculation with P. capsici, and the resistant species increased faster. The resistant plants retain their vascular structure in lignin revealed by histochemical analysis.

CONCLUSIONS

Our data provide critical information regarding target genes and a technological basis for future studies of black pepper genetic improvements, including transgenic breeding.

Transcriptome- Assisted Label-Free Quantitative Proteomics Analysis Reveals Novel Insights into Piper nigrum-Phytophthora capsici Phytopathosystem

Black pepper (Piper nigrum L.), a tropical spice crop of global acclaim, is susceptible to Phytophthora capsici, an oomycete pathogen which causes the highly destructive foot rot disease. A systematic understanding of this phytopathosystem has not been possible owing to lack of genome or proteome information. In this study, we explain an integrated transcriptome-assisted label-free quantitative proteomics pipeline to study the basal immune components of black pepper when challenged with P. capsici. We report a global identification of 532 novel leaf proteins from black pepper, of which 518 proteins were functionally annotated using BLAST2GO tool. A label-free quantitation of the protein datasets revealed 194 proteins common to diseased and control protein datasets of which 22 proteins showed significant up-regulation and 134 showed significant down-regulation. Ninety-three proteins were identified exclusively on P. capsici infected leaf tissues and 245 were expressed only in mock (control) infected samples. In-depth analysis of our data gives novel insights into the regulatory pathways of black pepper which are compromised during the infection. Differential down-regulation was observed in a number of critical pathways like carbon fixation in photosynthetic organism, cyano-amino acid metabolism, fructose, and mannose metabolism, glutathione metabolism, and phenylpropanoid biosynthesis. The proteomics results were validated with real-time qRT-PCR analysis. We were also able to identify the complete coding sequences for all the proteins of which few selected genes were cloned and sequence characterized for further confirmation. Our study is the first report of a quantitative proteomics dataset in black pepper which provides convincing evidence on the effectiveness of a transcriptome-based label-free proteomics approach for elucidating the host response to biotic stress in a non-model spice crop like P. nigrum, for which genome information is unavailable. Our dataset will serve as a useful resource for future studies in this plant. Data are available via ProteomeXchange with identifier PXD003887.

Comparison of De Novo Transcriptome Assemblers and k-mer Strategies Using the Killifish, Fundulus heteroclitus

Background

De novo assembly of non-model organism’s transcriptomes has recently been on the rise in concert with the number of de novo transcriptome assembly software programs. There is a knowledge gap as to what assembler software or k-mer strategy is best for construction of an optimal de novo assembly. Additionally, there is a lack of consensus on which evaluation metrics should be used to assess the quality of de novo transcriptome assemblies.

Result

Six different assembly strategies were evaluated from four different assemblers. The Trinity assembly was used in its default 25 single k-mer value while Bridger, Oases, and SOAPdenovo-Trans were performed with multiple k-mer strategies. Bridger, Oases, and SOAPdenovo-Trans used a small multiple k-mer (SMK) strategy consisting of the k-mer lengths of 21, 25, 27, 29, 31, and 33. Additionally, Oases and SOAPdenovo-Trans were performed using a large multiple k-mer (LMK) strategy consisting of k-mer lengths of 25, 35, 45, 55, 65, 75, and 85. Eleven metrics were used to evaluate each assembly strategy including three genome related evaluation metrics (contig number, N50 length, Contigs >1 kb, reads) and eight transcriptome evaluation metrics (mapped back to transcripts (RMBT), number of full length transcripts, number of open reading frames, Detonate RSEM-EVAL score, and percent alignment to the southern platyfish, Amazon molly, BUSCO and CEGMA databases). The assembly strategy that performed the best, that is it was within the top three of each evaluation metric, was the Bridger assembly (10 of 11) followed by the Oases SMK assembly (8 of 11), the Oases LMK assembly (6 of 11), the Trinity assembly (4 of 11), the SOAP LMK assembly (4 of 11), and the SOAP SMK assembly (3 of 11).

Conclusion

This study provides an in-depth multi k-mer strategy investigation concluding that the assembler itself had a greater impact than k-mer size regardless of the strategy employed. Additionally, the comprehensive performance transcriptome evaluation metrics utilized in this study identified the need for choosing metrics centered on user defined research goals. Based on the evaluation metrics performed, the Bridger assembly was able to construct the best assembly of the testis transcriptome in Fundulus heteroclitus.

Transcriptome Analysis of Differentially Expressed Genes Involved in Proanthocyanidin Accumulation in the Rhizomes of Fagopyrum dibotrys and an Irradiation-Induced Mutant

The rhizome of Fagopyrum dibotrys is a traditional Chinese medicine that has recently gained attention due to substantial findings regarding its bioactive proanthocyanidin (PA) compounds. However, the molecular mechanism underlying PA accumulation in F. dibotrys remains elusive. We previously obtained an irradiation-induced mutant (RM_R) of F. dibotrys that had a higher PA content compared to that of the wild-type (CK_R). The present study aimed to elucidate the molecular mechanism underlying PA accumulation in F. dibotrys by comparing the rhizome transcriptomes of the irradiation-induced mutant and wild-type using RNA-seq analysis. A total of 53,540 unigenes were obtained, of which 29,901 (55.84%) were annotated based on BLAST searches against public databases, and 501 unique sequences were differentially expressed between the two samples, which consisted of 204 up-regulated and 297 down-regulated unigenes. Further analysis showed that the expression patterns of some unigenes encoding enzymes involved in PAs biosynthesis in F. dibotrys rhizomes differed between RM_R and CK_R. In addition, we identified transcription factor families and several cytochrome P450s that may be involved in PA regulation in F. dibotrys. Finally, 12 unigenes that encode PA biosynthetic enzymes were confirmed by qRT-PCR analysis. This study sheds light on the molecular mechanism underlying radiation-mediated flavonoid accumulation and regulation in F. dibotrys rhizomes. These results will also provide a platform for further functional genomic research on this particular species.

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.

Large-scale transcriptional profiling of lignified tissues in Tectona grandis

BACKGROUND

Currently, Tectona grandis is one of the most valuable trees in the world and no transcript dataset related to secondary xylem is available. Considering how important the secondary xylem and sapwood transition from young to mature trees is, little is known about the expression differences between those successional processes and which transcription factors could regulate lignin biosynthesis in this tropical tree. Although MYB transcription factors are one of the largest superfamilies in plants related to secondary metabolism, it has not yet been characterized in teak. These results will open new perspectives for studies of diversity, ecology, breeding and genomic programs aiming to understand deeply the biology of this species.

RESULTS

We present a widely expressed gene catalog for T. grandis using Illumina technology and the de novo assembly. A total of 462,260 transcripts were obtained, with 1,502 and 931 genes differentially expressed for stem and branch secondary xylem, respectively, during age transition. Analysis of stem and branch secondary xylem indicates substantial similarity in gene ontologies including carbohydrate enzymes, response to stress, protein binding, and allowed us to find transcription factors and heat-shock proteins differentially expressed. TgMYB1 displays a MYB domain and a predicted coiled-coil (CC) domain, while TgMYB2, TgMYB3 and TgMYB4 showed R2R3-MYB domain and grouped with MYBs from several gymnosperms and flowering plants. TgMYB1, TgMYB4 and TgCES presented higher expression in mature secondary xylem, in contrast with TgMYB2, TgHsp1, TgHsp2, TgHsp3, and TgBi whose expression is higher in young lignified tissues. TgMYB3 is expressed at lower level in secondary xylem.

CONCLUSIONS

Expression patterns of MYB transcription factors and heat-shock proteins in lignified tissues are dissimilar when tree development was evaluated, obtaining more expression of TgMYB1 and TgMYB4 in lignified tissues of 60-year-old trees, and more expression in TgHsp1, TgHsp2, TgHsp3 and TgBi in stem secondary xylem of 12-year-old trees. We are opening a door for further functional characterization by reverse genetics and marker-assisted selection with those genes. Investigation of some of the key regulators of lignin biosynthesis in teak, however, could be a valuable step towards understanding how rigidity of teak wood and extractives content are different from most other woods. The obtained transcriptome data represents new sequences of T. grandis deposited in public databases, representing an unprecedented opportunity to discover several related-genes associated with secondary xylem such as transcription factors and stress-related genes in a tropical tree.

De Novo Assembly and Characterization of Fruit Transcriptome in Black Pepper (Piper nigrum)

Black pepper is one of the most popular and oldest spices in the world and valued for its pungent constituent alkaloids. Pinerine is the main bioactive compound in pepper alkaloids, which perform unique physiological functions. However, the mechanisms of piperine synthesis are poorly understood. This study is the first to describe the fruit transcriptome of black pepper by sequencing on Illumina HiSeq 2000 platform. A total of 56,281,710 raw reads were obtained and assembled. From these raw reads, 44,061 unigenes with an average length of 1,345 nt were generated. During functional annotation, 40,537 unigenes were annotated in Gene Ontology categories, Kyoto Encyclopedia of Genes and Genomes pathways, Swiss-Prot database, and Nucleotide Collection (NR/NT) database. In addition, 8,196 simple sequence repeats (SSRs) were detected. In a detailed analysis of the transcriptome, housekeeping genes for quantitative polymerase chain reaction internal control, polymorphic SSRs, and lysine/ornithine metabolism-related genes were identified. These results validated the availability of our database. Our study could provide useful data for further research on piperine synthesis in black pepper.

Transcriptomic screening for cyclotides and other cysteine-rich proteins in the metallophyte Viola baoshanensis

Cysteine (Cys)-rich proteins (CRPs) are frequently associated with plant defense and stress resistance. Viola baoshanensis is a cadmium (Cd) hyper-accumulating plant whose CRPs-based defense systems are so far poorly understood. Next generation sequencing (NGS) techniques and a specialist searching tool, CrpExcel, were employed for identifying CRPs in V. baoshanensis. The transcriptome sequences of V. baoshanensis were assembled primarily from 454FLX/Hiseq2000 reads of plant cDNA sequencing libraries. CrpExcel was then used to search the ORFs and 9687 CRPs were identified, and included zinc finger (ZF) proteins, lipid transfer proteins, thaumatins and cyclotide precursors. Real-time PCR results showed that all CRP genes tested are constitutively expressed, but the genes of defensive peptides showed greater up-regulated expression than those of ZF-proteins in Cd- and/or wounding (Wd) treatments of V. baoshanensis seedlings. The NGS-derived sequences of cyclotide precursor genes were verified by RT-PCR and ABI3730 sequencing studies, and 32 novel cyclotides were identified in V. baoshanensis. In general, the metal-binding sites of ZF-containing CRPs also represented the potential vulnerable targets of toxic metals. This study provides broad insights into CRPs-based defense systems and stress-vulnerable targets in V. baoshanensis. It now brings the number of cyclotide sequences in V. baoshanensis to 53 and based on projections from this work, the number of cyclotides in the Violaceae is now conservatively estimated to be >30000.

Comparison of root transcriptomes and expressions of genes involved in main medicinal secondary metabolites from Bupleurum chinense and Bupleurum scorzonerifolium, the two Chinese official Radix bupleuri source species

Radix bupleuri, roots of Bupleurum species, is a widely used traditional Chinese medicine. Here, we compared the root transcriptomes of both Bupleurum chinense DC. and Bupleurum scorzonerifolium Willd. A total of 313 483 and 342 263 high quality expressed sequence tags were obtained, respectively. In addition, 17 117 (59.2%) and 19 416 (62.8%) unigenes for B. chinense and B. scorzonerifolium had homologous genes in the opposite dataset. For B. chinense, Kyoto Encyclopedia of Genes and Genomes database (KEGG) annotation identified carbohydrate metabolism, energy metabolism and amino acid metabolism as the three highest groups in the metabolism category. For B. scorzonerifolium, the lipid metabolism group had the most unigenes. Genes that may participate in the biosynthesis of terpenoid, triterpenoid, sterol, lignan and flavonoids were identified according to their annotations. (Tri)terpenoid-related genes were predominantly found in B. chinense. The expressions of certain genes were analyzed by quantitative reverse transcription-polymerase chain reaction (qRT-PCR) in the roots of the two species. A total of 558 putative transcription factors (TFs) and 137 transcriptional regulators (TRs) among 1364 TFs and 327 TRs, and 610 TFs and 129 TRs among 1600 TFs and 323 TRs were specific for B. chinense and B. scorzonerifolium, respectively. Our transcriptome comparison reflects the different types and proportions of metabolites synthesized by the two species. The data, especially, those genes involved in the biosynthesis of particular types of metabolites, will provide the basis for further investigations of the secondary metabolite repertoire of the two Bupleurum species, as well as other species from the genus of Bupleurum.

De novo assembly of Eugenia uniflora L. transcriptome and identification of genes from the terpenoid biosynthesis pathway

Pitanga (Eugenia uniflora L.) is a member of the Myrtaceae family and is of particular interest due to its medicinal properties that are attributed to specialized metabolites with known biological activities. Among these molecules, terpenoids are the most abundant in essential oils that are found in the leaves and represent compounds with potential pharmacological benefits. The terpene diversity observed in Myrtaceae is determined by the activity of different members of the terpene synthase and oxidosqualene cyclase families. Therefore, the aim of this study was to perform a de novo assembly of transcripts from E. uniflora leaves and to annotation to identify the genes potentially involved in the terpenoid biosynthesis pathway and terpene diversity. In total, 72,742 unigenes with a mean length of 1048bp were identified. Of these, 43,631 and 36,289 were annotated with the NCBI non-redundant protein and Swiss-Prot databases, respectively. The gene ontology categorized the sequences into 53 functional groups. A metabolic pathway analysis with KEGG revealed 8,625 unigenes assigned to 141 metabolic pathways and 40 unigenes predicted to be associated with the biosynthesis of terpenoids. Furthermore, we identified four putative full-length terpene synthase genes involved in sesquiterpenes and monoterpenes biosynthesis, and three putative full-length oxidosqualene cyclase genes involved in the triterpenes biosynthesis. The expression of these genes was validated in different E. uniflora tissues.

De novo assembly of the perennial ryegrass transcriptome using an RNA-Seq strategy

BACKGROUND

Perennial ryegrass is a highly heterozygous outbreeding grass species used for turf and forage production. Heterozygosity can affect de-Bruijn graph assembly making de novo transcriptome assembly of species such as perennial ryegrass challenging. Creating a reference transcriptome from a homozygous perennial ryegrass genotype can circumvent the challenge of heterozygosity. The goals of this study were to perform RNA-sequencing on multiple tissues from a highly inbred genotype to develop a reference transcriptome. This was complemented with RNA-sequencing of a highly heterozygous genotype for SNP calling.

RESULT

De novo transcriptome assembly of the inbred genotype created 185,833 transcripts with an average length of 830 base pairs. Within the inbred reference transcriptome 78,560 predicted open reading frames were found of which 24,434 were predicted as complete. Functional annotation found 50,890 transcripts with a BLASTp hit from the Swiss-Prot non-redundant database, 58,941 transcripts with a Pfam protein domain and 1,151 transcripts encoding putative secreted peptides. To evaluate the reference transcriptome we targeted the high-affinity K+ transporter gene family and found multiple orthologs. Using the longest unique open reading frames as the reference sequence, 64,242 single nucleotide polymorphisms were found. One thousand sixty one open reading frames from the inbred genotype contained heterozygous sites, confirming the high degree of homozygosity.

CONCLUSION

Our study has developed an annotated, comprehensive transcriptome reference for perennial ryegrass that can aid in determining genetic variation, expression analysis, genome annotation, and gene mapping.

Comparison of the transcriptomes of ginger (Zingiber officinale Rosc.) and mango ginger (Curcuma amada Roxb.) in response to the bacterial wilt infection

Bacterial wilt in ginger (Zingiber officinale Rosc.) caused by Ralstonia solanacearum is one of the most important production constraints in tropical, sub-tropical and warm temperature regions of the world. Lack of resistant genotype adds constraints to the crop management. However, mango ginger (Curcuma amada Roxb.), which is resistant to R. solanacearum, is a potential donor, if the exact mechanism of resistance is understood. To identify genes involved in resistance to R. solanacearum, we have sequenced the transcriptome from wilt-sensitive ginger and wilt-resistant mango ginger using Illumina sequencing technology. A total of 26387032 and 22268804 paired-end reads were obtained after quality filtering for C. amada and Z. officinale, respectively. A total of 36359 and 32312 assembled transcript sequences were obtained from both the species. The functions of the unigenes cover a diverse set of molecular functions and biological processes, among which we identified a large number of genes associated with resistance to stresses and response to biotic stimuli. Large scale expression profiling showed that many of the disease resistance related genes were expressed more in C. amada. Comparative analysis also identified genes belonging to different pathways of plant defense against biotic stresses that are differentially expressed in either ginger or mango ginger. The identification of many defense related genes differentially expressed provides many insights to the resistance mechanism to R. solanacearum and for studying potential pathways involved in responses to pathogen. Also, several candidate genes that may underline the difference in resistance to R. solanacearum between ginger and mango ginger were identified. Finally, we have developed a web resource, ginger transcriptome database, which provides public access to the data. Our study is among the first to demonstrate the use of Illumina short read sequencing for de novo transcriptome assembly and comparison in non-model species of Zingiberaceae.

RNA sequencing analysis of the gametophyte transcriptome from the liverwort, Marchantia polymorpha

The liverwort Marchantia polymorpha is a member of the most basal lineage of land plants (embryophytes) and likely retains many ancestral morphological, physiological and molecular characteristics. Despite its phylogenetic importance and the availability of previous EST studies, M. polymorpha's lack of economic importance limits accessible genomic resources for this species. We employed Illumina RNA-Seq technology to sequence the gametophyte transcriptome of M. polymorpha. cDNA libraries from 6 different male and female developmental tissues were sequenced to delineate a global view of the M. polymorpha transcriptome. Approximately 80 million short reads were obtained and assembled into a non-redundant set of 46,533 transcripts (> =  200 bp) from 46,070 loci. The average length and the N50 length of the transcripts were 757 bp and 471 bp, respectively. Sequence comparison of assembled transcripts with non-redundant proteins from embryophytes resulted in the annotation of 43% of the transcripts. The transcripts were also compared with M. polymorpha expressed sequence tags (ESTs), and approximately 69.5% of the transcripts appeared to be novel. Twenty-one percent of the transcripts were assigned GO terms to improve annotation. In addition, 6,112 simple sequence repeats (SSRs) were identified as potential molecular markers, which may be useful in studies of genetic diversity. A comparative genomics approach revealed that a substantial proportion of the genes (35.5%) expressed in M. polymorpha were conserved across phylogenetically related species, such as Selaginella and Physcomitrella, and identified 580 genes that are potentially unique to liverworts. Our study presents an extensive amount of novel sequence information for M. polymorpha. This information will serve as a valuable genomics resource for further molecular, developmental and comparative evolutionary studies, as well as for the isolation and characterization of functional genes that are involved in sex differentiation and sexual reproduction in this liverwort.

RNA sequencing read depth requirement for optimal transcriptome coverage in Hevea brasiliensis

BACKGROUND

One of the concerns of assembling de novo transcriptomes is determining the amount of read sequences required to ensure a comprehensive coverage of genes expressed in a particular sample. In this report, we describe the use of Illumina paired-end RNA-Seq (PE RNA-Seq) reads from Hevea brasiliensis (rubber tree) bark to devise a transcript mapping approach for the estimation of the read amount needed for deep transcriptome coverage.

FINDINGS

We optimized the assembly of a Hevea bark transcriptome based on 16 Gb Illumina PE RNA-Seq reads using the Oases assembler across a range of k-mer sizes. We then assessed assembly quality based on transcript N50 length and transcript mapping statistics in relation to (a) known Hevea cDNAs with complete open reading frames, (b) a set of core eukaryotic genes and (c) Hevea genome scaffolds. This was followed by a systematic transcript mapping process where sub-assemblies from a series of incremental amounts of bark transcripts were aligned to transcripts from the entire bark transcriptome assembly. The exercise served to relate read amounts to the degree of transcript mapping level, the latter being an indicator of the coverage of gene transcripts expressed in the sample. As read amounts or datasize increased toward 16 Gb, the number of transcripts mapped to the entire bark assembly approached saturation. A colour matrix was subsequently generated to illustrate sequencing depth requirement in relation to the degree of coverage of total sample transcripts.

CONCLUSIONS

We devised a procedure, the "transcript mapping saturation test", to estimate the amount of RNA-Seq reads needed for deep coverage of transcriptomes. For Hevea de novo assembly, we propose generating between 5-8 Gb reads, whereby around 90% transcript coverage could be achieved with optimized k-mers and transcript N50 length. The principle behind this methodology may also be applied to other non-model plants, or with reads from other second generation sequencing platforms.

Conserved changes in the dynamics of metabolic processes during fruit development and ripening across species

Computational analyses of molecular phenotypes traditionally aim at identifying biochemical components that exhibit differential expression under various scenarios (e.g. environmental and internal perturbations) in a single species. High-throughput metabolomics technologies allow the quantification of (relative) metabolite levels across developmental stages in different tissues, organs, and species. Novel methods for analyzing the resulting multiple data tables could reveal preserved dynamics of metabolic processes across species. The problem we address in this study is 2-fold. (1) We derive a single data table, referred to as a compromise, which captures information common to the investigated set of multiple tables containing data on different fruit development and ripening stages in three climacteric (i.e. peach [Prunus persica] and two tomato [Solanum lycopersicum] cultivars, Ailsa Craig and M82) and two nonclimacteric (i.e. strawberry [Fragaria × ananassa] and pepper [Capsicum chilense]) fruits; in addition, we demonstrate the power of the method to discern similarities and differences between multiple tables by analyzing publicly available metabolomics data from three tomato ripening mutants together with two tomato cultivars. (2) We identify the conserved dynamics of metabolic processes, reflected in the data profiles of the corresponding metabolites that contribute most to the determined compromise. Our analysis is based on an extension to principal component analysis, called STATIS, in combination with pathway overenrichment analysis. Based on publicly available metabolic profiles for the investigated species, we demonstrate that STATIS can be used to identify the metabolic processes whose behavior is similarly affected during fruit development and ripening. These findings ultimately provide insights into the pathways that are essential during fruit development and ripening across species.

cDNA normalization by hydroxyapatite chromatography to enrich transcriptome diversity in RNA-seq applications

Second-generation sequencing (SGS) has become the preferred method for RNA transcriptome profiling of organisms and single cells. However, SGS analysis of transcriptome diversity (including protein-coding transcripts and regulatory non-coding RNAs) is inefficient unless the sample of interest is first depleted of nucleic acids derived from ribosomal RNA (rRNA), which typically account for up to 95% of total intracellular RNA content. Here we describe a novel microscale hydroxyapatite chromatography (HAC) normalization method to remove eukaryotic and prokaryotic high abundant rRNA species, thereby increasing sequence coverage depth and transcript diversity across non-rRNA populations. RNA-seq analysis of Escherichia coli K-12 and human intracellular total RNA showed that HAC-based normalization enriched for all non-ribosomal RNA species regardless of RNA transcript abundance or length when compared with untreated controls. Microcolumn HAC normalization generated rRNA-depleted cDNA libraries comparable to the well-established duplex specific nuclease (DSN) normalization and Ribo-Zero rRNA-depletion methods, thus establishing microscale HAC as an effective, cost saving, and non-destructive alternative normalization technique.

De novo transcriptome sequencing reveals a considerable bias in the incidence of simple sequence repeats towards the downstream of 'Pre-miRNAs' of black pepper

Next generation sequencing has an advantageon transformational development of species with limited available sequence data as it helps to decode the genome and transcriptome. We carried out the de novo sequencing using illuminaHiSeq™ 2000 to generate the first leaf transcriptome of black pepper (Piper nigrum L.), an important spice variety native to South India and also grown in other tropical regions. Despite the economic and biochemical importance of pepper, a scientifically rigorous study at the molecular level is far from complete due to lack of sufficient sequence information and cytological complexity of its genome. The 55 million raw reads obtained, when assembled using Trinity program generated 2,23,386 contigs and 1,28,157 unigenes. Reports suggest that the repeat-rich genomic regions give rise to small non-coding functional RNAs. MicroRNAs (miRNAs) are the most abundant type of non-coding regulatory RNAs. In spite of the widespread research on miRNAs, little is known about the hair-pin precursors of miRNAs bearing Simple Sequence Repeats (SSRs). We used the array of transcripts generated, for the in silico prediction and detection of ‘43 pre-miRNA candidates bearing different types of SSR motifs’. The analysis identified 3913 different types of SSR motifs with an average of one SSR per 3.04 MB of thetranscriptome. About 0.033% of the transcriptome constituted ‘pre-miRNA candidates bearing SSRs’. The abundance, type and distribution of SSR motifs studied across the hair-pin miRNA precursors, showed a significant bias in the position of SSRs towards the downstream of predicted ‘pre-miRNA candidates’. The catalogue of transcripts identified, together with the demonstration of reliable existence of SSRs in the miRNA precursors, permits future opportunities for understanding the genetic mechanism of black pepper and likely functions of ‘tandem repeats’ in miRNAs.