Transcriptomic Analysis of the Regulation of Rhizome Formation in Temperate and Tropical Lotus (Nelumbo nucifera)

Two plants improve stress response of a subterranean herbivore by downregulating amphetamine addiction pathways

Introduction

Captivity serves as the primary method for enhancing animal survival and productivity. However, the stress induced by confinement can hinder animal growth and reproduction. The administration of drugs to captive animals can effectively regulate their stress response and can also be used inartificial breeding, reproduction, and experimental animalization of wild species. The plateau zokor (Eospalax baileyi), a subterranean herbivore, experiences significant stress during the captive process owing to its unique habitat.

Methods

In our study, we utilized Radix astragali (RA) and Acanthopanax senticosus (AS) extracts to intervene in the stress response of plateau zokors.

Results

Our findings demonstrated that RA and AS treatment considerably improved food intake and reduced weight loss, stress-related behavior, and stress hormone levels in plateau zokors. Furthermore, the excitatory pathway of amphetamine addition in the hypothalamus was suppressed by RA and AS treatment, acting through the Grin and Prkc gene families. Notably, after RA treatment, the extracellular matrix-receptor interaction pathway, enriched by the Col1a1/3a1/1a2/6a1 gene, was significantly upregulated, potentially enhancing the immune function of captive plateau zokors.

Discussion

In conclusion, our research demonstrates that RA and AS treatment can effectively alleviate the stress response of plateau zokors in captive environments. The downregulation of the excitation pathway and upregulation of the immune pathway offer valuable insights into the response and potential mechanisms of plant-based drugs in mitigating animal stress.

Integrating Full-Length and Second-Generation Transcriptomics to Reveal Differentially Expressed Genes Associated with the Development of Corydalis yanhusuo Tuber

Corydalis yanhusuo is a medicinal herb in China that has been widely used to treat various kinds of pain. The tuber is the main organ of C. yanhusuo used for medicinal purposes, but changes in related genes during the development of the tuber have rarely been reported. To identify the differentially expressed genes during tuber development, C. yanhusuo full-length transcriptomic sequencing was performed using single-molecule real-time technology, and tubers at three development stages were selected for comparative transcriptome analysis. A total of 90,496 full-length non-chimeric transcripts were obtained, and 19,341 transcripts were annotated in at least one public database. A total of 9221 differentially expressed genes were identified during the swelling process of C. yanhusuo tuber. A Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis revealed that differentially expressed genes associated with a "starch and sucrose metabolism pathway", "phenylpropanoid biosynthesis pathway", "isoquinoline alkaloid biosynthesis pathway", "zeatin biosynthesis pathway", and "brassinosteroid biosynthesis pathway" were predominantly enriched. In addition, the genes involved in cell wall metabolism were potentially associated with tuber swelling. These processes regulated and were involved in C. yanhusuo tuber development. The results provide a foundation for further research on tuber formation in medicinal plants.

Petal size is controlled by the MYB73/TPL/HDA19-miR159-CKX6 module regulating cytokinin catabolism in Rosa hybrida

The size of plant lateral organs is determined by well-coordinated cell proliferation and cell expansion. Here, we report that miR159, an evolutionarily conserved microRNA, plays an essential role in regulating cell division in rose (Rosa hybrida) petals by modulating cytokinin catabolism. We uncover that Cytokinin Oxidase/Dehydrogenase6 (CKX6) is a target of miR159 in petals. Knocking down miR159 levels results in the accumulation of CKX6 transcripts and earlier cytokinin clearance, leading to a shortened cell division period and smaller petals. Conversely, knocking down CKX6 causes cytokinin accumulation and a prolonged developmental cell division period, mimicking the effects of exogenous cytokinin application. MYB73, a R2R3-type MYB transcription repressor, recruits a co-repressor (TOPLESS) and a histone deacetylase (HDA19) to form a suppression complex, which regulates MIR159 expression by modulating histone H3 lysine 9 acetylation levels at the MIR159 promoter. Our work sheds light on mechanisms for ensuring the correct timing of the exit from the cell division phase and thus organ size regulation by controlling cytokinin catabolism.

Histology, physiology, and transcriptomic and metabolomic profiling reveal the developmental dynamics of annual shoots in tree peonies (Paeonia suffruticosa Andr.)

The development of tree peony annual shoots is characterized by "withering", which is related to whether there are bud points in the leaf axillaries of annual shoots. However, the mechanism of "withering" in tree peony is still unclear. In this study, Paeonia ostii 'Fengdan' and P. suffruticosa 'Luoyanghong' were used to investigate dynamic changes of annual shoots through anatomy, physiology, transcriptome, and metabolome. The results demonstrated that the developmental dynamics of annual shoots of the two cultivars were comparable. The withering degree of P. suffruticosa 'Luoyanghong' was higher than that of P. ostii 'Fengdan', and their upper internodes of annual flowering shoots had a lower degree of lignin deposition, cellulose, C/N ratio, showing no obvious sclerenchyma, than the bottom ones and the whole internodes of vegetative shoot, which resulted in the "withering" of upper internodes. A total of 36 phytohormone metabolites were detected, of which 33 and 31 were detected in P. ostii 'Fengdan' and P. suffruticosa 'Luoyanghong', respectively. In addition, 302 and 240 differentially expressed genes related to lignin biosynthesis, carbon and nitrogen metabolism, plant hormone signal transduction, and zeatin biosynthesis were screened from the two cultivars. Furtherly, 36 structural genes and 40 transcription factors associated with the development of annual shoots were highly co-expressed, and eight hub genes involved in this developmental process were identified. Consequently, this study explained the developmental dynamic on the varied annual shoots through multi-omics, providing a theoretical foundation for germplasm innovation and the mechanized harvesting of tree peony annual shoots.

Screening of Early Flowering Lotus (Nelumbo nucifera Gaertn.) Cultivars and Effects of Different Cultivars on Flowering Period

Flowering time is an important trait that determines the breeding process of ornamental plants. The flowering period of lotus (Nelumbo nucifera Gaertn.) is mainly concentrated in June-August. During this period, the weather is hot and there are few tourists, which made many lotus scenic spots difficult to operate. People have a strong demand for early flowering lotus cultivars. In this paper, 30 lotus cultivars with high ornamental value were selected as materials and their phenological periods were observed for two consecutive years in 2019 and 2020. A number of cultivars with early flowering potential and stable flowering periods, such as 'Fenyanzi', 'Chengshanqiuyue', 'Xianghumingyue' and 'Wuzhilian', were screened by K-Means clustering method. The relationship between accumulated temperature and flowering time of 19 lotus cultivars at different growth stages was analyzed. It was found that lotus cultivars with early flowering traits could adapt well to the changes of early environmental temperature and were not affected by low temperature. On the other hand, by analyzing the relationship between different traits and flowering time of three typical cultivars, such as rhizome weight, phenological period, etc., it shows that the nutrient content of the rhizome and the early morphology of plants will affect the flowering time. These results provide a reference for the formation of a systematic lotus early flowering cultivar breeding mechanism and the establishment of a perfect flowering regulation technology system, which can further improve the ornamental value of lotus and promote industrial development.

Can the Wild Perennial, Rhizomatous Rice Species Oryza longistaminata be a Candidate for De Novo Domestication?

As climate change intensifies, the development of resilient rice that can tolerate abiotic stresses is urgently needed. In nature, many wild plants have evolved a variety of mechanisms to protect themselves from environmental stresses. Wild relatives of rice may have abundant and virtually untapped genetic diversity and are an essential source of germplasm for the improvement of abiotic stress tolerance in cultivated rice. Unfortunately, the barriers of traditional breeding approaches, such as backcrossing and transgenesis, make it challenging and complex to transfer the underlying resilience traits between plants. However, de novo domestication via genome editing is a quick approach to produce rice with high yields from orphans or wild relatives. African wild rice, Oryza longistaminata, which is part of the AA-genome Oryza species has two types of propagation strategies viz. vegetative propagation via rhizome and seed propagation. It also shows tolerance to multiple types of abiotic stress, and therefore O. longistaminata is considered a key candidate of wild rice for heat, drought, and salinity tolerance, and it is also resistant to lodging. Importantly, O. longistaminata is perennial and propagates also via rhizomes both of which are traits that are highly valuable for the sustainable production of rice. Therefore, O. longistaminata may be a good candidate for de novo domestication through genome editing to obtain rice that is more climate resilient than modern elite cultivars of O. sativa.

Transcriptomic and Metabolomic Analyses Reveal That Fullerol Improves Drought Tolerance in Brassica napus L

Carbon nanoparticles have potential threats to plant growth and stress tolerance. The polyhydroxy fullerene-fullerol (one of the carbon nanoparticles) could increase biomass accumulation in several plants subjected to drought; however, the underlying molecular and metabolic mechanisms governed by fullerol in improving drought tolerance in Brassica napus remain unclear. In the present study, exogenous fullerol was applied to the leaves of B. napus seedlings under drought conditions. The results of transcriptomic and metabolomic analyses revealed changes in the molecular and metabolic profiles of B. napus. The differentially expressed genes and the differentially accumulated metabolites, induced by drought or fullerol treatment, were mainly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to carbohydrate metabolism (e.g., "carbon metabolism" and "galactose metabolism"), amino acid metabolism (e.g., "biosynthesis of amino acids" and "arginine and proline metabolism"), and secondary metabolite metabolism (e.g., "biosynthesis of secondary metabolites"). For carbohydrate metabolism, the accumulation of oligosaccharides (e.g., sucrose) was decreased, whereas that of monosaccharides (e.g., mannose and myo-inositol) was increased by drought. With regard to amino acid metabolism, under drought stress, the accumulation of amino acids such as phenylalanine and tryptophan decreased, whereas that of glutamate and proline increased. Further, for secondary metabolite metabolism, B. napus subjected to soil drying showed a reduction in phenolics and flavonoids, such as hyperoside and trans-3-coumaric acid. However, the accumulation of carbohydrates was almost unchanged in fullerol-treated B. napus subjected to drought. When exposed to water shortage, the accumulation of amino acids, such as proline, was decreased upon fullerol treatment. However, that of phenolics and flavonoids, such as luteolin and trans-3-coumaric acid, was enhanced. Our findings suggest that fullerol can alleviate the inhibitory effects of drought on phenolics and flavonoids to enhance drought tolerance in B. napus.

Identification of LsPIN1 gene and its potential functions in rhizome turning of Leymus secalinus

Background

Continuous tilling and the lateral growth of rhizomes confer rhizomatous grasses with the unique ability to laterally expand, migrate and resist disturbances. They play key roles especially in degraded grasslands, deserts, sand dunes, and other fragile ecological system. The rhizomatous plant Leymus secalinus has both rhizome buds and tiller buds that grow horizontally and upward at the ends of rhizome differentiation and elongation, respectively. The mechanisms of rhizome formation and differentiation in L. secalinus have not yet been clarified.

Results

In this study, we found that the content of gibberellin A3 (GA3) and indole-3-acetic acid (IAA) were significantly higher in upward rhizome tips than in horizontal rhizome tips; by contrast, the content of methyl jasmonate and brassinolide were significantly higher in horizontal rhizome tips than in upward rhizome tips. GA3 and IAA could stimulate the formation and turning of rhizomes. An auxin efflux carrier gene, LsPIN1, was identified from L. secalinus based on previous transcriptome data. The conserved domains of LsPIN1 and the relationship of LsPIN1 with PIN1 genes from other plants were analyzed. Subcellular localization analysis revealed that LsPIN1 was localized to the plasma membrane. The length of the primary roots (PRs) and the number of lateral roots (LRs) were higher in Arabidopsis thaliana plants overexpressing LsPIN1 than in wild-type (Col-0) plants. Auxin transport was altered and the gravitropic response and phototropic response were stronger in 35S:LsPIN1 transgenic plants compared with Col-0 plants. It also promoted auxin accumulation in root tips.

Conclusion

Our findings indicated that LsPIN1 plays key roles in auxin transport and root development. Generally, our results provide new insights into the regulatory mechanisms underlying rhizome development in L. secalinus.

Nontarget Metabolites of Rhizomes of Edible Sacred Lotus Provide New Insights into Rhizome Browning

The “edible rhizome” variant of Nelumbo nucifera with various cultivars has a long history of use as a food in East Asia. In this study, 48 target metabolites were untargeted and identified in 212 rhizome cultivars (tropical and temperate types) using ultraperformance liquid chromatography-electrospray ionization quadrupole time-of-flight high-resolution mass spectrometry; among these, 32 compounds were newly reported in the rhizome. Combined with the browning phenotype of 212 lotus rhizomes, (epi) catechin, norarmepavine, and N-feruloyl-3-methoxytyramine were used as predominant chemical markers to separate different degrees of lotus rhizome browning. p-Coumaroyltyramine and N-trans-feruloyltyramine were selected as predominant chemical markers to investigate the differential expression between tropical and temperate lotus using principal component analysis and orthogonal partial least squares discriminant analysis. Shared and unique structure plots were used to compare the outcomes of the ecotype and browning OPLS model, showing that variation in tropical lotus rhizome browning is not obvious; this will be of great importance for genetic improvement by providing a hereditary basis.

Transcriptomic analysis of succulent stem development of Chinese kale (Brassica oleracea var. alboglabra Bailey) and its synthetic allotetraploid via RNA sequencing

Chinese kale (Brassica oleracea var. alboglabra Bailey, CC) is a succulent stem vegetable in the Brassica family. Its allotetraploid (AACC) vegetable germplasm, which was synthesized via distant hybridization with the colloquially named 'yellow turnip' (B. rapa L. ssp. rapifera Matzg., AA), has a swelling stem similar to CC. To address the molecular mechanism of stem development for CC and AACC, RNA sequencing (RNA-seq) was used to investigate transcriptional regulation of their stem development at three key stages including 28 days, 42 days and the bolting stage (BS) after sowing. As a result, 32,642, 32,665, 33,816, 32,147, 32,293 and 32,275 genes were identified in six corresponding cDNA libraries. Among them, 25,459 genes were co-expressed, while 7,183, 7,206, 8,357, 6,688, 6,834 and 6,814 genes were specifically expressed. Additionally, a total of 29,222 differentially expressed genes (DEGs) were found for functional enrichment as well as many genes involved in plant hormones including gibberellin (GA), abscisic acid (ABA), cytokinin (CTK) and auxin (AUX). Based on gene expression consistency between CC and AACC, the gene families including DELLA, GID, PYR/PYL, PP2C, A-ARR and AUX/IAA might be related to stem development. Among these, eight genes including Bo00834s040, Bo5g093140, Bo6g086770, Bo9g070200, Bo7g116570, Bo3g054410, Bo7g093470 and Bo5g136600 may play important roles in stem development based on their remarkable expression levels as confirmed by qRT-PCR. These findings provide a new theoretical basis for understanding the molecular mechanism of stem development in Brassica vegetable stem breeding.

Genome-wide association study of traits in sacred lotus uncovers MITE-associated variants underlying stamen petaloid and petal number variations

Understanding the genetic variants responsible for floral trait diversity is important for the molecular breeding of ornamental flowers. Widely used in water gardening for thousands of years, the sacred lotus exhibits a wide range of diversity in floral organs. Nevertheless, the genetic variations underlying various morphological characteristics in lotus remain largely unclear. Here, we performed a genome-wide association study of sacred lotus for 12 well-recorded ornamental traits. Given a moderate linkage disequilibrium level of 32.9 kb, we successfully identified 149 candidate genes responsible for seven flower traits and plant size variations, including many pleiotropic genes affecting multiple floral-organ-related traits, such as NnKUP2. Notably, we found a 2.75-kb presence-and-absence genomic fragment significantly associated with stamen petaloid and petal number variations, which was further confirmed by re-examining another independent population dataset with petal number records. Intriguingly, this fragment carries MITE transposons bound by siRNAs and is related to the expression differentiation of a nearby candidate gene between few-petalled and double-petalled lotuses. Overall, these genetic variations and candidate genes responsible for diverse lotus traits revealed by our GWAS highlight the role of transposon variations, particularly MITEs, in shaping floral trait diversity.

Identification of QTLs and a putative candidate gene involved in rhizome enlargement of Asian lotus (Nelumbo nucifera)

QTL mapping studies identified three reliable QTLs of rhizome enlargement in lotus. NnBEL6 located within the confidence interval of the major QTL cqREI-LG2 is a key candidate gene enhancing rhizome enlargement. Lotus (Nelumbo) is perennial aquatic plant with nutritional, pharmacological, and ornamental significance. Rhizome is an underground lotus stem that acts as a storage organ and as a reproductive tissue for asexual production. The enlargement of lotus rhizome is an important adaptive strategy for surviving the cold winter. The aims of this study were to identify quantitative trait loci (QTLs) for rhizome enlargement traits including rhizome enlargement index (REI) and number of enlarged rhizome (NER), and to uncover their associated candidate genes. A high-density genetic linkage map was constructed, consisting of 2935 markers binned from 236,840 SNPs. A total of 14 significant QTLs were detected for REI and NER, which explained 6.7-22.3% of trait variance. Three QTL regions were repeatedly identified in at least 2 years, and a major QTL, designated cqREI-LG2, with a rhizome-enlargement effect and about 20% of the phenotypic contribution was identified across the 3 climatic years. A candidate NnBEL6 gene located within the confidence interval of cqREI-LG2 was considered to be putatively involved in lotus rhizome enlargement. The expression of NnBEL6 was exclusively induced by rhizome swelling. Sequence comparison of NnBEL6 among lotus cultivars revealed a functional Indel site in its promoter that likely initiates the rhizome enlargement process. Transgenic potato assay was used to confirm the role of NnBEL6 in inducing tuberization. The successful identification QTLs and functional validation of NnBEL6 gene reported in this study will enrich our knowledge on the genetic basis of rhizome enlargement in lotus.

Genome-Wide Identification and Expression Analysis of Dof Transcription Factors in Lotus (Nelumbo nucifera Gaertn.)

Lotus (Nelumbo nucifera Gaertn.) is a traditional Chinese aquatic flower with high ornamental and economic value, but water salinity seriously affects lotus cultivation and distribution. The Dof transcription factors (TFs) play a crucial function in the regulatory network of growth and defense in plants. However, no systematic investigations of the Dof TFs in lotus have been performed. In this study, comprehensive searches of the lotus genome yielded 29 potential NnDofs. We carried out a series of standardized analyses, which include physical properties, multiple sequence alignment, phylogenetic analysis, gene structure, motif composition, cis-acting element prediction, chromosome distribution, and synteny analysis. The results showed that segment duplication probably caused the NnDofs gene family expansion. The potential functions of NnDofs in lotus development and stress conditions are speculated by promoter analysis. Furthermore, a complete expression investigation of NnDofs utilizing an RNA-seq atlas and quantitative real-time polymerase chain reaction (qRT-PCR) was performed. The majority of the NnDofs exhibit tissue-specific expression patterns, and many genes have been identified as being extremely sensitive to salt stressors. Overall, this study is the first to report a genome-wide assessment of the Dof family in lotus, and the findings offer vital insights for prospective functional studies on lotus salinity stress.

Plateau zokors (Eospalax baileyi) respond to secondary metabolites from the roots of Stellera chamaejasme by enhancing hepatic inflammatory factors and metabolic pathway genes

Herbivores rarely consume toxic plants. An increase in the proportion of toxic plant secondary metabolites (PSMs) in poisonous plants can promote detoxification and related metabolic capacity of animals. Poisonous plants with thick taproots like Stellera chamaejasme (SC) are important stored food for the plateau zokor (Eospalax baileyi) during the winter and promote the development of detoxification mechanisms in this animal. In this study, plateau zokors were administered gavages of 0.2, 1.05, and 2.10 ml/kg SC water extracts. Serum samples were collected from plateau zokors to measure the levels of transaminases and oxidative stress. Transcriptome analysis was conducted to evaluate the differential genes of multiple metabolic pathways to investigate the relationship between the physiological processes and metabolic adaptation capacity of these animals in response to SC. After SC administration, plateau zokors showed significant hepatic granular degeneration and inflammatory reactions in the liver and aspartate aminotransferase, alanine aminotransferase, and malondialdehyde levels increased in a dose-dependent manner. Further, differential expression was also found in the plateau zokor livers, with most enrichment in inflammation and detoxification metabolism pathways. The metabolic adaptation responses in P450 xenobiotic clearance, bile secretion, and pancreatic secretion (Gusb, Hmgcr, Gstm1, Gstp1, and Eobag004630005095) were verified by mRNA network analysis as key factors related to the mechanism. Plateau zokors respond to SC PSMs through changes in liver physiology, biochemistry, and genes in multiple metabolic pathways, validating our hypothesis that plateau zokors can metabolize PSMs when they ingest toxic plants.

Integrated mRNA and Small RNA Sequencing Reveals a microRNA Regulatory Network Associated with Starch Biosynthesis in Lotus (Nelumbo nucifera Gaertn.) Rhizomes

Internode starch biosynthesis is one of the most important traits in lotus rhizome because of its relation to crop productivity. Understanding the microRNA (miRNA) and mRNA expression profiles related to lotus internode starch biosynthesis would help develop molecular improvement strategies, but they are not yet well-investigated. To identify genes and miRNAs involved in internode starch biosynthesis, the cDNA and small RNA libraries of Z6-1, Z6-2, and Z6-3 were sequenced, and their expression were further studied. Through combined analyses of transcriptome data and small RNA sequencing data, a complex co-expression regulatory network was constructed, in which 20 miRNAs could modulate starch biosynthesis in different internodes by tuning the expression of 10 target genes. QRT-PCR analysis, transient co-expression experiment and dual luciferase assay comprehensively confirmed that NnumiR396a down-regulated the expression of NnSS2 and ultimately prevents the synthesis of amylopectin, and NnumiR396b down-regulated the expression of NnPGM2 and ultimately prevents the synthesis of total starch. Our results suggest that miRNAs play a critical role in starch biosynthesis in lotus rhizome, and that miRNA-mediated networks could modulate starch biosynthesis in this tissue. These results have provided important insights into the molecular mechanism of starch biosynthesis in developing lotus rhizome.

Studies on Lotus Genomics and the Contribution to Its Breeding

Lotus (Nelumbo nucifera), under the Nelumbonaceae family, is one of the relict plants possessing important scientific research and economic values. Because of this, much attention has been paid to this species on both its biology and breeding among the scientific community. In the last decade, the genome of lotus has been sequenced, and several high-quality genome assemblies are available, which have significantly facilitated functional genomics studies in lotus. Meanwhile, re-sequencing of the natural and genetic populations along with different levels of omics studies have not only helped to classify the germplasm resources but also to identify the domestication of selected regions and genes controlling different horticultural traits. This review summarizes the latest progress of all these studies on lotus and discusses their potential application in lotus breeding.

Genome-Wide Characterization and Comprehensive Analysis of NAC Transcription Factor Family in Nelumbo nucifera

NAC (NAM, ATAF, and CUC) is a ubiquitously expressed plant-specific transcription factor (TF) family which is involved in the regulation of various biological processes. However, a systematic characterization of NAC gene family is yet to be reported in lotus. Here, 82 NnNAC genes which included five predicted membrane-bound NAC proteins were identified in the lotus genome. Phylogenetic analysis revealed seven-subfamily clusters (I–VII) of NnNAC proteins, with homologous gene pairs displaying similar conserved motifs and gene structure characteristics. Transactivation assay of NnNAC proteins revealed an extensive transcriptional activation capacity which is mediated by the highly divergent C-terminal activation domain (AD). Expression analysis of NnNAC genes in lotus tissues showed high transcript levels in root, stamen, petal and seed coat. In addition, 30 and 29 differentially expressed NnNAC candidate genes putatively involved in lotus seed development and response to complete submergence stress, respectively, were identified. Overall, our study provides potentially useful candidate gene resources for future molecular breeding of lotus varieties with novel agronomic traits.

Transcriptomic analysis of tuberous root in two sweet potato varieties reveals the important genes and regulatory pathways in tuberous root development

Background

Tuberous root formation and development is a complex process in sweet potato, which is regulated by multiple genes and environmental factors. However, the regulatory mechanism of tuberous root development is unclear.

Results

In this study, the transcriptome of fibrous roots (R0) and tuberous roots in three developmental stages (Rl, R2, R3) were analyzed in two sweet potato varieties, GJS-8 and XGH. A total of 22,914 and 24,446 differentially expressed genes (DEGs) were identified in GJS-8 and XGH respectively, 15,920 differential genes were shared by GJS-8 and XGH. KEGG pathway enrichment analysis showed that the DEGs shared by GJS-8 and XGH were mainly involved in “plant hormone signal transduction” “starch and sucrose metabolism” and “MAPK signal transduction”. Trihelix transcription factor (Tai6.25300) was found to be closely related to tuberous root enlargement by the comprehensive analysis of these DEGs and weighted gene co-expression network analysis (WGCNA).

Conclusion

A hypothetical model of genetic regulatory network for tuberous root development of sweet potato is proposed, which emphasizes that some specific signal transduction pathways like “plant hormone signal transduction” “Ca²⁺signal” “MAPK signal transduction” and metabolic processes including “starch and sucrose metabolism” and “cell cycle and cell wall metabolism” are related to tuberous root development in sweet potato. These results provide new insights into the molecular mechanism of tuberous root development in sweet potato.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08670-x.

Comparative analyses of American and Asian lotus genomes reveal insights into petal color, carpel thermogenesis and domestication

Nelumbo lutea (American lotus), which differs from Nelumbo nucifera (Asian lotus) morphologically, is one of the two remaining species in the basal eudicot family Nelumbonaceae. Here, we assembled the 843-Mb genome of American lotus into eight pseudochromosomes containing 31 382 protein-coding genes. Comparative analyses revealed conserved synteny without large chromosomal rearrangements between the genomes of American and Asian lotus and identified 29 533 structural variants (SVs). Carotenoid and anthocyanin pigments determine the yellow and red petal colors of American and Asian lotus, respectively. The structural genes encoding enzymes of the carotenoid and anthocyanin biosynthesis pathways were conserved between two species but differed in expression. We detected SVs caused by repetitive sequence expansion or contraction among the anthocyanin biosynthesis regulatory MYB genes. Further transient overexpression of candidate NnMYB5 induced anthocyanin accumulation in lotus petals. Alternative oxidase (AOX), uncoupling proteins (UCPs), and sugar metabolism and transportation contributed to carpel thermogenesis. Carpels produce heat with sugars transported from leaves as the main substrates, because there was weak tonoplast sugar transporter (TST) activity, and with SWEETs were highly expressed during thermogenesis. Cell proliferation-related activities were particularly enhanced in the warmer carpels compared with stamens during the cold night before blooming, which suggested that thermogenesis plays an important role in flower protogyny. Population genomic analyses revealed deep divergence between American and Asian lotus, and independent domestication affecting seed, rhizome, and flower traits. Our findings provide a high-quality reference genome of American lotus for exploring the genetic divergence and variation between two species and revealed possible genomic bases for petal color, carpel thermogenesis and domestication in lotus.

Expression rewiring and methylation of non-coding RNAs involved in rhizome phenotypic variations of lotus ecotypes

Non-coding RNAs (ncRNAs), including miRNAs, lncRNAs, and circRNAs, emerge as crucial components for gene regulation. Nelumbo nucifera (lotus), a horticulturally important plant, differentiates into a temperate ecotype of enlarged rhizomes and a tropical ecotype of thin rhizomes. Nevertheless, whether and how ncRNAs can be rewired in expression and differentially methylated contributing to adaptive divergence of this storage organ in lotus ecotypes is unclear. Herein, we study the expression behaviors and DNA methylation patterns of ncRNAs in temperate and tropical lotus rhizomes. By whole transcriptome sequencing, we found both mRNAs and lncRNAs have divergent expression patterns between ecotypes, whereas miRNAs and circRNAs tended to be accession-specific or noisier in expression. The differentially expressed ncRNAs are involved in phenotypic differentiation of lotus rhizome between ecotypes, as the genes that interacted with them in the competing endogenous RNA network are enriched in functions including carbohydrate metabolism and plant hormone signaling, being critical to rhizome enlargement. Intriguingly, ncRNA-targeted genes are less prone to show positive selection or differential expression during ecotypic divergence due to constraints from ncRNA-mRNA interactions. The methylation levels of ncRNAs generally tend to be higher in temperate lotus than in tropical lotus, and differential methylation of lncRNAs also tends to have expression changes. Overall, our study of ncRNAs and their targets highlights the role of ncRNAs in rhizome growth variation between lotus ecotypes through expression rewiring and methylation modification.

Transcriptome-Wide Characterization of Alkaloids and Chlorophyll Biosynthesis in Lotus Plumule

Lotus plumule is a green tissue in the middle of seeds that predominantly accumulates bisbenzylisoquinoline alkaloids (bis-BIAs) and chlorophyll (Chl). However, the biosynthetic mechanisms of these two metabolites remain largely unknown in lotus. This study used physiological and RNA sequencing (RNA-Seq) approaches to characterize the development and molecular mechanisms of bis-BIAs and Chl biosynthesis in lotus plumule. Physiological analysis revealed that exponential plumule growth occurred between 9 and 15 days after pollination (DAP), which coincided with the onset of bis-BIAs biosynthesis and its subsequent rapid accumulation. Transcriptome analysis of lotus plumule identified a total of 8,725 differentially expressed genes (DEGs), representing ~27.7% of all transcripts in the lotus genome. Sixteen structural DEGs, potentially associated with bis-BIAs biosynthesis, were identified. Of these, 12 encoded O-methyltransferases (OMTs) are likely involved in the methylation and bis-BIAs diversity in lotus. In addition, functionally divergent paralogous and redundant homologous gene members of the BIAs biosynthesis pathway, as well as transcription factors co-expressed with bis-BIAs and Chl biosynthesis genes, were identified. Twenty-two genes encoding 16 conserved enzymes of the Chl biosynthesis pathway were identified, with the majority being significantly upregulated by Chl biosynthesis. Photosynthesis and Chl biosynthesis pathways were simultaneously activated during lotus plumule development. Moreover, our results showed that light-driven Pchlide reduction is essential for Chl biosynthesis in the lotus plumule. These results will be useful for enhancing our understanding of alkaloids and Chl biosynthesis in plants.

Roles of Abscisic Acid and Gibberellins in Stem/Root Tuber Development

Root and tuber crops are of great importance. They not only contribute to feeding the population but also provide raw material for medicine and small-scale industries. The yield of the root and tuber crops is subject to the development of stem/root tubers, which involves the initiation, expansion, and maturation of storage organs. The formation of the storage organ is a highly intricate process, regulated by multiple phytohormones. Gibberellins (GAs) and abscisic acid (ABA), as antagonists, are essential regulators during stem/root tuber development. This review summarizes the current knowledge of the roles of GA and ABA during stem/root tuber development in various tuber crops.

De novo and comparative transcriptomic analysis explain morphological differences in Panax notoginseng taproots

Background

Panax notoginseng (Burk.) F. H. Chen (PN) belonging to the genus Panax of family Araliaceae is widely used in traditional Chinese medicine to treat various diseases. PN taproot, as the most vital organ for the accumulation of bioactive components, presents a variable morphology (oval or long), even within the same environment. However, no related studies have yet explained the molecular mechanism of phenotypic differences. To investigate the cause of differences in the taproot phenotype, de novo and comparative transcriptomic analysis on PN taproot was performed.

Results

A total of 133,730,886 and 114,761,595 paired-end clean reads were obtained based on high-throughput sequencing from oval and long taproot samples, respectively. 121,955 unigenes with contig N50 = 1,774 bp were generated by using the de novo assembly transcriptome, 63,133 annotations were obtained with the BLAST. And then, 42 genes belong to class III peroxidase (PRX) gene family, 8 genes belong to L-Ascorbate peroxidase (APX) gene family, and 55 genes belong to a series of mitogen-activated protein kinase (MAPK) gene family were identified based on integrated annotation results. Differentially expressed genes analysis indicated substantial up-regulation of PnAPX3 and PnPRX45, which are related to reactive oxygen species metabolism, and the PnMPK3 gene, which is related to cell proliferation and plant root development, in long taproots compared with that in oval taproots. Furthermore, the determination results of real-time quantitative PCR, enzyme activity, and H₂O₂ content verified transcriptomic analysis results.

Conclusion

These results collectively demonstrate that reactive oxygen species (ROS) metabolism and the PnMPK3 gene may play vital roles in regulating the taproot phenotype of PN. This study provides further insights into the genetic mechanisms of phenotypic differences in other species of the genus Panax.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08283-w.

RhWRKY33 Positively Regulates Onset of Floral Senescence by Responding to Wounding- and Ethylene-Signaling in Rose Plants

Rose plants are one of the most important horticultural crops, whose commercial value mainly depends on long-distance transportation, and wounding and ethylene are the main factors leading to their quality decline and accelerated senescence in the process. However, underlying molecular mechanisms of crosstalk between wounding and ethylene in the regulation of flower senescence remain poorly understood. In relation to this, transcriptome analysis was performed on rose flowers subjected to various treatments, including control, wounding, ethylene, and wounding- and ethylene- (EW) dual treatment. A large number of differentially expressed genes (DEGs) were identified, ranging from 2,442 between the ethylene- and control-treated groups to 4,055 between the EW- and control-treated groups. Using weighted gene co-expression network analysis (WGCNA), we identified a hub gene RhWRKY33 (rchiobhmchr5g0071811), accumulated in the nucleus, where it may function as a transcription factor. Moreover, quantitative reverse transcription PCR (RT-qPCR) results showed that the expression of RhWRKY33 was higher in the wounding-, ethylene, and EW-treated petals than in the control-treated petals. We also functionally characterized the RhWRKY33 gene through virus-induced gene silencing (VIGS). The silencing of RhWRKY33 significantly delayed the senescence process in the different treatments (control, wounding, ethylene, and EW). Meanwhile, we found that the effect of RhWRKY33-silenced petals under ethylene and EW dual-treatment were stronger than those under wounding treatment in delaying the petal senescence process, implying that RhWRKY33 is closely involved with ethylene and wounding mediated petal senescence. Overall, the results indicate that RhWRKY33 positively regulates the onset of floral senescence mediated by both ethylene and wounding signaling, but relies heavily on ethylene signaling.

Molecular cloning, characterization and expression analysis of three key starch synthesis-related genes from the bulb of a rare lily germplasm, Lilium brownii var. giganteum

Starch is the predominant compound in bulb scales, and previous studies have shown that bulblet development is closely associated with starch enrichment. However, how starch synthesis affects bulbification at the molecular level is unclear. In this study, we demonstrate that Lilium brownii var. giganteum, a wild lily with a giant bulb in nature, and L. brownii, the native species, have different starch levels and characteristics according to cytological and ultra-structural observations. We cloned the complete sequence of three key gene-encoding enzymes (LbgAGPS, LbgGBSS, andLbgSSIII) during starch synthesis by rapid amplification of 5' and 3' complementary DNA (cDNA) ends (RACE) technology. Bioinformatics analysis revealed that the proteins deduced by these genes contain the canonical conserved domains. Constructed phylogenetic trees confirmed the evolutionary relationships with proteins from other species, including monocotyledons and dicotyledons. The transcript levels of various tissues and time course samples obtained during bulblet development uncovered relatively high expression levels in bulblets and gradual increase expression accompanying bulblet growth. Moreover, a set of single nucleotide polymorphisms (SNPs) was discovered in the AGPS genes of four lily genotypes, and a purifying selection fashion was predicted according to the non-synonymous/synonymous (Ka/Ks) values. Taken together, our results suggested that key starch-synthesizing genes might play important roles in bulblet development and lead to distinctive phenotypes in bulblet size.

Distinct methylome patterns contribute to ecotypic differentiation in the growth of the storage organ of a flowering plant (sacred lotus)

DNA methylation is an epigenetic modification involved in phenotypic diversity, plant development, and environmental responses. However, the mechanisms of DNA methylation underpinning the adaption of lotus (Nelumbo nucifera) ecotypes to high and low latitudes remain unsolved, especially adaptive evolution of their storage organs. Tropical and temperate lotus ecotypes have thin and enlarged rhizomes which are adapted to low and high latitudes, respectively. Here, we investigated the DNA methylomes and transcriptomes of rhizomes of the temperate and tropical lotus to address this issue. Compared with that of the tropical lotus, the DNA of the temperate lotus was significantly more hypermethylated, indicating an increase in global DNA methylation in the lotus, with rhizome enlargement. Meanwhile, genes associated with differentially methylated regions in their promoters tended to be differentially expressed between the two ecotypes. Interestingly, the genes with their expression negatively correlated with methylation levels in their promoters and genomic regions displayed significantly higher transposon coverage, while the genes showing a significant positive correlation between expression and methylation showed lesser transposon coverage. Further, we identified that DNA methylation, especially in the promoter region, was significantly correlated with the expression of many starch-biosynthetic, gibberellin-, and brassinosteroid-signalling genes associated with rhizome differentiation. Overall, our study unveiled that distinct global and local methylation patterns between the two lotus ecotypes contribute to their expression differences and adaptive phenotypic divergence of their storage organs, highlighting the role of DNA methylation in shaping the ecotypic differentiation of lotus.

Biased allelic expression in tissues of F1 hybrids between tropical and temperate lotus (Nelumbo nuicfera)

The genome-wide allele-specific expression in F1 hybrids from the cross of tropical and temperate lotus unveils how cis-regulatory divergences affect genes in key pathways related to ecotypic divergence. Genetic variation, particularly cis-regulatory variation, plays a crucial role in phenotypic variation and adaptive evolution in plants. Temperate and tropical lotus, the two ecotypes of Nelumbo nucifera, show distinction in the degree of rhizome enlargement, which is associated with winter dormancy. To understand the roles of genome-wide cis-regulatory divergences on adaptive evolution of temperate and tropical lotus (Nelumbo nucifera), here we performed allele-specific expression (ASE) analyses on the tissues including flowers, leaves and rhizome from F1 hybrids of tropical and temperate lotus. For all investigated tissues in F1s, about 36% of genes showed ASE and about 3% of genes showed strong consistent ASE. Most of ASEs were biased towards the tropical parent in all surveyed samples, indicating that the tropical genome might be dominant over the temperate genome in gene expression of tissues from their F1 hybrids. We found that promoter sequences with similar allelic expression are more conserved than genes with significant or conditional ASE, suggesting the cis-regulatory sequence divergence underlie the allelic expression bias. We further uncovered biased genes being related to phenotypic differentiation between two lotus ecotypes, especially metabolic and phytohormone-related pathways in the rhizome. Overall, our study provides a global landscape of cis-regulatory variations between two lotus ecotypes and highlights their roles in rhizome growth variation for the climatic adaptation.

Proteome analysis of bermudagrass stolons and rhizomes provides new insights into the adaptation of plant stems to aboveground and underground growth

As an important perennial warm-season turfgrass species, bermudagrass (Cynodon dactylon L.) forms underground-growing rhizomes and aboveground-growing stolons simultaneously, making it a fast propagating clonal plant with strong regeneration ability. In the current study, we compared the internode proteomes of rhizomes and stolons at the same developmental stage in the bermudagrass cultivar Yangjiang using iTRAQ. The results indicated that 228 protein species were differentially accumulated in the two specialized stems. In agreement with the different contents of starch, chlorophyll, anthocyanin and H₂O₂ in the two types of stems, photosynthesis and flavonoid biosynthesis were enriched with differentially accumulated protein species (DAPs) in stolons, whereas starch and sucrose metabolism, glycolysis, and H₂O₂ metabolism were enriched with DAPs in rhizomes. Burying stolons in the soil resulted in the gradual degradation of chlorophyll and anthocyanin, accumulation of starch, and increment of H₂O₂, which is similar to the physiological characteristics of rhizomes. These results collectively revealed that stolons and rhizomes of bermudagrass have significant differences at the proteome level and light might play important regulatory roles in the discrepancy of the proteome profiles and specialization of the two stems, providing new insights into the adaptation of plant stems to aboveground and underground growth. BIOLOGICAL SIGNIFICANCE: As two types of specialized stems that grow underground and aboveground respectively, rhizomes and stolons play important roles in overwintering and ecological invasion of many perennial and clonal plants. However, because rhizomes and stolons rarely coexist in single plant species, the differences between the two stems remain unclear at the molecular level. In this study, through an iTRAQ comparative proteomic analysis, we reported the identification of 228 differentially accumulated protein species (DAPs) in rhizomes and stolons of bermudagrass for the first time. We found that the 228 DAPs were interconnected to form protein networks in regulating diverse cellular activities and biochemical reactions. We also observed that stolons growing underground showed similar physiological activities and DAP expression as those of underground-growing rhizomes, suggesting that light might play important regulatory roles in the specialization of stolons and rhizomes. These results expanded our understanding of the mysterious adaption of plant stems to different growth conditions.

Genome-wide analysis of AP2/ERF superfamily in lotus (Nelumbo nucifera) and the association between NnADAP and rhizome morphology

Background

The AP2/ERF family is widely present in plants and plays a crucial regulatory role in plant growth and development. As an essential aquatic horticultural model plant, lotus has an increasingly prominent economic and research value.

Results

We have identified and analysed the AP2/ERF gene family in the lotus. Initially, 121 AP2/ERF family genes were identified. By analysing their gene distribution and protein structure, and their expression patterns during the development of lotus rhizome, combined with previous studies, we obtained an SNP (megascaffold_20:3578539) associated with lotus rhizome phenotype. This SNP was in the NnADAP gene of the AP2 subfamily, and the changes in SNP (C/T) caused amino acid conversion (proline/leucine). We constructed a population of 95 lotus varieties for SNP verification. Through population typing experiments, we found that the group with SNP CC had significantly larger lotus rhizome and higher soluble sugar content among the population.

Conclusions

In conclusion, we speculate that the alteration of the SNP in the NnADAP can affect the size and sugar content of the lotus rhizome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07473-w.

Genetic dissection of rhizome yield-related traits in Nelumbo nucifera through genetic linkage map construction and QTL mapping

Lotus (Nelumbo nucifera) is a perennial aquatic plant with great value in ornamentation, nutrition, and medicine. Being a storage organ, lotus rhizome is not only used for vegetative reproduction, but also as a popular vegetable in Southeast Asia. Rhizome development, especially enlargement, largely determines its yield and hence becomes one of the major concerns in rhizome lotus breeding and cultivation. To obtain the genetic characteristic of this trait, and discover markers or genes associated with this trait, an F₂ population was generated by crossing between temperate and tropical cultivars with contrasting rhizome enlargement. Based on this F₂ population and Genotyping-by-Sequencing (GBS) technique, a genetic map was constructed with 1475 bin markers containing 12,113 SNP markers. Six traits associated with rhizome yield were observed over 3 years. Quantitative trait locus (QTL) mapping analysis identified 22 QTLs that are associated with at least one of these traits, among which 9 were linked with 3 different intervals. Comparison of the genes located in these three intervals with our previous transcriptomic data showed that light and phytohormone signaling might contribute to the development and enlargement of lotus rhizome. The QTLs obtained here could also be used for marker-assisted breeding of rhizome lotus.

A CONSTANS-LIKE gene of Nelumbo nucifera could promote potato tuberization

CONSTANS-LIKE 5 of Nelumbo nucifera is capable of promoting potato tuberization through CONSTANS-FLOWERING LOCUS T and gibberellin signaling pathways with a probable association with lotus rhizome enlargement. Lotus (Nelumbo nucifera) is an aquatic plant that is affiliated to the Nelumbonaceace family. It is widely used as an ornamental, vegetable, and medicinal herb with its rhizome being a popular vegetable. To explore the molecular mechanism underlying its rhizome enlargement, we conducted a systematic analysis on the CONSTANS-LIKE (COL) gene family, with the results, indicating that this gene plays a role in regulating potato tuber expansion. These analyses included phylogenetic relationships, gene structure, and expressional patterns of lotus COL family genes. Based on these analyses, NnCOL5 was selected for further study on its potential function in lotus rhizome formation. NnCOL5 was shown to be located in the nucleus, and its expression was positively associated with the enlargement of lotus rhizome. Besides, the overexpression of NnCOL5 in potato led to increased tuber weight and starch content under short-day conditions without changing the number of tubers. Further analysis suggested that the observed tuber changes might be mediated by affecting the expression of genes in CO-FT and GA signaling pathways. These results provide valuable insight in understanding the functions of COL gene as well as the enlargement of lotus rhizome.

Developmental regulation of stolon and rhizome

Stolons and rhizomes are modified stems for vegetative reproduction. While stolons grow above the ground, rhizomes grow beneath the ground. Stolons and rhizomes maintain the genotypes of hybrids and hence are invaluable for agricultural propagation. Diploid strawberry is a model for studying stolon development. At the axillary meristems, gibberellins and MADS box gene SOC1 promote stolon formation, while the DELLA repressor inhibits stolon development. Photoperiod regulates stolon formation through regulating GA biosynthesis or balancing asexual with sexual mode of reproduction in the axillary meristems. In rhizomatous wild rice, the BLADE-ON-PETIOLE gene promotes sheath-to-blade ratio to confer rhizome tip stiffness and support underground growth. Together, this review aims to encourage further investigations into stolon and rhizome to benefit agriculture and environment.

New insight into the rapid growth of the Mikania micrantha stem based on DIA proteomic and RNA-Seq analysis

Mikania micrantha is one of the world's most invasive plants, which causes severe damage to natural ecosystems and agroforestry systems due to its rapid stem growth. This work investigated the proteomic and transcriptomic profiles of M. micrantha in different stem tissues (pre-internode, post-internode, and internode), as well as in adventitious roots and primary roots with the final goal of elucidating differentially expressed genes and proteins responsible for the rapid growth of stem. The objective was approached by using DIA-based proteomic and RNA-Seq technologies. More than seven giga-transcriptome clean reads were sequenced, and 5196 protein species were identified. Differentially expressed genes identified in all stem tissues were significantly enriched in photosynthesis and carbon fixation, suggesting that the stem possesses a strong photosynthetic capacity in order to maintain the energy supply for this species. Analysis of differentially expressed proteins showed that proteins related to photosystem I/II and the cytochrome b6/f complex, such as D1, D2, and cp43, were also highly accumulated in the adventitious roots, corroborating the transcriptome analysis results. These results provided basic proteomic and transcriptional expression information about the M. micrantha stem and adventitious root, thereby improving our understanding of the molecular mechanism underlying rapid growth in this species. SIGNIFICANCE: This is the first study to investigate the proteomic and transcriptomic profiles of Mikania micrantha, a highly invasive plant, in different stem tissues (pre-internode, post-internode, and internode), as well as in adventitious and primary roots, using the latest DIA-based (data-independent acquisition mode) proteomic and RNA-Seq technologies. A comprehensive study was carried out, and differentially expressed genes and differentially expressed proteins identified in the pre-internode, post-internode, and internode tissues were significantly enriched during photosynthesis and carbon fixation, suggesting that the M. micrantha stem possesses a strong photosynthetic capacity that allows the plant to maintain a high energy supply. Enriched plant hormone signal transduction pathway analysis revealed an interaction between auxin and other phytohormones involved in adventitious root development. The study provided basic data on the molecular mechanism of M. micrantha vegetative propagation and the rapid growth of its stem. The novel scientific content of this study successfully builds upon the limited information currently available on the subject, therefore warranting publication.

Nelumbo genome database, an integrative resource for gene expression and variants of Nelumbo nucifera

Sacred lotus (Nelumbo nucifera, or lotus) is one of the most widely grown aquatic plant species with important uses, such as in water gardening and in vegetable and herbal medicine. A public genomic database of lotus would facilitate studies of lotus and other aquatic plant species. Here, we constructed an integrative database: the Nelumbo Genome Database (NGD, http://nelumbo.biocloud.net). This database is a collection of the most updated lotus genome assembly and contains information on both gene expression in different tissues and coexpression networks. In the NGD, we also integrated genetic variants and key traits from our 62 newly sequenced lotus cultivars and 26 previously reported cultivars, which are valuable for lotus germplasm studies. As applications including BLAST, BLAT, Primer, Annotation Search, Variant and Trait Search are deployed, users can perform sequence analyses and gene searches via the NGD. Overall, the valuable genomic resources provided in the NGD will facilitate future studies on population genetics and molecular breeding of lotus.

Measurement(s)	reference genome data • whole genome sequencing • transcriptome
Technology Type(s)	Hi-C • PacBio Sequel System • Illumina sequencing • RNA sequencing • DNA sequencing
Sample Characteristic - Organism	Nelumbo nucifera

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.13487271

The Establishment of an Efficient Callus Induction System for Lotus (Nelumbo nucifera)

The lotus (Nelumbo nucifera) is one of the most popular aquatic plants in Asia, and has emerged as a novel model for studying flower and rhizome development, and primary and secondary metabolite accumulation. Here, we developed a highly efficient callus induction system for the lotus by optimizing a series of key factors that affect callus formation. The highest efficient callus production was induced on immature cotyledon and embryo explants grown on Murashige and Skoog (MS) basal medium containing an optimized combination of 3 mg/L 2,4-dichlorophenoxyacetic acid (2,4-D) and 0.5 mg/L 6-benzylaminopurine (6-BA). In addition, lotus callus induction was proven to be influenced by lotus genotypes, light conditions, the developmental stages of explants and the time of explant sampling. Collecting immature cotyledons from seeds of the genotype “Shilihe 1”, at 9 days post pollination, and to culture the explants in darkness, are proposed as the optimum conditions for lotus callus induction. Interestingly, highly efficient callus induction was also observed in explants of immature embryo derived aseptic seedlings; and a small amount of lotus benzylisoquinoline alkaloid (BIA) and obvious expression of BIA biosynthetic genes were detected in lotus callus.

Comprehensive analysis of AGPase genes uncovers their potential roles in starch biosynthesis in lotus seed

BACKGROUND

Starch in the lotus seed contains a high proportion of amylose, which endows lotus seed a promising property in the development of hypoglycemic and low-glycemic index functional food. Currently, improving starch content is one of the major goals for seed-lotus breeding. ADP-glucose pyrophosphorylase (AGPase) plays an essential role in regulating starch biosynthesis in plants, but little is known about its characterization in lotus.

RESULTS

We describe the nutritional compositions of lotus seed among 30 varieties with starch as a major component. Comparative transcriptome analysis showed that AGPase genes were differentially expressed in two varieties (CA and JX) with significant different starch content. Seven putative AGPase genes were identified in the lotus genome (Nelumbo nucifera Gaertn.), which could be grouped into two subfamilies. Selective pressure analysis indicated that purifying selection acted as a vital force in the evolution of AGPase genes. Expression analysis revealed that lotus AGPase genes have varying expression patterns, with NnAGPL2a and NnAGPS1a as the most predominantly expressed, especially in seed and rhizome. NnAGPL2a and NnAGPS1a were co-expressed with a number of starch and sucrose metabolism pathway related genes, and their expressions were accompanied by increased AGPase activity and starch content in lotus seed.

CONCLUSIONS

Seven AGPase genes were characterized in lotus, with NnAGPL2a and NnAGPS1a, as the key genes involved in starch biosynthesis in lotus seed. These results considerably extend our understanding on lotus AGPase genes and provide theoretical basis for breeding new lotus varieties with high-starch content.

Revealing the full-length transcriptome of caucasian clover rhizome development

BACKGROUND

Caucasian clover (Trifolium ambiguum M. Bieb.) is a strongly rhizomatous, low-crowned perennial leguminous and ground-covering grass. The species may be used as an ornamental plant and is resistant to cold, arid temperatures and grazing due to a well-developed underground rhizome system and a strong clonal reproduction capacity. However, the posttranscriptional mechanism of the development of the rhizome system in caucasian clover has not been comprehensively studied. Additionally, a reference genome for this species has not yet been published, which limits further exploration of many important biological processes in this plant.

RESULT

We adopted PacBio sequencing and Illumina sequencing to identify differentially expressed genes (DEGs) in five tissues, including taproot (T1), horizontal rhizome (T2), swelling of taproot (T3), rhizome bud (T4) and rhizome bud tip (T5) tissues, in the caucasian clover rhizome. In total, we obtained 19.82 GB clean data and 80,654 nonredundant transcripts were analysed. Additionally, we identified 78,209 open reading frames (ORFs), 65,227 coding sequences (CDSs), 58,276 simple sequence repeats (SSRs), 6821 alternative splicing (AS) events, 2429 long noncoding RNAs (lncRNAs) and 4501 putative transcription factors (TFs) from 64 different families. Compared with other tissues, T5 exhibited more DEGs, and co-upregulated genes in T5 are mainly annotated as involved in phenylpropanoid biosynthesis. We also identified betaine aldehyde dehydrogenase (BADH) as a highly expressed gene-specific to T5. A weighted gene co-expression network analysis (WGCNA) of transcription factors and physiological indicators were combined to reveal 11 hub genes (MEgreen-GA3), three of which belong to the HB-KNOX family, that are up-regulated in T3. We analysed 276 DEGs involved in hormone signalling and transduction, and the largest number of genes are associated with the auxin (IAA) signalling pathway, with significant up-regulation in T2 and T5.

CONCLUSIONS

This study contributes to our understanding of gene expression across five different tissues and provides preliminary insight into rhizome growth and development in caucasian clover.

Characterization of rhizome transcriptome and identification of a rhizomatous ER body in the clonal plant Cardamine leucantha

The rhizome is a plant organ that develops from a shoot apical meristem but penetrates into belowground environments. To characterize the gene expression profile of rhizomes, we compared the rhizome transcriptome with those of the leaves, shoots and roots of a rhizomatous Brassicaceae plant, Cardamine leucantha. Overall, rhizome transcriptomes were characterized by the absence of genes that show rhizome-specific expression and expression profiles intermediate between those of shoots and roots. Our results suggest that both endogenous developmental factors and external environmental factors are important for controlling the rhizome transcriptome. Genes that showed relatively high expression in the rhizome compared to shoots and roots included those related to belowground defense, control of reactive oxygen species and cell elongation under dark conditions. A comparison of transcriptomes further allowed us to identify the presence of an ER body, a defense-related belowground organelle, in epidermal cells of the C. leucantha rhizome, which is the first report of ER bodies in rhizome tissue.

Biochemical and Comparative Transcriptome Analyses Reveal Key Genes Involved in Major Metabolic Regulation Related to Colored Leaf Formation in Osmanthus fragrans 'Yinbi Shuanghui' during Development

Osmanthus fragrans 'Yinbi Shuanghui' not only has a beautiful shape and fresh floral fragrance, but also rich leaf colors that change, making the tree useful for landscaping. In order to study the mechanisms of color formation in O. fragrans 'Yinbi Shuanghui' leaves, we analyzed the colored and green leaves at different developmental stages in terms of leaf pigment content, cell structure, and transcriptome data. We found that the chlorophyll content in the colored leaves was lower than that of green leaves throughout development. By analyzing the structure of chloroplasts, the colored leaves demonstrated more stromal lamellae and low numbers of granum thylakoid. However, there was a large number of plastoglobuli. Using transcriptome sequencing, we demonstrated that the expression of differentially expressed genes (DEGs) involved in chlorophyll degradation was upregulated, i.e., heme oxygennase-1 (HO1), pheophorbide a oxidase (PAO), and chlorophyllase-2 (CLH2), affecting the synthesis of chlorophyll in colored leaves. The stay-green gene (SGR) was upregulated in colored leaves. Genes involved in carotenoid synthesis, i.e., phytoene synthase 1 (PSY1) and 1-Deoxyxylulose-5-phosphate synthase (DXS), were downregulated in colored leaves, impeding the synthesis of carotenoids. In the later stage of leaf development, the downregulated expression of Golden2-Like (GLK) inhibited chloroplast development in colored leaves. Using weighted gene co-expression network analysis (WGCNA) to investigate the correlation between physiological indicators and DEGs, we chose the modules with the highest degree of relevance to chlorophyll degradation and carotenoid metabolism. A total of five genes (HSFA2, NFYC9, TCP20, WRKY3, and WRKY4) were identified as hub genes. These analyses provide new insights into color formation mechanisms in O. fragrans 'Yinbi Shuanghui' leaves at the transcriptional level.

The Control of Developmental Phase Transitions by microRNAs and Their Targets in Seed Plants

Seed plants usually undergo various developmental phase transitions throughout their lifespan, mainly including juvenile-to-adult and vegetative-to-reproductive transitions, as well as developmental transitions within organ/tissue formation. MicroRNAs (miRNAs), as a class of small endogenous non-coding RNAs, are involved in the developmental phase transitions in plants by negatively regulating the expression of their target genes at the post-transcriptional level. In recent years, cumulative evidence has revealed that five miRNAs, miR156, miR159, miR166, miR172, and miR396, are key regulators of developmental phase transitions in plants. In this review, the advanced progress of the five miRNAs and their targets in regulating plant developmental transitions, especially in storage organ formation, are summarized and discussed, combining our own findings with the literature. In general, the functions of the five miRNAs and their targets are relatively conserved, but their functional divergences also emerge to some extent. In addition, potential research directions of miRNAs in regulating plant developmental phase transitions are prospected.

Integrated mRNA and miRNA transcriptome analysis reveals a regulatory network for tuber expansion in Chinese yam (Dioscorea opposita)

BACKGROUND

Yam tuber is a storage organ, derived from the modified stem. Tuber expansion is a complex process, and depends on the expressions of genes that can be influenced by environmental and endogenous factors. However, little is known about the regulatory mechanism of tuber expansion. In order to identify the genes and miRNAs involved in tuber expansion, we examined the mRNAs and small RNAs in Dioscorea opposita (Chinese yam) cv. Guihuai 16 tuber during its initiation and expansion stages.

RESULTS

A total of 14,238 differentially expressed genes in yam tuber at its expansion stage were identified by using RNA sequencing technology. Among them, 5723 genes were up-regulated, and 8515 genes were down-regulated. Functional analysis revealed the coordination of tuber plant involved in processes of cell events, metabolism, biosynthesis, and signal transduction pathways at transcriptional level, suggesting that these differentially expressed genes are somehow involved in response to tuber expansion, including CDPK, CaM, CDL, SAUR, DELLA, SuSy, and expansin. In addition, 541 transcription factor genes showed differential expression during the expansion stage at transcriptional level. MADS, bHLH, and GRAS were involved in cell differentiation, division, and expansion, which may relate to tuber expansion. Noteworthy, data analysis revealed that 22 known tuber miRNAs belong to 10 miRNA families, and 50 novel miRNAs were identified. The integrated analysis of miRNA-mRNA showed that 4 known miRNAs and 11 genes formed 14 miRNA-target mRNA pairs were co-expressed in expansion stage. miRNA160, miRNA396, miRNA535 and miRNA5021 may be involved in complex network to regulate cell division and differentiation in yam during its expansion stage.

CONCLUSION

The mRNA and miRNA datasets presented here identified a subset of candidate genes and miRNAs that are putatively associated with tuber expansion in yam, a hypothetical model of genetic regulatory network associated with tuber expansion in yam was put forward, which may provide a foundation for molecular regulatory mechanism researching on tuber expansion in Dioscorea species.

Integrated mRNA and miRNA transcriptome analysis reveals a regulatory network for tuber expansion in Chinese yam (Dioscorea opposita)

Background

Yam tuber is a storage organ, derived from the modified stem. Tuber expansion is a complex process, and depends on the expressions of genes that can be influenced by environmental and endogenous factors. However, little is known about the regulatory mechanism of tuber expansion. In order to identify the genes and miRNAs involved in tuber expansion, we examined the mRNAs and small RNAs in Dioscorea opposita (Chinese yam) cv. Guihuai 16 tuber during its initiation and expansion stages.

Results

A total of 14,238 differentially expressed genes in yam tuber at its expansion stage were identified by using RNA sequencing technology. Among them, 5723 genes were up-regulated, and 8515 genes were down-regulated. Functional analysis revealed the coordination of tuber plant involved in processes of cell events, metabolism, biosynthesis, and signal transduction pathways at transcriptional level, suggesting that these differentially expressed genes are somehow involved in response to tuber expansion, including CDPK, CaM, CDL, SAUR, DELLA, SuSy, and expansin. In addition, 541 transcription factor genes showed differential expression during the expansion stage at transcriptional level. MADS, bHLH, and GRAS were involved in cell differentiation, division, and expansion, which may relate to tuber expansion. Noteworthy, data analysis revealed that 22 known tuber miRNAs belong to 10 miRNA families, and 50 novel miRNAs were identified. The integrated analysis of miRNA-mRNA showed that 4 known miRNAs and 11 genes formed 14 miRNA-target mRNA pairs were co-expressed in expansion stage. miRNA160, miRNA396, miRNA535 and miRNA5021 may be involved in complex network to regulate cell division and differentiation in yam during its expansion stage.

Conclusion

The mRNA and miRNA datasets presented here identified a subset of candidate genes and miRNAs that are putatively associated with tuber expansion in yam, a hypothetical model of genetic regulatory network associated with tuber expansion in yam was put forward, which may provide a foundation for molecular regulatory mechanism researching on tuber expansion in Dioscorea species.

Towards Improved Molecular Identification Tools in Fine Fescue (Festuca L., Poaceae) Turfgrasses: Nuclear Genome Size, Ploidy, and Chloroplast Genome Sequencing

Fine fescues (Festuca L., Poaceae) are turfgrass species that perform well in low-input environments. Based on morphological characteristics, the most commonly-utilized fine fescues are divided into five taxa: three are subspecies within F. rubra L. and the remaining two are treated as species within the F. ovina L. complex. Morphologically, these five taxa are very similar; both identification and classification of fine fescues remain challenging. In an effort to develop identification methods for fescues, we used flow cytometry to estimate genome size and ploidy level and sequenced the chloroplast genome of all five taxa. Fine fescue chloroplast genome sizes ranged from 133,331 to 133,841 bp and contained 113-114 genes. Phylogenetic relationship reconstruction using whole chloroplast genome sequences agreed with previous work based on morphology. Comparative genomics suggested unique repeat signatures for each fine fescue taxon that could potentially be used for marker development for taxon identification.

Integrated transcriptome and miRNA analysis uncovers molecular regulators of aerial stem-to-rhizome transition in the medical herb Gynostemma pentaphyllum

BACKGROUND

Gynostemma pentaphyllum is an important perennial medicinal herb belonging to the family Cucurbitaceae. Aerial stem-to-rhizome transition before entering the winter is an adaptive regenerative strategy in G. pentaphyllum that enables it to survive during winter. However, the molecular regulation of aerial stem-to-rhizome transition is unknown in plants. Here, integrated transcriptome and miRNA analysis was conducted to investigate the regulatory network of stem-to-rhizome transition.

RESULTS

Nine transcriptome libraries prepared from stem/rhizome samples collected at three stages of developmental stem-to-rhizome transition were sequenced and a total of 5428 differentially expressed genes (DEGs) were identified. DEGs associated with gravitropism, cell wall biosynthesis, photoperiod, hormone signaling, and carbohydrate metabolism were found to regulate stem-to-rhizome transition. Nine small RNA libraries were parallelly sequenced, and seven significantly differentially expressed miRNAs (DEMs) were identified, including four known and three novel miRNAs. The seven DEMs targeted 123 mRNAs, and six pairs of miRNA-target showed significantly opposite expression trends. The GpmiR166b-GpECH2 module involved in stem-to-rhizome transition probably promotes cell expansion by IBA-to-IAA conversion, and the GpmiR166e-GpSGT-like module probably protects IAA from degradation, thereby promoting rhizome formation. GpmiR156a was found to be involved in stem-to-rhizome transition by inhibiting the expression of GpSPL13A/GpSPL6, which are believed to negatively regulate vegetative phase transition. GpmiR156a and a novel miRNA Co.47071 co-repressed the expression of growth inhibitor GpRAV-like during stem-to-rhizome transition. These miRNAs and their targets were first reported to be involved in the formation of rhizomes. In this study, the expression patterns of DEGs, DEMs and their targets were further validated by quantitative real-time PCR, supporting the reliability of sequencing data.

CONCLUSIONS

Our study revealed a comprehensive molecular network regulating the transition of aerial stem to rhizome in G. pentaphyllum. These results broaden our understanding of developmental phase transitions in plants.

Comparative transcriptome and metabolome analyses provide new insights into the molecular mechanisms underlying taproot thickening in Panax notoginseng

BACKGROUND

Taproot thickening is a complex biological process that is dependent on the coordinated expression of genes controlled by both environmental and developmental factors. Panax notoginseng is an important Chinese medicinal herb that is characterized by an enlarged taproot as the main organ of saponin accumulation. However, the molecular mechanisms of taproot enlargement are poorly understood.

RESULTS

A total of 29,957 differentially expressed genes (DEGs) were identified during the thickening process in the taproots of P. notoginseng. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment revealed that DEGs associated with "plant hormone signal transduction," "starch and sucrose metabolism," and "phenylpropanoid biosynthesis" were predominantly enriched. Further analysis identified some critical genes (e.g., RNase-like major storage protein, DA1-related protein, and Starch branching enzyme I) and metabolites (e.g., sucrose, glucose, fructose, malate, and arginine) that potentially control taproot thickening. Several aspects including hormone crosstalk, transcriptional regulation, homeostatic regulation between sugar and starch, and cell wall metabolism, were identified as important for the thickening process in the taproot of P. notoginseng.

CONCLUSION

The results provide a molecular regulatory network of taproot thickening in P. notoginseng and facilitate the further characterization of the genes responsible for taproot formation in root medicinal plants or crops.

Proteomic analysis showing the signaling pathways involved in the rhizome enlargement process in Nelumbo nucifera

Background

Rhizome is the storage underground stem of lotus (Nelumbo nucifera), which is enlarged before winter season and could be used for asexual propagation. In addition, the enlarged rhizome is a nutritional vegetable with abundant starch, proteins, and vitamins. Enlargement of lotus rhizome is not only significance for itself to survive from the cold winter, but also important for its economic value.

Results

To explore the mechanism underlying its enlargement, integrative analyses of morphology, physiology and proteomics were conducted on the rhizome at stolon, middle, and enlarged stages. Morphological observation and physiological analyses showed that rhizomes were gradually enlarged during this process, in which the starch accumulation was also initiated. Quantitative proteomic analysis on the rhizomes at these three stages identified 302 stage-specific proteins (SSPs) and 172 differently expressed proteins (DEPs), based on which GO and KEGG enrichment analyses were conducted. The results indicated that light and auxin signal might be transduced through secondary messenger Ca²⁺, and play important roles in lotus rhizome enlargement.

Conclusion

These results will provide new insights into understanding the mechanism of lotus rhizome enlargement. Meanwhile, some candidate genes might be useful for further studies on this process, as well as breeding of rhizome lotus.

The Latest Studies on Lotus (Nelumbo nucifera)-an Emerging Horticultural Model Plant

Lotus (Nelumbo nucifera) is a perennial aquatic basal eudicot belonging to a small family Nelumbonaceace, which contains only one genus with two species. It is an important horticultural plant, with its uses ranging from ornamental, nutritional to medicinal values, and has been widely used, especially in Southeast Asia. Recently, the lotus obtained a lot of attention from the scientific community. An increasing number of research papers focusing on it have been published, which have shed light on the mysteries of this species. Here, we comprehensively reviewed the latest advancement of studies on the lotus, including phylogeny, genomics and the molecular mechanisms underlying its unique properties, its economic important traits, and so on. Meanwhile, current limitations in the research of the lotus were addressed, and the potential prospective were proposed as well. We believe that the lotus will be an important model plant in horticulture with the generation of germplasm suitable for laboratory operation and the establishment of a regeneration and transformation system.