High throughput deep sequencing reveals the important roles of microRNAs during sweetpotato storage at chilling temperature

De novo assembly of the complete mitochondrial genome of sweet potato (Ipomoea batatas [L.] Lam) revealed the existence of homologous conformations generated by the repeat-mediated recombination

Sweet potato (Ipomoea batatas [L.] Lam) is an important food crop, an excellent fodder crop, and a new type of industrial raw material crop. The lack of genomic resources could affect the process of industrialization of sweet potato. Few detailed reports have been completed on the mitochondrial genome of sweet potato. In this research, we sequenced and assembled the mitochondrial genome of sweet potato and investigated its substructure. The mitochondrial genome of sweet potato is 270,304 bp with 23 unique core genes and 12 variable genes. We detected 279 pairs of repeat sequences and found that three pairs of direct repeats could mediate the homologous recombination into four independent circular molecules. We identified 70 SSRs in the whole mitochondrial genome of sweet potato. The longest dispersed repeat in mitochondrial genome was a palindromic repeat with a length of 915 bp. The homologous fragments between the chloroplast and mitochondrial genome account for 7.35% of the mitochondrial genome. We also predicted 597 RNA editing sites and found that the rps3 gene was edited 54 times, which occurred most frequently. This study further demonstrates the existence of multiple conformations in sweet potato mitochondrial genomes and provides a theoretical basis for the evolution of higher plants and cytoplasmic male sterility breeding.

Blocking IbmiR319a Impacts Plant Architecture and Reduces Drought Tolerance in Sweet Potato

MicroRNA319 (miR319) plays a key role in plant growth, development, and multiple resistance by repressing the expression of targeted TEOSINTE BRANCHED/CYCLOIDEA/PCF (TCP) genes. Two members, IbmiR319a and IbmiR319c, were discovered in the miR319 gene family in sweet potato (Ipomoea batatas [L.] Lam). Here, we focused on the biological function and potential molecular mechanism of the response of IbmiR319a to drought stress in sweet potato. Blocking IbmiR319a in transgenic sweet potato (MIM319) resulted in a slim and tender phenotype and greater sensitivity to drought stress. Microscopic observations revealed that blocking IbmiR319a decreased the cell width and increased the stomatal distribution in the adaxial leaf epidermis, and also increased the intercellular space in the leaf and petiole. We also found that the lignin content was reduced, which led to increased brittleness in MIM319. Quantitative real-time PCR showed that the expression levels of key genes in the lignin biosynthesis pathway were much lower in the MIM319 lines than in the wild type. Ectopic expression of IbmiR319a-targeted genes IbTCP11 and IbTCP17 in Arabidopsis resulted in similar phenotypes to MIM319. We also showed that the expression of IbTCP11 and IbTCP17 was largely induced by drought stress. Transcriptome analysis indicated that cell growth-related pathways, such as plant hormonal signaling, were significantly downregulated with the blocking of IbmiR319a. Taken together, our findings suggest that IbmiR319a affects plant architecture by targeting IbTCP11/17 to control the response to drought stress in sweet potato.

Molecular Characterization and Target Prediction of Candidate miRNAs Related to Abiotic Stress Responses and/or Storage Root Development in Sweet Potato

Sweet potato is a tuberous root crop with strong environmental stress resistance. It is beneficial to study its storage root formation and stress responses to identify sweet potato stress- and storage-root-thickening-related regulators. Here, six conserved miRNAs (miR156g, miR157d, miR158a-3p, miR161.1, miR167d and miR397a) and six novel miRNAs (novel 104, novel 120, novel 140, novel 214, novel 359 and novel 522) were isolated and characterized in sweet potato. Tissue-specific expression patterns suggested that miR156g, miR157d, miR158a-3p, miR167d, novel 359 and novel 522 exhibited high expression in fibrous roots or storage roots and were all upregulated in response to storage-root-related hormones (indole acetic acid, IAA; zeaxanthin, ZT; abscisic acid, ABA; and gibberellin, GAs). The expression of miR156g, miR158a-3p, miR167d, novel 120 and novel 214 was induced or reduced dramatically by salt, dehydration and cold or heat stresses. Moreover, these miRNAs were all upregulated by ABA, a crucial hormone modulator in regulating abiotic stresses. Additionally, the potential targets of the twelve miRNAs were predicted and analyzed. Above all, these results indicated that these miRNAs might play roles in storage root development and/or stress responses in sweet potato as well as provided valuable information for the further investigation of the roles of miRNA in storage root development and stress responses.

Genome-Wide Identification of TCP Transcription Factors Family in Sweet Potato Reveals Significant Roles of miR319-Targeted TCPs in Leaf Anatomical Morphology

Plant-specific TCP transcription factors play vital roles in the controlling of growth, development, and the stress response processes. Extensive researches have been carried out in numerous species, however, there hasn't been any information available about TCP genes in sweet potato (Ipomoea batatas L.). In this study, a genome-wide analysis of TCP genes was carried out to explore the evolution and function in sweet potato. Altogether, 18 IbTCPs were identified and cloned. The expression profiles of the IbTCPs differed dramatically in different organs or different stages of leaf development. Furthermore, four CIN-clade IbTCP genes contained miR319-binding sites. Blocking IbmiR319 significantly increased the expression level of IbTCP11/17 and resulted in a decreased photosynthetic rate due to the change in leaf submicroscopic structure, indicating the significance of IbmiR319-targeted IbTCPs in leaf anatomical morphology. A systematic analyzation on the characterization of the IbTCPs together with the primary functions in leaf anatomical morphology were conducted to afford a basis for further study of the IbmiR319/IbTCP module in association with leaf anatomical morphology in sweet potato.

Transcriptome Analysis Reveals Key Genes Involved in Weevil Resistance in the Hexaploid Sweetpotato

Because weevils are the most damaging pests of sweetpotato, the development of cultivars resistant to weevil species is considered the most important aspect in sweetpotato breeding. However, the genes and the underlying molecular mechanisms related to weevil resistance are yet to be elucidated. In this study, we performed an RNA sequencing-based transcriptome analysis using the resistant Kyushu No. 166 (K166) and susceptible Tamayutaka cultivars. The weevil resistance test showed a significant difference between the two cultivars at 30 days after the inoculation, specifically in the weevil growth stage and the suppressed weevil pupation that was only observed in K166. Differential expression and gene ontology analyses revealed that the genes upregulated after inoculation in K166 were related to phosphorylation, metabolic, and cellular processes. Because the weevil resistance was considered to be related to the suppression of larval pupation, we investigated the juvenile hormone (JH)-related genes involved in the inhibition of insect metamorphosis. We found that the expression of some terpenoid-related genes, which are classified as plant-derived JHs, was significantly increased in K166. This is the first study involving a comprehensive gene expression analysis that provides new insights about the genes and mechanisms associated with weevil resistance in sweetpotato.

Integration of Transcriptional and Post-transcriptional Analysis Revealed the Early Response Mechanism of Sugarcane to Cold Stress

Cold stress causes major losses to sugarcane production, yet the precise molecular mechanisms that cause losses due to cold stress are not well-understood. To survey miRNAs and genes involved in cold tolerance, RNA-seq, miRNA-seq, and integration analyses were performed on Saccharum spontaneum. Results showed that a total of 118,015 genes and 6,034 of these differentially expressed genes (DEGs) were screened. Protein–protein interaction (PPI) analyses revealed that ABA signaling via protein phosphatase 2Cs was the most important signal transduction pathway and late embryogenesis abundant protein was the hub protein associated with adaptation to cold stress. Furthermore, a total of 856 miRNAs were identified in this study and 109 of them were differentially expressed in sugarcane responding to cold stress. Most importantly, the miRNA–gene regulatory networks suggested the complex post-transcriptional regulation in sugarcane under cold stress, including 10 miRNAs−42 genes, 16 miRNAs−70 genes, and three miRNAs−18 genes in CT vs. LT0.5, CT vs. LT1, and CT0.5 vs. LT1, respectively. Specifically, key regulators from 16 genes encoding laccase were targeted by novel-Chr4C_47059 and Novel-Chr4A_40498, while five LRR-RLK genes were targeted by Novel-Chr6B_65233 and Novel-Chr5D_60023, 19 PPR repeat proteins by Novel-Chr5C_57213 and Novel-Chr5D_58065. Our findings suggested that these miRNAs and cell wall-related genes played vital regulatory roles in the responses of sugarcane to cold stress. Overall, the results of this study provide insights into the transcriptional and post-transcriptional regulatory network underlying the responses of sugarcane to cold stress.

Integrated transcriptome, small RNA and degradome sequencing approaches proffer insights into chlorogenic acid biosynthesis in leafy sweet potato

Leafy sweet potato is rich in total phenolics (TP) which play key roles in health protection, the chlorogenic acid (CGA) constitutes the major components of phenolic compounds in leafy sweet potato. Unfortunately, the mechanism of CGA biosynthesis in leafy sweet potato is unclear. To dissect the mechanisms of CGA biosynthesis, we performed transcriptome, small RNA (sRNA) and degradome sequencing of one low-CGA content and one high-CGA content genotype at two stages. A total of 2,333 common differentially expressed genes (DEGs) were identified, and the enriched DEGs were related to photosynthesis, starch and sucrose metabolism and phenylpropanoid biosynthesis. The functional genes, such as CCR, CCoAOMT and HCT in the CGA biosynthetic pathway were down-regulated, indicating that the way to lignin was altered, and two possible CGA biosynthetic routes were hypothesized. A total of 38 DE miRNAs were identified, and 1,799 targets were predicated for 38 DE miRNAs by using in silico approaches. The target genes were enriched in lignin and phenylpropanoid catabolic processes. Transcription factors (TFs) such as apetala2/ethylene response factor (AP2/ERF) and Squamosa promoter binding protein-like (SPL) predicated in silico were validated by degradome sequencing. Association analysis of the DE miRNAs and transcriptome datasets identified that miR156 family negatively targeted AP2/ERF and SPL. Six mRNAs and six miRNAs were validated by qRT-PCR, and the results showed that the expression levels of the mRNAs and miRNAs were consistent with the sequencing data. This study established comprehensive functional genomic resources for the CGA biosynthesis, and provided insights into the molecular mechanisms involving in this process. The results also enabled the first perceptions of the regulatory roles of mRNAs and miRNAs, and offered candidate genes for leafy sweet potato improvements.

Optimization of reference genes for qRT-PCR analysis of microRNA expression under abiotic stress conditions in sweetpotato

Sweetpotato (Ipomoea batatas. L) is an important food crop, harvested for its nutrient-rich tuberous roots. Drought and salt stresses are two major factors limiting the sweetpotato production. Since microRNAs (miRNAs) are well known to play crucial roles in regulation of plant stress responses, quantitative profiling of miRNA expression under stress conditions will facilitate identification and genetic manipulation of novel miRNAs to improve stress tolerance. Real-time quantitative reverse transcription PCR (qRT-PCR) is a commonly used tool for this purpose, but not without challenges. Although stem-loop and poly(A)-tail modified qRT-PCR methods were developed for characterizing miRNA expression, accurate profiling of miRNAs is still difficult in many plant species because of a lack of reliable reference genes for normalizing miRNA transcripts. To identify reference genes that are suitable for normalizing miRNA expression in sweetpotato, the expression stability of eight candidate miRNAs and two commonly used reference genes were tested in 96 samples involving four tissues and two cultivars under drought and salt stress treatments. Data analysis using the geNorm, NormFinder and Bestkeeper algorithms demonstrated that miRn60, miR482, and their combination were reliable references. We further validated the reference genes by expression analysis of the well-characterized miR319 and miR156 that regulate drought and salt stress responses, respectively. The reference genes identified in this study will facilitate future miRNA analysis under abiotic stress conditions in sweetpotato.

Identification of candidate miRNAs related in storage root development of sweet potato by high throughput sequencing

Sweet potato (Ipomoea batatas L.) is a food consumed worldwide, an industrial raw material and new energy crop. The storage root is the most economical part of the crop. However, the mechanism of storage root initiation and development is still unclear. In this study, conserved and novel miRNAs during storage root development were identified by high-throughput sequencing technology by constructing small RNA libraries from sweet potato fibrous roots (F) and storage roots at four different developmental stages (storage roots with different diameters: 1 cm, D1; 3 cm, D3; 5 cm, D5 and 10 cm, D10). A total of 61 known miRNAs and 471 novel miRNAs were identified. In addition, 145 differentially expressed miRNAs were identified in the F library compared with the four storage root libraries, with 30 known miRNAs and 115 novel miRNAs. Moreover, the targets of the differentially expressed miRNAs were predicted and their network was further investigated by GO analysis using our previous transcriptome data. The GO analysis revealed that antioxidant activity and binding process were the most enriched terms of the target genes. The secondary structure and expression of six candidate miRNAs including three conserved miRNAs and three novel miRNAs were investigated and their predicted targets were validated by qRT-PCR. The results showed that the expression levels of the miRNAs were all consistent with the sequencing data. Most of the miRNAs and their corresponding targets had obvious negative correlations. This study contributed to elucidating the potential miRNA mediated regulatory mechanism of storage root development in sweet potato. The specific differentially expressed miRNAs in sweet potato storage roots can be used to breed high-yield sweet potatoes and other tuberous root crops.

Selection of stable reference genes for gene expression analysis in sweet potato (Ipomoea batatas L.)

Gene expression analysis is one of the most common and important studies in biology and biomedicine. No matter for traditional blotting analysis or currently commonly used PCR strategy, all need a stable reference gene for normalizing the gene expression. To screen and select housekeeping genes as the most stable reference genes, quantitative real-time PCR (qRT-PCR) was employed to analyze the expression of sixteen commonly used reference genes (IbelF, Ibα-tubulin, IbHIS, IbCOX, IbGAPDH, IbH2B1, IbARF, IbCYC, Ibβ-tubulin, IbACT, IbEFl-a, IbG14, IbPLD, IbRPL2, IbUBQ, IbUBI) in five different tissues under two different temperature stresses in sweet potato. Data analysis by the Delta CT, geNorm, NormFinder, and BestKeeper programs revealed that IbelF is the most stable gene and IbUBI is the least stable gene as reference. Our study also shows that combination of two or more genes as reference is a better choice, rendering more substantiated expression data for comparison. This study provides evidence for selecting reference genes in sweet potato gene expression analysis.

sRIS: A Small RNA Illustration System for Plant Next-Generation Sequencing Data Analysis

Small RNA (sRNA), such as microRNA (miRNA) and short interfering RNA, are well-known to control gene expression based on degradation of target mRNA in plants. A considerable amount of research has applied next-generation sequencing (NGS) to reveal the regulatory pathways of plant sRNAs. Consequently, numerous bioinformatics tools have been developed for the purpose of analyzing sRNA NGS data. However, most methods focus on the study of sRNA expression profiles or novel miRNAs predictions. The analysis of sRNA target genes is usually not integrated into their pipelines. As a result, there is still no means available for identifying the interaction mechanisms between host and virus or the synergistic effects between two viruses. For the present study, a comprehensive system, called the Small RNA Illustration System (sRIS), has been developed. This system contains two main components. The first is for sRNA overview analysis and can be used not only to identify miRNA but also to investigate virus-derived small interfering RNA. The second component is for sRNA target prediction, and it employs both bioinformatics calculations and degradome sequencing data to enhance the accuracy of target prediction. In addition, this system has been designed so that figures and tables for the outputs of each analysis can be easily retrieved and accessed, making it easier for users to quickly identify and quantify their results. sRIS is available at http://sris.itps.ncku.edu.tw/.

Chilling and Heat Stress-Induced Physiological Changes and MicroRNA-Related Mechanism in Sweetpotato (Ipomoea batatas L.)

Sweetpotato (Ipomoea batatas (L.) Lam.) is an important industrial and food crop. Both chilling and heat stress inhibits sweetpotato growth and development and then affects yield. However, the physiological and molecular mechanisms of sweetpotato response to chilling and heat stress is unclear. In this study, we investigated the effect of extreme temperature on sweetpotato physiological response, with a focus on oxidative stress and the potential microRNA (miRNA)-mediated molecular mechanism. Our results showed that both chilling and heat stress resulted in accumulation of reactive oxygen species (ROS), including H₂O₂ and O₂ ^-, and caused oxidative stress in sweetpotato. This further affected the activities of oxidative stress-related enzymes and products, including SOD, POD, and MDA. Both chilling and heat stress inhibited POD activities but induced the enzyme activities of SOD and MDA. This suggests that sweetpotato cells initiated its own defense mechanism to handle extreme temperature-caused oxidative damage. Oxidative damage and repair are one mechanism that sweetpotato plants respond to extreme temperatures. Another potential mechanism is miRNA-mediated gene response. Chilling and heat stress altered the expression of stress-responsive miRNAs in sweetpotato seedlings. These miRNAs regulate sweetpotato response to extreme stress through targeting individual protein-coding genes.

Antioxidative capacity is highly associated with the storage property of tuberous roots in different sweetpotato cultivars

The activities and gene expression of antioxidative enzymes and the ROS content were analyzed in two typical storage-tolerant cultivars (Xushu 32 and Shangshu 19) and another two storage-sensitive cultivars (Yanshu 25 and Sushu 16) to explore the association between the storage capacity of sweetpotato (Ipomoea batatas (L.) Lam) and ROS scavenging capability. The storage roots of the storage-tolerant cultivars maintained higher activities and expression levels of antioxidative enzymes, including ascorbate peroxidase (APX), peroxidase (POD), catalase (CAT), and superoxide dismutase (SOD); lower activity and expression of lipoxygenase (LOX); and lower accumulation of ROS metabolites compared with the storage-sensitive cultivars. The antioxidative capability and ROS parameters of leaves were positively correlated with those of storage roots. Our results provide valuable insight for evaluating the storability of sweetpotato cultivars by analyzing the capabilities of the antioxidative system and the contents of ROS metabolites.

Uncovering anthocyanin biosynthesis related microRNAs and their target genes by small RNA and degradome sequencing in tuberous roots of sweetpotato

BACKGROUND

Compared with white-fleshed sweetpotato (WFSP), purple-fleshed sweetpotato (PFSP) is a desirable resource for functional food development because of the abundant anthocyanin accumulation in its tuberous roots. Some studies have shown that the expression regulation mediated by miRNA plays an important role in anthocyanin biosynthesis in plants. However, few miRNAs and their corresponding functions related to anthocyanin biosynthesis in tuberous roots of sweetpotato have been known.

RESULTS

In this study, small RNA (sRNA) and degradome libraries from the tuberous roots of WFSP (Xushu-18) and PFSP (Xuzishu-3) were constructed, respectively. Totally, 191 known and 33 novel miRNAs were identified by sRNA sequencing, and 180 target genes cleaved by 115 known ib-miRNAs and 5 novel ib-miRNAs were identified by degradome sequencing. Of these, 121 miRNAs were differently expressed between Xushu-18 and Xuzishu-3. Integrated analysis of sRNA, degradome sequencing, GO, KEGG and qRT-PCR revealed that 26 differentially expressed miRNAs and 36 corresponding targets were potentially involved in the anthocyanin biosynthesis. Of which, an inverse correlation between the expression of ib-miR156 and its target ibSPL in WFSP and PFSP was revealed by both qRT-PCR and sRNA sequencing. Subsequently, ib-miR156 was over-expressed in Arabidopsis. Interestingly, the ib-miR156 over-expressing plants showed suppressed abundance of SPL and a purplish phenotype. Concomitantly, upregulated expression of four anthocyanin pathway genes was detected in transgenic Arabidopsis plants. Finally, a putative ib-miRNA-target model involved in anthocyanin biosynthesis in sweetpotato was proposed.

CONCLUSIONS

The results represented a comprehensive expression profiling of miRNAs related to anthocyanin accumulation in sweetpotato and provided important clues for understanding the regulatory network of anthocyanin biosynthesis mediated by miRNA in tuberous crops.

Antioxidative system in sweet potato root is activated by low-temperature storage

BACKGROUND

Sweet potato is susceptible to chilling injury during low-temperature storage. To explore the correlation between chilling injury and reactive oxygen species (ROS) metabolism, the content of ROS and the activities and gene expression of antioxidant enzymes were analyzed in the typical storage-tolerant cultivar Xushu 32 and storage-sensitive cultivar Yanshu 25.

RESULTS

The activities of antioxidant enzymes including ascorbate peroxidase (APX), superoxide dismutase (SOD), catalase (CAT) and glutathione reductase (GR) were enhanced rapidly in the early period of storage in response to chilling stress. Thereafter, the content of ROS metabolites increased consistently due to gradual decrease in ROS scavenging enzymes. Storage-tolerant cultivar Xushu 32 had higher antioxidant enzyme activities and gene expressions as well as higher content of antioxidant metabolites and lower content of ROS metabolites compared with storage-sensitive cultivar Yanshu 25, suggesting that the capacity of ROS scavenging by antioxidant enzymes and antioxidants is highly associated with the tolerance of sweet potato to chilling stress.

CONCLUSION

These results indicated that the antioxidative system is activated in the storage root of sweet potato and the antioxidative capacity is positively associated with better storage performance in the storage-tolerant cultivar. © 2019 Society of Chemical Industry.