Comparative analysis of the root transcriptomes of cultivated sweetpotato (Ipomoea batatas [L.] Lam) and its wild ancestor (Ipomoea trifida [Kunth] G. Don)

Regulation of root development in nitrogen-susceptible and nitrogen-tolerant sweet potato cultivars under different nitrogen and soil moisture conditions

BACKGROUND

Due to unreasonable nitrogen (N) application and water supply, sweet potato vines tend to grow excessively. Early development of storage roots is conducive to inhibiting vine overgrowth. Hence, we investigated how N and soil moisture affect early root growth and development.

RESULTS

A pot experiment was conducted using the sweet potato cultivars Jishu26 (J26, N-susceptible) and Xushu32 (X32, N-tolerant). Two N application rates of 50 (N1) and 150 mg kg- 1 (N2) and two water regimes, drought stress (DS) (W1) and normal moisture (W2), were applied to each cultivar. For J26, the lowest expansion root weight was observed in the N2W2 treatment, while for X32, the N1W2 and N2W2 treatments resulted in higher root weights compared to other treatments. The interaction between N rates and water regimes significantly affected root surface area and volume in J26. Root cross-sections revealed that N2W2 increased the percentage of root area covered by xylem vessels and decreased the amount of secondary xylem vessels (SXV) in J26. However, in X32, it increased the number of SXV. A high N rate reduced the 13 C distribution ratio in J26 expansion roots, but had no significant effect on X32. In J26, N2W2 inhibited starch synthesis in roots by downregulating the expression of AGPa, AGPb, GBSS I, and SBE I.

CONCLUSION

The observed effects were more pronounced in J26. For X32, relatively high N and moisture levels did not significantly impact storage root development. Therefore, special attention should be paid to N supply and soil moisture for N-susceptible cultivars during the early growth stage.

Form and contour: breeding and genetics of organ shape from wild relatives to modern vegetable crops

Shape is a primary determinant of consumer preference for many horticultural crops and it is also associated with many aspects of marketing, harvest mechanics, and postharvest handling. Perceptions of quality and preference often map to specific shapes of fruits, tubers, leaves, flowers, roots, and other plant organs. As a result, humans have greatly expanded the palette of shapes available for horticultural crops, in many cases creating a series of market classes where particular shapes predominate. Crop wild relatives possess organs shaped by natural selection, while domesticated species possess organs shaped by human desires. Selection for visually-pleasing shapes in vegetable crops resulted from a number of opportunistic factors, including modification of supernumerary cambia, allelic variation at loci that control fundamental processes such as cell division, cell elongation, transposon-mediated variation, and partitioning of photosynthate. Genes that control cell division patterning may be universal shape regulators in horticultural crops, influencing the form of fruits, tubers, and grains in disparate species. Crop wild relatives are often considered less relevant for modern breeding efforts when it comes to characteristics such as shape, however this view may be unnecessarily limiting. Useful allelic variation in wild species may not have been examined or exploited with respect to shape modifications, and newly emergent information on key genes and proteins may provide additional opportunities to regulate the form and contour of vegetable crops.

mRNA-Seq and miRNA-Seq Analyses Provide Insights into the Mechanism of Pinellia ternata Bulbil Initiation Induced by Phytohormones

Pinellia ternata (Thunb.) Breit (abbreviated as P. ternata) is a plant with an important medicinal value whose yield is restricted by many factors, such as low reproductive efficiency and continuous cropping obstacles. As an essential breeding material for P. ternata growth and production, the bulbils have significant advantages such as a high survival rate and short breeding cycles. However, the location effect, influencing factors, and molecular mechanism of bulbil occurrence and formation have not been fully explored. In this study, exogenously applied phytohormones were used to induce in vitro petiole of P. ternata to produce bulbil structure. Transcriptome sequencing of mRNA and miRNA were performed in the induced petiole (TCp) and the induced bulbil (TCb). Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed for the identification of key genes and pathways involved in bulbil development. A total of 58,019 differentially expressed genes (DEGs) were identified. The GO and KEGG analysis indicated that DEGs were mainly enriched in plant hormone signal transduction and the starch and sucrose metabolism pathway. The expression profiles of miR167a, miR171a, and miR156a during bulbil induction were verified by qRT-PCR, indicating that these three miRNAs and their target genes may be involved in the process of bulbil induction and play an important role. However, further molecular biological experiments are required to confirm the functions of the identified bulbil development-related miRNAs and targets.

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

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

Comparative transcriptomics analysis pipeline for the meta-analysis of phylogenetically divergent datasets (CoRMAP)

Background

Transcriptional regulation is a fundamental mechanism underlying biological functions. In recent years, a broad array of RNA-Seq tools have been used to measure transcription levels in biological experiments, in whole organisms, tissues, and at the single cell level. Collectively, this is a vast comparative dataset on transcriptional processes across organisms. Yet, due to technical differences between the studies (sequencing, experimental design, and analysis) extracting usable comparative information and conducting meta-analyses remains challenging.

Results

We introduce Comparative RNA-Seq Metadata Analysis Pipeline (CoRMAP), a meta-analysis tool to retrieve comparative gene expression data from any RNA-Seq dataset using de novo assembly, standardized gene expression tools and the implementation of OrthoMCL, a gene orthology search algorithm. It employs the use of orthogroup assignments to ensure the accurate comparison of gene expression levels between experiments and species. Here we demonstrate the use of CoRMAP on two mouse brain transcriptomes with similar scope, that were collected several years from each other using different sequencing technologies and analysis methods. We also compare the performance of CoRMAP with a functional mapping tool, previously published.

Conclusion

CoRMAP provides a framework for the meta-analysis of RNA-Seq data from divergent taxonomic groups. This method facilitates the retrieval and comparison of gene expression levels from published data sets using standardized assembly and analysis. CoRMAP does not rely on reference genomes and consequently facilitates direct comparison between diverse studies on a range of organisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04972-9.

Comparative Transcriptome Analyses of Different Rheum officinale Tissues Reveal Differentially Expressed Genes Associated with Anthraquinone, Catechin, and Gallic Acid Biosynthesis

Rheum officinale Baill. is an important traditional Chinese medicinal herb, its dried roots and rhizomes being widely utilized to cure diverse diseases. However, previous studies mainly focused on the active compounds and their pharmacological effects, and the molecular mechanism underlying the biosynthesis of these ingredients in R. officinale is still elusive. Here, we performed comparative transcriptome analyses to elucidate the differentially expressed genes (DEGs) in the root, stem, and leaf of R. officinale. A total of 236,031 unigenes with N50 of 769 bp was generated, 136,329 (57.76%) of which were annotated. A total of 5884 DEGs was identified after the comparative analyses of different tissues; 175 and 126 key enzyme genes with tissue-specific expression were found in the anthraquinone, catechin/gallic acid biosynthetic pathway, respectively, and some of these key enzyme genes were verified by qRT-PCR. The phylogeny of the PKS III family in Polygonaceae indicated that probably only PL_741 PKSIII1, PL_11549 PKSIII5, and PL_101745 PKSIII6 encoded PKSIII in the polyketide pathway. These results will shed light on the molecular basis of the tissue-specific accumulation and regulation of secondary metabolites in R. officinale, and lay a foundation for the future genetic diversity, molecular assisted breeding, and germplasm resource improvement of this essential medicinal plant.

Comparative Transcriptome Profiling Reveals the Genes Involved in Storage Root Expansion in Sweetpotato (Ipomoea batatas (L.) Lam.)

Sweetpotato (Ipomoea batatas (L.) Lam.) is recognized as one of the most important root crops in the world by the Food and Agriculture Organization of the United Nations. The yield of sweetpotato is closely correlated with the rate of storage root (SR) formation and expansion. At present, most of the studies on sweetpotato SR expansion are focused on the physiological mechanism. To explore the SR expansion mechanism of sweetpotato, we performed transcriptome sequencing of SR harvested at 60, 90, 120, and 150 days after planting (DAP) to analyze two sweetpotato lines, Xuzishu 8 and its crossing progenies named Xu 18-192, which were selected from an F1 segregation population of Xuzishu 8 and Meiguohong, in which SR expansion was delayed significantly. A total of 57,043 genes were produced using transcriptome sequencing, of which 1312 were differentially expressed genes (DEGs) in four SR growth periods of the sweetpotato lines. The combination of the KEGG and trend analysis revealed several key candidate genes involved in SR expansion. The SBEI gene involved in starch metabolism, and transcription factors ARF6, NF-YB3 and NF-YB10 were all significantly up-regulated during SR expansion. The data from this study provide insights into the complex mechanisms of SR formation and expansion in sweetpotato and identify new candidate genes for increasing the yield of sweetpotato.

Genome-Wide Identification and Expression Analysis of Expansin Gene Family in the Storage Root Development of Diploid Wild Sweetpotato Ipomoea trifida

Expansins play important roles in root growth and development, but investigation of the expansin gene family has not yet been reported in Ipomoea trifida, and little is known regarding storage root (SR) development. In this work, we identified a total of 37 expansins (ItrEXPs) in our previously reported SR-forming I. trifida strain Y22 genome, which included 23 ItrEXPAs, 4 ItrEXPBs, 2 ItrEXLAs and 8 ItrEXLBs. The phylogenetic relationship, genome localization, subcellular localization, gene and protein structure, promoter cis-regulating elements, and protein interaction network were systematically analyzed to reveal the possible roles of ItrEXPs in the SR development of I. trifida. The gene expression profiling in Y22 SR development revealed that ItrEXPAs and ItrEXLBs were down-regulated, and ItrEXPBs were up-regulated while ItrEXLAs were not obviously changed during the critical period of SR expansion, and might be beneficial to SR development. Combining the tissue-specific expression in young SR transverse sections of Y22 and sweetpotato tissue, we deduced that ItrEXLB05, ItrEXLB07 and ItrEXLB08 might be the key genes for initial SR formation and enlargement, and ItrEXLA02 might be the key gene for root growth and development. This work provides new insights into the functions of the expansin gene family members in I. trifida, especially for EXLA and EXLB subfamilies genes in SR development.

Differential response of physiology and metabolic response to drought stress in different sweetpotato cultivars

Sweetpotato (Ipomoea batatas [L.] Lam) is a widely cultivated food crop with generally good adaptability. However, drought stress can cause a significant decline in yield. To reveal the response mechanism of sweetpotato to drought stress, an integrated physiological, proteomic and metabolomic investigation was conducted in leaves of two sweetpotato varieties with differing responses to drought stress, drought-resistant Wanzishu56 (WZ56) and a more sensitive variety, Ningzishu2(NZ2). Physiological analysis showed that the variety with better drought tolerance had superior performance in water retention capacity and photosynthetic efficiency under drought stress. A total of 1140 proteins were identified within the two varieties. Among them, 192 differentially expressed proteins were detected under drought conditions, including 97 that were up-regulated. Functional analysis showed that these up-regulated proteins were primarily involved in photosynthesis, reactive oxygen species metabolism, organonitrogen compound metabolism, and precursor metabolite catabolism and energy generation. All differentially expressed proteins in WZ56 that were involved in photosynthetic and glutathione metabolic processes were up-regulated. Enzyme activity assays were carried out to validate the proteomics data. Moreover, 75 metabolites were found to have a higher expression level in WZ56 than NZ2 under drought stress. The higher concentration of carbohydrates, amino acids, flavonoids and organic acids found in drought-stressed leaves of WZ56 suggested that these metabolites may improve the drought resistance of sweetpotato. This study uncovered specific-proteins and metabolites associated with drought resistance, providing new insights into the molecular mechanisms of drought tolerance in sweetpotato.

Lignin Synthesis, Affected by Sucrose in Lotus (Nelumbo nucifera) Seedlings, Was Involved in Regulation of Root Formation in the Arabidopsis thanliana

Adventitious roots (ARs) have an unmatched status in plant growth and metabolism due to the degeneration of primary roots in lotuses. In the present study, we sought to assess the effect of sucrose on ARs formation and observed that lignin synthesis was involved in ARs development. We found that the lignification degree of the ARs primordium was weaker in plants treated with 20 g/L sucrose than in 50 g/L sucrose treatment and control plants. The contents of lignin were lower in plants treated with 20 g/L sucrose and higher in plants treated with 50 g/L sucrose. The precursors of monomer lignin, including p-coumaric acid, caffeate, sinapinal aldehyde, and ferulic acid, were lower in the GL50 library than in the GL20 library. Further analysis revealed that the gene expression of these four metabolites had no novel difference in the GL50/GL20 libraries. However, a laccase17 gene (NnLAC17), involved in polymer lignin synthesis, had a higher expression in the GL50 library than in the GL20 library. Therefore, NnLAC17 was cloned and the overexpression of NnLAC17 was found to directly result in a decrease in the root number in transgenic Arabidopsis plants. These findings suggest that lignin synthesis is probably involved in ARs formation in lotus seedlings.

Nitrogen utilization characteristics and early storage root development in nitrogen-tolerant and nitrogen-susceptible sweet potato

In recent years, sweet potato has been cultivated not only in marginal lands but also in fertile plains in northern China. The fertile nitrogen (N)-rich soil may inhibit storage root formation. Cultivars with different N tolerances and split application of reduced N rates should be considered. To investigate the effects of N on the N utilization, root differentiation, and storage root formation of cultivars with different N tolerances, the cultivars Jishu26 (J26) and Xushu32 (X32) were treated with three N levels supplied by urea: 0 (N0), 200 (N1) and 400 mg kg^-1 (N2). With increasing N rates, "X32" absorbed less N in plants and distributed more N to developing storage roots than "J26." The storage root development of "J26" was sensitive to both N1 and N2, while that of "X32" was only sensitive to N2. High N nutrition upregulated the expression of certain genes during storage root formation, such as PAL, CHI, F3H, C₄ H, 4CL, CAD, α-amylase, and β-amylase. Under N1 and N2, "X32" led to an increased sugar supply in sink organs and downregulated the expression of genes related to lignin and flavonoid synthesis, which promoted the C flux toward starch metabolism, thus reducing lignification and promoting starch accumulation during storage root formation. These results provide evidence for the effects of N on the C distribution in different metabolic pathways by regulating the expression of related key genes. N-tolerant cultivars are suitable in fertile plain areas because of the earlier formation of storage roots under high N conditions.

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.

Full-length transcriptome analysis provides new insights into the early bolting occurrence in medicinal Angelica sinensis

Angelica sinensis (Oliv.) Diels root part is an integral component of traditional Chinese medicine, widely prescribed to improve blood circulation and blood stasis. However, early bolting of A. sinensis compromises the quality of the roots and hence is a major limitation for yield of medicinal materials. To date, little information about the molecular mechanisms underlying bolting is available for this important medicinal plant. To identify genes putatively involved in early bolting, we have conducted the transcriptome analysis of the shoot tips of the early-bolting plants and non-bolting (normal) plants of A. sinensis, respectively, using a combination of third-generation sequencing and next-generation sequencing. A total of 43,438 non-redundant transcripts were collected and 475 unique differentially expressed genes (DEGs) were identified. Gene annotation and functional analyses revealed that DEGs were highly involved in plant hormone signaling and biosynthesis pathways, three main flowering pathways, pollen formation, and very-long-chain fatty acids biosynthesis pathways. The levels of endogenous hormones were also changed significantly in the early bolting stage of A. sinensis. This study provided new insights into the transcriptomic control of early bolting in A. sinensis, which could be further applied to enhance the yield of medicinally important raw materials.

Seasonal Variation in Transcriptomic Profiling of Tetrastigma hemsleyanum Fully Developed Tuberous Roots Enriches Candidate Genes in Essential Metabolic Pathways and Phytohormone Signaling

Tetrastigma hemsleyanum Diels et Gilg (Sanyeqing, SYQ) is a perennial climbing liana and an endemic plant to southern China. Its tuberous roots (TRs) are used in traditional Chinese medicine for treating some diseases such as high fever, pneumonia, asthma, hepatitis, and cancers. However, the mechanisms underlying the development of TR and the content of flavonoids and phenylpropanoids (FPs) are not well-understood. In this study, we performed a transcriptomic analysis of 12 fully developed TR (FD-TR) samples harvested in four seasons [spring (Sp), summer (Su), autumn (Au), and winter (Wi)] using the RNA-Sequencing (RNA-Seq). We obtained a total of 78.54 Gb raw data and 65,578 unigenes. Then, the unigenes were annotated by using six databases such as non-redundant protein database (NR), Pfam, eggNOG, SWISSProt, Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene ontology (GO). The transcriptomic profiling showed closer relationships between the samples obtained in Su and Au than those obtained in Sp and Wi based on the results of both total unigenes and differentially expressed genes (DEGs). Three pathways, including the biosynthesis of FPs, metabolism of starch and sucrose, and signaling of phytohormones, were highly enriched, suggesting a gene-level seasonal variation. Based on the numbers of DEGs, brassinosteroid (BR) signal transduction factors appeared to play a key role in modulating the development of TRs while most of the auxin signaling genes were mainly activated in Wi and Sp FD-TRs. Most genes in the biosynthesis and biodegradation of starch and biodegradation of cellulose were activated in Wi FD-TRs. As determined by the high performance liquid chromatography (HPLC) and aluminum nitrate colorimetric method, the contents of total flavonoids and most detected FP components increased from Sp to Au but decreased in Wi. Enhanced expression levels of some genes in the biosynthetic pathways of FPs were detected in Su and Au samples, which corroborated well with metabolite content. Our findings provide the first transcriptomic and biochemical data on a seasonal variation in the composition of medically important metabolites in SYQ FD-TRs.

A high-density genetic map and multiple environmental tests reveal novel quantitative trait loci and candidate genes for fibre quality and yield in cotton

A high-density linkage map of an intraspecific RIL population was constructed using 6187 bins to identify QTLs for fibre quality- and yield-related traits in upland cotton by whole-genome resequencing. Good fibre quality and high yield are important production goals in cotton (Gossypium hirsutum L.), which is a leading natural fibre crop worldwide. However, a greater understanding of the genetic variants underlying fibre quality- and yield-related traits is still required. In this study, a large-scale population including 588 F₇ recombinant inbred lines, derived from an intraspecific cross between the upland cotton cv. Nongdamian13, which exhibits high quality, and Nongda601, which exhibits a high yield, was genotyped by using 232,946 polymorphic single-nucleotide polymorphisms obtained via a whole-genome resequencing strategy with 4.3-fold genome coverage. We constructed a high-density bin linkage map containing 6187 bin markers spanning 4478.98 cM with an average distance of 0.72 cM. We identified 58 individual quantitative trait loci (QTLs) and 25 QTL clusters harbouring 94 QTLs, and 119 previously undescribed QTLs controlling 13 fibre quality and yield traits across eight environments. Importantly, the QTL counts for fibre quality in the Dt subgenome were more than two times that in the At subgenome, and chromosome D02 harboured the greatest number of QTLs and clusters. Furthermore, we discovered 24 stable QTLs for fibre quality and 12 stable QTLs for yield traits. Four novel major stable QTLs related to fibre length, fibre strength and lint percentage, and seven previously unreported candidate genes with significantly differential expression between the two parents were identified and validated by RNA-seq. Our research provides valuable information for improving the fibre quality and yield in cotton breeding.

Integrative transcriptomics reveals genotypic impact on sugar beet storability

An integrative comparative transcriptomic approach on six sugar beet varieties showing different amount of sucrose loss during storage revealed genotype-specific main driver genes and pathways characterizing storability. Sugar beet is next to sugar cane one of the most important sugar crops accounting for about 15% of the sucrose produced worldwide. Since its processing is increasingly centralized, storage of beet roots over an extended time has become necessary. Sucrose loss during storage is a major concern for the sugar industry because the accumulation of invert sugar and byproducts severely affect sucrose manufacturing. This loss is mainly due to ongoing respiration, but changes in cell wall composition and pathogen infestation also contribute. While some varieties can cope better during storage, the underlying molecular mechanisms are currently undiscovered. We applied integrative transcriptomics on six varieties exhibiting different levels of sucrose loss during storage. Already prior to storage, well storable varieties were characterized by a higher number of parenchyma cells, a smaller cell area, and a thinner periderm. Supporting these findings, transcriptomics identified changes in genes involved in cell wall modifications. After 13 weeks of storage, over 900 differentially expressed genes were detected between well and badly storable varieties, mainly in the category of defense response but also in carbohydrate metabolism and the phenylpropanoid pathway. These findings were confirmed by gene co-expression network analysis where hub genes were identified as main drivers of invert sugar accumulation and sucrose loss. Our data provide insight into transcriptional changes in sugar beet roots during storage resulting in the characterization of key pathways and hub genes that might be further used as markers to improve pathogen resistance and storage properties.

Multiple QTL Mapping in Autopolyploids: A Random-Effect Model Approach with Application in a Hexaploid Sweetpotato Full-Sib Population

In developing countries, the sweetpotato, Ipomoea batatas (L.) Lam. [Formula: see text], is an important autopolyploid species, both socially and economically. However, quantitative trait loci (QTL) mapping has remained limited due to its genetic complexity. Current fixed-effect models can fit only a single QTL and are generally hard to interpret. Here, we report the use of a random-effect model approach to map multiple QTL based on score statistics in a sweetpotato biparental population ('Beauregard' × 'Tanzania') with 315 full-sibs. Phenotypic data were collected for eight yield component traits in six environments in Peru, and jointly adjusted means were obtained using mixed-effect models. An integrated linkage map consisting of 30,684 markers distributed along 15 linkage groups (LGs) was used to obtain the genotype conditional probabilities of putative QTL at every centiMorgan position. Multiple interval mapping was performed using our R package QTLpoly and detected a total of 13 QTL, ranging from none to four QTL per trait, which explained up to 55% of the total variance. Some regions, such as those on LGs 3 and 15, were consistently detected among root number and yield traits, and provided a basis for candidate gene search. In addition, some QTL were found to affect commercial and noncommercial root traits distinctly. Further best linear unbiased predictions were decomposed into additive allele effects and were used to compute multiple QTL-based breeding values for selection. Together with quantitative genotyping and its appropriate usage in linkage analyses, this QTL mapping methodology will facilitate the use of genomic tools in sweetpotato breeding as well as in other autopolyploids.

Low Nitrogen Fertilization Alter Rhizosphere Microorganism Community and Improve Sweetpotato Yield in a Nitrogen-Deficient Rocky Soil

Sweetpotato can be cultivated in the reclaimed rocky soil in Sichuan Basin, China, which benefits from the release of mineral nutrients in the rocky soil by microorganisms. Shortage of nitrogen (N) in the rocky soil limits sweetpotato yield, which can be compensated through N fertilization. Whereas high N fertilization inhibits biological N fixation and induces unintended environmental consequences. However, the effect of low N fertilization on microorganism community and sweetpotato yield in the N-deficient rocky soil is still unclear. We added a low level of 1.5 g urea/m² to a rocky soil cultivated with sweetpotato, and measured rocky soil physiological and biochemical properties, rhizosphere microbial diversity, sweetpotato physiological properties and transcriptome. When cultivating sweetpotato in the rocky soil, low N fertilization (1.5 g urea/m²) not only improved total N (TN) and available N (AN) in the rocky soil, but also increased available phosphorus (AP), available potassium (AK), and nitrogenase and urease activity. Interestingly, although low N fertilization could reduce bacterial diversity through affecting sweetpotato root exudates and rocky soil properties, the relative abundance of P and K-solubilizing bacteria, N-fixing and urease-producing bacteria increased under low N fertilization, and the relative abundance of plant pathogens decreased. Furthermore, low N fertilization increased the phytohormones, such as zeatin riboside, abscisic acid, and methyl jasmonate contents in sweetpotato root. Those increases were consistent with our transcriptome findings: the inhibition of the lignin synthesis, the promotion of the starch synthesis, and the upregulated expression of Expansin, thus resulting in promoting the formation of tuberous roots and further increasing the sweetpotato yield by half, up to 3.3 kg/m². This study indicated that low N fertilization in the N-deficient rocky soil improved this soil quality through affecting microorganism community, and further increased sweetpotato yield under regulation of phytohormones pathway.

Functional genomics by integrated analysis of transcriptome of sweet potato (Ipomoea batatas (L.) Lam.) during root formation

BACKGROUND

Sweet potato is easily propagated by cuttings. But the molecular biological mechanism of adventitious root formation are not yet clear.

OBJECTIVE

To understand the molecular mechanisms of adventitious root formation from stem cuttings in sweet potato.

METHODS

RNA-seq analysis was performed using un-rooted stem (0 day) and rooted stem (3 days). Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, comparison with Arabidopsis transcription factors (TFs) of DEGs were conducted to investigate the characteristics of genes and TFs involved in root formation. In addition, qRT-PCR analysis using roots at 0, 3, 6, 9, and 12 days after planting was performed to confirm RNA-seq reliability and related genes expression.

RESULTS

42,459 representative transcripts and 2092 DEGs were obtained through the RNA-seq analysis. The DEGs indicated the GO terms related to the single-organism metabolic process and cell periphery, and involved in the biosynthesis of secondary metabolites, and phenylpropanoid biosynthesis in KEGG pathways. The comparison with Arabidopsis thaliana TF database showed that 3 TFs (WRKY, NAC, bHLH) involved in root formation of sweet potato. qRT-PCR analysis, which was conducted to confirm the reliability of RNA-seq analysis, indicated that some metabolisms including oxidative stress and wounding, transport, hormone may be involved in adventitious root formation.

CONCLUSIONS

The detected genes related to secondary metabolism, some hormone (auxin, gibberellin), transports, etc. and 3 TFs (WRKY, NAC, bHLH) may have functions in adventitious roots formation. This results provide valuable resources for future research on the adventitious root formation of sweet potato.

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.

Full-length transcriptome analysis of shade-induced promotion of tuber production in Pinellia ternata

BACKGROUND

Pinellia ternata is native to China and has been used as a traditional herb due to its antiemetic, antitussive, analgesic, and anxiolytic effects. When exposed to strong light intensity and high temperature during the reproductive growth process, P. ternata withers in a phenomenon known as "sprout tumble", which largely limits tuber production. Shade was previously found to delay sprout tumble formation (STF); however, no information exists regarding this process at the molecular level. Hence, we determined the genes involved in tuber development and STF in P. ternata.

RESULTS

Compared to that with natural sun-light (control), shade significantly induced chlorophyll accumulation, increased chlorophyll fluorescence parameters including initial fluorescence, maximal fluorescence, and qP, and dramatically repressed chlorophyll a:b and NPQ. Catalase (CAT) activity was largely induced by shade, and tuber products were largely increased in this environment. Transcriptome profiles of P. ternata grown in natural sun-light and shaded environments were analyzed by a combination of next generation sequencing (NGS) and third generation single-molecule real-time (SMRT) sequencing. Corrections of SMRT long reads based on NGS short reads yielded 136,163 non-redundant transcripts, with an average N50 length of 2578 bp. In total, 6738 deferentially-expressed genes (DEGs) were obtained from the comparisons, specifically D5S vs D5CK, D20S vs D20CK, D20S vs D5S, and D20CK vs D5CK, of which, 6384 DEGs (94.8%) were generated from the D20S vs D20CK comparison. Gene annotation and functional analyses revealed that these genes were related to auxin signal transduction, polysaccharide and sugar metabolism, phenylpropanoid biosynthesis, and photosynthesis. Moreover, the expression of genes enriched in photosynthesis appeared to be significantly altered by shade. The expression patterns of 16 candidate genes were consistent with changes in their transcript abundance as identified by RNA-Seq, and these might contribute to STF and tuber production.

CONCLUSION

The full-length transcripts identified in this study have provided a more accurate depiction of P. ternata gene transcription. Further, we identified potential genes involved in STF and tuber growth. Such data could serve as a genetic resource and a foundation for further research on this important traditional herb.

Genome-wide genetic diversity detection and population structure analysis in sweetpotato (Ipomoea batatas) using RAD-seq

Sweetpotato (Ipomoea batatas L.) is one of the most important food and grain-forage crops globally. It has been planted in >100 countries. Due to the complexity of the sweetpotato genome, its research is far behind other major food crops. At present, limited information about the sweetpotato genome is available. Thus, it is central to find an efficient approach for the investigation of sweetpotato genome. In this study, RAD-seq (Restriction site-associated DNA sequencing) was used to evaluate sweetpotato genetic structure diversity and to develop relevant SSR markers. The study yielded >128 Gb reliable sequence data from 81 sweetpotato accessions. By analyzing polymorphic tags from each accession, a total of 55,622 restriction-site associated DNA sequencing tags (RAD-seq) were found, containing 907,010 SNP. Genetic analysis divided 81 accessions into five major clusters based on their SNP genotype, which matches the results of genetic analysis and the genetic family tree. In addition, 18,320 SSRs loci were detected and 9336 SSR primer pairs were developed. Eighty-three primer pairs were amplified in different sweetpotato genotypes, 76 of which successfully amplified polymorphism bands. These results provide significant information about sweetpotato genome, which can be used to identify novel gene and to further develop the gene chip. And more significant, clustering results based on the SNP genotype provide an essential reference for breeders to match parent plants in breeding program. Additionally, SSR markers developed in this study will supply a wealth of markers for marker-assisted selection in sweetpotato breeding.

Single-molecule real-time sequencing identifies massive full-length cDNAs and alternative-splicing events that facilitate comparative and functional genomics study in the hexaploid crop sweet potato

Background

Sweet potato (Ipomoea batatas (L.) Lam.) is one of the most important crops in many developing countries and provides a candidate source of bioenergy. However, neither a complete reference genome nor large-scale full-length cDNA sequences for this outcrossing hexaploid crop are available, which in turn impedes progress in research studies in I. batatas functional genomics and molecular breeding.

Methods

In this study, we sequenced full-length transcriptomes in I. batatas and its diploid ancestor I. trifida by single-molecule real-time sequencing and Illumina second-generation sequencing technologies. With the generated datasets, we conducted comprehensive intraspecific and interspecific sequence analyses and experimental characterization.

Results

A total of 53,861/51,184 high-quality long-read transcripts were obtained, which covered about 10,439/10,452 loci in the I. batatas/I. trifida genome. These datasets enabled us to predict open reading frames successfully in 96.83%/96.82% of transcripts and identify 34,963/33,637 full-length cDNA sequences, 1,401/1,457 transcription factors, 25,315/27,090 simple sequence repeats, 1,656/1,389 long non-coding RNAs, and 5,251/8,901 alternative splicing events. Approximately, 32.34%/38.54% of transcripts and 46.22%/51.18% multi-exon transcripts underwent alternative splicing in I. batatas/I. trifida. Moreover, we validated one alternative splicing event in each of 10 genes and identified tuberous-root-specific expressed isoforms from a starch-branching enzyme, an alpha-glucan phosphorylase, a neutral invertase, and several ABC transporters. Overall, the collection and analysis of large-scale long-read transcripts generated in this study will serve as a valuable resource for the I. batatas research community, which may accelerate the progress in its structural, functional, and comparative genomics studies.

Molecular Dating of the Emergence of Anaerobic Rumen Fungi and the Impact of Laterally Acquired Genes

The anaerobic gut fungi (AGF), or Neocallimastigomycota, inhabit the rumen and alimentary tract of herbivorous mammals, where they play important roles in the degradation of plant fiber. Comparative genomic and phylogenomic analyses of the AGF have long been hampered by their fastidious growth condition, as well as their large (up to 200 Mb) and AT-biased (78 to 84%) genomes. We sequenced 21 AGF transcriptomes and combined them with 5 available AGF genome sequences to explore their evolutionary relationships, time their divergence, and characterize gene gain/loss patterns associated with their evolution. We estimate that the most recent common ancestor of the AGF diverged 66 (±10) million years ago, a time frame that coincides with the evolution of grasses (Poaceae), as well as the mammalian transition from insectivory to herbivory. The concordance of independent estimations suggests that AGF have been important in shaping the success of mammalian herbivory transition by improving the efficiency of energy acquisition from recalcitrant plant materials. Comparative genomics identified multiple lineage-specific genes in the AGF, two of which were acquired from rumen gut bacteria and animal hosts via horizontal gene transfer (HGT). A third AGF domain, plant-like polysaccharide lyase, represents a novel gene in fungi that potentially aids AGF to degrade pectin. Analysis of genomic and transcriptomic sequences confirmed both the presence and expression of these lineage-specific genes in nearly all AGF clades. These genetic elements may contribute to the exceptional abilities of AGF to degrade plant biomass and enable metabolism of the rumen microbes and animal hosts.IMPORTANCE Anaerobic fungi living in the rumen of herbivorous mammals possess an extraordinary ability to degrade plant biomass. We examined the origin and genomic composition of these poorly characterized anaerobic gut fungi using both transcriptome and genomic data. Phylogenomics and molecular dating analyses found remarkable concurrence of the divergence times of the rumen fungi, the forage grasses, and the dietary shift of ancestral mammals from primarily insectivory to herbivory. Comparative genomics identified unique machinery in these fungi to utilize plant polysaccharides. The rumen fungi were also identified with the ability to code for three protein domains with putative functions in plant pectin degradation and microbial defense, which were absent from all other fungal organisms (examined over 1,000 fungal genomes). Two of these domains were likely acquired from rumen gut bacteria and animal hosts separately via horizontal gene transfer. The third one is a plant-like polysaccharide lyase, representing a unique fungal enzyme with potential pectin breakdown abilities.

Transcriptome analysis and codominant markers development in caper, a drought tolerant orphan crop with medicinal value

Caper (Capparis spinosa L.) is a xerophytic shrub cultivated for its flower buds and fruits, used as food and for their medicinal properties. Breeding programs and even proper taxonomic classification of the genus Capparis has been hampered so far by the lack of reliable genetic information and molecular markers. Here, we present the first genomic resource for C. spinosa, generated by transcriptomic approach and de novo assembly. The sequencing effort produced nearly 80 million clean reads assembled into 124,723 unitranscripts. Careful annotation and comparison with public databases revealed homologs to genes with a key role in important metabolic pathways linked to abiotic stress tolerance and bio-compounds production, such purine, thiamine and phenylpropanoid biosynthesis, α-linolenic acid and lipid metabolism. Additionally, a panel of genes involved in stomatal development/distribution and encoding for Stress Associated Proteins (SAPs) was also identified. We also used the transcriptomic data to uncover novel molecular markers for caper. Out of 50 SSRs tested, 14 proved polymorphic and represent the first set of SSR markers for the genus Capparis. This transcriptome will be an important contribution to future studies and breeding programs for this orphan crop, aiding to the development of improved varieties to sustain agriculture in arid conditions.

Transcriptome analysis reveals the effects of grafting on sweetpotato scions during the full blooming stages

BACKGROUND

Sweetpotato (Ipomoea batatas) is a hexaploid plant and generally most genotypes do not flower at all in sub-tropics. Heterografting was carried out between sweetpotato cultivar 'Xushu 18' and Japanese morning glory (Ipomoea nil). With sweetpotato as 'scion' and I. nil as 'rootstock', sweetpotato was induced flowering in the autumn. However, little is known about the molecular mechanisms underlying sweetpotato responses to grafting, especially during the full blooming stages.

OBJECTIVES

To investigate the poorly understood molecular responses underlying the grafting-induced phenotypic processes in sweetpotato at full anthesis.

METHODS

In this study, to explore the transcriptome diversity and complexity of sweetpotato, PacBio Iso-Seq and Illumina RNA-seq analysis were combined to obtain full-length transcripts and to profile the changes in gene expression of five tissues: scion flowers (SF), scion leaves (SL), scion stems (SS), own-rooted leaves (OL) and own-rooted stems (OS).

RESULTS

A total of 138,151 transcripts were generated with an average length of 2255 bp, and more than 72% (100,396) of the transcripts were full-length. During full blooming, to examine the difference in gene expression of sweetpotato under grafting and natural growth conditions, 7905, 7795 and 15,707 differentially expressed genes were detected in pairwise comparisons of OS versus SS, OL versus SL and SL versus SF, respectively. Moreover, differential transcription of genes associated with anthocyanin biosynthesis, light pathway and photosynthesis, ethylene signal transduction pathway was observed in scion responses to grafting.

CONCLUSION

Our study is useful in understanding the molecular basis of grafting-induced flowering in grafted sweetpotatoes, and will lay a foundation for further research on sweetpotato breeding in the future.

Comparative transcriptome analysis of roots, stems and leaves of Isodon amethystoides reveals candidate genes involved in Wangzaozins biosynthesis

BACKGROUND

Isodon amethystoides (Ben-th) Cy Wu et Hsuan is an important traditional medicinal plant endowed with pharmacological properties effective in the treatment of various diseases, including pulmonary tuberculosis. The tetracyclic diterpenoids, Wangzaozins (Wangzaozin A, glaucocalyxin A, glaucocalyxin B), are the major bioactive compounds of I. amethystoides. However, the molecular information about the biosynthesis of these compounds still remains unclear.

RESULTS

An examination of the accumulated levels of Wangzaozins in I. amethystoides revealed considerable variations in the root, stem, and leaf tissues of this plant, indicating possible differences in metabolite biosynthesis and accumulation among various tissues. To better elucidate the tetracyclic diterpenoid biosynthesis pathway, we generated transcriptome sequences from the root, stem, and leaf tissues, and performed de novo sequence assembly, yielding 230,974 transcripts and 114,488 unigenes, with average N50 lengths of 1914 and 1241 bp, respectively. Putative functions could be assigned to 73,693 transcripts (31.9%) based on BLAST searches against annotation databases, including GO, KEGG, Swiss-Prot, NR, and Pfam. Moreover, the candidate genes involving in the diterpenoid biosynthesis, such as CPS, KSL, were also analyzed. The expression profiles of eight transcripts, involving the tetracyclic diterpenoid biosynthesis, were validated in different I. amethystoides tissues by qRT-PCR, unraveling the gene expression profile of the pathway. The differential expressions of ISPD, ISPF and ISPH (MEP pathway), and IaCPS and IaKSL (diterpenoid pathway) candidate genes in leaves and roots, may contribute to the high accumulation of Wangzaozins in I. amethystoides leaves.

CONCLUSION

The genomic dataset and analyses reported here lay the foundations for further research on this important medicinal plant.

Genetic Diversity and Population Structure of the USDA Sweetpotato (Ipomoea batatas) Germplasm Collections Using GBSpoly

Sweetpotato (Ipomoea batatas) plays a critical role in food security and is the most important root crop worldwide following potatoes and cassava. In the United States (US), it is valued at over $700 million USD. There are two sweetpotato germplasm collections (Plant Genetic Resources Conservation Unit and US Vegetable Laboratory) maintained by the USDA, ARS for sweetpotato crop improvement. To date, no genome-wide assessment of genetic diversity within these collections has been reported in the published literature. In our study, population structure and genetic diversity of 417 USDA sweetpotato accessions originating from 8 broad geographical regions (Africa, Australia, Caribbean, Central America, Far East, North America, Pacific Islands, and South America) were determined using single nucleotide polymorphisms (SNPs) identified with a genotyping-by-sequencing (GBS) protocol, GBSpoly, optimized for highly heterozygous and polyploid species. Population structure using Bayesian clustering analyses (STRUCTURE) with 32,784 segregating SNPs grouped the accessions into four genetic groups and indicated a high degree of mixed ancestry. A neighbor-joining cladogram and principal components analysis based on a pairwise genetic distance matrix of the accessions supported the population structure analysis. Pairwise FST values between broad geographical regions based on the origin of accessions ranged from 0.017 (Far East - Pacific Islands) to 0.110 (Australia - South America) and supported the clustering of accessions based on genetic distance. The markers developed for use with this collection of accessions provide an important genomic resource for the sweetpotato community, and contribute to our understanding of the genetic diversity present within the US sweetpotato collection and the species.

Genetic Diversity and Population Structure of the USDA Sweetpotato (Ipomoea batatas) Germplasm Collections Using GBSpoly

Sweetpotato (Ipomoea batatas) plays a critical role in food security and is the most important root crop worldwide following potatoes and cassava. In the United States (US), it is valued at over $700 million USD. There are two sweetpotato germplasm collections (Plant Genetic Resources Conservation Unit and US Vegetable Laboratory) maintained by the USDA, ARS for sweetpotato crop improvement. To date, no genome-wide assessment of genetic diversity within these collections has been reported in the published literature. In our study, population structure and genetic diversity of 417 USDA sweetpotato accessions originating from 8 broad geographical regions (Africa, Australia, Caribbean, Central America, Far East, North America, Pacific Islands, and South America) were determined using single nucleotide polymorphisms (SNPs) identified with a genotyping-by-sequencing (GBS) protocol, GBSpoly, optimized for highly heterozygous and polyploid species. Population structure using Bayesian clustering analyses (STRUCTURE) with 32,784 segregating SNPs grouped the accessions into four genetic groups and indicated a high degree of mixed ancestry. A neighbor-joining cladogram and principal components analysis based on a pairwise genetic distance matrix of the accessions supported the population structure analysis. Pairwise F_ST values between broad geographical regions based on the origin of accessions ranged from 0.017 (Far East - Pacific Islands) to 0.110 (Australia - South America) and supported the clustering of accessions based on genetic distance. The markers developed for use with this collection of accessions provide an important genomic resource for the sweetpotato community, and contribute to our understanding of the genetic diversity present within the US sweetpotato collection and the species.

Starch content differences between two sweet potato accessions are associated with specific changes in gene expression

Sweet potato (Ipomoea batatas (L.) Lam.) is one of the most important root crops in the world. Initial formation and development of storage roots (SRs) are key factors affecting its yields. In order to study the molecular mechanism and regulatory networks of the SRs development process, we have analyzed root transcriptomes between the high and low starch content sweet potato accessions at three different developmental stages. In this study, we assembled 46,840 unigenes using Illumina paired-end sequencing reads and identified differentially expressed genes (DEGs) between two accessions. The numbers of DEGs were increased with the development of SRs, indicating that the difference between two accessions is enlarging with the maturation. DEGs were mainly enriched in starch biosynthesis, plant hormones regulatory, and genetic information processing pathways. Then, expression patterns of DEGs that are most significant and starch biosynthesis related were validated using qRT-PCR. Our results provide valuable resources to future study on molecular mechanisms of SRs development and candidate genes for starch content improvement in sweet potato.