Genome-wide transcriptome analysis reveals small RNA profiles involved in early stages of stolon-to-tuber transitions in potato under photoperiodic conditions

The phased short-interfering RNA siRD29(-) regulates GIBBERELLIN 3-OXIDASE 3 during stolon-to-tuber transitions in potato

Phased short-interfering RNAs (phasiRNAs) fine tune various stages of growth, development, and stress responses in plants. Potato (Solanum tuberosum) tuberization is a complex process, wherein a belowground modified stem (stolon) passes through developmental stages like swollen stolon and minituber before it matures to a potato. Previously, we identified several phasiRNA-producing loci (PHAS) from stolon-to-tuber transition stages. However, whether phasiRNAs mediate tuber development remains unknown. Here, we show that a gene encoding NB-ARC DOMAIN-CONTAINING DISEASE RESISTANCE PROTEIN (StRGA4; a PHAS locus) is targeted by Stu-microRNA482c to generate phasiRNAs. Interestingly, we observed that one of the phasiRNAs, referred as short-interfering RNA D29(-), i.e. siRD29(-), targets the gibberellin (GA) biosynthesis gene GIBBERELLIN 3-OXIDASE 3 (StGA3ox3). Since regulation of bioactive GA levels in stolons controls tuber development, we hypothesized that a gene regulatory module, Stu-miR482c-StRGA4-siRD29(-)-StGA3ox3, could govern tuber development. Through transient expression assays and small RNA sequencing, generation of siRD29(-) and its phase was confirmed in planta. Notably, the expression of StGA3ox3 was higher in swollen stolon compared to stolon, whereas siRD29(-) showed a negative association with StGA3ox3 expression. Antisense (AS) lines of StGA3ox3 produced more tubers compared to wild type. As expected, StRGA4 overexpression (OE) lines had high levels of siRD29(-) and mimicked the phenotypes of StGA3ox3-AS lines, indicating the functionality of this module in potato. In vitro tuberization assays (with or without a GA inhibitor) using StGA3ox3 antisense lines and overexpression lines of StGA3ox3 or StRGA4 revealed that StGA3ox3 controls the tuber stalk development. Taken together, our findings suggest that a phasiRNA, siRD29(-), mediates the regulation of StGA3ox3 during stolon-to-tuber transitions in potato.

Identification of miRNAs Involving Potato-Phytophthora infestans Interaction

sRNAs (small RNAs) play an important role in regulation of plant immunity against a variety of pathogens. In this study, sRNA sequencing analysis was performed to identify miRNAs (microRNAs) during the interaction of potato and Phytophthora infestans. Totally, 171 potato miRNAs were identified, 43 of which were annotated in the miRNA database and 128 were assigned as novel miRNAs in this study. Those potato miRNAs may target 878 potato genes and half of them encode resistance proteins. Fifty-three potato miRNAs may target 194 P. infestans genes. Three potato miRNAs (novel 72, 133, and 140) were predicted to have targets only in the P. infestans genome. miRNAs transient expression and P. infestans inoculation assay showed that miR396, miR166, miR6149-5P, novel133, or novel140 promoted P. infestans colonization, while miR394 inhibited colonization on Nicotiana benthamiana leaves. An artificial miRNA target (amiRNA) degradation experiment demonstrated that miR394 could target both potato gene (PGSC0003DMG400034305) and P. infestans genes. miR396 targets the multicystatin gene (PGSC0003DMG400026899) and miR6149-5p could shear the galactose oxidase F-box protein gene CPR30 (PGSC0003DMG400021641). This study provides new information on the aspect of cross-kingdom immune regulation in potato-P. infestans interaction at the sRNAs regulation level.

Solanum tuberosum Microtuber Development under Darkness Unveiled through RNAseq Transcriptomic Analysis

Potato microtuber (MT) development through in vitro techniques are ideal propagules for producing high quality potato plants. MT formation is influenced by several factors, i.e., photoperiod, sucrose, hormones, and osmotic stress. We have previously developed a protocol of MT induction in medium with sucrose (8% w/v), gelrite (6g/L), and 2iP as cytokinin under darkness. To understand the molecular mechanisms involved, we performed a transcriptome-wide analysis. Here we show that 1715 up- and 1624 down-regulated genes were involved in this biological process. Through the protein-protein interaction (PPI) network analyses performed in the STRING database (v11.5), we found 299 genes tightly associated in 14 clusters. Two major clusters of up-regulated proteins fundamental for life growth and development were found: 29 ribosomal proteins (RPs) interacting with 6 PEBP family members and 117 cell cycle (CC) proteins. The PPI network of up-regulated transcription factors (TFs) revealed that at least six TFs-MYB43, TSF, bZIP27, bZIP43, HAT4 and WOX9-may be involved during MTs development. The PPI network of down-regulated genes revealed a cluster of 83 proteins involved in light and photosynthesis, 110 in response to hormone, 74 in hormone mediate signaling pathway and 22 related to aging.

Screening of differentially expressed microRNAs and target genes in two potato varieties under nitrogen stress

BACKGROUND

A reasonable supply of nitrogen (N) fertilizer is essential for obtaining high-quality, high-level, and stable potato yields, and an improvement in the N utilization efficiency can effectively reduce N fertilizer use. It is important to use accurate, straightforward, and efficient transgenic breeding techniques for the identification of genes that can improve nitrogen use efficiency, thus enabling us to achieve the ultimate goal of breeding N-efficient potato varieties. In recent years, some of the mechanisms of miRNAs have been elucidated via the analysis of the correlation between the expression levels of potato miRNA target genes and regulated genes under conditions of stress, but the role of miRNAs in the inhibition/expression of key genes regulating N metabolism under N stress is still unclear. Our study aimed to identify the role played by specific enzymes and miRNAs in the responses of plants to N stress.

RESULTS

The roots and leaves of the N-efficient potato variety, Yanshu4 ("Y"), and N-inefficient potato variety, Atlantic ("D"), were collected at the seedling and budding stages after they were exposed to different N fertilizer treatments. The miRNAs expressed differentially under the two types of N stress and their corresponding target genes were first predicted using miRNA and degradome analysis. Then, quantitative polymerase chain reaction (qRT-PCR) was performed to verify the expression of differential miRNAs that were closely related to N metabolism. Finally, the shearing relationship between stu-miR396-5p and its target gene StNiR was determined by analyzing luciferase activity levels. The results showed that NiR activity increased significantly with an increase in the applied N levels from the seedling stage to the budding stage, and NiR responded significantly to different N treatments. miRNA sequencing enabled us to predict 48 families with conserved miRNAs that were mainly involved in N metabolism, carbon metabolism, and amino acid biosynthesis. The differences in the expression of the following miRNAs were identified via screening (high expression levels and P < 0.05): stu-miR396-5p, stu-miR408b-3p_R-1, stu-miR3627-3p, stu-miR482a-3p, stu-miR8036-3p, stu-miR482a-5p, stu-miR827-5p, stu-miR156a_L-1, stu-miR827-3p, stu-miR172b-5p, stu-miR6022-p3_7, stu-miR398a-5p, and stu-miR166c-5p_L-3. Degradome analysis showed that most miRNAs had many-to-many relationships with target genes. The main target genes involved in N metabolism were NiR, NiR1, NRT2.5, and NRT2.7. qRT-PCR analysis showed that there were significant differences in the expression levels of stu-miR396-5p, stu-miR8036-3p, and stu-miR482a-3p in the leaves and roots of the Yanshu4 and Atlantic varieties at the seedling and budding stages under conditions that involved no N and excessive N application; the expression of these miRNAs was induced in response to N stress. The correlation between the differential expression of stu-miR396-5p and its corresponding target gene NiR was further verified by determining the luciferase activity level and was found to be strongly negative.

CONCLUSION

The activity of NiR was significantly positively correlated with N application from the seedling to the budding stage. Differential miRNAs and target genes showed a many-to-many relationship with each other. The expression of stu-miR396-5p, stu-miR482a-3p, and stu-miR8036-3p in the roots and leaves of the Yanshu4 and Atlantic varieties at the seedling and budding stages was notably different under two types of N stress. Under two types of N stress, stu-miR396-5p was down-regulated in Yanshu4 in the seedling-stage and shoot-stage roots, and up-regulated in seedling-stage roots and shoot-stage leaves; stu-miR482a-3p was up-regulated in the seedling and shoot stages. The expression of stu-miR8036-3p was up-regulated in the leaves and roots at the seedling and budding stages, and down-regulated in roots under both types of N stress. The gene expressing the key enzyme involved in N metabolism, StNiR, and the stu-miR396-5p luciferase assay reporter gene had a strong regulatory relationship with each other. This study provides candidate miRNAs related to nitrogen metabolism and highlights that differential miRNAs play a key role in nitrogen stress in potato, providing insights for future research on miRNAs and their target genes in nitrogen metabolic pathways and breeding nitrogen-efficient potatoes.

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.

Starch metabolism in potato <i>Solanum tuberosum</i> L

Starch is a major storage carbohydrate in plants. It is an important source of calories in the human and animal diet. Also, it is widely used in various industries. Native starch consists of water-insoluble semicrystalline granules formed by natural glucose polymers amylose and amylopectin. The physicochemical properties of starch are determined by the amylose:amylopectin ratio in the granule and degrees of their polymerization and phosphorylation. Potato Solanum tuberosum L. is one of the main starch-producing crops. Growing industrial needs necessitate the breeding of plant varieties with increased starch content and specified starch properties. This task demands detailed information on starch metabolism in the producing plant. It is a complex process, requiring the orchestrated work of many enzymes, transporter and targeting proteins, transcription factors, and other regulators. Two types of starch are recognized with regard to their biological functions. Transitory starch is synthesized in chloroplasts of photosynthetic organs and degraded in the absence of light, providing carbohydrates for cell needs. Storage starch is synthesized and stored in amyloplasts of storage organs: grains and tubers. The main enzymatic reactions of starch biosynthesis and degradation, as well as carbohydrate transport and metabolism, are well known in the case of transitory starch of the model plant Arabidopsis thaliana. Less is known about features of starch metabolism in storage organs, in particular, potato tubers. Several issues remain obscure: the roles of enzyme isoforms and different regulatory factors in tissues at various plant developmental stages and under different environmental conditions; alternative enzymatic processes; targeting and transport proteins. In this review, the key enzymatic reactions of plant carbohydrate metabolism, transitory and storage starch biosynthesis, and starch degradation are discussed, and features specific for potato are outlined. Attention is also paid to the known regulatory factors affecting starch metabolism

Development of aerial and belowground tubers in potato is governed by photoperiod and epigenetic mechanism

Plants exhibit diverse developmental plasticity and modulate growth responses under various environmental conditions. Potato (Solanum tuberosum), a modified stem and an important food crop, serves as a substantial portion of the world's subsistence food supply. In the past two decades, crucial molecular signals have been identified that govern the tuberization (potato development) mechanism. Interestingly, microRNA156 overexpression in potato provided the first evidence for induction of profuse aerial stolons and tubers from axillary meristems under short-day (SD) photoperiod. A similar phenotype was noticed for overexpression of epigenetic modifiers-MUTICOPY SUPRESSOR OF IRA1 (StMSI1) or ENAHNCER OF ZESTE 2 (StE[z]2), and knockdown of B-CELL-SPECIFIC MOLONEY MURINE LEUKEMIA VIRUS INTEGRATION SITE 1 (StBMI1). This striking phenotype represents a classic example of modulation of plant architecture and developmental plasticity. Differentiation of a stolon to a tuber or a shoot under in vitro or in vivo conditions symbolizes another example of organ-level plasticity and dual fate acquisition in potato. Stolon-to-tuber transition is governed by SD photoperiod, mobile RNAs/proteins, phytohormones, a plethora of small RNAs and their targets. Recent studies show that polycomb group proteins control microRNA156, phytohormone metabolism/transport/signaling and key tuberization genes through histone modifications to govern tuber development. Our comparative analysis of differentially expressed genes between the overexpression lines of StMSI1, StBEL5 (BEL1-LIKE transcription factor [TF]), and POTATO HOMEOBOX 15 TF revealed more than 1,000 common genes, indicative of a mutual gene regulatory network potentially involved in the formation of aerial and belowground tubers. In this review, in addition to key tuberization factors, we highlight the role of photoperiod and epigenetic mechanism that regulates the development of aerial and belowground tubers in potato.

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.

smartPARE: An R Package for Efficient Identification of True mRNA Cleavage Sites

Degradome sequencing is commonly used to generate high-throughput information on mRNA cleavage sites mediated by small RNAs (sRNA). In our datasets of potato (Solanum tuberosum, St) and Phytophthora infestans (Pi), initial predictions generated high numbers of cleavage site predictions, which highlighted the need of improved analytic tools. Here, we present an R package based on a deep learning convolutional neural network (CNN) in a machine learning environment to optimize discrimination of false from true cleavage sites. When applying smartPARE to our datasets on potato during the infection process by the late blight pathogen, 7.3% of all cleavage windows represented true cleavages distributed on 214 sites in P. infestans and 444 sites in potato. The sRNA landscape of the two organisms is complex with uneven sRNA production and cleavage regions widespread in the two genomes. Multiple targets and several cases of complex regulatory cascades, particularly in potato, was revealed. We conclude that our new analytic approach is useful for anyone working on complex biological systems and with the interest of identifying cleavage sites particularly inferred by sRNA classes beyond miRNAs.

Hormonal and transcriptional analyses provides new insights into the molecular mechanisms underlying root thickening and isoflavonoid biosynthesis in Callerya speciosa (Champ. ex Benth.) Schot

Callerya speciosa (Champ. ex Benth.) Schot is a traditional Chinese medicine characterized by tuberous roots as the main organ of isoflavonoid accumulation. Root thickening and isoflavonoid accumulation are two major factors for yield and quality of C. speciosa. However, the underlying mechanisms of root thickening and isoflavonoid biosynthesis have not yet been elucidated. Here, integrated morphological, hormonal and transcriptomic analyses of C. speciosa tuberous roots at four different ages (6, 12, 18, 30 months after germination) were performed. The growth cycle of C. speciosa could be divided into three stages: initiation, rapid-thickening and stable-thickening stage, which cued by the activity of vascular cambia. Endogenous changes in phytohormones were associated with developmental changes during root thickening. Jasmonic acid might be linked to the initial development of tuberous roots. Abscisic acid seemed to be essential for tuber maturation, whereas IAA, cis-zeatin and gibberellin 3 were considered essential for rapid thickening of tuberous roots. A total of 4337 differentially expressed genes (DEGs) were identified during root thickening, including 15 DEGs participated in isoflavonoid biosynthesis, and 153 DEGs involved in starch/sucrose metabolism, hormonal signaling, transcriptional regulation and cell wall metabolism. A hypothetical model of genetic regulation associated with root thickening and isoflavonoid biosynthesis in C. speciosa is proposed, which will help in understanding the underlying mechanisms of tuberous root formation and isoflavonoid biosynthesis.

Invertase inhibitors in potato: towards a biochemical and molecular understanding of cold-induced sweetening

Invertase inhibitors classified as cell wall/apoplastic and vacuolar belonging to the pectin methylesterase family, play a major role in cold-induced sweetening (CIS) process of potato tubers. The CIS process is controlled at the post-translational level via an interaction between invertase (cell wall/apoplastic and vacuolar) by their compartment-specific inhibitors (cell wall/apoplastic and vacuolar). Invertase inhibitors have been cloned, sequenced and functionally characterized from potato cultivars differing in their CIS ability. The secondary structure of the invertase inhibitors consisted of seven alpha-helices and four conserved cysteine residues. The well-conserved three amino acids i.e. Pro-Lys-Phe are known to interact with invertase. Location of the genes encoding cell wall/apoplastic and vacuolar invertase inhibitors on potato chromosome number twelve in a tandem orientation without any intervening genes suggest their divergence into the cell wall and vacuole forms following the event of gene duplication. Under cold storage conditions, the vacuolar invertase inhibitor gene showed developmentally regulated alternative splicing and produce hybrid mRNAs which were the result of mRNA splicing of an upstream region of vacuolar invertase inhibitor gene to a downstream region of the apoplastic invertase inhibitor gene. Transgenic potato tubers overexpressing invertase inhibitors resulted in decreased invertase activity, low reducing sugars and improved processing quality making invertase inhibitors highly potential candidate genes for overcoming CIS. Clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) gene-editing technology offers transgene-free breeding for developing CIS resistant potato cultivars. Moreover, the post-transcriptional regulation of invertase inhibitors during cold storage can be warranted. This review summarizes progress and current knowledge on biochemical and molecular approaches used for the understanding of invertase inhibitors with special reference to key findings in potato.

Genome-wide identification and characterization of microRNAs by small RNA sequencing for low nitrogen stress in potato

Nitrogen is an important nutrient for plant growth and tuber quality of potato. Since potato crop requires high dose of N, improving nitrogen use efficiency (NUE) of plant is an inevitable approach to minimize N fertilization. The aim of this study was to identify and characterize microRNAs (miRNAs) by small RNA sequencing in potato plants grown in aeroponic under two contrasting N (high and low) regimes. A total of 119 conserved miRNAs belonging to 41 miRNAs families, and 1002 putative novel miRNAs were identified. From total, 52 and 54 conserved miRNAs, and 404 and 628 putative novel miRNAs were differentially expressed in roots and shoots, respectively under low N stress. Of total 34,135 predicted targets, the gene ontology (GO) analysis indicated that maximum targets belong to biological process followed by molecular function and cellular component. Eexpression levels of the selected miRNAs and targets were validated by real time-quantitative polymerase chain reaction (RT-qPCR) analysis. Two predicted targets of potential miRNAs (miR397 and miR398) were validated by 5' RLM-RACE (RNA ligase mediated rapid amplification of cDNA ends). In general, predicted targets are associated with stress-related, kinase, transporters and transcription factors such as universal stress protein, heat shock protein, salt-tolerance protein, calmodulin binding protein, serine-threonine protein kinsae, Cdk10/11- cyclin dependent kinase, amino acid transporter, nitrate transporter, sugar transporter, transcription factor, F-box family protein, and zinc finger protein etc. Our study highlights that miR397 and miR398 play crucial role in potato during low N stress management. Moreover, study provides insights to modulate miRNAs and their predicted targets to develop N-use efficient potato using transgenic/genome-editing tools in future.

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.

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.

Overexpression of Solanum habrochaites microRNA319d (sha-miR319d) confers chilling and heat stress tolerance in tomato (S. lycopersicum)

BACKGROUND

MicroRNA319 (miR319) acts as an essential regulator of gene expression during plant development and under stress conditions. Although the role of miR319a in regulating leaf development has been well studied in tomato (Solanum lycopersicum), the function of the recently discovered wild tomato Solanum habrochaites miRNA319d (sha-miR319d) remains poorly understood. In this study, we overexpressed sha-miR319d in cultivated tomato 'Micro-Tom' to further investigate its role in tomato temperature stress responses.

RESULTS

Under chilling or heat stress, sha-miR319d-overexpressing plants showed enhanced stress tolerance, including lower relative electrolyte leakage (REL), malondialdehyde (MDA) concentration, O₂^- generation and H₂O₂ concentration and higher chlorophyll contents and Fv/Fm values than wild-type (WT) plants. Overexpression of sha-miR319d enhanced the activities of superoxide dismutase (SOD) and catalase (CAT), with possible correlation with elevated expression levels of the genes FeSOD, CuZnSOD and CAT. Moreover, different expression levels of key genes involved in chilling (MYB83 and CBF1), heat (HsfA1a, HsfA1b and Hsp90), and reactive oxygen species (ROS) (ZAT12 and ZAT10) signaling in transgenic plants and WT were determined, suggesting a role for sha-miR319d in regulating tomato temperature stress via chilling, heat and ROS signaling. Silencing GAMYB-like1 increased tomato chilling tolerance as well as the expression levels of CBF1, CuZnSOD, CAT, APX1, APX2, ZAT12 and ZAT10. Additionally, overexpression of sha-miR319d in tomato caused plant leaf crinkling and reduced height.

CONCLUSIONS

Overexpression of sha-miR319d confers chilling and heat stress tolerance in tomato. Sha-miR319d regulates tomato chilling tolerance, possibly by inhibiting expression of GAMYB-like1 and further alters chilling, heat and ROS signal transduction. Our research provides insight for further study of the role of sha-miR319d in tomato growth and stress regulation and lays a foundation for the genetic improvement of tomato.

Post-transcriptional Regulation of FLOWERING LOCUS T Modulates Heat-Dependent Source-Sink Development in Potato

Understanding tuberization in the major crop plant potato (Solanum tuberosum L.) is of importance to secure yield even under changing environmental conditions. Tuber formation is controlled by a homolog of the floral inductor FLOWERING LOCUS T, referred to as SP6A. To gain deeper insights into its function, we created transgenic potato plants overexpressing a codon-optimized version of SP6A, SP6A^cop, to avoid silencing effects. These plants exhibited extremely early tuberization at the juvenile stage, hindering green biomass development and indicating a tremendous shift in the source sink balance. The meristem identity was altered in dormant buds of transgenic tubers. This strong phenotype, not being reported so far for plants overexpressing an unmodified SP6A, could be due to post-transcriptional regulation. In fact, a putative SP6A-specific small regulatory RNA was identified in potato. It was effectively repressing SP6A mRNA accumulation in transient assays as well as in leaves of young potato plants prior to tuber formation. SP6A expression is downregulated under heat, preventing tuberization. The molecular mechanism has not been elucidated yet. We showed that this small RNA is strongly upregulated under heat. The importance of the small RNA was demonstrated by overexpression of a target mimicry construct, which led to an increased SP6A expression, enabling tuberization even under continuous heat conditions, which abolished tuber formation in the wild-type. Thus, our study describes an additional regulatory mechanism for SP6A besides the well-known pathway that integrates both developmental and environmental signals to control tuberization and is therefore a promising target for breeding of heat-tolerant potato.