An Improved Genome Assembly of Azadirachta indica A. Juss

Fruit-In-Sight: A deep learning-based framework for secondary metabolite class prediction using fruit and leaf images

Fruits produce a wide variety of secondary metabolites of great economic value. Analytical measurement of the metabolites is tedious, time-consuming, and expensive. Additionally, metabolite concentrations vary greatly from tree to tree, making it difficult to choose trees for fruit collection. The current study tested whether deep learning-based models can be developed using fruit and leaf images alone to predict a metabolite's concentration class (high or low). We collected fruits and leaves (n = 1045) from neem trees grown in the wild across 0.6 million sq km, imaged them, and measured concentration of five metabolites (azadirachtin, deacetyl-salannin, salannin, nimbin and nimbolide) using high-performance liquid chromatography. We used the data to train deep learning models for metabolite class prediction. The best model out of the seven tested (YOLOv5, GoogLeNet, InceptionNet, EfficientNet_B0, Resnext_50, Resnet18, and SqueezeNet) provided a validation F1 score of 0.93 and a test F1 score of 0.88. The sensitivity and specificity of the fruit model alone in the test set were 83.52 ± 6.19 and 82.35 ± 5.96, and 79.40 ± 8.50 and 85.64 ± 6.21, for the low and the high classes, respectively. The sensitivity was further boosted to 92.67± 5.25 for the low class and 88.11 ± 9.17 for the high class, and the specificity to 100% for both classes, using a multi-analyte framework. We incorporated the multi-analyte model in an Android mobile App Fruit-In-Sight that uses fruit and leaf images to decide whether to 'pick' or 'not pick' the fruits from a specific tree based on the metabolite concentration class. Our study provides evidence that images of fruits and leaves alone can predict the concentration class of a secondary metabolite without using expensive laboratory equipment and cumbersome analytical procedures, thus simplifying the process of choosing the right tree for fruit collection.

Deciphering the key pathway for triterpenoid biosynthesis in Azadirachta indica A. Juss.: a comprehensive review of omics studies in nature's pharmacy

Since ancient times, Azadirachta indica, or Neem, has been a well-known species of plant that produces a broad range of bioactive terpenoid chemicals that are involved in a variety of biological functions. Understanding the molecular mechanisms that are responsible for the biosynthesis and control of terpenoid synthesis is majorly dependent on successfully identifying the genes that are involved in their production. This review provides an overview of the recent developments concerning the identification of genes in A. indica that are responsible for the production of terpenoids. Numerous candidate genes encoding enzymes that are involved in the terpenoid biosynthesis pathway have been found through the use of transcriptomic and genomic techniques. These candidate genes include those that are responsible for the precursor synthesis, cyclization, and modification of terpenoid molecules. In addition, cutting-edge omics technologies, such as metabolomics and proteomics, have helped to shed light on the intricate regulatory networks that govern terpenoid biosynthesis. These networks are responsible for the production of terpenoids. The identification and characterization of genes involved in terpenoid biosynthesis in A. indica presents potential opportunities for genetic engineering and metabolic engineering strategies targeted at boosting terpenoid production as well as discovering novel bioactive chemicals.

De novo transcriptome analysis of Dysoxylum binectariferum to unravel the biosynthesis of pharmaceutically relevant specialized metabolites

The tropical tree, D. binectariferum, is a prominent source of chromone alkaloid rohitukine, which is used in the semi-syntheses of anticancer molecules such as flavopiridol and P-276-00. The biosynthetic pathway of rohitukine or its derivatives is currently unknown in plants. Here, we explored chromone alkaloid biosynthesis in D. binectariferum through targeted transcriptome sequencing. Illumina sequencing of leaves and roots of a year-old D. binectariferum seedling generated, 42.43 and 38.74 million paired-end short reads, respectively. Quality filtering and de novo assembly of the transcriptome generated 274,970 contigs and 126,788 unigenes with an N50 contig length of 1560 bp. The assembly generated 117,619 translated unigene protein sequences and 51,598 non-redundant sequences. Nearly 80% of these non-redundant sequences were annotated to publicly available protein and nucleotide databases, suggesting the completeness and effectiveness of the transcriptome assembly. Using the assembly, we identified a chalcone synthase (CHS) and three type III polyketide synthases (PKS-III; non-CHS type) that are likely to be involved in the biosynthesis of chromone ring/noreugenin moiety of rohitukine. We also identified key enzymes like lysine decarboxylase in the piperidine pathway that make the piperidine moiety of rohitukine. Besides these, the upstream enzymes in flavonoid biosynthesis like phenylalanine ammonia-lyase (PAL), trans-cinnamate 4-hydroxylase (C4H),4-coumarate-CoA ligase (4CL), and chalcone isomerase (CHI) have also been identified. Also, terpene synthases that are likely to be involved in the biosynthesis of various terpenoid scaffolds have been identified. Together, the D. binectariferum transcriptome resource forms a basis for further exploration of biosynthetic pathways of these valuable compounds through functional validation of the candidate genes and metabolic engineering in heterologous hosts. Additionally, the transcriptome dataset generated will serve as an important resource for research on functional genomics and enzyme discovery in D. binectariferum and comparative analysis with other Meliaceae family members.

Meliaceae genomes provide insights into wood development and limonoids biosynthesis

Meliaceae is a useful plant family owing to its high-quality timber and its many limonoids that have pharmacological and biological activities. Although some genomes of Meliaceae species have been reported, many questions regarding their unique family features, namely wood quality and natural products, have not been answered. In this study, we provide the whole-genome sequence of Melia azedarach comprising 237.16 Mb with a contig N50 of 8.07 Mb, and an improved genome sequence of Azadirachta indica comprising 223.66 Mb with a contig N50 of 8.91 Mb. Moreover, genome skimming data, transcriptomes and other published genomes were comprehensively analysed to determine the genes and proteins that produce superior wood and valuable limonoids. Phylogenetic analysis of chloroplast genomes, single-copy gene families and single-nucleotide polymorphisms revealed that Meliaceae should be classified into two subfamilies: Cedreloideae and Melioideae. Although the Meliaceae species did not undergo additional whole-genome duplication events, the secondary wall biosynthetic genes of the woody Cedreloideae species, Toona sinensis, expanded significantly compared to those of A. indica and M. azedarach, especially in downstream transcription factors and cellulose/hemicellulose biosynthesis-related genes. Moreover, expanded special oxidosqualene cyclase catalogues can help diversify Sapindales skeletons, and the clustered genes that regulate terpene chain elongation, cyclization and modification would support their roles in limonoid biosynthesis. The expanded clans of terpene synthase, O-methyltransferase and cytochrome P450, which are mainly derived from tandem duplication, are responsible for the different limonoid classes among the species. These results are beneficial for further investigations of wood development and limonoid biosynthesis.

Iridoid glycosides from Picrorhiza genus endemic to the Himalayan region: phytochemistry, biosynthesis, pharmacological potential and biotechnological intercessions to boost production

Iridoid glycosides are monoterpenoids synthesized in several plant species known to exhibit a diverse range of pharmacological activities. They are used as important bioactive ingredients in many commercially available drug formulations and as lead compounds in pharmaceutical research. The genus Picrorhiza comprises two medicinally important herbs endemic to the Himalayan region viz. Picrorhiza kurrooa Royle and Picrorhiza scrophulariiflora Hong. The medicinal properties of these two species are mainly due to iridoid glycosides present in their root, rhizome, and leaves. Unregulated harvesting from the wild, habitat specificity, narrow distribution range, small population size and lack of organized cultivation led to the enrolling of these species in the endangered category by the International Union for Conservation of Nature and Natural Resources (IUCN). Therefore, there is a need for immediate biotechnological and molecular interventions. Such intercessions will open up new vistas for large-scale propagation, development of genomic/transcriptomic resources for understanding the biosynthetic pathway, the possibility of genetic/metabolic manipulations, and possible commercialization of iridoid glycosides. The current review article elucidates the phytochemistry and pharmacological importance of iridoid glycosides from the genus Picrorhiza. In addition, the role of biotechnological approaches and opportunities offered by next-generation sequencing technologies in overcoming challenges associated with the genetic engineering of these species are also discussed.

Achievements and perspectives of synthetic biology in botanical insecticides

Botanical insecticides are the origin of all insecticidal compounds. They have been widely used to control pests in crops for a long time. Currently, the commercial production of botanical insecticides extracted from plants is limited because of insufficient raw material supply. Synthetic biology is a promising and effective approach for addressing the current problems of the production of botanical insecticides. It is an emerging biological research hotspot in the field of botanical insecticides. However, the biosynthetic pathways of many botanical insecticides are not completely elucidated. On the other hand, the cytotoxicity of botanical pesticides and low efficiency of these biosynthetic enzymes in new hosts make it still challenging for their heterologous production. In the present review, we summarized the recent developments in the heterologous production of botanical insecticides, analyzed the current challenges, and discussed the feasible production strategies, focusing on elucidating biosynthetic pathways, enzyme engineering, host engineering, and cytotoxicity engineering. Looking to the future, synthetic biology promises to further advance heterologous production of more botanical pesticides.

Genome structure and evolutionary history of frankincense producing Boswellia sacra

Boswellia sacra Flueck (family Burseraceae) tree is wounded to produce frankincense. We report its de novo assembled genome (667.8 Mb) comprising 18,564 high-confidence protein-encoding genes. Comparing conserved single-copy genes across eudicots suggest >97% gene space assembly of B. sacra genome. Evolutionary history shows B. sacra gene-duplications derived from recent paralogous events and retained from ancient hexaploidy shared with other eudicots. The genome indicated a major expansion of Gypsy retroelements in last 2 million years. The B. sacra genetic diversity showed four clades intermixed with a primary genotype—dominating most resin-productive trees. Further, the stem transcriptome revealed that wounding concurrently activates phytohormones signaling, cell wall fortification, and resin terpenoid biosynthesis pathways leading to the synthesis of boswellic acid—a key chemotaxonomic marker of Boswellia. The sequence datasets reported here will serve as a foundation to investigate the genetic determinants of frankincense and other resin-producing species in Burseraceae.

•

Assembly and architecture of frankincense producing Boswellia sacra Flueck

•

Comparative genomics and evolutionary history of frankincense tree within orders

•

Transcriptome of stem part and gene expression patterns of wounding to the tree

•

Resin biosynthesis pathway and related CYP450 enzymes and gene families

Research progress of meliaceous limonoids from 2011 to 2021

Covering: July 2010 to December 2021Limonoids, a kind of natural tetranortriterpenoids with diverse skeletons and valuable insecticidal and medicinal bioactivities, are the characteristic metabolites of most plants of the Meliaceae family. The chemistry and bioactivities of meliaceous limonoids are a continuing hot area of natural products research; to date, about 2700 meliaceous limonoids have been identified. In particular, more than 1600, including thirty kinds of novel rearranged skeletons, have been isolated and identified in the past decade due to their wide distribution and abundant content in Meliaceae plants and active biosynthetic pathways. In addition to the discovery of new structures, many positive medicinal bioactivities of meliaceous limonoids have been investigated, and extensive achievements regarding the chemical and biological synthesis have been made. This review summarizes the recent research progress in the discovery of new structures, medicinal and agricultural bioactivities, and chem/biosynthesis of limonoids from the plants of the Meliaceae family during the past decade, with an emphasis on the discovery of limonoids with novel skeletons, the medicinal bioactivities and mechanisms, and chemical synthesis. The structures, origins, and bioactivities of other new limonoids were provided as ESI. Studies published from July 2010 to December 2021 are reviewed, and 482 references are cited.

Genomic Analysis Based on Chromosome-Level Genome Assembly Reveals an Expansion of Terpene Biosynthesis of Azadirachta indica

Azadirachta indica (neem), an evergreen tree of the Meliaceae family, is a source of the potent biopesticide azadirachtin. The lack of a chromosome-level assembly impedes an in-depth understanding of its genome architecture and the comparative genomic analysis of A. indica. Here, a high-quality genome assembly of A. indica was constructed using a combination of data from Illumina, PacBio, and Hi-C technology, which is the first chromosome-scale genome assembly of A. indica. Based on the length of our assembly, the genome size of A. indica is estimated to be 281 Mb anchored to 14 chromosomes (contig N50 = 6 Mb and scaffold N50 = 19 Mb). The genome assembly contained 115 Mb repetitive elements and 25,767 protein-coding genes. Evolutional analysis revealed that A. indica didn't experience any whole-genome duplication (WGD) event after the core eudicot γ event, but some genes and genome segment might likely experienced recent duplications. The secondary metabolite clusters, TPS genes, and CYP genes were also identified. Comparative genomic analysis revealed that most of the A. indica-specific TPS genes and CYP genes were located on the terpene-related clusters on chromosome 13. It is suggested that chromosome 13 may play an important role in the specific terpene biosynthesis of A. indica. The gene duplication events may be responsible for the terpene biosynthesis expansion in A. indica. The genomic dataset and genomic analysis created for A. indica will shed light on terpene biosynthesis in A. indica and facilitate comparative genomic research of the family Meliaceae.

Building Biofoundry India: challenges and path forward

Biofoundry is a place where biomanufacturing meets automation. The highly modular structure of a biofoundry helps accelerate the design–build–test–learn workflow to deliver products fast and in a streamlined fashion. In this perspective, we describe our efforts to build Biofoundry India, where we see the facility add a substantial value in supporting research, innovation and entrepreneurship. We describe three key areas of our focus, harnessing the potential of non-expressing parts of the sequenced genomes, using deep learning in pathway reconstruction and synthesising enzymes and metabolites. Toward the end, we describe specific challenges in building such facility in India and the path to mitigate some of those working with the other biofoundries worldwide.

Chloroplast Genome Sequencing, Comparative Analysis, and Discovery of Unique Cytoplasmic Variants in Pomegranate (Punica granatum L.)

Here we report on comprehensive chloroplast (cp) genome analysis of 16 pomegranate (Punica granatum L.) genotypes representing commercial cultivars, ornamental and wild types, through large-scale sequencing and assembling using next-generation sequencing (NGS) technology. Comparative genome analysis revealed that the size of cp genomes varied from 158,593 bp (in wild, “1201” and “1181”) to 158,662 bp (cultivar, “Gul-e-Shah Red”) among the genotypes, with characteristic quadripartite structures separated by a pair of inverted repeats (IRs). The higher conservation for the total number of coding and non-coding genes (rRNA and tRNA) and their sizes, and IRs (IR-A and IR-B) were observed across all the cp genomes. Interestingly, high variations were observed in sizes of large single copy (LSC, 88,976 to 89,044 bp) and small single copy (SSC, 18,682 to 18,684 bp) regions. Although, the structural organization of newly assembled cp genomes were comparable to that of previously reported cp genomes of pomegranate (“Helow,” “Tunisia,” and “Bhagawa”), the striking differences were observed with the Lagerstroemia lines, viz., Lagerstroemia intermedia (NC_0346620) and Lagerstroemia speciosa (NC_031414), which clearly confirmed previous findings. Furthermore, phylogenetic analysis also revealed that members outside the genus Punica were clubbed into a separate clade. The contraction and expansion analysis revealed that the structural variations in IRs, LSC, and SSC have significantly accounted for the evolution of cp genomes of Punica and L. intermedia over the periods. Microsatellite survey across cp genomes resulted in the identification of a total of 233 to 234 SSRs, with majority of them being mono- (A/T or C/G, 164–165 numbers), followed by di- (AT/AT or AG/CT, 54), tri- (6), tetra- (8), and pentanucleotides (1). Furthermore, the comparative structural variant analyses across cp genomes resulted in the identification of many varietal specific SNP/indel markers. In summary, our study has offered a successful development of large-scale cp genomics resources to leverage future genetic, taxonomical, and phylogenetic studies in pomegranate.

Azadirachta indica MicroRNAs: Genome-Wide Identification, Target Transcript Prediction, and Expression Analyses

MicroRNAs are short, endogenous, non-coding RNAs, liable for essential regulatory function. Numerous miRNAs have been identified and studied in plants with known genomic or small RNA resources. Despite the availability of genomic and transcriptomic resources, the miRNAs have not been reported in the medicinal tree Azadirachta indica (Neem) till date. Here for the first time, we report extensive identification of miRNAs and their possible targets in A. indica which might help to unravel their therapeutic potential. A comprehensive search of miRNAs in the A. indica genome by C-mii tool was performed. Overall, 123 miRNAs classified into 63 families and their stem-loop hairpin structures were predicted. The size of the A. indica (ain)-miRNAs ranged between 19 and 23 nt in length, and their corresponding ain-miRNA precursor sequence MFEI value averaged as -1.147 kcal/mol. The targets of ain-miRNAs were predicted in A. indica as well as Arabidopsis thaliana plant. The gene ontology (GO) annotation revealed the involvement of ain-miRNA targets in developmental processes, transport, stress, and metabolic processes including secondary metabolism. Stem-loop qRT-PCR was carried out for 25 randomly selected ain-miRNAs and differential expression patterns were observed in different A. indica tissues. Expression of miRNAs and its targets shows negative correlation in a dependent manner.

Multi-tissue transcriptome analysis using hybrid-sequencing reveals potential genes and biological pathways associated with azadirachtin A biosynthesis in neem (azadirachta indica)

BACKGROUND

Azadirachtin A is a triterpenoid from neem tree exhibiting excellent activities against over 600 insect species in agriculture. The production of azadirachtin A depends on extraction from neem tissues, which is not an eco-friendly and sustainable process. The low yield and discontinuous supply of azadirachtin A impedes further applications. The biosynthetic pathway of azadirachtin A is still unknown and is the focus of our study.

RESULTS

We attempted to explore azadirachtin A biosynthetic pathway and identified the key genes involved by analyzing transcriptome data from five neem tissues through the hybrid-sequencing (Illumina HiSeq and Pacific Biosciences Single Molecule Real-Time (SMRT)) approach. Candidates were first screened by comparing the expression levels between the five tissues. After phylogenetic analysis, domain prediction, and molecular docking studies, 22 candidates encoding 2,3-oxidosqualene cyclase (OSC), alcohol dehydrogenase, cytochrome P450 (CYP450), acyltransferase, and esterase were proposed to be potential genes involved in azadirachtin A biosynthesis. Among them, two unigenes encoding homologs of MaOSC1 and MaCYP71CD2 were identified. A unigene encoding the complete homolog of MaCYP71BQ5 was reported. Accuracy of the assembly was verified by quantitative real-time PCR (qRT-PCR) and full-length PCR cloning.

CONCLUSIONS

By integrating and analyzing transcriptome data from hybrid-seq technology, 22 differentially expressed genes (DEGs) were finally selected as candidates involved in azadirachtin A pathway. The obtained reliable and accurate sequencing data provided important novel information for understanding neem genome. Our data shed new light on understanding the biosynthesis of other triterpenoids in neem trees and provides a reference for exploring other valuable natural product biosynthesis in plants.

Identification of key enzymes responsible for protolimonoid biosynthesis in plants: Opening the door to azadirachtin production

Limonoids are natural products made by members of the Meliaceae and Rutaceae families. Some limonoids (e.g., azadirachtin) are toxic to insects yet harmless to mammals. The use of limonoids in crop protection and other applications currently depends on extraction from limonoid-producing plants. Metabolic engineering offers opportunities to generate crop plants with enhanced insect resistance and also to produce high-value limonoids (e.g., for pharmaceutical use) by expression in heterologous hosts. However, to achieve this the enzymes responsible for limonoid biosynthesis must first be characterized. Here we identify 3 conserved enzymes responsible for the biosynthesis of the protolimonoid melianol, a precursor to limonoids, from Melia azedarach and Citrus sinensis, so paving the way for limonoid metabolic engineering and diversification.

Limonoids are natural products made by plants belonging to the Meliaceae (Mahogany) and Rutaceae (Citrus) families. They are well known for their insecticidal activity, contribution to bitterness in citrus fruits, and potential pharmaceutical properties. The best known limonoid insecticide is azadirachtin, produced by the neem tree (Azadirachta indica). Despite intensive investigation of limonoids over the last half century, the route of limonoid biosynthesis remains unknown. Limonoids are classified as tetranortriterpenes because the prototypical 26-carbon limonoid scaffold is postulated to be formed from a 30-carbon triterpene scaffold by loss of 4 carbons with associated furan ring formation, by an as yet unknown mechanism. Here we have mined genome and transcriptome sequence resources for 3 diverse limonoid-producing species (A. indica, Melia azedarach, and Citrus sinensis) to elucidate the early steps in limonoid biosynthesis. We identify an oxidosqualene cyclase able to produce the potential 30-carbon triterpene scaffold precursor tirucalla-7,24-dien-3β-ol from each of the 3 species. We further identify coexpressed cytochrome P450 enzymes from M. azedarach (MaCYP71CD2 and MaCYP71BQ5) and C. sinensis (CsCYP71CD1 and CsCYP71BQ4) that are capable of 3 oxidations of tirucalla-7,24-dien-3β-ol, resulting in spontaneous hemiacetal ring formation and the production of the protolimonoid melianol. Our work reports the characterization of protolimonoid biosynthetic enzymes from different plant species and supports the notion of pathway conservation between both plant families. It further paves the way for engineering crop plants with enhanced insect resistance and producing high-value limonoids for pharmaceutical and other applications by expression in heterologous hosts.

First complete chloroplast genomics and comparative phylogenetic analysis of Commiphora gileadensis and C. foliacea: Myrrh producing trees

Commiphora gileadensis and C. foliacea (family Burseraceae) are pantropical in nature and known for producing fragrant resin (myrrh). Both the tree species are economically and medicinally important however, least genomic understanding is available for this genus. Herein, we report the complete chloroplast genome sequences of C. gileadensis and C. foliacea and comparative analysis with related species (C. wightii and Boswellia sacra). A modified chloroplast DNA extraction method was adopted, followed with next generation sequencing, detailed bioinformatics and PCR analyses. The results revealed that the cp genome sizes of C. gileadensis and C. foliacea, are 160,268 and 160,249 bp, respectively, with classic quadripartite structures that comprises of inverted repeat's pair. Overall, the organization of these cp genomes, GC contents, gene order, and codon usage were comparable to other cp genomes in angiosperm. Approximately, 198 and 175 perfect simple sequence repeats were detected in C. gileadensis and C. foliacea genomes, respectively. Similarly, 30 and 25 palindromic, 15 and 25 forward, and 20 and 25 tandem repeats were determined in both the cp genomes, respectively. Comparison of these complete cp genomes with C. wightii and B. sacra revealed significant sequence resemblance and comparatively highest deviation in intergenic spacers. The phylo-genomic comparison showed that C. gileadensis and C. foliacea form a single clade with previously reported C. wightii and B. sacra from family Burseraceae. Current study reports for the first time the cp genomics of species from Commiphora, which could be helpful in understanding genetic diversity and phylogeny of this myrrh producing species.

An improved genome assembly for Larimichthys crocea reveals hepcidin gene expansion with diversified regulation and function

Larimichthys crocea (large yellow croaker) is a type of perciform fish well known for its peculiar physiological properties and economic value. Here, we constructed an improved version of the L. crocea genome assembly, which contained 26,100 protein-coding genes. Twenty-four pseudo-chromosomes of L. crocea were also reconstructed, comprising 90% of the genome assembly. This improved assembly revealed several expansions in gene families associated with olfactory detection, detoxification, and innate immunity. Specifically, six hepcidin genes (LcHamps) were identified in L. crocea, possibly resulting from lineage-specific gene duplication. All LcHamps possessed similar genomic structures and functional domains, but varied substantially with respect to expression pattern, transcriptional regulation, and biological function. LcHamp1 was associated specifically with iron metabolism, while LcHamp2s were functionally diverse, involving in antibacterial activity, antiviral activity, and regulation of intracellular iron metabolism. This functional diversity among gene copies may have allowed L. crocea to adapt to diverse environmental conditions.

Comparison of three assembly strategies for a heterozygous seedless grapevine genome assembly

BACKGROUND

De novo heterozygous assembly is an ongoing challenge requiring improved assembly approaches. In this study, three strategies were used to develop de novo Vitis vinifera 'Sultanina' genome assemblies for comparison with the inbred V. vinifera (PN40024 12X.v2) reference genome and a published Sultanina ALLPATHS-LG assembly (AP). The strategies were: 1) a default PLATANUS assembly (PLAT_d) for direct comparison with AP assembly, 2) an iterative merging strategy using METASSEMBLER to combine PLAT_d and AP assemblies (MERGE) and 3) PLATANUS parameter modifications plus GapCloser (PLAT*_GC).

RESULTS

The three new assemblies were greater in size than the AP assembly. PLAT*_GC had the greatest number of scaffolds aligning with a minimum of 95% identity and ≥1000 bp alignment length to V. vinifera (PN40024 12X.v2) reference genome. SNP analysis also identified additional high quality SNPs. A greater number of sequence reads mapped back with zero-mismatch to the PLAT_d, MERGE, and PLAT*_GC (>94%) than was found in the AP assembly (87%) indicating a greater fidelity to the original sequence data in the new assemblies than in AP assembly. A de novo gene prediction conducted using seedless RNA-seq data predicted > 30,000 coding sequences for the three new de novo assemblies, with the greatest number (30,544) in PLAT*_GC and only 26,515 for the AP assembly. Transcription factor analysis indicated good family coverage, but some genes found in the VCOST.v3 annotation were not identified in any of the de novo assemblies, particularly some from  the MYB and ERF families.

CONCLUSIONS

The PLAT_d and PLAT*_GC had a greater number of synteny blocks with the V. vinifera (PN40024 12X.v2) reference genome than AP or MERGE. PLAT*_GC provided the most contiguous assembly with only 1.2% scaffold N, in contrast to AP (10.7% N), PLAT_d (6.6% N) and Merge (6.4% N). A PLAT*_GC pseudo-chromosome assembly with chromosome alignment to the reference genome V. vinifera, (PN40024 12X.v2) provides new information for use in seedless grape genetic mapping studies. An annotated de novo gene prediction for the PLAT*_GC assembly, aligned with VitisNet pathways provides new seedless grapevine specific transcriptomic resource that has excellent fidelity with the seedless short read sequence data.

Transcriptome and metabolite analyses in Azadirachta indica: identification of genes involved in biosynthesis of bioactive triterpenoids

Azadirachta indica A. Juss, commonly known as Neem, is the reservoir of triterpenoids of economic importance. Metabolite analysis of different developmental stages of leaf and fruit suggests tissue-specific accumulation of the major triterpenoids in this important tree. Though biosynthesis of these complex molecules requires substrate flux from the isoprenoid pathway, enzymes involved in late biosynthetic steps remain uncharacterized. We established and analyzed transcriptome datasets from leaf and fruit and identified members of gene families involved in intermediate steps of terpenoid backbone biosynthesis and those related to secondary transformation leading to the tissue-specific triterpenoid biosynthesis. Expression analysis suggests differential expression of number of genes between leaf and fruit and probable participation in the biosynthesis of fruit-specific triterpenoids. Genome-wide analysis also identified members of gene families putatively involved in secondary modifications in late biosynthetic steps leading to the synthesis of highly oxygenated triterpenoids. Expression and molecular docking analyses suggest involvement of specific members of CYP450 family in secondary modifications for the biosynthesis of bioactive triterpenoids. This study generated rich genomic resource and identified genes involved in biosynthesis of important molecules, which will aid in the advancement of tools for functional genomics and elucidation of the biosynthesis of triterpenoid from this important tree.