A genome assembly of ginger (Zingiber officinale Roscoe) provides insights into genome evolution and 6-gingerol biosynthesis.

Ginger is cultivated in tropical and subtropical regions and is one of the most crucial spices worldwide owing to its special taste and scent. Here, we present a high-quality genome assembly for 'Small Laiwu Ginger', a famous cultivated ginger in northern China. The ginger genome was phased into two haplotypes, haplotype A (1.55Gb), and haplotype B (1.44Gb). Analysis of Ty1/Copia and Ty3/Gypsy LTR retrotransposon families revealed that both have undergone multiple retrotransposon bursts about 0-1 million years ago. In addition to a recent whole-genome duplication event, there has been a lineage-specific expansion of genes involved in stilbenoid, diarylheptanoid, and gingerol biosynthesis, thereby enhancing 6-gingerol biosynthesis. Furthermore, we focused on the biosynthesis of 6-gingerol, the most important gingerol, and screened key transcription factors ZoMYB106 and ZobHLH148 that regulate 6-gingerol synthesis by transcriptomic and metabolomic analysis in the ginger rhizome at four growth stages. The results of yeast one-hybrid, electrophoretic mobility shift, and dual-luciferase reporter gene assays showed that both ZoMYB106 and ZobHLH148 bind to the promoters of the key rate-limiting enzyme genes ZoCCOMT1 and ZoCCOMT2 in the 6-gingerol synthesis pathway and promote their transcriptional activities. The reference genome, transcriptome, and metabolome data pave the way for further research on the molecular mechanism underlying the biosynthesis of 6-gingerol. Furthermore, it provides precious new resources for the study on the biology and molecular breeding of ginger.

Integrative analyses of metabolome and transcriptome reveal the dynamic accumulation and regulatory network in rhizomes and fruits of Polygonatum cyrtonema Hua.

BACKGROUND: According to Chinese ancient books, both fruits and rhizomes of Polygonatum cyrtonema Hua have medicinal and edible values. Up to now, there is no report about the metabolite profiles and regulatory network in fruits and different year-old rhizomes of P. cyrtonema. RESULTS: In this study, we performed integrative analyses of metabolome and transcriptome to reveal the dynamic accumulation and regulatory network of fruits and different year-old rhizomes in P. cyrtonema. The relative content of phenolic acids, lignans and coumarins, flavonoids and alkaloids increased with growth years, while steroids and lipids decreased with it. In addition, the relative content of nucleotides and derivatives, flavonoids, organic acids, steroids and lipids in fruits were higher than rhizomes. Genes that might relate to the biosynthesis of polysaccharides, flavonoids, triterpene saponins and alkaloids biosynthesis were further analyzed by transcriptome analysis, including sacA, GMPP, PMM, CCoAOMT, CHI, ANR, CHS, DXS, GGPS, ZEP, CYP72A219 and so on, for their expressions were positively correlated with the relative content of the metabolites. Additionally, the correlation network in sugar and aromatic amino acids metabolites were constructed to further illustrate the biosynthesis of polysaccharides, flavonoids and alkaloids in P. cyrtonema, and some transcription factors (TFs) were screened, such as C2C2, MYB, bZIP, GRAS and NAC. CONCLUSIONS: This study can deepen our understanding of the accumulation patterns and molecular mechanism of the main compounds in P. cyrtonema, and provide reference for the standardize production of P. cyrtonema.

Spatiotemporal transcriptomic atlas of rhizome formation in Oryza longistaminata.

Rhizomes are modified stems that grow underground and produce new individuals genetically identical to the mother plant. Recently, a breakthrough has been made in efforts to convert annual grains into perennial ones by utilizing wild rhizomatous species as donors, yet the developmental biology of this organ is rarely studied. Oryza longistaminata, a wild rice species featuring strong rhizomes, provides a valuable model for exploration of rhizome development. Here, we first assembled a double-haplotype genome of O. longistaminata, which displays a 48-fold improvement in contiguity compared to the previously published assembly. Furthermore, spatiotemporal transcriptomics was performed to obtain the expression profiles of different tissues in O. longistaminata rhizomes and tillers. Two spatially reciprocal cell clusters, the vascular bundle 2 cluster and the parenchyma 2 cluster, were determined to be the primary distinctions between the rhizomes and tillers. We also captured meristem initiation cells in the sunken area of parenchyma located at the base of internodes, which is the starting point for rhizome initiation. Trajectory analysis further indicated that the rhizome is regenerated through de novo generation. Collectively, these analyses revealed a spatiotemporal transcriptional transition underlying the rhizome initiation, providing a valuable resource for future perennial crop breeding.

Bud shapes dictate tiller-rhizome transition in African perennial rice (Oryza longistaminata).

KEY MESSAGE: Rhizome formation of Oryza longistaminata was dependent on the bud shape. The loci qBS3.1, qBS3.2 and qBS3.3 for controlling rhizome formation were functional redundant under Oryza longistaminata background. The rhizome, a root-like underground stem, is the key organ for grasses to achieve perennial growth. Oryza longistaminata, the only rhizomatous wild Oryza species with the same AA genome as cultivated rice, is an important germplasm for developing perennial rice. Our study found that the rhizome formation of O. longistaminata was dependent on the bud shape: the dome-like axillary bud (dome bud) usually penetrated through the leaf sheaths, developing into rhizome (extravaginal branching), but the flat axillary bud (flat bud) wrapped by the leaf sheaths only developed into tiller (intravaginal branching). The genetic loci (QTL) controlling the bud shape (BS) were mapped by entire population genotyping method (F2 population from crossing O. longistaminata with Balilla (Oryza sativa) and selective genotyping mapping method (BC1F2 population from backcrossing F1 with Balilla). A total of twelve loci were identified, including four major-effect QTL: qBS2, qBS3.1, qBS3.2 and qBS3.3, and the genetic network of these twelve loci was established. The dome bud lost the potential to develop into rhizome with the increase in backcross generations under Balilla background. Considering the rapid loss of rhizome under Balilla background, the near-isogenic lines under O. longistaminata background were used to identify the effect of major-effect loci. According to the BC3F2, BC4F2 and BC5F2 under O. longistaminata background, there was some functional redundancy among qBS3.1, qBS3.2 and qBS3.3. Our results provided a new perspective for analyzing the genetic basis of perenniality and laid the foundation for fine mapping and verification of related genes.

[Difference of effective component content and key enzyme gene expression of biosynthesis in Atractylodes chinensis with different leaf shapes].

In this study, four Atractylodes chinensis(A. chinensis) with different leaf shapes, such as the split leaf, long and narrow leaf, oval leaf, and large round leaf, were used as experimental materials to establish a method for simultaneously determining atractylodin, atractylenolide Ⅰ, ß-eudesmol, and atractylon in the rhizome of A. chinensis. The expression of key enzyme genes for biosynthesis of acetyl-CoA carboxylase(ACC), 3-hydroxy-3-methylglutaryl-CoA reductase(HMGR), and farnesyl pyrophosphate synthase(FPPS) was detected by real-time fluorescence quantitative polymerase chain reaction(qRT-PCR). High performance liquid chromatography(HPLC) was used to compare the difference in the content of four active components in A. chinensis with different leaf shapes, and the correlation between the content of active components and the expression of key enzyme genes in biosynthesis was discussed. The results show that there was good linearity among atractylodin, atractylenolide Ⅰ, ß-eudesmol, and atractylon in the range of 3.30-33.00 µg·mL~(-1)(r =0.999 7), 12.04-120.40 µg·mL~(-1)(r =0.999 5), 29.16-291.60 µg·mL~(-1)(r =0.999 5), and 14.20-142.00 µg·mL~(-1)(r =0.999 5), respectively. The average recoveries were 99.77%(RSD=2.1%), 98.56%(RSD=1.2%), 103.0%(RSD=1.2%), and 100.6%(RSD=1.5%), respectively. The method was accurate and had good reproducibility, which could be used to simultaneously detect atractylodin, atractylenolide Ⅰ, ß-eudesmol, and atractylon. The results showed that there were significant differences in the content of four active components in A. chinensis with different leaf shapes. The content of atractylodin, atractylenolide Ⅰ, and ß-eudesmol in A. chinensis with split leaves was the highest, which were 1.341 9, 5.237 2, and 12.084 3 mg·g~(-1), respectively. The content of atractylon in A. chinensis with long and narrow leaves was the highest(5.470 1 mg·g~(-1)). The content of atractylodin, atractylenolide Ⅰ, ß-eudesmol, and atractylon in A. chinensis with oval leaves was the lowest. The total content of the four effective components in descending order was A. chinensis with split leaves > A. chinensis with long and narrow leaves > A. chinensis with large round leaves > A. chinensis with oval leaves. The gene expression levels of key enzymes ACC, HMGR, and FPPS in A. chinensis with split leaves were the highest(P < 0.05), and the gene expression levels of key enzymes ACC and HMGR in A. chinensis with oval leaves were the lowest(P < 0.05). The gene expression level of key enzyme FPPS in A. chinensis with large round leaves was the lowest. In A. chinensis with different leaf shapes, the key enzyme gene ACC was significantly positively correlated with the polyacetylene component, namely atractylodin(P < 0.01), and the key enzyme genes HMGR and FPPS were positively correlated with the sesquiterpene components, namely atractylenolide Ⅰ, ß-eudesmol, and atractylon. In summary, the quality of A. chinensis with split leaves is the best, and the biosynthesis of atractylodin is significantly correlated with the gene expression of key enzyme ACC, which provides a theoretical basis for screening and optimizing the germplasm resources of A. chinensis and improving the quality of medicinal materials.

Combined analysis of microRNA and mRNA profiles provides insights into the pathogenic resistant mechanisms of the lotus rhizome rot.

Lotus rhizome rot caused by Fusarium oxysporum is a common vascular fungal disease in plants that significantly impacts the yield. However, only a few studies have studied the mechanism of Nelumbo nucifera responding to lotus rhizome rot. Here, we investigated the pathogenic genes and miRNAs in lotus rhizome rot to uncover the pathogenic resistant mechanisms by transcriptome and small RNA sequencing of lotus roots after inoculation with Fusarium oxysporum. GO and KEGG functional enrichment analysis showed that differential miRNAs were mostly enriched in starch and sucrose metabolism, biosynthesis of secondary metabolites, glutathione metabolism, brassinosteroid biosynthesis and flavonoid biosynthesis pathways. Twenty-seven upregulated miRNAs, 19 downregulated miRNAs and their target genes were identified. Correlation analysis found that miRNAs negatively regulate target genes, which were also enriched in starch and sucrose metabolism and glutathione metabolism pathways. Their expression was measured by reverse transcription quantitative PCR (qRT-PCR), and the results were consistent with the transcriptome analysis, thus verifying the reliability of transcriptome data. We selected three miRNAs (miRNA858-y, miRNA171-z and a novel miRNA novel-m0005-5p) to test the relationship between miRNAs and their target genes. The activity of the GUS testing assay indicated that miRNA could decrease the GUS activity by inhibiting the expression of their target genes. Collectively, this study provides a comprehensive analysis of transcriptome and small RNA sequencing of lotus root after inoculation with Fusarium oxysporum, and we identified candidate miRNAs and their target genes for breeding strategies of Nelumbo nucifera.

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

KEY MESSAGE: 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.

The chrysanthemum DEAD-box RNA helicase CmRH56 regulates rhizome outgrowth in response to drought stress.

Plants have evolved complex mechanisms to reprogram growth in response to drought stress. In herbaceous perennial plant species, the rhizome, which is normally an organ for propagation and food storage, can also support plant growth in stressful environments, and allows the plant to perennate and survive stress damage. However, the mechanisms that regulate rhizome growth in perennial herbs during abiotic stresses are unknown. Here, we identified a chrysanthemum (Chrysanthemum morifolium) DEAD-box RNA helicase gene, CmRH56, that is specifically expressed in the rhizome shoot apex. Knock down of CmRH56 transcript levels decreased the number of rhizomes and enhanced drought stress tolerance. We determined that CmRH56 represses the expression of a putative gibberellin (GA) catabolic gene, GA2 oxidase6 (CmGA2ox6). Exogenous GA treatment and silencing of CmGA2ox6 resulted in more rhizomes. These results demonstrate that CmRH56 suppresses rhizome outgrowth under drought stress conditions by blocking GA biosynthesis.

Hormone and carbohydrate metabolism associated genes play important roles in rhizome bud full-year germination of Cephalostachyum pingbianense.

Cephalostachyum pingbianense is the only woody bamboo species that can produce bamboo shoots in four seasons under natural conditions. So far, the regulatory mechanism of shoot bud differentiation and development is unknown. In the present study, indole-3-acetic acid (IAA), zeatin riboside (ZR), gibberellin A3 (GA3 ) and abscisic acid (ABA) contents determination, RNA sequencing and differentially expressed gene analysis were performed on dormant rhizome bud (DR), growing rhizome bud (GR), and germinative bud (GB) in each season. The results showed that the contents of IAA and ZR increased while ABA content decreased, and GA3 content was stable during bud transition from dormancy to germination in each season. Moreover, rhizome bud germination was cooperatively regulated by multiple pathways such as carbohydrate metabolism, hormone signal transduction, cell wall biogenesis, temperature response, and water transport. The inferred hub genes among these candidates were identified by protein-protein interaction network analyses, most of which were involved in hormone and carbohydrate metabolism, such as HK and BGLU4 in spring, IDH and GH3 in winter, GPI and talA/talB in summer and autumn. It is speculated that dynamic phytohormone changes and differential expression of these genes promote the release of rhizome bud dormancy and contribute to the phenological characteristics of full-year shooting. Moreover, the rhizome buds of C. pingbianense may not suffer from ecodormancy in winter. These findings would help accumulate knowledge on shooting mechanisms in woody bamboos and provide a physiological insight into germplasm conservation and forest management of C. pingbianense.

A nucleotide signature for the identification of Pinelliae Rhizoma (Banxia) and its products.

BACKGROUND: Ensuring the authenticity of raw materials is a key step prior to producing Chinese patent medicines. Pinellia ternata (Thunb.) Breit. is the botanical origin of Pinelliae Rhizoma (Banxia), a traditional Chinese medicine used to treat cough, insomnia, nausea, inflammation, epilepsy, and so on. Unfortunately, authentic Pinelliae Rhizoma is often adulterated by morphologically indistinguishable plant material due to the insufficient regulatory procedures of processed medicinal plant products. Thus, it is important to develop a molecular assay based on species-specific nucleotide signatures and primers to efficiently distinguish authentic Pinelliae Rhizoma from its adulterants. METHODS AND RESULTS: The ITS2 region of 67 Pinelliae Rhizoma and its common adulterants were sequenced. Eight single nucleotide polymorphisms within a 28-43 bp stretch of ITS2 were used to develop six primer pairs to amplify these species-specific regions. We assayed 56 Pinelliae Rhizoma products sold on the Chinese market, including medicinal slices, powder and Chinese patent medicines, which revealed that about 66% of products were adulterated. The most common adulterants were Pinellia pedatisecta (found in 57% of the assayed products), Arisaema erubescens (9%), Typhonium giganteum (2%) and Typhonium flagelliforme (2%). CONCLUSIONS: A severe adulteration condition was revealed in the traditional medicine market. The species-specific nucleotide assays developed in this study can be applied to reliably identify Pinelliae Rhizoma and its adulterants, aiding in the authentication and quality control of processed products on the herbal market.

Differences in gene expression and endophytic bacterial diversity in Atractylodes macrocephala Koidz. rhizomes from different growth years.

Atractylodes macrocephala Koidz. (AMK) is widely used in traditional Chinese medicine owing to its pharmacological activity. Here, we aimed to characterize the differentially expressed genes (DEGs) of one- and three-year growth (OYG and TYG) rhizomes of AMK, combined with endophytic bacterial diversity analysis using high-throughput RNA sequencing. A total of 114 572 unigenes were annotated using six public databases. In all, 3570 DEGs revealed a clear difference, of which 936 and 2634 genes were upregulated and downregulated, respectively. The results of KEGG pathway analysis indicated that DEGs corresponding to terpenoid synthesis gene were downregulated in TYG rhizomes. In addition, 414 424 sequences corresponding to the 16S rRNA gene were divided into 1267 operational taxonomic units (OTUs). Moreover, the diversity of endophytic bacteria changed with species in the OYG (773) and TYG (1201) rhizomes at the OTU level, and Proteobacteria, Actinobacteria, and Bacteroidetes were the dominant phyla. A comparison of species differences among different growth years revealed that some species were significantly different, such as Actinomycetes, Variovorax, and Cloacibacterium. Interestingly, the decrease in the function-related metabolism of terpenoids and polyketides was correlated with the low expression of terpene synthesis genes in TYG rhizomes, as assessed using PICRUSt2. These data provide a scientific basis for elucidating the mechanisms underlying metabolite accumulation and endophytic bacterial diversity in relation to the growth years in AMK.

Transcriptome Analysis of Different Sections of Rhizome in Polygonatum sibiricum Red. and Mining Putative Genes Participate in Polysaccharide Biosynthesis.

As a traditional Chinese herb, the rhizomes of Polygonatum sibiricum Red. are rich in various compounds which have plenty of pharmacological applications and biological activities. Among them, Polygonatum sibiricum polysaccharides (PSP) are the main active ingredients and exhibit a broad range of pharmacological. Based on previous researches, identifying genes involved in PSP biosynthesis will help delineate such pathway at the molecular level. In that case, we performed RNA sequencing analysis for two sections of P. sibiricum Red.'s rhizomes significantly different in PSP content. A total of 435,858 unigenes were obtained by assembling transcripts from both sections and 29,548 (6.77%) ones were annotated in all seven public databases. Analyzing count data of RNA-seq, 13,460 differential expression genes (DEGs) between two sections of rhizomes were acquired. After DEGs were mapped to KEGG databases, twelve represented KEGG pathways related to PSP biosynthesis were summed up. And most DEGs were assigned to the pathway of "Starch and sucrose metabolism". Finally, seventeen candidate genes whose expression levels were related to the polysaccharide content, were considered involving PSP biosynthesis in P. sibiricum Red. The present study lays a foundation of researching the molecular mechanisms of PSP biosynthesis.

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.

[Comparison of transcriptome of Atractylodes lancea rhizome and exploration of genes for sesquiterpenoid biosynthesis].

This study compared the transcriptome of Atractylodes lancea rhizome at different development stages and explored genes encoding the key enzymes of the sesquiterpenoid biosynthesis pathway. Specifically, Illumina NovaSeq 6000 was employed for sequencing the cDNA libraries of A. lancea rhizome samples at the growth stage(SZ), flowering stage(KH), and harvesting stage(CS), respectively. Finally, a total of 388 201 748 clean reads were obtained, and 16 925, 8 616, and 13 702 differentially expressed genes(DEGs) were identified between SZ and KH, KH and CS, and SZ and CS, separately. Among them, 53 genes were involved in the sesquiterpenoid biosynthesis pathways: 9 encoding 6 enzymes of the mevalonic acid(MVA) pathway, 15 encoding 7 enzymes of the 2-C-methyl-D-erythritol-4-phosphate(MEP) pathway, and 29 of sesquiterpenoid and triterpenoid biosynthesis pathway. Weighted gene co-expression network analysis(WGCNA) yielded 12 genes related to sesquiterpenoid biosynthesis for the SZ, 1 gene for the KH, and 1 gene for CS, and several candidate genes for sesquiterpenoid biosynthesis were discovered based on the co-expression network. This study laid a solid foundation for further research on the sesquiterpenoid biosynthesis pathway, analysis of the regulation mechanism, and mechanism for the accumulation of sesquiterpenoids in A. lancea.

Resequencing of 296 cultivated and wild lotus accessions unravels its evolution and breeding history.

Lotus (family: Nelumbonaceae) are perennial aquatic plants that represent one of the most ancient basal dicots. In the present study, we resequenced 296 lotus accessions from various geographical locations and germplasms to explore their genomic diversity and population structure. This germplasm set consisted of four accessions of American wild lotus and 292 accessions of Asian lotus, which were divided into four subgroups: wild, rhizome, flower and seed. Total single nucleotide polymorphisms (SNPs) suggested that the wild lotus had the highest variant number (7 191 010). Population structure and genome diversity analysis indicated that the American wild lotus demonstrated a distant genetic relationship with the Asian lotus. Furthermore, the seed and rhizome lotus groups had not originated from a single source but rather had a more complex multisource origin. Besides that, the seed lotus showed higher genetic diversity, which might have been due to the gene flow from the flower lotus to seed lotus by artificial crossing, and the rhizome lotus showed a much lower genetic diversity than the other groups. The present study provides SNP markers for lotus genomic diversity analysis, which will be useful for guiding lotus breeding.

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.

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

KEY MESSAGE: 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.

The complete plastid genome sequence of the enigmatic moss, Takakia lepidozioides (Takakiopsida, Bryophyta): evolutionary perspectives on the largest collection of genes in mosses and the intensive RNA editing.

KEY MESSAGE: Complete chloroplast genome sequence of a moss, Takakia lepidozioides (Takakiopsida) is reported. The largest collection of genes in mosses and the intensive RNA editing were discussed from evolutionary perspectives. We assembled the entire plastid genome sequence of Takakia lepidozioides (Takakiopsida), emerging from the first phylogenetic split among extant mosses. The genome sequences were assembled into a circular molecule 149,016 bp in length, with a quadripartite structure comprising a large and a small single-copy region separated by inverted repeats. It contained 88 genes coding for proteins, 32 for tRNA, four for rRNA, two open reading frames, and at least one pseudogene (tufA). This is the largest number of genes of all sequenced plastid genomes in mosses and Takakia is the only moss that retains the seven coding genes ccsA, cysA, cysT, petN rpoA, rps16 and trnPGGG. Parsimonious interpretation of gene loss suggests that the last common ancestor of bryophytes had all seven genes and that mosses lost at least three of them during their diversification. Analyses of the plastid transcriptome identified the extraordinary frequency of RNA editing with more than 1100 sites. We indicated a close correlation between the monoplastidy of vegetative tissue and the intensive RNA editing sites in the plastid genome in land plant lineages. Here, we proposed a hypothesis that the small population size of plastids in each vegetative cell of some early diverging land plants, including Takakia, might cause the frequent fixation of mutations in plastid genome through the intracellular genetic drift and that deleterious mutations might be continuously compensated by RNA editing during or following transcription.

An evolutionarily conserved NIMA-related kinase directs rhizoid tip growth in the basal land plant Marchantia polymorpha.

Tip growth is driven by turgor pressure and mediated by the polarized accumulation of cellular materials. How a single polarized growth site is established and maintained is unclear. Here, we analyzed the function of NIMA-related protein kinase 1 (MpNEK1) in the liverwort Marchantia polymorpha In the wild type, rhizoid cells differentiate from the ventral epidermis and elongate through tip growth to form hair-like protrusions. In Mpnek1 knockout mutants, rhizoids underwent frequent changes in growth direction, resulting in a twisted and/or spiral morphology. The functional MpNEK1-Citrine protein fusion localized to microtubule foci in the apical growing region of rhizoids. Mpnek1 knockouts exhibited increases in both microtubule density and bundling in the apical dome of rhizoids. Treatment with the microtubule-stabilizing drug taxol phenocopied the Mpnek1 knockout. These results suggest that MpNEK1 directs tip growth in rhizoids through microtubule organization. Furthermore, MpNEK1 expression rescued ectopic outgrowth of epidermal cells in the Arabidopsis thaliana nek6 mutant, strongly supporting an evolutionarily conserved NEK-dependent mechanism of directional growth. It is possible that such a mechanism contributed to the evolution of the early rooting system in land plants.

Association among starch storage, metabolism, related genes and growth of Moso bamboo (Phyllostachys heterocycla) shoots.

BACKGROUND: Both underground rhizomes/buds and above-ground Moso bamboo (Phyllostachys heterocycla) shoots/culms/branches are connected together into a close inter-connecting system in which nutrients are transported and shared among each organ. However, the starch storage and utilization mechanisms during bamboo shoot growth remain unclear. This study aimed to reveal in which organs starch was stored, how carbohydrates were transformed among each organ, and how the expression of key genes was regulated during bamboo shoot growth and developmental stages which should lay a foundation for developing new theoretical techniques for bamboo cultivation. RESULTS: Based on changes of the NSC content, starch metabolism-related enzyme activity and gene expression from S0 to S3, we observed that starch grains were mainly elliptical in shape and proliferated through budding and constriction. Content of both soluble sugar and starch in bamboo shoot peaked at S0, in which the former decreased gradually, and the latter initially decreased and then increased as shoots grew. Starch synthesis-related enzymes (AGPase, GBSS and SBE) and starch hydrolase (α-amylase and ß-amylase) activities exhibited the same dynamic change patterns as those of the starch content. From S0 to S3, the activity of starch synthesis-related enzyme and starch amylase in bamboo rhizome was significantly higher than that in bamboo shoot, while the NSC content in rhizomes was obviously lower than that in bamboo shoots. It was revealed by the comparative transcriptome analysis that the expression of starch synthesis-related enzyme-encoding genes were increased at S0, but reduced thereafter, with almost the same dynamic change tendency as the starch content and metabolism-related enzymes, especially during S0 and S1. It was revealed by the gene interaction analysis that AGPase and SBE were core genes for the starch and sucrose metabolism pathway. CONCLUSIONS: Bamboo shoots were the main organ in which starch was stored, while bamboo rhizome should be mainly functioned as a carbohydrate transportation channel and the second carbohydrate sink. Starch metabolism-related genes were expressed at the transcriptional level during underground growth, but at the post-transcriptional level during above-ground growth. It may be possible to enhance edible bamboo shoot quality for an alternative starch source through genetic engineering.