Comprehensive analyses of genomes, transcriptomes and metabolites of neem tree

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.

Comprehensive Dissipation of Azadirachtin in Grapes and Tomatoes: The Effect of Bacillus thuringiensis and Tentative Identification of Unknown Metabolites

Neem oil is a biopesticide normally applied together with Bacillus thuringiensis (Bt). However, neither its dissipation nor the influence of Bt has been previously evaluated. In this study, dissipation of neem oil was investigated when it was applied alone or together with Bt at 3 and 22 °C. A methodology involving solid-liquid extraction and liquid chromatography-high-resolution mass spectrometry was developed for that purpose. The method was validated obtaining recoveries from 87 to 103%, with relative standard deviations lower than 19% and limits of quantification from 5 to 10 μg/kg. Azadirachtin A (AzA) dissipation was fit to a single first order, being faster when neem oil was applied together with Bt and at 22 °C (RL₅₀ = 12-21 days) than alone and at 3 °C (RL₅₀ = 14-25 days). Eight related compounds were found in real samples with similar dissipation curves compared to AzA, and five unknown metabolites were identified in degraded samples, with increasing concentrations during parent compound degradation.

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.

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.

De Novo Assembly and Species-Specific Marker Development as a Useful Tool for the Identification of Scutellaria L. Species

Scutellaria L. (family Lamiaceae) includes approximately 470 species found in most parts of the world and is commonly known as skullcaps. Scutellaria L. is a medicinal herb used as a folk remedy in Korea and East Asia, but it is difficult to identify and classify various subspecies by morphological methods. Since Scutellaria L. has not been studied genetically, to expand the knowledge of species in the genus Scutellaria L., de novo whole-genome assembly was performed in Scutellaria indica var. tsusimensis (H. Hara) Ohwi using the Illumina sequencing platform. We aimed to develop a molecular method that could be used to classify S.indica var. tsusimensis (H. Hara) Ohwi, S. indica L. and three other Scutellaria L. species. The assembly results for S.indica var. tsusimensis (H. Hara) Ohwi revealed a genome size of 318,741,328 bp and a scaffold N50 of 78,430. The assembly contained 92.08% of the conserved BUSCO core gene set and was estimated to cover 94.65% of the genome. The obtained genes were compared with previously registered Scutellaria nucleotide sequences and similar regions using the NCBI BLAST service, and a total of 279 similar nucleotide sequences were detected. By selecting the 279 similar nucleotide sequences and nine chloroplast DNA barcode genes, primers were prepared so that the size of the PCR product was 100 to 1000 bp. As a result, a species-specific primer set capable of distinguishing five species of Scutellaria L. was developed.

A complete sequence of mitochondrial genome of Neolamarckia cadamba and its use for systematic analysis

Neolamarckia cadamba is an important tropical and subtropical tree for timber industry in southern China and is also a medicinal plant because of the secondary product cadambine. N. cadamba belongs to Rubiaceae family and its taxonomic relationships with other species are not fully evaluated based on genome sequences. Here, we report the complete sequences of mitochondrial genome of N. cadamba, which is 414,980 bp in length and successfully assembled in two genome circles (109,836 bp and 305,144 bp). The mtDNA harbors 83 genes in total, including 40 protein-coding genes (PCGs), 31 transfer RNA genes, 6 ribosomal RNA genes, and 6 other genes. The base composition of the whole genome is estimated as 27.26% for base A, 22.63% for C, 22.53% for G, and 27.56% for T, with the A + T content of 54.82% (54.45% in the small circle and 54.79% in the large circle). Repetitive sequences account for ~ 0.14% of the whole genome. A maximum likelihood (ML) tree based on DNA sequences of 24 PCGs supports that N. cadamba belongs to order Gentianales. A ML tree based on rps3 gene of 60 species in family Rubiaceae shows that N. cadamba is more related to Cephalanthus accidentalis and Hymenodictyon parvifolium and belongs to the Cinchonoideae subfamily. The result indicates that N. cadamba is genetically distant from the species and genera of Rubiaceae in systematic position. As the first sequence of mitochondrial genome of N. cadamba, it will provide a useful resource to investigate genetic variation and develop molecular markers for genetic breeding in the future.

Development, characterization, functional annotation and validation of genomic and genic-SSR markers using de novo next generation sequencing in Melia dubia Cav

Melia dubia Cav. (Meliaceae), a fast-growing tropical tree finds use in plywood, pulp and high-value solid wood products. To increase its productivity, we must essentially capture genetic diversity and identify genotypes with superior wood properties. This study aimed to develop novel microsatellite markers from genomic data and validate the markers in M. dubia. Direct Seq-to-SSR approach was adopted and using an in-house Perl script, 426,390 SSR markers identified. For validation, selected 151 markers, of which 50 were genomic markers chosen randomly, and 101 were genic markers identified through BLAST2GO. Amplification was observed in all loci, and 81.4% generated high-quality, reproducible amplicons of the expected size. Out of 50 genomic markers, we used ten highly polymorphic markers to assess genetic diversity among 75 genotypes from three populations. One hundred fourteen alleles were recorded, with a moderate level of diversity and a positive fixation index. Twenty-nine genic markers representing 13 enzymes showing polymorphism for wood stiffness were selected for diversity assessment of 24 genotypes (12 genotypes each with high and low-stress wave velocity). The product size ranged from 87 to 279, covering the majority of the genome. Cluster and structure analysis segregated ~ 80% of the genotypes based on the trait. This is the first report of the development of genic markers from a genomic survey and has proved efficient in differentiating genotypes based on the trait. The markers developed in this study will be useful for genetic mapping, diversity estimation, marker-assisted selection for desired traits and breeding for wood traits in M. dubia.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02858-w.

Characterization of microRNAs from neem (Azadirachta indica) and their tissue-specific expression study in leaves and stem

Neem (Azadirachta indica) is a very popular traditional medicinal plant used since ancient times to treat numerous ailments. MicroRNAs (miRNAs) are highly conserved, non-coding, short RNA molecules that play important regulatory roles in plant development and metabolism. In this study, deploying a high stringent genome-wide computational-based approach and following a set of strict filtering norms a total of 44 potential conserved neem miRNAs belonging to 21 families and their corresponding 48 potential target transcripts were identified. Important targets include Squamosa promoter binding protein-like proteins, NAC, Scarecrow proteins, Auxin response factor, and F-box proteins. A biological network has also been developed to understand the miRNA-mediated gene regulation using the minimum free energy (MFE) values of the miRNA-target interaction. Moreover, six selected miRNAs were reported to be involved in secondary metabolism in other plant species (miR156a, miR156l, miR160, miR164, miR171, miR395) were validated by qPCR and their tissue-specific differential expression pattern was observed in leaves and stem. Except for ain-miR395, all the other miRNAs were found overexpressed in the stem as compared to leaves. To the best of our knowledge, this is the first report of neem miRNAs and we believe the finding of the present study will be useful for the functional genomic study of medicinal plants.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02839-z.

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.

The 3-hydroxy-3-methylglutaryl-coenzyme A reductase 5 gene from Malus domestica enhances oxidative stress tolerance in Arabidopsis thaliana

3-Hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) is the first rate-limiting enzyme regulating the synthesis of terpenoids upstream of the mevalonate (MVA) pathway. In higher plants, members of the HMGR genes families play an important role in plant growth and development and in response to various environmental stresses. In the present study, a novel HMGR gene, designated MdHMGR5, was isolated from apple (Malus domestica L.) and characterized. Expression of MdHMGR5 enhanced the activity of HMGR enzyme in transgenic Arabidopsis thaliana L. plants. Under oxidative stress, transgenic A. thaliana plants over-expressing MdHMGR5 had a higher germination rate, a longer main root length, higher chlorophyll and proline content, and higher activities of antioxidant enzymes. On the other hand, malondialdehyde (MDA) content, relative conductivity and reactive oxygen species (ROS) production rate were significantly lower than in wild type plants. These results indicated that over-expression of MdHMGR5 enhanced plant tolerance to oxidative stress by scavenging ROS in transgenic plants. Over-expression of MdHMGR5 also affected the expression levels of genes in mevalonic acid and 2C-methyl-D-erythritol 4-phosphate (MVA and MEP) pathways of A. thaliana plants. These results indicate that over-expression of MdHMGR5 enhances tolerance to oxidative stress by maintaining photosynthesis and scavenging ROS in transgenic A. thaliana 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.

Genomic resources for the Neotropical tree genus Cedrela (Meliaceae) and its relatives

BACKGROUND

Tree species in the genus Cedrela P. Browne are threatened by timber overexploitation across the Neotropics. Genetic identification of processed timber can be used to supplement wood anatomy to assist in the taxonomic and source validation of protected species and populations of Cedrela. However, few genetic resources exist that enable both species and source identification of Cedrela timber products. We developed several 'omic resources including a leaf transcriptome, organelle genome (cpDNA), and diagnostic single nucleotide polymorphisms (SNPs) that may assist the classification of Cedrela specimens to species and geographic origin and enable future research on this widespread Neotropical tree genus.

RESULTS

We designed hybridization capture probes to enrich for thousands of genes from both freshly preserved leaf tissue and from herbarium specimens across eight Meliaceae species. We first assembled a draft de novo transcriptome for C. odorata, and then identified putatively low-copy genes. Hybridization probes for 10,001 transcript models successfully enriched 9795 (98%) of these targets, and analysis of target capture efficiency showed that probes worked effectively for five Cedrela species, with each species showing similar mean on-target sequence yield and depth. The probes showed greater enrichment efficiency for Cedrela species relative to the other three distantly related Meliaceae species. We provide a set of candidate SNPs for species identification of four of the Cedrela species included in this analysis, and present draft chloroplast genomes for multiple individuals of eight species from four genera in the Meliaceae.

CONCLUSIONS

Deforestation and illegal logging threaten forest biodiversity globally, and wood screening tools offer enforcement agencies new approaches to identify illegally harvested timber. The genomic resources described here provide the foundation required to develop genetic screening methods for Cedrela species identification and source validation. Due to their transferability across the genus and family as well as demonstrated applicability for both fresh leaves and herbarium specimens, the genomic resources described here provide additional tools for studies examining the ecology and evolutionary history of Cedrela and related species in the Meliaceae.

De novo transcriptome sequencing and SSR markers development for Cedrela balansae C.DC., a native tree species of northwest Argentina

The endangered Cedrela balansae C.DC. (Meliaceae) is a high-value timber species with great potential for forest plantations that inhabits the tropical forests in Northwestern Argentina.Research on this species is scarce because of the limited genetic and genomic information available. Here, we explored the transcriptome of C. balansae using 454 GS FLX Titanium next-generation sequencing (NGS) technology. Following de novo assembling, we identified 27,111 non-redundant unigenes longer than 200 bp, and considered these transcripts for further downstream analysis. The functional annotation was performed searching the 27,111 unigenes against the NR-Protein and the Interproscan databases. This analysis revealed 26,977 genes with homology in at least one of the Database analyzed. Furthermore, 7,774 unigenes in 142 different active biological pathways in C. balansae were identified with the KEGG database. Moreover, after in silico analyses, we detected 2,663 simple sequence repeats (SSRs) markers. A subset of 70 SSRs related to important "stress tolerance" traits based on functional annotation evidence, were selected for wet PCR-validation in C. balansae and other Cedrela species inhabiting in northwest and northeast of Argentina (C. fissilis, C. saltensis and C. angustifolia). Successful transferability was between 77% and 93% and thanks to this study, 32 polymorphic functional SSRs for all analyzed Cedrela species are now available. The gene catalog and molecular markers obtained here represent a starting point for further research, which will assist genetic breeding programs in the Cedrela genus and will contribute to identifying key populations for its preservation.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

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.

Induction of Autophagy and Apoptosis via PI3K/AKT/TOR Pathways by Azadirachtin A in Spodoptera litura Cells

Azadirachtin is one of the most effective botanical insecticides and has been widely used in pest control. Toxicological reports show that azadirachtin can induce apoptosis in various insect cell lines. However, studies of azadirachtin-induced autophagy in cultured insect cells are lacking. This study reports that azadirachtin A significantly inhibits cell proliferation by inducing autophagic and apoptotic cell death in Spodoptera litura cultured cell line (SL-1 cell). Characteristic autophagolysosome and Atg8-PE (phosphatidylethanolamine) accumulation were observed by electron microscopy and western blotting, indicating that azadirachtin triggered autophagy in SL-1 cell. Furthermore, azadirachtin inhibited survival signaling by blocking the activation of PI3K, AKT and the down-stream target of rapamycin. Similar to the positive control of starvation, azadirachtin induced the activation of insulin receptor (InR) via a cellular feedback mechanism. In addition, the autophagy-related 5 (Atg5), a molecular switch of autophagy and apoptosis, was truncated (tAtg5) to trigger cytochrome c release into the cytoplasm under azadirachtin stress, which indicated that azadirachtin induced apoptosis through autophagy. Our findings suggest that azadirachtin primarily induced autophagy in SL-1 cell by dysregulating InR- and PI3K/AKT/TOR pathways, then stimulated apoptosis by activating tAtg5.

Gene expression profiling analysis reveals a crucial gene regulating metabolism in adventitious roots of neem (Azadirachta indica)

Neem is a widely used traditional plant containing bioactive secondary metabolites, especially azadirachtin.