Next Generation Sequencing Technologies: The Doorway to the Unexplored Genomics of Non-Model Plants

Identification and expression profiling of microRNAs in leaf tissues of Foeniculum vulgare Mill. under salinity stress

Foeniculum vulgare Mill. commonly known as fennel, is a globally recognized aromatic medicinal plant and culinary herb with widespread popularity due to its antimicrobial, antioxidant, carminative, and diuretic properties, among others. Although the phenotypic effects of salinity stress have been previously explored in fennel, the molecular mechanisms underlying responses to elevated salinity in this plant remain elusive. MicroRNAs (miRNAs) are tiny, endogenous, and extensively conserved non-coding RNAs (ncRNAs) typically ranging from 20 to 24 nucleotides (nt) in length that play a major role in a myriad of biological functions. In fact, a number of miRNAs have been extensively associated with responses to abiotic stress in plants. Consequently, employing computational methodologies and rigorous filtering criteria, 40 putative miRNAs belonging to 25 different families were characterized from fennel in this study. Subsequently, employing the psRNATarget tool, a total of 67 different candidate target transcripts for the characterized fennel miRNAs were predicted. Additionally, the expression patterns of six selected fennel miRNAs (i.e. fvu-miR156a, fvu-miR162a-3p, fvu-miR166a-3p, fvu-miR167a-5p, fvu-miR171a-3p, and fvu-miR408-3p) were analyzed under salinity stress conditions via qPCR. This article holds notable significance as it identifies not only 40 putative miRNAs in fennel, a non-model plant, but also pioneers the analysis of their expression under salinity stress conditions.

PABLOG: a Primer Analysis tool using a Bee-Like approach on Orthologous Genes

RNA-seq data is currently generated in numerous non-model organisms that lack a reference genome. Nevertheless, the confirmation of gene expression levels using RT-qPCR remains necessary, and the existing techniques do not seamlessly interface with the omics pipeline workflow. Developing primers for many targets by utilising orthologous genes can be a laborious, imprecise, and subjective process, particularly for plant species that are not commonly studied and do not have a known genome. We have developed a primer design tool, named PABLOG, that analyses the alignments generated from long or short RNA-seq reads and a reference orthologous gene. PABLOG scans, much like a bee searching several flowers for pollen, and presents a sorted list of potential exon-exon junction locations, ranked according to their reliability. Through computational analysis across the whole genomes of several non-model species, we demonstrate that PABLOG performs more effectively than other methods in identifying exon-exon junctions since it generates significantly fewer false-positive results. Examination of candidate regions at the gene level, in conjunction with laboratory studies, shows that the suggested primers successfully amplified particular targets in non-model plants without any presence of genomic contamination. Our tool includes a consensus sequence feature that enables the complete process of primer design, from aligning with the target gene to determining amplification parameters. The utility can be accessed via the GitHub repository located at: https://github.com/tools4plant-omics/PABLOG.

MCPtaggR: R package for accurate genotype calling in reduced representation sequencing data by eliminating error-prone markers based on genome comparison

Reduced representation sequencing (RRS) offers cost-effective, high-throughput genotyping platforms such as genotyping-by-sequencing (GBS). RRS reads are typically mapped onto a reference genome. However, mapping reads harbouring mismatches against the reference can potentially result in mismapping and biased mapping, leading to the detection of error-prone markers that provide incorrect genotype information. We established a genotype-calling pipeline named mappable collinear polymorphic tag genotyping (MCPtagg) to achieve accurate genotyping by eliminating error-prone markers. MCPtagg was designed for the RRS-based genotyping of a population derived from a biparental cross. The MCPtagg pipeline filters out error-prone markers prior to genotype calling based on marker collinearity information obtained by comparing the genome sequences of the parents of a population to be genotyped. A performance evaluation on real GBS data from a rice F2 population confirmed its effectiveness. Furthermore, our performance test using a genome assembly that was obtained by genome sequence polishing on an available genome assembly suggests that our pipeline performs well with converted genomes, rather than necessitating de novo assembly. This demonstrates its flexibility and scalability. The R package, MCPtaggR, was developed to provide functions for the pipeline and is available at https://github.com/tomoyukif/MCPtaggR.

Diosgenin biosynthesis investigation in medicinal herb (Tribulus terrestris) by transcriptome analysis

Next-generation sequencing (NGS) has revolutionized the analysis of specific genes, pathways, and their regulation in various species. Tribulus terrestris L., an annual medicinal herb of Zygophyllaceae family, has gained significant attention due to its diverse medicinal properties, including anti-inflammatory, antimicrobial, and anti-cancer effects. Diosgenin, a steroidal saponin, is the major bioactive compound responsible for the medicinal importance of T. terrestris. However, there is a paucity of information regarding the genes involved in the diosgenin biosynthetic pathway in T. terrestris. To address this gap, this study aimed to identify candidate genes associated with diosgenin biosynthesis through whole transcriptome profiling. A total of ∼7.9 GB of data, comprising 482 million reads, was obtained and assembled into 148,871 unigenes. Subsequently, functional annotations were assigned to 50 % of the unigenes using sequence similarity searches against the NCBI non-redundant (NR), Uniprot, KEGG, Pfam, GO, and COG databases, primarily based on Gene Ontology and KEGG-KAAS pathways. The majority of unigenes associated with the biosynthesis of the steroidal diosgenin backbone exhibited up-regulation in the fruit, leaf, and root tissues, except the SQE gene in root. The differential expression of selected genes was further validated through quantitative real-time polymerase chain reaction (qRT-PCR). Additionally, the study identified 21,026 unigenes related to transcription factors and 15,551 unigenes containing simple sequence repeats (SSR). Notably, di-nucleotide SSR motifs exhibited a high repeat frequency. These findings greatly enhance our understanding of the diosgenin biosynthesis pathway and provide a basis for future research in molecular investigation and metabolic engineering, specifically for boosting diosgenin content.

Transcriptomic and metabolomic profiling reveals the drought tolerance mechanism of Illicium difengpi (Schisandraceae)

Illicium difengpi (Schisandraceae), an endangered medicinal plant endemic to karst areas, is highly tolerant to drought and thus can be used as an ideal material for investigating adaptive mechanism to drought stress. The understanding of the drought tolerance of I. difengpi, especially at the molecular level, is lacking. In the present study, we aimed to clarify the molecular mechanism underlying drought tolerance in endemic I. difengpi plant in karst regions. The response characteristics of transcripts and changes in metabolite abundance of I. difengpi subjected to drought and rehydration were analyzed, the genes and key metabolites responsive to drought and rehydration were screened, and some important biosynthetic and secondary metabolic pathways were identified. A total of 231,784 genes and 632 metabolites were obtained from transcriptome and metabolome analyses, and most of the physiological metabolism in drought-treated I. difengpi plants recovered after rehydration. There were more upregulated genes than downregulated genes under drought and rehydration treatments, and rehydration treatment induced stable expression of 65.25% of genes, indicating that rehydration alleviated drought stress to some extent. Drought and rehydration treatment generated flavonoids, phenolic acids, flavonols, amino acids and their derivatives, as well as metabolites such as saccharides and alcohols in the leaves of I. difengpi plants, which alleviated the injury caused by excessive reactive oxygen species. The integration of transcriptome and metabolome analyses showed that, under drought stress, I. difengpi increased glutathione, flavonoids, polyamines, soluble sugars and amino acids, contributing to cell osmotic potential and antioxidant activity. The results show that the high drought tolerance and recovery after rehydration are the reasons for the normal growth of I. difengpi in karst mountain areas.

Recent Advances and Techniques for Identifying Novel Antibacterial Targets

BACKGROUND

With the emergence of drug-resistant bacteria, the development of new antibiotics is urgently required. Target-based drug discovery is the most frequently employed approach for the drug development process. However, traditional drug target identification techniques are costly and time-consuming. As research continues, innovative approaches for antibacterial target identification have been developed which enabled us to discover drug targets more easily and quickly.

METHODS

In this review, methods for finding drug targets from omics databases have been discussed in detail including principles, procedures, advantages, and potential limitations. The role of phage-driven and bacterial cytological profiling approaches is also discussed. Moreover, current article demonstrates the advancements being made in the establishment of computational tools, machine learning algorithms, and databases for antibacterial target identification.

RESULTS

Bacterial drug targets successfully identified by employing these aforementioned techniques are described as well.

CONCLUSION

The goal of this review is to attract the interest of synthetic chemists, biologists, and computational researchers to discuss and improve these methods for easier and quicker development of new drugs.

Genome-wide transcriptional profiling and physiological investigation elucidating the molecular mechanism of multiple abiotic stress response in Stevia rebaudiana Bertoni

Considering the major source of plant-derived low/non-calorie steviol glycosides (SGs), comprehensive physiological, biochemical, and deep transcriptional investigations were conducted to explicit deeper insight into multiple abiotic stress responses in Stevia rebaudiana. The physiological indicators including photosynthesis, chlorophyll, relative water content, shoot growth, electrolyte leakage, and SG biosynthesis were negatively impacted under drought (DS), followed by salinity (SS) and waterlogging (WS). Global transcriptional analysis revealed significant upregulated expression of the genes encoding for ROS detoxification (GST, SOD, APX, glutathione peroxidase), osmotic adjustment (alpha-trehalose-phosphate and S-adenosylmethionine decarboxylase), ion transporters (CAX, NHX, CNGS, VPPase, VATPase), water channel (PIP1, TIP) and abiotic stress-responsive candidate genes (LEA, HSPs, and Dehydrins) regulating abiotic stress response in S. rebaudiana. These inferences were complemented with predicted interactome network that revealed regulation of energy metabolism by key stress-responsive genes (GST, HKT1, MAPKs, P5CSs, PIP), transcription factors (HSFA2, DREB1A, DREB2A), and abiotic stress responsive pathways (ABA, ethylene, ion stress). This is the first detailed study to comprehend the molecular regulation of stress response and their interplay under DS, SS, and WS. The key genes and regulators can be functionally validated, and will facilitate targeted gene editing for genetic improvement of crop sustainability under changing environmental conditions in S. rebaudiana.

Recent advances in triterpenoid pathway elucidation and engineering

Triterpenoids are among the most assorted class of specialized metabolites found in all the taxa of living organisms. Triterpenoids are the leading active ingredients sourced from plant species and are utilized in pharmaceutical and cosmetic industries. The triterpenoid precursor 2,3-oxidosqualene, which is biosynthesized via the mevalonate (MVA) pathway is structurally diversified by the oxidosqualene cyclases (OSCs) and other scaffold-decorating enzymes such as cytochrome P450 monooxygenases (P450s), UDP-glycosyltransferases (UGTs) and acyltransferases (ATs). A majority of the bioactive triterpenoids are harvested from the native hosts using the traditional methods of extraction and occasionally semi-synthesized. These methods of supply are time-consuming and do not often align with sustainability goals. Recent advancements in metabolic engineering and synthetic biology have shown prospects for the green routes of triterpenoid pathway reconstruction in heterologous hosts such as Escherichia coli, Saccharomyces cerevisiae and Nicotiana benthamiana, which appear to be quite promising and might lead to the development of alternative source of triterpenoids. The present review describes the biotechnological strategies used to elucidate complex biosynthetic pathways and to understand their regulation and also discusses how the advances in triterpenoid pathway engineering might aid in the scale-up of triterpenoid production in engineered hosts.

Genetic Enhancement of Cereals Using Genomic Resources for Nutritional Food Security

Advances in genomics resources have facilitated the evolution of cereal crops with enhanced yield, improved nutritional values, and heightened resistance to various biotic and abiotic stresses. Genomic approaches present a promising avenue for the development of high-yielding varieties, thereby ensuring food and nutritional security. Significant improvements have been made within the omics domain, specifically in genomics, transcriptomics, and proteomics. The advent of Next-Generation Sequencing (NGS) techniques has yielded an immense volume of data, accompanied by substantial progress in bioinformatic tools for proficient analysis. The synergy between genomics and computational tools has been acknowledged as pivotal for unravelling the intricate mechanisms governing genome-wide gene regulation. Within this review, the essential genomic resources are delineated, and their harmonization in the enhancement of cereal crop varieties is expounded upon, with a paramount focus on fulfilling the nutritional requisites of humankind. Furthermore, an encompassing compendium of the available genomic resources for cereal crops is presented, accompanied by an elucidation of their judicious utilization in the advancement of crop attributes.

Development of EST-SSRs based on the transcriptome of Castanopsis carlesii and cross-species transferability in other Castanopsis species

Castanopsis carlesii (Hemsl.) Hay. is a widely distributed and dominant tree species native to subtropical China with significant ecological and economic value. Due to serious human-related disturbance, its wild resources have been increasingly reduced, and whether may result in the loss of genetic diversity. However, no population genetics studies of natural C. carlesii have been reported to date. Microsatellite markers have been a useful tool in population genetics. Therefore, we developed EST-SSR markers based on the transcriptome sequencing of C. carlesii leaves. A total of 149,380,224 clean reads were obtained, and 63,012 nonredundant unigenes with a mean length of 1,034 bp were assembled and annotated based on sequence similarity searches in the Nr, Nt, KO, SwissProt, PFAM, KOG, and GO databases. The results showed that only 5,559 (8.82%) unigenes were annotated in all seven databases, but 46,338 (73.53%) could be annotated in at least one database. A total of 31,459 potential EST-SSRs were identified in 18,690 unigenes, with an average frequency of one SSR approximately 2 kb. Among the 100 EST-SSR primer pairs designed, 49 primer pairs successfully produced the expected product by amplification, with a success rate of 49%, but only 20 primer pairs showed abundant polymorphisms. Polymorphisms were verified using 25 samples from C. carlesii in Qimen, Anhui. A total of 119 alleles were detected, with a mean number of alleles (Na) of 5.95 per locus and a mean polymorphism information content (PIC) of 0.6125. All the 20 newly developed EST-SSR markers were verified in other Castanopsis species (C. sclerophylla, C. lamontii, C. fargesii, C. eyrei and C. jucunda). Sixteen primer pairs showed successful amplification in all five Castanopsis species (80%), and the transferability ratios ranged from 90% to 100%. These developed EST-SSR markers can be applied to population genetic and germplasm evaluations of C. carlesii and related species.

Diversity analyses in two ornamental and large-genome Ranunculaceae species based on a low-cost Klenow NGS-based protocol

Persian buttercup (Ranunculus asiaticus L.) and poppy anemone (Anemone coronaria L.) are ornamental, outcrossing, perennial species belonging to the Ranunculaceae family, characterized by large and highly repetitive genomes. We applied K-seq protocol in both species to generate high-throughput sequencing data and produce a large number of genetic polymorphisms. The technique entails the application of Klenow polymerase-based PCR using short primers designed by analyzing k-mer sets in the genome sequence. To date the genome sequence of both species has not been released, thus we designed primer sets based on the reference the genome sequence of the related species Aquilegia oxysepala var. kansuensis (Brühl). A whole of 11,542 SNPs were selected for assessing genetic diversity of eighteen commercial varieties of R. asiaticus, while 1,752 SNPs for assessing genetic diversity in six cultivars of A. coronaria. UPGMA dendrograms were constructed and in R. asiaticus integrated in with PCA analysis. This study reports the first molecular fingerprinting within Persian buttercup, while the results obtained in poppy anemone were compared with a previously published SSR-based fingerprinting, proving K-seq to be an efficient protocol for the genotyping of complex genetic backgrounds.

A review on efforts for improvement in medicinally important chemical constituents inAconitum through biotechnological interventions

The genus Aconitum belongs to the family Ranunculaceae, is endowed with more than 350 species on the earth. Medicinally important aconitine type of diterpenoid alkaloids are the characteristic compounds in most of the Aconitum species. The present review endeavored the major research carried out in the field of genetic resource characterization, pharmacological properties, phytochemistry, major factors influencing quantity, biosynthetic pathways and processing methods for recovery of active ingredients, variety improvement, propagation methods, and important metabolite production through cell/organ culture of various Aconitum species. More than 450 derivatives of aconitine-type C19 and C20-diterpenoid alkaloids along with a few other non-alkaloidal compounds, such as phenylpropanoids, flavonoids, terpenoids, and fatty acids, have been identified in the genus. A few Aconitum species and their common diterpenoid alkaloid compounds are also well characterized for analgesic, inflammatory and cytotoxic properties. However, the different isolated compound needs to be validated for supporting other traditional therapeutical uses of the plant species. Aconitine alkaloids shared common biosynthesis pathway, but their diversification mechanism remains unexplored in the genus. Furthermore, the process needs to be developed on secondary metabolite recovery, mass-scale propagation methods, and agro-technologies for maintaining the quality of products. Many species are losing their existence in nature due to over-exploitation or anthropogenic factors; thus, temporal monitoring of the population status in its habitat, and suitable management programs for ascertaining conservation needs to be developed.

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.

Transcriptomic and functional analyses reveal the molecular mechanisms underlying Fe-mediated tobacco resistance to potato virus Y infection

Potato virus Y (PVY) mainly infects Solanaceous crops, resulting in considerable losses in the yield and quality. Iron (Fe) is involved in various biological processes in plants, but its roles in resistance to PVY infection has not been reported. In this study, foliar application of Fe could effectively inhibit early infection of PVY, and a full-length transcriptome and Illumina RNA sequencing was performed to investigate its modes of action in PVY-infected Nicotiana tabacum. The results showed that 18,074 alternative splicing variants, 3,654 fusion transcripts, 3,086 long non-coding RNAs and 14,403 differentially expressed genes (DEGs) were identified. Specifically, Fe application down-regulated the expression levels of the DEGs related to phospholipid hydrolysis, phospholipid signal, cell wall biosynthesis, transcription factors (TFs) and photosystem I composition, while those involved with photosynthetic electron transport chain (PETC) were up-regulated at 1 day post inoculation (dpi). At 3 dpi, these DEGs related to photosystem II composition, PETC, molecular chaperones, protein degradation and some TFs were up-regulated, while those associated with light-harvesting, phospholipid hydrolysis, cell wall biosynthesis were down-regulated. At 9 dpi, Fe application had little effects on resistance to PVY infection and transcript profiles. Functional analysis of these potentially critical DEGs was thereafter performed using virus-induced gene silencing approaches and the results showed that NbCat-6A positively regulates PVY infection, while the reduced expressions of NbWRKY26, NbnsLTP, NbFAD3 and NbHSP90 significantly promote PVY infection in N. benthamiana. Our results elucidated the regulatory network of Fe-mediated resistance to PVY infection in plants, and the functional candidate genes also provide important theoretical bases to further improve host resistance against PVY infection.

Data Science and Plant Metabolomics

The study of plant metabolism is one of the most complex tasks, mainly due to the huge amount and structural diversity of metabolites, as well as the fact that they react to changes in the environment and ultimately influence each other. Metabolic profiling is most often carried out using tools that include mass spectrometry (MS), which is one of the most powerful analytical methods. All this means that even when analyzing a single sample, we can obtain thousands of data. Data science has the potential to revolutionize our understanding of plant metabolism. This review demonstrates that machine learning, network analysis, and statistical modeling are some techniques being used to analyze large quantities of complex data that provide insights into plant development, growth, and how they interact with their environment. These findings could be key to improving crop yields, developing new forms of plant biotechnology, and understanding the relationship between plants and microbes. It is also necessary to consider the constraints that come with data science such as quality and availability of data, model complexity, and the need for deep knowledge of the subject in order to achieve reliable outcomes.

Molecular Mechanisms of Plant Regeneration from Differentiated Cells: Approaches from Historical Tissue Culture Systems

Plants can exert remarkable capacity for cell reprogramming even from differentiated cells. This ability allows plants to regenerate tissues/organs and even individuals in nature and in vitro. In recent decades, Arabidopsis research has uncovered molecular mechanisms of plant regeneration; however, our understanding of how plant cells retain both differentiated status and developmental plasticity is still obscure. In this review, we first provide a brief outlook of the representative modes of plant regeneration and key factors revealed by Arabidopsis research. We then re-examine historical tissue culture systems that enable us to investigate the molecular details of cell reprogramming in differentiated cells and discuss the different approaches, specifically highlighting our recent progress in shoot regeneration from the epidermal cell of Torenia fournieri.

Weak genetic differentiation but strong climate-induced selective pressure toward the rear edge of mountain pine in north-eastern Spain

Local differentiation at distribution limits may influence species' adaptive capacity to environmental changes. However, drivers, such gene flow and local selection, are still poorly understood. We focus on the role played by range limits in mountain forests to test the hypothesis that relict tree populations are subjected to genetic differentiation and local adaptation. Two alpine treelines of mountain pine (Pinus uncinata Ram. ex DC) were investigated in the Spanish Pyrenees. Further, an isolated relict population forming the species' southernmost distribution limit in north-eastern Spain was also investigated. Using genotyping by sequencing, a genetic matrix conformed by single nucleotide polymorphisms (SNPs) was obtained. This matrix was used to perform genotype-environment and genotype-phenotype associations, as well as to model risk of non-adaptedness. Increasing climate seasonality appears as an essential element in the interpretation of SNPs subjected to selective pressures. Genetic differentiations were overall weak. The differences in leaf mass area and radial growth rate, as well as the identification of several SNPs subjected to selective pressures, exceeded neutral predictions of differentiation among populations. Despite genetic drift might prevail in the isolated population, the Fst values (0.060 and 0.066) showed a moderate genetic drift and Nm values (3.939 and 3.555) indicate the presence of gene flow between the relict population and both treelines. Nonetheless, the SNPs subjected to selection pressures provide evidences of possible selection in treeline ecotones. Persistence in range boundaries seems to involve several selective pressures in species' traits, which were significantly related to enhanced drought seasonality at the limit of P. uncinata distribution range. We conclude that gene flow is unlikely to constrain adaptation in the P. uncinata rear edge, although this species shows vulnerability to future climate change scenarios involving warmer and drier conditions.

Integrated omic techniques and their genomic features for invasive weeds

Many emerging invasive weeds display rapid adaptation against different stressful environments compared to their natives. Rapid adaptation and dispersal habits helped invasive populations have strong diversity within the population compared to their natives. Advances in molecular marker techniques may lead to an in-depth understanding of the genetic diversity of invasive weeds. The use of molecular techniques is rapidly growing, and their implications in invasive weed studies are considered powerful tools for genome purposes. Here, we review different approach used multi-omics by invasive weed studies to understand the functional structural and genomic changes in these species under different environmental fluctuations, particularly, to check the accessibility of advance-sequencing techniques used by researchers in genome sequence projects. In this review-based study, we also examine the importance and efficiency of different molecular techniques in identifying and characterizing different genes, associated markers, proteins, metabolites, and key metabolic pathways in invasive and native weeds. Use of these techniques could help weed scientists to further reduce the knowledge gaps in understanding invasive weeds traits. Although these techniques can provide robust insights about the molecular functioning, employing a single omics platform can rarely elucidate the gene-level regulation and the associated real-time expression of weedy traits due to the complex and overlapping nature of biological interactions. We conclude that different multi-omic techniques will provide long-term benefits in launching new genome projects to enhance the understanding of invasive weeds' invasion process.

Development and Complex Application of Methods for the Identification of Mutations in the FAD3A and FAD3B Genes Resulting in the Reduced Content of Linolenic Acid in Flax Oil

Flax is grown worldwide for seed and fiber production. Linseed varieties differ in their oil composition and are used in pharmaceutical, food, feed, and industrial production. The field of application primarily depends on the content of linolenic (LIN) and linoleic (LIO) fatty acids. Inactivating mutations in the FAD3A and FAD3B genes lead to a decrease in the LIN content and an increase in the LIO content. For the identification of the three most common low-LIN mutations in flax varieties (G-to-A in exon 1 of FAD3A substituting tryptophan with a stop codon, C-to-T in exon 5 of FAD3A leading to arginine to a stop codon substitution, and C-to-T in exon 2 of FAD3B resulting in histidine to tyrosine substitution), three approaches were proposed: (1) targeted deep sequencing, (2) high resolution melting (HRM) analysis, (3) cleaved amplified polymorphic sequences (CAPS) markers. They were tested on more than a thousand flax samples of various types and showed promising results. The proposed approaches can be used in marker-assisted selection to choose parent pairs for crosses, separate heterogeneous varieties into biotypes, and select genotypes with desired homozygous alleles of the FAD3A and FAD3B genes at the early stages of breeding for the effective development of varieties with a particular LIN and LIO content, as well as in basic studies of the molecular mechanisms of fatty acid synthesis in flax seeds to select genotypes adequate to the tasks.

Genetic Diversity and Population Structure of Chinese Chestnut (Castanea mollissima Blume) Cultivars Revealed by GBS Resequencing

Chinese chestnut (Castanea mollissima Bl.) is one of the earliest domesticated and cultivated fruit trees, and it is widely distributed in China. Because of the high quality of its nuts and its high resistance to abiotic and biotic stresses, Chinese chestnut could be used to improve edible chestnut varieties worldwide. However, the unclear domestication history and highly complex genetic background of Chinese chestnut have prevented the efficiency of breeding efforts. To explore the genetic diversity and structure of Chinese chestnut populations and generate new insights that could aid chestnut breeding, heterozygosity statistics, molecular variance analysis, ADMIXTURE analysis, principal component analysis, and phylogenetic analysis were conducted to analyze single nucleotide polymorphism data from 185 Chinese chestnut landraces from five geographical regions in China via genotyping by sequencing. Results showed that the genetic diversity level of the five populations from different regions was relatively high, with an observed heterozygosity of 0.2796-0.3427. The genetic diversity level of the population in the mid-western regions was the highest, while the population north of the Yellow River was the lowest. Molecular variance analysis showed that the variation among different populations was only 2.07%, while the intra-group variation reached 97.93%. The Chinese chestnut samples could be divided into two groups: a northern and southern population, separated by the Yellow River; however, some samples from the southern population were genetically closer to samples from the northern population. We speculate that this might be related to the migration of humans during the Han dynasty due to the frequent wars that took place during this period, which might have led to the introduction of chestnut to southern regions. Some samples from Shandong Province and Beijing City were outliers that did not cluster with their respective groups, and this might be caused by the special geographical, political, and economic significance of these two regions. The findings of our study showed the complex genetic relationships among Chinese chestnut landraces and the high genetic diversity of these resources.

Genome Survey Sequencing and Genetic Background Characterization of Ilex chinensis Sims (Aquifoliaceae) Based on Next-Generation Sequencing

Ilex chinensis Sims. is an evergreen arbor species with high ornamental and medicinal value that is widely distributed in China. However, there is a lack of molecular and genomic data for this plant, which severely restricts the development of its relevant research. To obtain the whole reference genome, we first conducted a genome survey of I. chinensis by next-generation sequencing (NGS) to perform de novo whole-genome sequencing. As a result, our estimates using k-mer and flow cytometric analysis suggested the genome size of I. chinensis to be around 618-655 Mb, with the GC content, heterozygous rate, and repeat sequence rate of 37.52%, 1.1%, and 38%, respectively. A total of 334,649 microsatellite motifs were detected from the I. chinensis genome data, which will provide basic molecular markers for germplasm characterization, genetic diversity, and QTL mapping studies for I. chinensis. In summary, the I. chinensis genome is complex with high heterozygosity and few repeated sequences. Overall, this is the first report on the genome features of I. chinensis, and the information may lay a strong groundwork for future whole-genome sequencing and molecular breeding studies of this species.

Next-generation sequencing in the biodiversity conservation of endangered medicinal plants

Medicinal plants have been used as traditional herbal medicines in the treatment of various types of diseases. However, the increased demand for these plants highlights the importance of conservation specifically for endangered species. Significant advancements in next-generation sequencing (NGS) technologies have accelerated medicinal plant research while reducing costs and time demands. NGS systems enable high-throughput whole genome sequencing as well as direct RNA sequencing and transcriptome analysis. The sequence data sets created can be used in a variety of areas of study, including biodiversity conservation, comparative genomics, transcriptomic analysis, single cell mining, metagenomics, epigenetics, molecular marker discovery, multi genome sequencing, and so on. Commercial sequencing service providers are constantly working to improve technologies to address bioinformatics problems in NGS data analysis. Several genome sequencing projects on medicinal plants have been completed recently and a few more are in the works. In some medicinal plants, massive NGS-based data has been developed. In the present review, we have attempted to briefly discuss advancements in NGS technology on medicinally essential plants in India. The review will also provide ideas for applying NGS technologies for exploring genomes of various endangered medicinal plants whose genome sequences are not normally available and thus provides valuable insights for the conservation of these vulnerable species.

De novo transcriptome assembly, gene annotation, and EST-SSR marker development of an important medicinal and edible crop, Amomum tsaoko (Zingiberaceae)

BACKGROUND

Amomum tsaoko is a medicinal and food dual-use crop that belongs to the Zingiberaceae family. However, the lack of transcriptomic and genomic information has limited the understanding of the genetic basis of this species. Here, we performed transcriptome sequencing of samples from different A. tsaoko tissues, and identified and characterized the expressed sequence tag-simple sequence repeat (EST-SSR) markers.

RESULTS

A total of 58,278,226 high-quality clean reads were obtained and de novo assembled to generate 146,911 unigenes with an N50 length of 2002 bp. A total of 128,174 unigenes were successfully annotated by searching seven protein databases, and 496 unigenes were identified as annotated as putative terpenoid biosynthesis-related genes. Furthermore, a total of 55,590 EST-SSR loci were detected, and 42,333 primer pairs were successfully designed. We randomly selected 80 primer pairs to validate their polymorphism in A. tsaoko; 18 of these primer pairs produced distinct, clear, and reproducible polymorphisms. A total of 98 bands and 96 polymorphic bands were amplified by 18 pairs of EST-SSR primers for the 72 A. tsaoko accessions. The Shannon's information index (I) ranged from 0.477 (AM208) to 1.701 (AM242) with an average of 1.183, and the polymorphism information content (PIC) ranged from 0.223 (AM208) to 0.779 (AM247) with an average of 0.580, indicating that these markers had a high level of polymorphism. Analysis of molecular variance (AMOVA) indicated relatively low genetic differentiation among the six A. tsaoko populations. Cross-species amplification showed that 14 of the 18 EST-SSR primer pairs have transferability between 11 Zingiberaceae species.

CONCLUSIONS

Our study is the first to provide transcriptome data of this important medicinal and edible crop, and these newly developed EST-SSR markers are a very efficient tool for germplasm evaluation, genetic diversity, and molecular marker-assisted selection in A. tsaoko.

Development of Genomic Resources in Mexican Bursera (Section: Bullockia: Burseraceae): Genome Assembly, Annotation, and Marker Discovery for Three Copal Species

Bursera comprises ~100 tropical shrub and tree species, with the center of the species diversification in Mexico. The genomic resources developed for the genus are scarce, and this has limited the study of the gene flow, local adaptation, and hybridization dynamics. In this study, based on ~155 million Illumina paired-end reads per species, we performed a de novo genome assembly and annotation of three Bursera species of the Bullockia section: Bursera bipinnata, Bursera cuneata, and Bursera palmeri. The total lengths of the genome assemblies were 253, 237, and 229 Mb for B. cuneata, B. palmeri, and B. bipinnata, respectively. The assembly of B. palmeri retrieved the most complete and single-copy BUSCOs (87.3%) relative to B. cuneata (86.5%) and B. bipinnata (76.6%). The ab initio gene prediction recognized between 21,000 and 32,000 protein-coding genes. Other genomic features, such as simple sequence repeats (SSRs), were also detected. Using the de novo genome assemblies as a reference, we identified single-nucleotide polymorphisms (SNPs) for a set of 43 Bursera individuals. Moreover, we mapped the filtered reads of each Bursera species against the chloroplast genomes of five Burseraceae species, obtaining consensus sequences ranging from 156 to 160 kb in length. Our work contributes to the generation of genomic resources for an important but understudied genus of tropical-dry-forest species.

An Integrative Transcriptional Network Revealed Spatial Molecular Interplay Underlying Alantolactone and Inulin Biosynthesis in Inula racemosa Hook f

Inula racemosa Hook. f. (Pushkarmula), a perennial Himalayan herb known for its aromatic and phytopharmaceutical attributes, is not yet explored at genomic/transcriptomic scale. In this study, efforts were made to unveil the global transcriptional atlas underlying organ-specific specialized metabolite biosynthesis by integrating RNA-Seq analysis of 433 million sequenced reads with the phytochemical analysis of leaf, stem, and root tissues. Overall, 7242 of 83,772 assembled nonredundant unigenes were identified exhibiting spatial expression in leaf (3761), root (2748), and stem (733). Subsequently, integration of the predicted transcriptional interactome network of 2541 unigenes (71,841 edges) with gene ontology and KEGG pathway enrichment analysis revealed isoprenoid, terpenoid, diterpenoid, and gibberellin biosynthesis with antimicrobial activities in root tissue. Interestingly, the root-specific expression of germacrene-mediated alantolactone biosynthesis (GAS, GAO, G8H, IPP, DMAP, and KAO) and antimicrobial activities (BZR1, DEFL, LTP) well-supported with both quantitative expression profiling and phytochemical accumulation of alantolactones (726.08 μg/10 mg) and isoalantolactones (988.59 μg/10 mg), which suggests “roots” as the site of alantolactone biosynthesis. A significant interaction of leaf-specific carbohydrate metabolism with root-specific inulin biosynthesis indicates source (leaf) to sink (root) regulation of inulin. Our findings comprehensively demonstrate the source-sink transcriptional regulation of alantolactone and inulin biosynthesis, which can be further extended for upscaling the targeted specialized metabolites. Nevertheless, the genomic resource created in this study can also be utilized for development of genome-wide functionally relevant molecular markers to expedite the breeding strategies for genetic improvement of I. racemosa.

Spatial Genomic Resource Reveals Molecular Insights into Key Bioactive-Metabolite Biosynthesis in Endangered Angelica glauca Edgew

Angelica glauca Edgew, which is an endangered medicinal and aromatic herb, is a rich source of numerous industrially important bioactive metabolites, including terpenoids, phenolics, and phthalides. Nevertheless, genomic interventions for the sustainable utilization and restoration of its genetic resources are greatly offset due to the scarcity of the genomic resources and key regulators of the underlying specialized metabolism. To unravel the global atlas of the specialized metabolism, the first spatial transcriptome sequencing of the leaf, stem, and root generated 109 million high-quality paired-end reads, assembled de novo into 81,162 unigenes, which exhibit a 61.53% significant homology with the six public protein databases. The organ-specific clustering grouped 1136 differentially expressed unigenes into four subclusters differentially enriched in the leaf, stem, and root tissues. The prediction of the transcriptional-interactome network by integrating enriched gene ontology (GO) and the KEGG metabolic pathways identified the key regulatory unigenes that correspond to terpenoid, flavonoid, and carotenoid biosynthesis in the leaf tissue, followed by the stem and root tissues. Furthermore, the stem and root-specific significant enrichments of phenylalanine ammonia lyase (PAL), cinnamate-4-hydroxylase (C4H), and caffeic acid 3-O-methyltransferase (COMT) indicate that phenylalanine mediated the ferulic acid biosynthesis in the stem and root. However, the root-specific expressions of NADPH-dependent alkenal/one oxidoreductase (NADPH-AOR), S-adenosyl-L-methionine-dependent methyltransferases (SDMs), polyketide cyclase (PKC), and CYP72A15 suggest the “root” as the primary site of phthalide biosynthesis. Additionally, the GC-MS and UPLC analyses corresponded to the organ-specific gene expressions, with higher contents of limonene and phthalide compounds in the roots, while there was a higher accumulation of ferulic acid in the stem, followed by in the root and leaf tissues. The first comprehensive genomic resource with an array of candidate genes of the key metabolic pathways can be potentially utilized for the targeted upscaling of aromatic and pharmaceutically important bioactive metabolites. This will also expedite genomic-assisted conservation and breeding strategies for the revival of the endangered A. glauca.

Metagenomics: An Approach for Unraveling the Community Structure and Functional Potential of Activated Sludge of a Common Effluent Treatment Plant

The common effluent treatment plant (CETP) located at Baddi treats the industrial effluent from various industries, leading to the pooling of a diverse range of substrates and metabolites. The nutrient loading and its availability decide the balance of the microbial community and its diversity. The samples thus collected from the activated sludge (BS14) of CETP and Sirsa river (SR1) from the vicinity of CETP effluent discharge were processed for the whole metagenome analysis to reveal the microbial community and its functional potential. The taxonomic classification of the BS14 sample showed the dominance of the bacterial community with 96% of abundance, whereas the SR1 was populated by eukaryotes representing 50.4% of the community of SR1. The bacterial community of SR1 was constituted of 47.2%. The functional analysis of BS14 and SR1 with GhostKOALA against the KEGG database assigned 43.7% and 27.8% of the open reading frames (ORFs) with functions. It revealed the xenobiotic degradation modules with complete pathways along with resistance against the beta-lactams. The analysis with the comprehensive antibiotic resistance database (CARD) revealed 33 and 32 unique types of antimicrobial resistance in BS14 and SR1, respectively. Both the samples were dominated by the beta-lactam resistance genes. The carbohydrate-active enzyme (CAZy) database assigned a total of 6,611 and 2,941 active enzymes to BS14 and SR1, respectively. In contrast, the glycosyl hydrolases (GH) and glycosyltransferases (GT) class of enzymes were found to be abundant in both the samples as compared with polysaccharide lyases (PL), auxiliary activities (AA), carbohydrate esterases (CE), and carbohydrate-binding module (CBM).

ddRAD Sequencing-Based Scanning of Genetic Variants in Sargassum fusiforme

Sargassum fusiforme is a commercially important brown seaweed that has experienced significant population reduction both from heavy exploitation and degradation of the environment. Cultivated breed strains are also in a state of population mixing. These population stressors make it necessary to investigate the population genetics to discover best practices to conserve and breed this seaweed. In this study, the genetic diversity and population structure of S. fusiforme were investigated by the genome-wide SNP data acquired from double digest restriction site-associated DNA sequencing (ddRAD-seq). We found a low genetic diversity and a slight population differentiation within and between wild and cultivated populations, and the effective population size of S. fusiforme had experienced a continuous decline. Tajima’s D analysis showed the population contraction in wild populations may be related to copper pollution which showed a consistent trend with the increase of the sea surface temperature. The potential selection signatures may change the timing or level of gene expression, and further experiments are needed to investigate the effect of the mutation on relevant pathways. These results suggest an urgent need to manage and conserve S. fusiforme resources and biodiversity considering the accelerating change of the environment.

Plant Flavonoid Production in Bacteria and Yeasts

Flavonoids, a major group of secondary metabolites in plants, are promising for use as pharmaceuticals and food supplements due to their health-promoting biological activities. Industrial flavonoid production primarily depends on isolation from plants or organic synthesis, but neither is a cost-effective or sustainable process. In contrast, recombinant microorganisms have significant potential for the cost-effective, sustainable, environmentally friendly, and selective industrial production of flavonoids, making this an attractive alternative to plant-based production or chemical synthesis. Structurally and functionally diverse flavonoids are derived from flavanones such as naringenin, pinocembrin and eriodictyol, the major basic skeletons for flavonoids, by various modifications. The establishment of flavanone-producing microorganisms can therefore be used as a platform for producing various flavonoids. This review summarizes metabolic engineering and synthetic biology strategies for the microbial production of flavanones. In addition, we describe directed evolution strategies based on recently-developed high-throughput screening technologies for the further improvement of flavanone production. We also describe recent progress in the microbial production of structurally and functionally complicated flavonoids via the flavanone modifications. Strategies based on synthetic biology will aid more sophisticated and controlled microbial production of various flavonoids.

RNA-Seq Analysis of Developing Grains of Wheat to Intrigue Into the Complex Molecular Mechanism of the Heat Stress Response

Heat stress is one of the significant constraints affecting wheat production worldwide. To ensure food security for ever-increasing world population, improving wheat for heat stress tolerance is needed in the presently drifting climatic conditions. At the molecular level, heat stress tolerance in wheat is governed by a complex interplay of various heat stress-associated genes. We used a comparative transcriptome sequencing approach to study the effect of heat stress (5°C above ambient threshold temperature of 20°C) during grain filling stages in wheat genotype K7903 (Halna). At 7 DPA (days post-anthesis), heat stress treatment was given at four stages: 0, 24, 48, and 120 h. In total, 115,656 wheat genes were identified, including 309 differentially expressed genes (DEGs) involved in many critical processes, such as signal transduction, starch synthetic pathway, antioxidant pathway, and heat stress-responsive conserved and uncharacterized putative genes that play an essential role in maintaining the grain filling rate at the high temperature. A total of 98,412 Simple Sequences Repeats (SSR) were identified from de novo transcriptome assembly of wheat and validated. The miRNA target prediction from differential expressed genes was performed by psRNATarget server against 119 mature miRNA. Further, 107,107 variants including 80,936 Single nucleotide polymorphism (SNPs) and 26,171 insertion/deletion (Indels) were also identified in de novo transcriptome assembly of wheat and wheat genome Ensembl version 31. The present study enriches our understanding of known heat response mechanisms during the grain filling stage supported by discovery of novel transcripts, microsatellite markers, putative miRNA targets, and genetic variant. This enhances gene functions and regulators, paving the way for improved heat tolerance in wheat varieties, making them more suitable for production in the current climate change scenario.

Plant GATA Factors: Their Biology, Phylogeny, and Phylogenomics

GATA factors are evolutionarily conserved transcription factors that are found in animals, fungi, and plants. Compared to that of animals, the size of the plant GATA family is increased. In angiosperms, four main GATA classes and seven structural subfamilies can be defined. In recent years, knowledge about the biological role and regulation of plant GATAs has substantially improved. Individual family members have been implicated in the regulation of photomorphogenic growth, chlorophyll biosynthesis, chloroplast development, photosynthesis, and stomata formation, as well as root, leaf, and flower development. In this review, we summarize the current knowledge of plant GATA factors. Using phylogenomic analysis, we trace the evolutionary origin of the GATA classes in the green lineage and examine their relationship to animal and fungal GATAs. Finally, we speculate about a possible conservation of GATA-regulated functions across the animal, fungal, and plant kingdoms.

The phylogeny of Salix revealed by whole genome re-sequencing suggests different sex-determination systems in major groups of the genus

BACKGROUND AND AIMS

The largest genus of Salicaceae sensu lato, Salix, has been shown to consist of two main clades: clade Salix, in which species have XY sex-determination systems (SDSs) on chromosome 7, and clade Vetrix including species with ZW SDSs on chromosome 15. Here, we test the utility of whole genome re-sequencing (WGR) for phylogenomic reconstructions of willows to infer changes between different SDSs.

METHODS

We used more than 1 TB of WGR data from 70 Salix taxa to ascertain single nucleotide polymorphisms on the autosomes, the sex-linked regions (SLRs) and the chloroplast genomes, for phylogenetic and species tree analyses. To avoid bias, we chose reference genomes from both groups, Salix dunnii from clade Salix and S. purpurea from clade Vetrix.

KEY RESULTS

Two main largely congruent groups were recovered: the paraphyletic Salix grade and the Vetrix clade. The autosome dataset trees resolved four subclades (C1-C4) in Vetrix. C1 and C2 comprise species from the Hengduan Mountains and adjacent areas and from Eurasia, respectively. Section Longifoliae (C3) grouped within the Vetrix clade but fell into the Salix clade in trees based on the chloroplast dataset analysis. Salix triandra from Eurasia (C4) was revealed as sister to the remaining species of clade Vetrix. In Salix, the polyploid group C5 is paraphyletic to clade Vetrix and subclade C6 is consistent with Argus's subgenus Protitea. Chloroplast datasets separated both Vetrix and Salix as monophyletic, and yielded C5 embedded within Salix. Using only diploid species, both the SLR and autosomal datasets yielded trees with Vetrix and Salix as well-supported clades.

CONCLUSION

WGR data are useful for phylogenomic analyses of willows. The different SDSs may contribute to the isolation of the two major groups, but the reproductive barrier between them needs to be studied.

Novel Insights Into Refugia at the Southern Margin of the Distribution Range of the Endangered Species Ulmus laevis

Riparian ecosystems, in long-time developed regions, are among the most heavily impacted by human activities; therefore, the distribution of tree riparian species, such as Ulmus laevis, is highly affected. This phenomenon is particularly relevant at the margins of the natural habitat of the species, where populations are small and rare. In these cases, it is difficult to distinguish between relics or introductions, but it is relevant for the restoration of natural habitats and conservation strategies. The aim of this study was to study the phylogeography of the southern distribution of the species. We sequenced the entire chloroplast (cp) genomes of 54 individuals from five sampled populations across different European regions to highlight polymorphisms and analyze their distribution. Thirty-two haplotypes were identified. All the sampled populations showed private haplotypes that can be considered an indicator of long-term residency, given the low mutation rate of organellar DNA. The network of all haplotypes showed a star-like topology, and Serbian haplotypes were present in all branches. The Balkan population showed the highest level of nucleotide and genetic diversity. Low genetic differentiation between populations was observed but we found a significant differentiation among Serbia vs. other provenances. Our estimates of divergent time of U. laevis samples highlight the early split of above all Serbian individuals from other populations, emphasizing the reservoir role of white elm genetic diversity of Serbian population.

Distribution of simple sequence repeats, transcription factors, and differentially expressed genes in the NGS-based transcriptome of male and female seabuckthorn (Hippophae salicifolia)

Seabuckthorn (Hippophae salicifolia) is a perennial, multipurpose wonder plant, popular for its immense medicinal, nutritional, and therapeutic properties. However, due to the lack of whole-genome-based studies, the molecular mechanism governing distinct sexual phenotypes is still not clear. We employed the high-throughput NGS Illumina NovaSeq paired-end technology to generate whole transcriptome profiles of male and female plants of H. salicifolia. In total, 3.2 million raw short reads were generated with an average length of 150 bp, including 50911358 reads from the male leaf tissue samples and 45850364 reads from the female leaf tissue samples. Clustering of the high-quality reads yielded de novo short read assembly of 50259 transcripts of >100 bp length. The final transcripts were assigned Gene Ontology (GO) terms. The digital expression of genes was studied using the DESeq2 of R package that identified 7180 differentially expressed genes (DEGs) between the male and female plant samples. Further, 10,850 simple sequence repeats, and 8,351 transcription factors, distributed in more than 85 transcription families, were also mined from the final assembled transcriptome. Next, COG and KEGG pathway analyses were performed to assign biological functional terms to the DEGs. The findings of the present study will provide a valuable resource for gene expression discovery and other functional genomics studies aiming towards the selection of candidate genes for the development of sex-specific markers in seabuckthorn and other closely related species.Communicated by Ramaswamy H. Sarma.

ATP-binding cassette transporters in nonmodel plants

Knowledge about plant ATP-binding cassette (ABC) proteins is of great value for sustainable agriculture, economic yield, and the generation of high-quality products, especially under unfavorable growth conditions. We have learned much about ABC proteins in model organisms, notably Arabidopsis thaliana; however, the importance of research dedicated to these transporters extends far beyond Arabidopsis biology. Recent progress in genomic and transcriptomic approaches for nonmodel and noncanonical model plants allows us to look at ABC transporters from a wider perspective and consider chemodiversity and functionally driven adaptation as distinctive mechanisms during their evolution. Here, by considering several representatives from agriculturally important families and recent progress in functional characterization of nonArabidopsis ABC proteins, we aim to bring attention to understanding the evolutionary background, distribution among lineages and possible mechanisms underlying the adaptation of this versatile transport system for plant needs. Increasing the knowledge of ABC proteins in nonmodel plants will facilitate breeding and development of new varieties based on, for example, genetic variations of endogenous genes and/or genome editing, representing an alternative to transgenic approaches.

Pathogen resistance in Sphagneticola trilobata (Singapore daisy): molecular associations and differentially expressed genes in response to disease from a widespread fungus

Understanding the molecular associations underlying pathogen resistance in invasive plant species is likely to provide useful insights into the effective control of alien plants, thereby facilitating the conservation of native biodiversity. In the current study, we investigated pathogen resistance in an invasive clonal plant, Sphagneticola trilobata, at the molecular level. Sphagneticola trilobata (i.e., Singapore daisy) is a noxious weed that affects both terrestrial and aquatic ecosystems, and is less affected by pathogens in the wild than co-occurring native species. We used Illumina sequencing to investigate the transcriptome of S. trilobata following infection by a globally distributed generalist pathogen (Rhizoctonia solani). RNA was extracted from leaves of inoculated and un-inoculated control plants, and a draft transcriptome of S. trilobata was generated to examine the molecular response of this species following infection. We obtained a total of 49,961,014 (94.3%) clean reads for control (un-inoculated plants) and 54,182,844 (94.5%) for the infected treatment (inoculated with R. solani). Our analyses facilitated the discovery of 117,768 de novo assembled contigs and 78,916 unigenes. Of these, we identified 3506 differentially expressed genes and 60 hormones associated with pathogen resistance. Numerous genes, including candidate genes, were associated with plant-pathogen interactions and stress response in S. trilobata. Many recognitions, signaling, and defense genes were differentially regulated between treatments, which were confirmed by qRT-PCR. Overall, our findings improve our understanding of the genes and molecular associations involved in plant defense of a rapidly spreading invasive clonal weed, and serve as a valuable resource for further work on mechanism of disease resistance and managing invasive plants.

Entailing the Next-Generation Sequencing and Metabolome for Sustainable Agriculture by Improving Plant Tolerance

Crop production is a serious challenge to provide food for the 10 billion individuals forecasted to live across the globe in 2050. The scientists' emphasize establishing an equilibrium among diversity and quality of crops by enhancing yield to fulfill the increasing demand for food supply sustainably. The exploitation of genetic resources using genomics and metabolomics strategies can help generate resilient plants against stressors in the future. The innovation of the next-generation sequencing (NGS) strategies laid the foundation to unveil various plants' genetic potential and help us to understand the domestication process to unmask the genetic potential among wild-type plants to utilize for crop improvement. Nowadays, NGS is generating massive genomic resources using wild-type and domesticated plants grown under normal and harsh environments to explore the stress regulatory factors and determine the key metabolites. Improved food nutritional value is also the key to eradicating malnutrition problems around the globe, which could be attained by employing the knowledge gained through NGS and metabolomics to achieve suitability in crop yield. Advanced technologies can further enhance our understanding in defining the strategy to obtain a specific phenotype of a crop. Integration among bioinformatic tools and molecular techniques, such as marker-assisted, QTLs mapping, creation of reference genome, de novo genome assembly, pan- and/or super-pan-genomes, etc., will boost breeding programs. The current article provides sequential progress in NGS technologies, a broad application of NGS, enhancement of genetic manipulation resources, and understanding the crop response to stress by producing plant metabolites. The NGS and metabolomics utilization in generating stress-tolerant plants/crops without deteriorating a natural ecosystem is considered a sustainable way to improve agriculture production. This highlighted knowledge also provides useful research that explores the suitable resources for agriculture sustainability.

Genome-wide identification and characterization of functionally relevant microsatellite markers from transcription factor genes of Tea (Camellia sinensis (L.) O. Kuntze)

Tea, being one of the most popular beverages requires large set of molecular markers for genetic improvement of quality, yield and stress tolerance. Identification of functionally relevant microsatellite or simple sequence repeat (SSR) marker resources from regulatory “Transcription factor (TF) genes” can be potential targets to expedite molecular breeding efforts. In current study, 2776 transcripts encoding TFs harbouring 3687 SSR loci yielding 1843 flanking markers were identified from traits specific transcriptome resource of 20 popular tea cultivars. Of these, 689 functionally relevant SSR markers were successfully validated and assigned to 15 chromosomes (Chr) of CSS genome. Interestingly, 589 polymorphic markers including 403 core-set of TF-SSR markers amplified 2864 alleles in key TF families (bHLH, WRKY, MYB-related, C2H2, ERF, C3H, NAC, FAR1, MYB and G2-like). Their significant network interactions with key genes corresponding to aroma, quality and stress tolerance suggests their potential implications in traits dissection. Furthermore, single amino acid repeat reiteration in CDS revealed presence of favoured and hydrophobic amino acids. Successful deployment of markers for genetic diversity characterization of 135 popular tea cultivars and segregation in bi-parental population suggests their wider utility in high-throughput genotyping studies in tea.

Metabolic and Transcriptomic Profiling of Lilium Leaves Infected With Botrytis elliptica Reveals Different Stages of Plant Defense Mechanisms

Botrytis elliptica, the causal agent of gray mold disease, poses a major threat to commercial Lilium production, limiting its ornamental value and yield. The molecular and metabolic regulation mechanisms of Lilium's defense response to B. elliptica infection have not been completely elucidated. Here, we performed transcriptomic and metabolomic analyses of B. elliptica resistant Lilium oriental hybrid "Sorbonne" to understand the molecular basis of gray mold disease resistance in gray mold disease. A total of 115 differentially accumulated metabolites (DAMs) were detected by comparing the different temporal stages of pathogen infection. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed the differentially expressed genes (DEGs) and DAMs were enriched in the phenylpropanoid and flavonoid pathways at all stages of infection, demonstrating the prominence of these pathways in the defense response of "Sorbonne" to B. elliptica. Network analysis revealed high interconnectivity of the induced defense response. Furthermore, time-course analysis of the transcriptome and a weighted gene coexpression network analysis (WGCNA) led to the identification of a number of hub genes at different stages, revealing that jasmonic acid (JA), salicylic acid (SA), brassinolide (BR), and calcium ions (Ca2+) play a crucial role in the response of "Sorbonne" to fungal infection. Our work provides a comprehensive perspective on the defense response of Lilium to B. elliptica infection, along with a potential transcriptional regulatory network underlying the defense response, thereby offering gene candidates for resistance breeding and metabolic engineering of Lilium.

Progressive Genomic Approaches to Explore Drought- and Salt-Induced Oxidative Stress Responses in Plants under Changing Climate

Drought and salinity are the major environmental abiotic stresses that negatively impact crop development and yield. To improve yields under abiotic stress conditions, drought- and salinity-tolerant crops are key to support world crop production and mitigate the demand of the growing world population. Nevertheless, plant responses to abiotic stresses are highly complex and controlled by networks of genetic and ecological factors that are the main targets of crop breeding programs. Several genomics strategies are employed to improve crop productivity under abiotic stress conditions, but traditional techniques are not sufficient to prevent stress-related losses in productivity. Within the last decade, modern genomics studies have advanced our capabilities of improving crop genetics, especially those traits relevant to abiotic stress management. This review provided updated and comprehensive knowledge concerning all possible combinations of advanced genomics tools and the gene regulatory network of reactive oxygen species homeostasis for the appropriate planning of future breeding programs, which will assist sustainable crop production under salinity and drought conditions.

Can functional genomic diversity provide novel insights into mechanisms of community assembly? A pilot study from an invaded alpine streambed

An important focus of community ecology, including invasion biology, is to investigate functional trait diversity patterns to disentangle the effects of environmental and biotic interactions. However, a notable limitation is that studies usually rely on a small and easy-to-measure set of functional traits, which might not immediately reflect ongoing ecological responses to changing abiotic or biotic conditions, including those that occur at a molecular or physiological level. We explored the potential of using the diversity of expressed genes-functional genomic diversity (FGD)-to understand ecological dynamics of a recent and ongoing alpine invasion. We quantified FGD based on transcriptomic data measured for 26 plant species occurring along adjacent invaded and pristine streambeds. We used an RNA-seq approach to summarize the overall number of expressed transcripts and their annotations to functional categories, and contrasted this with functional trait diversity (FTD) measured from a suite of characters that have been traditionally considered in plant ecology. We found greater FGD and FTD in the invaded community, independent of differences in species richness. However, the magnitude of functional dispersion was greater from the perspective of FGD than from FTD. Comparing FGD between congeneric alien-native species pairs, we did not find many significant differences in the proportion of genes whose annotations matched functional categories. Still, native species with a greater relative abundance in the invaded community compared with the pristine tended to express a greater fraction of genes at significant levels in the invaded community, suggesting that changes in FGD may relate to shifts in community composition. Comparisons of diversity patterns from the community to the species level offer complementary insights into processes and mechanisms driving invasion dynamics. FGD has the potential to illuminate cryptic changes in ecological diversity, and we foresee promising avenues for future extensions across taxonomic levels and macro-ecosystems.

Next-Generation Breeding Strategies for Climate-Ready Crops

Climate change is a threat to global food security due to the reduction of crop productivity around the globe. Food security is a matter of concern for stakeholders and policymakers as the global population is predicted to bypass 10 billion in the coming years. Crop improvement via modern breeding techniques along with efficient agronomic practices innovations in microbiome applications, and exploiting the natural variations in underutilized crops is an excellent way forward to fulfill future food requirements. In this review, we describe the next-generation breeding tools that can be used to increase crop production by developing climate-resilient superior genotypes to cope with the future challenges of global food security. Recent innovations in genomic-assisted breeding (GAB) strategies allow the construction of highly annotated crop pan-genomes to give a snapshot of the full landscape of genetic diversity (GD) and recapture the lost gene repertoire of a species. Pan-genomes provide new platforms to exploit these unique genes or genetic variation for optimizing breeding programs. The advent of next-generation clustered regularly interspaced short palindromic repeat/CRISPR-associated (CRISPR/Cas) systems, such as prime editing, base editing, and de nova domestication, has institutionalized the idea that genome editing is revamped for crop improvement. Also, the availability of versatile Cas orthologs, including Cas9, Cas12, Cas13, and Cas14, improved the editing efficiency. Now, the CRISPR/Cas systems have numerous applications in crop research and successfully edit the major crop to develop resistance against abiotic and biotic stress. By adopting high-throughput phenotyping approaches and big data analytics tools like artificial intelligence (AI) and machine learning (ML), agriculture is heading toward automation or digitalization. The integration of speed breeding with genomic and phenomic tools can allow rapid gene identifications and ultimately accelerate crop improvement programs. In addition, the integration of next-generation multidisciplinary breeding platforms can open exciting avenues to develop climate-ready crops toward global food security.

Genome-wide transcriptional analysis unveils the molecular basis of organ-specific expression of isosteroidal alkaloids biosynthesis in critically endangered Fritillaria roylei Hook

Fritillaria roylei Hook. is a critically endangered high altitude Himalayan medicinal plant species with rich source of pharmaceutically active structurally diverse steroidal alkaloids. Nevertheless, except few marker compounds, the chemistry of the plant remains unexplored. Therefore, in the current study, transcriptome sequencing efforts were made to elucidate isosteroidal alkaloids biosynthesis by creating first organ-specific genomic resource using bulb, stem, and leaf tissues derived from natural populations of Indian Himalayan region. Overall, 349.9 million high quality paired-end reads obtained using NovaSeq 6000 platform were assembled (de novo) into 82,848 unigenes and 31,061 isoforms. Functional annotation and organ specific differential expression (DE) analysis identified 2488 significant DE transcripts, grouped into three potential sub-clusters (sub-cluster I: 728 transcripts; sub-cluster II: 446 transcripts and Sub-cluster III: 1314 transcripts). Subsequently, pathway enrichment (GO, KEGG) and protein-protein network analysis revealed significantly higher enrichment of phenyl-propanoid and steroid backbone including terpenoid, sesquiterpenoid and triterpenoid biosynthesis in bulb. Additionally, upregulated expression of cytochrome P450, UDP-dependent Glucuronosyltransferase families and key transcription factor families (FAR1, bHLH, GRAS, C2H2, TCP and MYB) suggests 'bulb' as a primary site of MVA mediated isosteroidal alkaloids biosynthesis. The comprehensive elucidation of molecular insights in this study is a first step towards the understanding of isosteroidal alkaloid biosynthesis pathway in F. roylei. Furthermore, key genes and regulators identified here can facilitate metabolic engineering of potential bioactive compounds at industrial scale.

Understanding Omics Driven Plant Improvement and de novo Crop Domestication: Some Examples

In the current era, one of biggest challenges is to shorten the breeding cycle for rapid generation of a new crop variety having high yield capacity, disease resistance, high nutrient content, etc. Advances in the "-omics" technology have revolutionized the discovery of genes and bio-molecules with remarkable precision, resulting in significant development of plant-focused metabolic databases and resources. Metabolomics has been widely used in several model plants and crop species to examine metabolic drift and changes in metabolic composition during various developmental stages and in response to stimuli. Over the last few decades, these efforts have resulted in a significantly improved understanding of the metabolic pathways of plants through identification of several unknown intermediates. This has assisted in developing several new metabolically engineered important crops with desirable agronomic traits, and has facilitated the de novo domestication of new crops for sustainable agriculture and food security. In this review, we discuss how "omics" technologies, particularly metabolomics, has enhanced our understanding of important traits and allowed speedy domestication of novel crop plants.

Deciphering the transcriptomic regulation of heat stress responses in Nothofagus pumilio

Global warming is predicted to exert negative impacts on plant growth due to the damaging effect of high temperatures on plant physiology. Revealing the genetic architecture underlying the heat stress response is therefore crucial for the development of conservation strategies, and for breeding heat-resistant plant genotypes. Here we investigated the transcriptional changes induced by heat in Nothofagus pumilio, an emblematic tree species of the sub-Antarctic forests of South America. Through the performance of RNA-seq of leaves of plants exposed to 20°C (control) or 34°C (heat shock), we generated the first transcriptomic resource for the species. We also studied the changes in protein-coding transcripts expression in response to heat. We found 5,214 contigs differentially expressed between temperatures. The heat treatment resulted in a down-regulation of genes related to photosynthesis and carbon metabolism, whereas secondary metabolism, protein re-folding and response to stress were up-regulated. Moreover, several transcription factor families like WRKY or ERF were promoted by heat, alongside spliceosome machinery and hormone signaling pathways. Through a comparative analysis of gene regulation in response to heat in Arabidopsis thaliana, Populus tomentosa and N. pumilio we provide evidence of the existence of shared molecular features of heat stress responses across angiosperms, and identify genes of potential biotechnological application.

Plant synthetic biology for producing potent phyto-antimicrobials to combat antimicrobial resistance

Inappropriate and injudicious use of antimicrobial drugs in human health, hygiene, agriculture, animal husbandry and food industries has contributed significantly to rapid emergence and persistence of antimicrobial resistance (AMR), one of the serious global public health threats. The crisis of AMR versus slower discovery of newer antibiotics put forth a daunting task to control these drug-resistant superbugs. Several phyto-antimicrobials have been identified in recent years with direct-killing (bactericidal) and/or drug-resistance reversal (re-sensitization of AMR phenotypes) potencies. Phyto-antimicrobials may hold the key in combating AMR owing to their abilities to target major microbial drug-resistance determinants including cell membrane, drug-efflux pumps, cell communication and biofilms. However, limited distribution, low intracellular concentrations, eco-geographical variations, beside other considerations like dynamic environments, climate change and over-exploitation of plant-resources are major blockades in full potential exploration phyto-antimicrobials. Synthetic biology (SynBio) strategies integrating metabolic engineering, RNA-interference, genome editing/engineering and/or systems biology approaches using plant chassis (as engineerable platforms) offer prospective tools for production of phyto-antimicrobials. With expanding SynBio toolkit, successful attempts towards introduction of entire gene cluster, reconstituting the metabolic pathway or transferring an entire metabolic (or synthetic) pathway into heterologous plant systems highlight the potential of this field. Through this perspective review, we are presenting herein the current situation and options for addressing AMR, emphasizing on the significance of phyto-antimicrobials in this apparently post-antibiotic era, and effective use of plant chassis for phyto-antimicrobial production at industrial scales along with major SynBio tools and useful databases. Current knowledge, recent success stories, associated challenges and prospects of translational success are also discussed.

Comparative Analysis of SNP Discovery and Genotyping in Fagus sylvatica L. and Quercus robur L. Using RADseq, GBS, and ddRAD Methods

Next-generation sequencing of reduced representation genomic libraries (RRL) is capable of providing large numbers of genetic markers for population genetic studies at relatively low costs. However, one major concern of these types of markers is the precision of genotyping, which is related to the common problem of missing data, which appears to be particularly important in association and genomic selection studies. We evaluated three RRL approaches (GBS, RADseq, ddRAD) and different SNP identification methods (de novo or based on a reference genome) to find the best solutions for future population genomics studies in two economically and ecologically important broadleaved tree species, namely F. sylvatica and Q. robur. We found that the use of ddRAD method coupled with SNP calling based on reference genomes provided the largest numbers of markers (28 k and 36 k for beech and oak, respectively), given standard filtering criteria. Using technical replicates of samples, we demonstrated that more than 80% of SNP loci should be considered as reliable markers in GBS and ddRAD, but not in RADseq data. According to the reference genomes’ annotations, more than 30% of the identified ddRAD loci appeared to be related to genes. Our findings provide a solid support for using ddRAD-based SNPs for future population genomics studies in beech and oak.

PacBio full-length transcriptome of wild apple (Malus sieversii) provides insights into canker disease dynamic response

BACKGROUND

Valsa canker is a serious disease in the stem of Malus sieversii, caused by Valsa mali. However, little is known about the global response mechanism in M. sieversii to V. mali infection.

RESULTS

Phytohormone jasmonic acid (JA) and salicylic acid (SA) profiles and transcriptome analysis were used to elaborate on the dynamic response mechanism. We determined that the JA was initially produced to respond to the necrotrophic pathogen V. mali infection at the early response stage, then get synergistically transduced with SA to respond at the late response stage. Furthermore, we adopted Pacific Biosciences (PacBio) full-length sequencing to identify differentially expressed transcripts (DETs) during the canker response stage. We obtained 52,538 full-length transcripts, of which 8139 were DETs. Total 1336 lncRNAs, 23,737 alternative polyadenylation (APA) sites and 3780 putative transcription factors (TFs) were identified. Additionally, functional annotation analysis of DETs indicated that the wild apple response to the infection of V. mali involves plant-pathogen interaction, plant hormone signal transduction, flavonoid biosynthesis, and phenylpropanoid biosynthesis. The co-expression network of the differentially expressed TFs revealed 264 candidate TF transcripts. Among these candidates, the WRKY family was the most abundant. The MsWRKY7 and MsWRKY33 were highly correlated at the early response stage, and MsWRKY6, MsWRKY7, MsWRKY19, MsWRKY33, MsWRKY40, MsWRKY45, MsWRKY51, MsWRKY61, MsWRKY75 were highly correlated at the late stage.

CONCLUSIONS

The full-length transcriptomic analysis revealed a series of immune responsive events in M. sieversii in response to V. mali infection. The phytohormone signal pathway regulatory played an important role in the response stage. Additionally, the enriched disease resistance pathways and differentially expressed TFs dynamics collectively contributed to the immune response. This study provides valuable insights into a dynamic response in M. sieversii upon the necrotrophic pathogen V. mali infection, facilitates understanding of response mechanisms to canker disease for apple, and provides supports in the identification of potential resistance genes in M. sieversii.