Rediscovering medicinal plants' potential with OMICS: microsatellite survey in expressed sequence tags of eleven traditional plants with potent antidiabetic properties

Genome assembly and microsatellite marker development using Illumina and PacBio sequencing in Persicaria maackiana (Polygonaceae) from Korea

BACKGROUND

Persicaria maackiana (Regel) is a potential medicinal plant that exerts anti-diabetic effects. However, the lack of genomic information on P. maackiana hinders research at the molecular level.

OBJECTIVE

Herein, we aimed to construct a draft genome assembly and obtain comprehensive genomic information on P. maackiana using high-throughput sequencing tools PacBio Sequel II and Illumina.

METHODS

Persicaria maackiana samples from three natural populations in Gaecheon, Gichi, and Uiryeong reservoirs in South Korea were used to generate genomic DNA libraries, perform genome de novo assembly, gene ontology analysis, phylogenetic tree analysis, genotyping, and identify microsatellite markers.

RESULTS

The assembled P. maackiana genome yielded 32,179 contigs. Assessment of assembly integrity revealed 1503 (93.12%) complete Benchmarking Universal Single-Copy Orthologs. A total of 64,712 protein-coding genes were predicted and annotated successfully in the protein database. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs, 13,778 genes were annotated into 18 categories. Genes that activated AMPK were identified in the KEGG pathway. A total of 316,992 microsatellite loci were identified, and primers targeting the flanking regions were developed for 292,059 microsatellite loci. Of these, 150 primer sets were randomly selected for amplification, and 30 of these primer sets were identified as polymorphic. These primers amplified 3-9 alleles. The mean observed and expected heterozygosity were 0.189 and 0.593, respectively. Polymorphism information content values of the markers were 0.361-0.754.

CONCLUSION

Collectively, our study provides a valuable resource for future comparative genomics, phylogeny, and population studies of P. maackiana.

Crosstalk of Multi-Omics Platforms with Plants of Therapeutic Importance

From time immemorial, humans have exploited plants as a source of food and medicines. The World Health Organization (WHO) has recorded 21,000 plants with medicinal value out of 300,000 species available worldwide. The promising modern "multi-omics" platforms and tools have been proven as functional platforms able to endow us with comprehensive knowledge of the proteome, genome, transcriptome, and metabolome of medicinal plant systems so as to reveal the novel connected genetic (gene) pathways, proteins, regulator sequences and secondary metabolite (molecule) biosynthetic pathways of various drug and protein molecules from a variety of plants with therapeutic significance. This review paper endeavors to abridge the contemporary advancements in research areas of multi-omics and the information involved in decoding its prospective relevance to the utilization of plants with medicinal value in the present global scenario. The crosstalk of medicinal plants with genomics, transcriptomics, proteomics, and metabolomics approaches will be discussed.

Medico-Religious Plants Employed in Mauritius: A Survey Among Hindu Priests

Mauritius is renowned for its rich floral and cultural diversity. The aim of this study was to identify and document medicinal plants exclusively employed in rituals and to promote wellness among the Hindu community. In addition, it also highlights poorly documented plants with respect to their safety, efficacy, and herb-drug interaction inducing potential. The informant group included Hindu priests (n = 14) solely, while the relative frequency of citation (RFC) and use value (UV) were utilised to analyse the data. A total of 33 plants were reported in this survey. Sixteen plants were reported by informants to have both medicinal and religious values, while the highest RFC and UV index was obtained by Ocimum tenuiflorum L. (RFC = 0.79; UV = 0.71). Nonetheless, despite known and utilised commonly for years by humanity, various plant species mentioned in this study have been poorly evaluated in randomised clinical trials and for any risks of herb-drug interactions. Hence, there are concerns surrounding the safe use of these plant species particularly among patients suffering from chronic diseases.

Carica papaya: Potential Implications in Human Health

Background:

Carica papaya, a tree-like herb, is cultivated in more than 50 tropical and subtropical countries worldwide. The parts [leaves, fruit (ripe and unripe), seeds and latex] are used as food as well as traditional medicine in several ways for a number of diseases. Papaya, with a variety of phytochemicals like carotenoids, polyphenols, benzyl isothiocynates, benzyl glucosinates, prunasin (cyanogenic substrate), papain and chymopapain, alkaloids, phenolic compounds, flavonoids, vitamins (A, C, E), carotenoids, cyanogenicglucosides, cystatin, and glucosinolates exhibits significant health benefits ranging from digestive to immune modulation. These compounds have antioxidant, chemoprotective, anti-diabetic, anti-bacterial, anti-plasmodial and anti-fungal activities. Papaya aqueous leaf extract was evaluated for immunomodulatory and anti-tumor activities through cytokine modulation with Th1 type shift in human immune system through molecular studies. The platelet augmenting potential of aqueous leaf extract has been reported in numerous clinical studies and deserves special mention.

Objective:

The scientific knowledge of carica papaya in the post-genomic era including molecular studies and clinical trials is discussed in the review.

Methods:

The published literature on botany, chemical composition, ethnopharmacology, and uses of papaya in food and medical industry was searched through databases like Pub- Med, Scopus, and Google scholar to comprehend the benefits of Carica papaya for human use with around more than 600 published peer- reviewed papers.

Results:

Many traditional and novel uses of Carica papaya for the human benefit are detailed in the review that significantly adds to the scientific knowledge of curious readers.

Conclusion:

Each component of the papaya plant is rich in phytochemicals and is economically important. Most of the phytochemicals are linked to biological functions and influence a variety of cellular processes; hence they have implications in refining human health.

Mining and comparative survey of EST-SSR markers among members of Euphorbiaceae family

Euphorbiaceae represents flowering plants family of tropical and sub-tropical region rich in secondary metabolites of economic importance. To understand and assess the genetic makeup among the members, this study was undertaken to characterize and compare SSR markers from publicly available ESTs and GSSs of nine selected species of the family. Mining of SSRs was performed by MISA, primer designing by Primer3, while functional annotation, gene ontology (GO) and enrichment analysis were performed by Blast2GO. A total 12,878 number of SSRs were detected from 101,701 number of EST sequences. SSR density ranged from 1 SSR/3.22 kb to 1 SSR/15.65 kb. A total of 1873 primer pairs were designed for the annotated SSR-Contigs. About 77.07% SSR-ESTs could be assigned a significant match to the protein database. 3037 unique SSR-FDM were assigned and IPR003657 (WRKY Domain) was found to be the most dominant FDM among the members. 1810 unique GO terms obtained were further subjected to enrichment analysis to obtain 513 statistically significant GO terms mapped to the SSR containing ESTs. Most frequent enriched GO terms were, GO:0003824 for molecular function, GO:0006350 for biological process and GO:0005886 for cellular component, justifying the richness of defensive secondary metabolites and phytomedicine within the family. The results from this study provides tangible insight to genetic make-up and distribution of SSRs. Functional annotation corresponded many genes of unknown functions which may be considered as novel genes or genes responsible for stress specific secondary metabolites. Further studies are required to understand stress specific genes accountable for leveraging the synthesis of secondary metabolites.

Repurposing L-Menthol for Systems Medicine and Cancer Therapeutics? L-Menthol Induces Apoptosis through Caspase 10 and by Suppressing HSP90

The objective of the present study was to repurpose L-menthol, which is frequently used in oral health and topical formulations, for cancer therapeutics. In this article, we argue that monoterpenes such as L-menthol might offer veritable potentials in systems medicine, for example, as cheaper anti-cancer compounds. Other monoterpenes such as limonene, perillyl alcohol, and geraniol have been shown to induce apoptosis in various cancer cell lines, but their mechanisms of action are yet to be completely elucidated. Earlier, we showed that L-menthol modulates tubulin polymerization and apoptosis to inhibit cancer cell proliferation. In the present report, we used an apoptosis-related gene microarray in conjunction with proteomics analyses, as well as in silico interpretations, to study gene expression modulation in human adenocarcinoma Caco-2 cell line in response to L-menthol treatment. The microarray analysis identified caspase 10 as the important initiator caspase, instead of caspase 8. The proteomics analyses showed downregulation of HSP90 protein (also corroborated by its low transcript abundance), which in turn indicated inhibition of AKT-mediated survival pathway, release of pro-apoptotic factor BAD from BAD and BCLxL complex, besides regulation of other factors related to apoptosis. Based on the combined microarray, proteomics, and in silico data, a signaling pathway for L-menthol-induced apoptosis is being presented for the first time here. These data and literature analysis have significant implications for "repurposing" L-menthol beyond oral medicine, and in understanding the mode of action of plant-derived monoterpenes towards development of cheaper anticancer drugs in future.

Designing New Kinase Inhibitor Derivatives as Therapeutics Against Common Complex Diseases: Structural Basis of Microtubule Affinity-Regulating Kinase 4 (MARK4) Inhibition

Drug development for common complex diseases is in need of new molecular entities and actionable drug targets. MAP/microtubule affinity-regulating kinase 4 (MARK4) is associated with numerous diseases such as neurodegenerative disorders, obesity, cancer, and type 2 diabetes. Understanding the structural basis of ligands' (inhibitors) and substrates' binding to MARK4 is crucial to design new kinase inhibitors for therapeutic purposes. This study reports new observations on docking three well-known kinase inhibitors in the kinase domain of MARK4 variants and the calculated binding affinity. These variants of MARK4 are named as MARK4-F1 (59 N-terminal residues along with kinase domain) and MARK4-F2 (kinase domain of MARK4). We additionally performed molecular dynamics (MD) simulation and fluorescence binding studies to calculate the actual binding affinity of kinase inhibitors, BX-912, BX-795, and OTSSP167 (hydrochloride) for the MARK4. Docking analyses revealed that ligands bind in the large hydrophobic cavity of the kinase domain of MARK4 through several hydrophobic and hydrogen-bonded interactions. Simulations suggested that OTSSP167 (hydrochloride) is forming a stable complex, and hence the best inhibitor of MARK4. Intrinsic fluorescence of MARK4 was significantly quenched by addition of ligands, indicating their potential binding to MARK4. A lower KD value of MARK4 with OTSSP167 (hydrochloride) suggested that it is a better interacting partner than BX-912 and BX-795. These data form a basis for designing novel and potent OTSSP167 (hydrochloride) derivatives as therapeutic candidates against common complex diseases. The inhibitors designed as such might possibly suppress the growth of tumor-forming cells and be potentially applied for treatment of a wide range of human cancers as well.

Tissue-Specific Molecular Biomarker Signatures of Type 2 Diabetes: An Integrative Analysis of Transcriptomics and Protein-Protein Interaction Data

Type 2 diabetes mellitus is a major global public health burden. A complex metabolic disease, type 2 diabetes affects multiple different tissues, demanding a "systems medicine" approach to biomarker and novel diagnostic discovery, not to mention data integration across omics-es. In the present study, transcriptomics data from different tissues including beta-cells, pancreatic islets, arterial tissue, peripheral blood mononuclear cells, liver, and skeletal muscle of 228 samples were integrated with protein-protein interaction data and genome scale metabolic models to unravel the molecular and tissue-specific biomarker signatures of type 2 diabetes mellitus. Classifying differentially expressed genes, reconstruction and topological analysis of active protein-protein interaction subnetworks indicated that genomic reprogramming depends on the type of tissue, whereas there are common signatures at different levels. Among all tissue and cell types, Mannosidase Alpha Class 1A Member 2 was the common signature at genome level, and activation-ppara reaction, which stimulates a nuclear receptor protein, was found out as the mutual reporter at metabolic level. Moreover, miR-335 and miR-16-5p came into prominence in regulation of transcription at different tissues. On the other hand, distinct signatures were observed for different tissues at the metabolome level. Various coenzyme-A derivatives were significantly enriched metabolites in pancreatic islets, whereas skeletal muscle was enriched for cholesterol, malate, L-carnitine, and several amino acids. Results have showed utmost importance concerning relations between T2D and cancer, blood coagulation, neurodegenerative diseases, and specific metabolic and signaling pathways.

Proteomics Evidence for the Activity of the Putative Antibacterial Plant Alkaloid (-)-Roemerine: Mainstreaming Omics-Guided Drug Discovery

Discovery of new antibacterials with novel mechanisms is important to counteract the ingenious resistance mechanisms of bacteria. In this connection, omics-guided drug discovery offers a rigorous method in the quest of new antibacterials. (-)-Roemerine is a plant alkaloid that has been reported to possess putative antibacterial activity against Escherichia coli, Bacillus subtilis, and Salmonella typhimurium. The aim of the present study was to characterize the activity of (-)-roemerine in Escherichia coli TB1 using proteomics tools. With (-)-roemerine treatment, we found limited permeability through the outer membrane and repression of transport proteins involved in carbohydrate metabolism, resulting in poor carbon source availability. The shortfall of intracellular carbon sources in turn led to impaired cell growth. The reduction in the abundance of proteins related to translational machinery, amino acid biosynthesis, and metabolism was accompanied by a nutrient-limited state. The latter finding could suggest a metabolic shutdown in E. coli cells. High osmolarity was clearly not one of the reasons of bacterial death by (-)-roemerine. These observations collectively attest to the promise of plant omics and profiling of putative drug candidates using proteomics tools. Omics-guided drug discovery deserves greater attention in mainstream pharmacology so as to better understand the plants' medicinal potentials.

Zebrafish as a Model for Systems Medicine R&D: Rethinking the Metabolic Effects of Carrier Solvents and Culture Buffers Determined by (1)H NMR Metabolomics

Zebrafish is a frequently employed model organism in systems medicine and biomarker discovery. A crosscutting fundamental question, and one that has been overlooked in the field, is the "system-wide" (omics) effects induced in zebrafish by metabolic solvents and culture buffers. Indeed, any bioactivity or toxicity test requires that the target compounds are dissolved in an appropriate nonpolar solvent or aqueous media. It is important to know whether the solvent or the buffer itself has an effect on the zebrafish model organism. We evaluated the effects of two organic carrier solvents used in research with zebrafish, as well as in drug screening: dimethyl sulfoxide (DMSO) and ethanol, and two commonly used aqueous buffers (egg water and Hank's balanced salt solution). The effects of three concentrations (0.01, 0.1, and 1%) of DMSO and ethanol were tested in the 5-day-old zebrafish embryo using proton nuclear magnetic resonance ((1)H NMR) based metabolomics. DMSO (1% and 0.1%, but not 0.01%) exposure significantly decreased the levels of adenosine triphosphate (ATP), betaine, alanine, histidine, lactate, acetate, and creatine (p < 0.05). By contrast, ethanol exposure did not alter the embryos' metabolome at any concentration tested. The two different aqueous media noted above impacted the zebrafish embryo metabolome as evidenced by changes in valine, alanine, lactate, acetate, betaine, glycine, glutamate, adenosine triphosphate, and histidine. These results show that DMSO has greater effects on the embryo metabolome than ethanol, and thus is used with caution as a carrier solvent in zebrafish biomarker research and oral medicine. Moreover, the DMSO concentration should not be higher than 0.01%. Careful attention is also warranted for the use of the buffers egg water and Hank's balanced salt solution in zebrafish. In conclusion, as zebrafish is widely used as a model organism in life sciences, metabolome changes induced by solvents and culture buffers warrant further attention for robust systems science, and precision biomarkers that will stand the test of time.

Novel Vaccine Candidates against Brucella melitensis Identified through Reverse Vaccinology Approach

Global health therapeutics is a rapidly emerging facet of postgenomics medicine. In this connection, Brucella melitensis is an intracellular bacterium that causes the zoonotic infectious disease, brucellosis. Presently, no licensed vaccines are available for human brucellosis. Here, we report the identification of potential vaccine candidates against B. melitensis using a reverse vaccinology approach. Based on a systematic screening of exoproteome and secretome of B. melitensis 16 M, we identified eight proteins as potential vaccine candidates, including LPS-assembly protein LptD, a polysaccharide export protein, a cell surface protein, heme transporter BhuA, flagellin FliC, 7-alpha-hydroxysteroid dehydrogenase, immunoglobulin-binding protein EIBE, and hemagglutinin. Among these, the roles of BhuA and hemagglutinin in the virulence of Brucella are essential to establish infection. Roles of other proteins in the virulence are yet to be studied. Prediction of protein-protein interactions revealed that these proteins can interact with other proteins involved in virulence, secretion system, metabolism, and transport. From these eight potential vaccine candidates, we predicted three surface exposed novel antigenic epitopes that can induce both B-cell and T-cell immune responses. These peptides can be used for the development of either exclusive peptide vaccines or multi-component vaccines against human brucellosis. Reverse vaccinology is an important strategy for discovery of novel global health therapeutics.

Genome-Wide Analysis of Oleosin Gene Family in 22 Tree Species: An Accelerator for Metabolic Engineering of BioFuel Crops and Agrigenomics Industrial Applications?

Trees contribute to enormous plant oil reserves because many trees contain 50%-80% of oil (triacylglycerols, TAGs) in the fruits and kernels. TAGs accumulate in subcellular structures called oil bodies/droplets, in which TAGs are covered by low-molecular-mass hydrophobic proteins called oleosins (OLEs). The OLEs/TAGs ratio determines the size and shape of intracellular oil bodies. There is a lack of comprehensive sequence analysis and structural information of OLEs among diverse trees. The objectives of this study were to identify OLEs from 22 tree species (e.g., tung tree, tea-oil tree, castor bean), perform genome-wide analysis of OLEs, classify OLEs, identify conserved sequence motifs and amino acid residues, and predict secondary and three-dimensional structures in tree OLEs and OLE subfamilies. Data mining identified 65 OLEs with perfect conservation of the "proline knot" motif (PX5SPX3P) from 19 trees. These OLEs contained >40% hydrophobic amino acid residues. They displayed similar properties and amino acid composition. Genome-wide phylogenetic analysis and multiple sequence alignment demonstrated that these proteins could be classified into five OLE subfamilies. There were distinct patterns of sequence conservation among the OLE subfamilies and within individual tree species. Computational modeling indicated that OLEs were composed of at least three α-helixes connected with short coils without any β-strand and that they exhibited distinct 3D structures and ligand binding sites. These analyses provide fundamental information in the similarity and specificity of diverse OLE isoforms within the same subfamily and among the different species, which should facilitate studying the structure-function relationship and identify critical amino acid residues in OLEs for metabolic engineering of tree TAGs.

The Uyghur population and genetic susceptibility to type 2 diabetes: potential role for variants in CDKAL1, JAZF1, and IGF1 genes

Substantial evidence suggests that type 2 diabetes mellitus (T2DM) is a multi-factorial disease with a strong genetic component. A list of genetic susceptibility loci in populations of European and Asian ancestry has been established in the literature. Little is known on the inter-ethnic contribution of such established functional polymorphic variants. We performed a case-control study to explore the genetic susceptibility of 16 selected T2DM-related SNPs in a cohort of 102 Uyghur objects (51 cases and 51 controls). Three of the 16 SNPs showed significant association with T2DM in the Uyghur population. There were significant differences between the T2DM and control groups in frequencies of the risk allelic distributions of rs7754840 (CDKAL1) (p=0.014), rs864745 (JAZF1) (p=0.032), and rs35767 (IGF1) (p=0.044). Carriers of rs7754840-C, rs35767-A, and rs864745-C risk alleles had a 2.32-fold [OR (95% CI): 1.19-4.54], 2.06-fold [OR (95% CI): 1.02-4.17], 0.48-fold [OR (95% CI): 0.24-0.94] increased risk for T2DM, respectively. The cumulative risk allelic scores of these 16 SNPs differed significantly between the T2DM patients and the controls [17.1±8.1 vs. 15.4±7.3; OR (95%CI): 1.27(1.07-1.50), p=0.007]. This is the first study to evaluate genomic variation at 16 SNPs in respective T2DM candidate genes for the Uyghur population compared with other ethnic groups. The SNP rs7754840 in CDKAL1, rs864745 in JAZF1, and rs35767 in IGF1 might serve as potential susceptibility loci for T2DM in Uyghurs. We suggest a broader capture and study of the world populations, including who that are hitherto understudied, are essential for a comprehensive understanding of the genetic/genomic basis of T2DM.