Transfer learning enables identification of multiple types of RNA modifications using nanopore direct RNA sequencing.

Nanopore direct RNA sequencing (DRS) has emerged as a powerful tool for RNA modification identification. However, concurrently detecting multiple types of modifications in a single DRS sample remains a challenge. Here, we develop TandemMod, a transferable deep learning framework capable of detecting multiple types of RNA modifications in single DRS data. To train high-performance TandemMod models, we generate in vitro epitranscriptome datasets from cDNA libraries, containing thousands of transcripts labeled with various types of RNA modifications. We validate the performance of TandemMod on both in vitro transcripts and in vivo human cell lines, confirming its high accuracy for profiling m6A and m5C modification sites. Furthermore, we perform transfer learning for identifying other modifications such as m7G, Ψ, and inosine, significantly reducing training data size and running time without compromising performance. Finally, we apply TandemMod to identify 3 types of RNA modifications in rice grown in different environments, demonstrating its applicability across species and conditions. In summary, we provide a resource with ground-truth labels that can serve as benchmark datasets for nanopore-based modification identification methods, and TandemMod for identifying diverse RNA modifications using a single DRS sample.

Comparative physiological and transcriptome analysis between potassium-deficiency tolerant and sensitive sweetpotato genotypes in response to potassium-deficiency stress.

BACKGROUND: Sweetpotato is a typical ''potassium (K+) favoring'' food crop, which root differentiation process needs a large supply of potassium fertilizer and determine the final root yield. To further understand the regulatory network of the response to low potassium stress, here we analyze physiological and biochemical characteristics, and investigated root transcriptional changes in two sweetpotato genotypes, namely, - K tolerant "Xu32" and - K susceptible"NZ1". RESULT: We found Xu32 had the higher capability of K+ absorption than NZ1 with better growth performance, higher net photosynthetic rate and higher chlorophyll contents under low potassium stress, and identified 889 differentially expressed genes (DEGs) in Xu32, 634 DEGs in NZ1, 256 common DEGs in both Xu32 and NZ1. The Gene Ontology (GO) term in molecular function enrichment analysis revealed that the DEGs under low K+ stress are predominately involved in catalytic activity, binding, transporter activity and antioxidant activity. Moreover, the more numbers of identified DEGs in Xu32 than that in NZ1 responded to K+-deficiency belong to the process of photosynthesis, carbohydrate metabolism, ion transport, hormone signaling, stress-related and antioxidant system may result in different ability to K+-deficiency tolerance. The unique genes in Xu32 may make a great contribution to enhance low K+ tolerance, and provide useful information for the molecular regulation mechanism of K+-deficiency tolerance in sweetpotato. CONCLUSIONS: The common and distinct expression pattern between the two sweetpotato genotypes illuminate a complex mechanism response to low potassium exist in sweetpotato. The study provides some candidate genes, which can be used in sweetpotato breeding program for improving low potassium stress tolerance.

Multiomics analyses of two Leonurus species illuminate leonurine biosynthesis and its evolution.

The Lamiaceae family is renowned for its terpenoid-based medicinal components, but Leonurus, which has traditional medicinal uses, stands out for its alkaloid-rich composition. Leonurine, the principal active compound found in Leonurus, has demonstrated promising effects in reducing blood lipids and treating strokes. However, the biosynthetic pathway of leonurine remains largely unexplored. Here, we present the chromosome-level genome sequence assemblies of Leonurus japonicus, known for its high leonurine production, and Leonurus sibiricus, characterized by very limited leonurine production. By integrating genomics, RNA sequencing, metabolomics, and enzyme activity assay data, we constructed the leonurine biosynthesis pathway and identified the arginine decarboxylase (ADC), uridine diphosphate glucosyltransferase (UGT), and serine carboxypeptidase-like (SCPL) acyltransferase enzymes that catalyze key reactions in this pathway. Further analyses revealed that the UGT-SCPL gene cluster evolved by gene duplication in the ancestor of Leonurus and neofunctionalization of SCPL in L. japonicus, which contributed to the accumulation of leonurine specifically in L. japonicus. Collectively, our comprehensive study illuminates leonurine biosynthesis and its evolution in Leonurus.

Haplotype-based phylogenetic analysis and population genomics uncover the origin and domestication of sweetpotato.

The hexaploid sweetpotato (Ipomoea batatas) is one of the most important root crops worldwide. However, its genetic origin remains controversial, and its domestication history remains unknown. In this study, we used a range of genetic evidence and a newly developed haplotype-based phylogenetic analysis to identify two probable progenitors of sweetpotato. The diploid progenitor was likely closely related to Ipomoea aequatoriensis and contributed the B1 subgenome, IbT-DNA2, and the lineage 1 type of chloroplast genome to sweetpotato. The tetraploid progenitor of sweetpotato was most likely I. batatas 4x, which donated the B2 subgenome, IbT-DNA1, and the lineage 2 type of chloroplast genome. Sweetpotato most likely originated from reciprocal crosses between the diploid and tetraploid progenitors, followed by a subsequent whole-genome duplication. In addition, we detected biased gene exchanges between the subgenomes; the rate of B1 to B2 subgenome conversions was nearly three times higher than that of B2 to B1 subgenome conversions. Our analyses revealed that genes involved in storage root formation, maintenance of genome stability, biotic resistance, sugar transport, and potassium uptake were selected during the speciation and domestication of sweetpotato. This study sheds light on the evolution of sweetpotato and paves the way for improvement of this crop.

Gap-free genome assembly and CYP450 gene family analysis reveal the biosynthesis of anthocyanins in Scutellaria baicalensis.

Scutellaria baicalensis Georgi, a member of the Lamiaceae family, is a widely utilized medicinal plant. The flavones extracted from S. baicalensis contribute to numerous health benefits, including anti-inflammatory, antiviral, and anti-tumor activities. However, the incomplete genome assembly hinders biological studies on S. baicalensis. This study presents the first telomere-to-telomere (T2T) gap-free genome assembly of S. baicalensis through the integration of Pacbio HiFi, Nanopore ultra-long and Hi-C technologies. A total of 384.59 Mb of genome size with a contig N50 of 42.44 Mb was obtained, and all sequences were anchored into nine pseudochromosomes without any gap or mismatch. In addition, we analysed the major cyanidin- and delphinidin-based anthocyanins involved in the determination of blue-purple flower using a widely-targeted metabolome approach. Based on the genome-wide identification of Cytochrome P450 (CYP450) gene family, three genes (SbFBH1, 2, and 5) encoding flavonoid 3'-hydroxylases (F3'Hs) and one gene (SbFBH7) encoding flavonoid 3'5'-hydroxylase (F3'5'H) were found to hydroxylate the B-ring of flavonoids. Our studies enrich the genomic information available for the Lamiaceae family and provide a toolkit for discovering CYP450 genes involved in the flavonoid decoration.

m6A RNA demethylase AtALKBH9B promotes mobilization of a heat-activated long terminal repeat retrotransposon in Arabidopsis.

Transposons are mobile and ubiquitous DNA molecules that can cause vast genomic alterations. In plants, it is well documented that transposon mobilization is strongly repressed by DNA methylation; however, its regulation at the posttranscriptional level remains relatively uninvestigated. Here, we suggest that transposon RNA is marked by m6A RNA methylation and can be localized in stress granules (SGs). Intriguingly, SG-localized AtALKBH9B selectively demethylates a heat-activated retroelement, Onsen, and thereby releases it from spatial confinement, allowing for its mobilization. In addition, we show evidence that m6A RNA methylation contributes to transpositional suppression by inhibiting virus-like particle assembly and extrachromosomal DNA production. In summary, this study unveils a previously unknown role for m6A in the suppression of transposon mobility and provides insight into how transposons counteract the m6A-mediated repression mechanism by hitchhiking the RNA demethylase of the host.

Sweet potato ADP-glucose pyrophosphorylase small subunit affects vegetative growth, starch content and storage root yield.

The development of storage roots is a key factor determining the yields of crop plants, including sweet potato. Here, using combined bioinformatic and genomic approaches, we identified a sweet potato yield-related gene, ADP-glucose pyrophosphorylase (AGP) small subunit (IbAPS). We found that IbAPS positively affects AGP activity, transitory starch biosynthesis, leaf development, chlorophyll metabolism, and photosynthesis, ultimately affecting the source strength. IbAPS overexpression in sweet potato led to increased vegetative biomass and storage root yield. RNAi of IbAPS resulted in reduced vegetative biomass, accompanied with a slender stature and stunted root development. In addition to the effects on root starch metabolism, we found that IbAPS affects other storage root development-associated events, including lignification, cell expansion, transcriptional regulation, and production of the storage protein sporamins. A combinatorial analysis based on transcriptomes, as well as morphological and physiological data, revealed that IbAPS affects several pathways that determine development of vegetative tissues and storage roots. Our work establishes an important role of IbAPS in concurrent control of carbohydrate metabolism, plant growth, and storage root yield. We showed that upregulation of IbAPS results in superior sweet potato with increased green biomass, starch content, and storage root yield. The findings expand our understanding of the functions of AGP enzymes and advances our ability to increase the yield of sweet potato and, perhaps, other crop plants.

Exploring and exploiting genetics and genomics for sweetpotato improvement: Status and perspectives.

Sweetpotato (Ipomoea batatas (L.) Lam.) is one of the most important root crops cultivated worldwide. Because of its adaptability, high yield potential, and nutritional value, sweetpotato has become an important food crop, particularly in developing countries. To ensure adequate crop yields to meet increasing demand, it is essential to enhance the tolerance of sweetpotato to environmental stresses and other yield-limiting factors. The highly heterozygous hexaploid genome of I. batatas complicates genetic studies and limits improvement of sweetpotato through traditional breeding. However, application of next-generation sequencing and high-throughput genotyping and phenotyping technologies to sweetpotato genetics and genomics research has provided new tools and resources for crop improvement. In this review, we discuss the genomics resources that are available for sweetpotato, including the current reference genome, databases, and available bioinformatics tools. We systematically review the current state of knowledge on the polyploid genetics of sweetpotato, including studies of its origin and germplasm diversity and the associated mapping of important agricultural traits. We then outline the conventional and molecular breeding approaches that have been applied to sweetpotato. Finally, we discuss future goals for genetic studies of sweetpotato and crop improvement via breeding in combination with state-of-the-art multi-omics approaches such as genomic selection and gene editing. These approaches will advance and accelerate genetic improvement of this important root crop and facilitate its sustainable global production.

Specific Flavonoids and Their Biosynthetic Pathway in Scutellaria baicalensis.

Scutellaria baicalensis, is one of the most traditional medicinal plants in the Lamiaceae family, and has been widely used to treat liver and lung complaints and as a complementary cancer treatment in traditional Chinese medicine. The preparation from its roots, called "Huang Qin," is rich in specialized flavones such as baicalein, wogonin, and their glycosides which lack a 4'-hydroxyl group on the B ring (4'-deoxyflavones), with anti-tumor, antioxidant, and antiviral activities. Baicalein has recently been reported to inhibit the replication of the COVID-19 virus. These 4'-deoxyflavones are found only in the order Lamiales and were discovered in the genus Scutellaria, suggesting that a new metabolic pathway synthesizing 4'-deoxyflavones evolved recently in this genus. In this review, we focus on the class of 4'-deoxyflavones in S. baicalensis and their pharmacological properties. We also describe the apparent evolutionary route taken by the genes encoding enzymes involved in the novel, root-specific, biosynthetic pathway for baicalein and wogonin, which provides insights into the evolution of specific flavone biosynthetic pathways in the mint family.

Characterization of UDP-glycosyltransferase family members reveals how major flavonoid glycoside accumulates in the roots of Scutellaria baicalensis.

BACKGROUND: Flavonoid glycosides extracted from roots of Scutellaria baicalensis exhibit strong pharmaceutical antitumor, antioxidative, anti-inflammatory, and antiviral activities. UDP glycosyltransferase (UGT) family members are responsible for the transfer of a glycosyl moiety from UDP sugars to a wide range of acceptor flavonoids. Baicalin is the major flavonoid glycoside found in S. baicalensis roots, and its aglycone baicalein is synthesized from a specially evolved pathway that has been elucidated. However, it is necessary to carry out a genome-wide study of genes involved in 7-O-glucuronidation, the final biosynthesis step of baicalin, which might elucidate the relationship between the enzymes and the metabolic accumulation patterns in this medicinal plant. RESULTS: We reported the phylogenetic analysis, tissue-specific expression, biochemical characterization and evolutionary analysis of glucosyltransferases (SbUGTs) and glucuronosyltransferases (SbUGATs) genes based on the recently released genome of S. baicalensis. A total of 124 UGTs were identified, and over one third of them were highly expressed in roots. In vitro enzyme assays showed that 6 SbUGTs could use UDP-glucose as a sugar donor and convert baicalein to oroxin A (baicalein 7-O-glucoside), while 4 SbUGATs used only UDP-glucuronic acid as the sugar donor and catalyzed baicalein to baicalin. SbUGAT4 and SbUGT2 are the most highly expressed SbUGAT and SbUGT genes in root tissues, respectively. Kinetic measurements revealed that SbUGAT4 had a lower Km value and higher Vmax/Km ratio to baicalein than those of SbUGT2. Furthermore, tandem duplication events were detected in SbUGTs and SbUGATs. CONCLUSIONS: This study demonstrated that glucosylation and glucuronidation are two major glycosylated decorations in the roots of S. baicalensis. Higher expression level and affinity to substrate of SbUGAT4, and expansion of this gene family contribute high accumulation of baicalin in the root of S. baicalensis.

Genome-wide identification of agronomically important genes in outcrossing crops using OutcrossSeq.

Many important crops (e.g., tuber, root, and tree crops) are cross-pollinating. For these crops, no inbred lines are available for genetic study and breeding because they are self-incompatible, clonally propagated, or have a long generation time, making the identification of agronomically important genes difficult, particularly in crops with a complex autopolyploid genome. In this study, we developed a method, OutcrossSeq, for mapping agronomically important loci in outcrossing crops based on whole-genome low-coverage resequencing of a large genetic population, and designed three computation algorithms in OutcrossSeq for different types of outcrossing populations. We applied OutcrossSeq to a tuberous root crop (sweet potato, autopolyploid), a tree crop (walnut tree, highly heterozygous diploid), and hybrid crops (double-cross populations) to generate high-density genotype maps for the outcrossing populations, which enable precise identification of genomic loci underlying important agronomic traits. Candidate causative genes at these loci were detected based on functional clues. Taken together, our results indicate that OutcrossSeq is a robust and powerful method for identifying agronomically important genes in heterozygous species, including polyploids, in a cost-efficient way. The OutcrossSeq software and its instruction manual are available for downloading at www.xhhuanglab.cn/tool/OutcrossSeq.html.

The genome of Tripterygium wilfordii and characterization of the celastrol biosynthesis pathway.

Tripterygium wilfordii is a vine from the Celastraceae family that is used in traditional Chinese medicine (TCM). The active ingredient, celastrol, is a friedelane-type pentacyclic triterpenoid with putative roles as an antitumor, immunosuppressive, and anti-obesity agent. Here, we report a reference genome assembly of T. wilfordii with high-quality annotation using a hybrid sequencing strategy. The total genome size obtained is 340.12 Mb, with a contig N50 value of 3.09 Mb. We successfully anchored 91.02% of sequences into 23 pseudochromosomes using high-throughput chromosome conformation capture (Hi-C) technology. The super-scaffold N50 value was 13.03 Mb. We also annotated 31,593 structural genes, with a repeat percentage of 44.31%. These data demonstrate that T. wilfordii diverged from Malpighiales species approximately 102.4 million years ago. By integrating genome, transcriptome and metabolite analyses, as well as in vivo and in vitro enzyme assays of two cytochrome P450 (CYP450) genes, TwCYP712K1 and TwCYP712K2, it is possible to investigate the second biosynthesis step of celastrol and demonstrate that this was derived from a common ancestor. These data provide insights and resources for further investigation of pathways related to celastrol, and valuable information to aid the conservation of resources, as well as understand the evolution of Celastrales.

Genome-wide characterization and expression analysis of HAK K+ transport family in Ipomoea.

The potassium transporter high-affinity K+ transporter/K+ uptake permease/K+ transporter (HAK/KUP/KT) family plays a vital role in potassium uptake, and potassium ion (K+)-mediated environmental stress. In the present study, we identified 22 IbHAK/KUP/KT (HAK) genes in sweet potato [Ipomoea batata (L.) Lam] and the same number of HAK genes from sweet potato wild relative Ipomoea trifida. Phylogeny analysis indicated that the HAKs can be divided into five clades. Chromosomal distribution and genome synteny analyses revealed two tandem-duplicated gene pairs IbHAK16/17 and IbHAK17/18 on chromosomes 13 and eight segmental-duplicated gene pairs on chromosomes 1, 3, 5, 8, 10, 12, 14 among the IbHAK gene family. Eleven orthologous HAK gene pairs between I. batata and I. trifida were involved in the duplication of genomic blocks based on comparative genomic analysis. The Ka/Ks ratios of these IbHAK genes ranged from 0.02 to 0.55(< 1), further indicated that purifying selection was the primary force driving the evolution of HAKs in Ipomoea. A heat map based on RNA-seq data showed that 13 HAKs in Xushu32 (a K+-tolerant sweet potato genotype) and 10 HAKs in Ningzi1 (a K+-sensitive sweet potato genotype) in response to K+ deficiency stress. Quantitative real-time PCR (qRT-PCR) analysis revealed IbHAK2, -3, -8, -10, -11, -18, -19, and -21 were induced in both Xushu32 and Ningzi1 under low K+ stress. Compared with other IbHAK genes, IbHAK8 showed more strongly upregulation after exposure to drought and salt stress. Furthermore, co-expression analysis showed that only IbHAK8 of 22 IbHAK genes involved in network interactions with 30 genes related to abiotic and biotic stresses. Taken together, these results are helpful for further functional studies on IbHAK and molecular breeding of sweet potato. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13205-020-02552-3.

Ancient events and climate adaptive capacity shaped distinct chloroplast genetic structure in the oak lineages.

BACKGROUND: Understanding the origin of genetic variation is the key to predict how species will respond to future climate change. The genus Quercus is a species-rich and ecologically diverse woody genus that dominates a wide range of forests and woodland communities of the Northern Hemisphere. Quercus thus offers a unique opportunity to investigate how adaptation to environmental changes has shaped the spatial genetic structure of closely related lineages. Furthermore, Quercus provides a deep insight into how tree species will respond to future climate change. This study investigated whether closely related Quercus lineages have similar spatial genetic structures and moreover, what roles have their geographic distribution, ecological tolerance, and historical environmental changes played in the similar or distinct genetic structures. RESULTS: Despite their close relationships, the three main oak lineages (Quercus sections Cyclobalanopsis, Ilex, and Quercus) have different spatial genetic patterns and occupy different climatic niches. The lowest level and most homogeneous pattern of genetic diversity was found in section Cyclobalanopsis, which is restricted to warm and humid climates. The highest genetic diversity and strongest geographic genetic structure were found in section Ilex, which is due to their long-term isolation and strong local adaptation. The widespread section Quercus is distributed across the most heterogeneous range of environments; however, it exhibited moderate haplotype diversity. This is likely due to regional extinction during Quaternary climatic fluctuation in Europe and North America. CONCLUSIONS: Genetic variations of sections Ilex and Quercus were significantly predicted by geographic and climate variations, while those of section Cyclobalanopsis were poorly predictable by geographic or climatic diversity. Apart from the different historical environmental changes experienced by different sections, variation of their ecological or climatic tolerances and physiological traits induced varying responses to similar environment changes, resulting in distinct spatial genetic patterns.

The Reference Genome Sequence of Scutellaria baicalensis Provides Insights into the Evolution of Wogonin Biosynthesis.

Scutellaria baicalensis Georgi is important in Chinese traditional medicine where preparations of dried roots, "Huang Qin," are used for liver and lung complaints and as complementary cancer treatments. We report a high-quality reference genome sequence for S. baicalensis where 93% of the 408.14-Mb genome has been assembled into nine pseudochromosomes with a super-N50 of 33.2 Mb. Comparison of this sequence with those of closely related species in the order Lamiales, Sesamum indicum and Salvia splendens, revealed that a specialized metabolic pathway for the synthesis of 4'-deoxyflavone bioactives evolved in the genus Scutellaria. We found that the gene encoding a specific cinnamate coenzyme A ligase likely obtained its new function following recent mutations, and that four genes encoding enzymes in the 4'-deoxyflavone pathway are present as tandem repeats in the genome of S. baicalensis. Further analyses revealed that gene duplications, segmental duplication, gene amplification, and point mutations coupled to gene neo- and subfunctionalizations were involved in the evolution of 4'-deoxyflavone synthesis in the genus Scutellaria. Our study not only provides significant insight into the evolution of specific flavone biosynthetic pathways in the mint family, Lamiaceae, but also will facilitate the development of tools for enhancing bioactive productivity by metabolic engineering in microbes or by molecular breeding in plants. The reference genome of S. baicalensis is also useful for improving the genome assemblies for other members of the mint family and offers an important foundation for decoding the synthetic pathways of bioactive compounds in medicinal plants.

The Application and Limitation of Universal Chloroplast Markers in Discriminating East Asian Evergreen Oaks.

The East Asian subtropics mostly occupied by evergreen broad-leaved forests (EBLFs), is one of the global diversity centers for evergreen oaks. Evergreen oaks are keystone canopy trees in EBLFs with important ecosystem function and crucial significance for regional biodiversity conservation. However, the species composition and diversity of Asian evergreen oaks are poorly understood. Here, we test whether the four chloroplast markers atpI-atpH, matK, psbA-trnH, and ycf1, can discriminate the two evergreen oak sections in Asia - Cyclobalanopsis and Ilex. Two hundred and seventy-two individuals representing 57 species were scanned and 17 species from other oaks sections were included for phylogenetic reconstruction. The genetic diversity of the Quercus sections was also compared. Overall, we found that universal chloroplast DNA (cpDNA) barcoding markers could resolve two clades in Quercus, i.e., subgenus Cerris (Old World Clade) and subgenus Quercus (New World Clade). The chloroplast markers distinguished the main sections, with few exceptions. Each cpDNA region showed no barcoding gap and none of them provided good resolution at the species level. The best species resolution (27.78%) was obtained when three or four markers were combined and analyzed using BLAST. The high conservation of the cpDNA and complicated evolutionary patterns, due to incomplete lineage sorting, interspecific hybridization and introgressions may hinder the ability of cpDNA markers to discriminate different species. When comparing diversification pattern across Quercus sections (Cyclobalanopsis, Ilex, Cerris, Quercus, and Protobalanus), we found that section Ilex was the most genetically diverse, and section Cyclobalanopsis was lower genetically diverse. This diversification pattern may have resulted from the interplay of the Eurasia Cenozoic tectonic movements, climate changes and different niches of their ancestral lineages.