Transcriptome analysis reveals novel genes potentially involved in photoperiodic tuberization in potato

Transcriptome Analysis Reveals Novel Genes Potentially Involved in Tuberization in Potato

The formation and development of tubers, the primary edible and economic organ of potatoes, directly affect their yield and quality. The regulatory network and mechanism of tuberization have been preliminarily revealed in recent years, but plenty of relevant genes remain to be discovered. A few candidate genes were provided due to the simplicity of sampling and result analysis of previous transcriptomes related to tuberization. We sequenced and thoroughly analyzed the transcriptomes of thirteen tissues from potato plants at the tuber proliferation phase to provide more reference information and gene resources. Among them, eight tissues were stolons and tubers at different developmental stages, which we focused on. Five critical periods of tuberization were selected to perform an analysis of differentially expressed genes (DEGs), according to the results of the tissue correlation. Compared with the unswollen stolons (Sto), 2751, 4897, 6635, and 9700 DEGs were detected in the slightly swollen stolons (Sto1), swollen stolons (Sto2), tubers of proliferation stage 1 (Tu1), and tubers of proliferation stage 4 (Tu4). A total of 854 transcription factors and 164 hormone pathway genes were identified in the DEGs. Furthermore, three co-expression networks associated with Sto-Sto1, Sto2-Tu1, and tubers of proliferation stages two to five (Tu2-Tu5) were built using the weighted gene co-expression network analysis (WGCNA). Thirty hub genes (HGs) and 30 hub transcription factors (HTFs) were screened and focalized in these networks. We found that five HGs were reported to regulate tuberization, and most of the remaining HGs and HTFs co-expressed with them. The orthologs of these HGs and HTFs were reported to regulate processes (e.g., flowering, cell division, hormone synthesis, metabolism and signal transduction, sucrose transport, and starch synthesis) that were also required for tuberization. Such results further support their potential to control tuberization. Our study provides insights and countless candidate genes of the regulatory network of tuberization, laying the foundation for further elucidating the genetic basis of tuber development.

Exploring the Potential Role of Ribosomal Proteins to Enhance Potato Resilience in the Face of Changing Climatic Conditions

Potatoes have emerged as a key non-grain crop for food security worldwide. However, the looming threat of climate change poses significant risks to this vital food source, particularly through the projected reduction in crop yields under warmer temperatures. To mitigate potential crises, the development of potato varieties through genome editing holds great promise. In this study, we performed a comprehensive transcriptomic analysis to investigate microtuber development and identified several differentially expressed genes, with a particular focus on ribosomal proteins-RPL11, RPL29, RPL40 and RPL17. Our results reveal, by protein-protein interaction (PPI) network analyses, performed with the highest confidence in the STRING database platform (v11.5), the critical involvement of these ribosomal proteins in microtuber development, and highlighted their interaction with PEBP family members as potential microtuber activators. The elucidation of the molecular biological mechanisms governing ribosomal proteins will help improve the resilience of potato crops in the face of today's changing climatic conditions.

Changes in Structural and Thermodynamic Properties of Starch during Potato Tuber Dormancy

The main reserve polysaccharide of plants-starch-is undoubtedly important for humans. One of the main sources of starch is the potato tuber, which is able to preserve starch for a long time during the so-called dormancy period. However, accumulated data show that this dormancy is only relative, which raises the question of the possibility of some kind of starch restructuring during dormancy periods. Here, the effect of long-term periods of tuber rest (at 2-4 °C) on main parameters of starches of potato tubers grown in vivo or in vitro were studied. Along with non-transgenic potatoes, Arabidopsis phytochrome B (AtPHYB) transformants were investigated. Distinct changes in starch micro and macro structures-an increase in proportion of amorphous lamellae and of large-sized and irregular-shaped granules, as well as shifts in thickness of the crystalline lamellae-were detected. The degree of such alterations, more pronounced in AtPHYB-transgenic tubers, increased with the longevity of tuber dormancy. By contrast, the polymorphic crystalline structure (B-type) of starch remained unchanged regardless of dormancy duration. Collectively, our data support the hypothesis that potato starch remains metabolically and structurally labile during the entire tuber life including the dormancy period. The revealed starch remodeling may be considered a process of tuber preadaptation to the upcoming sprouting stage.

Transcriptomic analysis of Chinese yam (Dioscorea polystachya Turcz.) variants indicates brassinosteroid involvement in tuber development

Dioscorea is an important but underutilized genus of flowering plants that grows predominantly in tropical and subtropical regions. Several species, known as yam, develop large underground tubers and aerial bulbils that are used as food. The Chinese yam (D. polystachya Turcz.) is one of the few Dioscorea species that grows well in temperate regions and has been proposed as a climate-resilient crop to enhance food security in Europe. However, the fragile, club-like tubers are unsuitable for mechanical harvesting, which is facilitated by shorter and thicker storage organs. Brassinosteroids (BRs) play a key role in plant cell division, cell elongation and proliferation, as well as in the gravitropic response. We collected RNA-Seq data from the head, middle and tip of two tuber shape variants: F60 (long, thin) and F2000 (short, thick). Comparative transcriptome analysis of F60 vs. F2000 revealed 30,229 differentially expressed genes (DEGs), 1,393 of which were differentially expressed in the growing tip. Several DEGs are involved in steroid/BR biosynthesis or signaling, or may be regulated by BRs. The quantification of endogenous BRs revealed higher levels of castasterone (CS), 28-norCS, 28-homoCS and brassinolide in F2000 compared to F60 tubers. The highest BR levels were detected in the growing tip, and CS was the most abundant (439.6 ± 196.41 pmol/g in F2000 and 365.6 ± 112.78 pmol/g in F60). Exogenous 24-epi-brassinolide (epi-BL) treatment (20 nM) in an aeroponic system significantly increased the width-to-length ratio (0.045 ± 0.002) compared to the mock-treated plants (0.03 ± 0.002) after 7 weeks, indicating that exogenous epi-BL produces shorter and thicker tubers. In this study we demonstrate the role of BRs in D. polystachya tuber shape, providing insight into the role of plant hormones in yam storage organ development. We found that BRs can influence tuber shape in Chinese yam by regulating the expression of genes involved cell expansion. Our data can help to improve the efficiency of Chinese yam cultivation, which could provide an alternative food source and thus contribute to future food security in Europe.

Biochemical defense mechanism associated with host-specific disease resistance pathways against Rhizoctonia solani AG3-PT potatoes canker disease

Screening for resistance in 40 potato genotypes to Rhizoctonia solani AG-3PT-stem-canker, antioxidant enzymes activity as well as total phenol compounds were documented. Rhizoctonia solani AG-3PT-stem-canker is one of the most devastating diseases that leads to severe economic losses in potatoes, Solanum tuberosum globally. Crop management and eugenic practices, especially the use of resistance can be effective in reducing the disease incidence. However, the information about potato-R. Solani interaction is still limited. This study explored screening for resistance in forty potato genotypes to R. solani, analyzing biomass growth parameters (BGPs), as well as antioxidant enzymes activity of which peroxidase/peroxide-reductases (POXs), superoxide dismutase (SOD), polyphenol oxidase (PPO), catalase (CAT), phenylalanine ammonia-lyase (PAL), β-1,3-glucanase (GLU) and total phenol compounds (TPCs) were taken into account. In addition, we analyzed up-regulation of two gene markers (PR-1 and Osmotin), using reverse transcription quantitative PCR (RT-qPCR). For which, the resistant 'Savalan', partially resistant 'Agria', partially susceptible 'Sagita' and susceptible 'Pashandi' were selected to explore the trails in their roots and leaves over the time courses of 1, 2 and 3-weeks post inoculation (wpi) following inoculation. Cluster analysis divided potatoes into four distinct groups, based on disease severity scales (0-100%) significance. The BGPs, shoot and root length, fresh and dry weight, and root volume were also significantly higher in infected potatoes compared to non-inoculated controls. Antioxidant enzymes activity also indicated the highest increased levels for POX (fourfold at 3wpi), CAT (1.5-fold at 3wpi), SOD (6.8-fold at 1wpi), and PAL (2.7-fold at 3wpi) in the resistant genotype, 'Savalan', whereas the highest activity was recorded in TPC (twofold at 1 wpi), PPO (threefold at 3wpi), and GLU (2.3-fold at 1wpi) in partially resistant genotypes. Although the defense-related enzymatic activities were sharply elevated in the resistant and partially resistant genotypes following inoculation, no significant correlations were between the activity trends of the related enzymes. The two related gene markers also showed comprehensive transcriptional responses up to 3.4-fold, predominantly in resistant genotypes. Surprisingly, the PR-1 gene marker, basically resistant to Wilting agent Verticillium dahlia was overexpressed in resistant 'Savalan' and 'Agria' against R. solani AG3-PT. Similar results were obtained on Osmotin gene marker resistant to late-blight P. infestans, and early-blight Alternaria solani that similarly modulates immunity against R. solani. Furthermore, there was a significant correlation between resistance, enzyme activity, and gene expression in the aforesaid cultivars. Studying the physiological metabolic pathways of antioxidant enzymes activity appears to be an important direction in research to elucidate resistance to R. solani in potatoes.

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.

Physiological and protein profiling analysis provides insight into the underlying molecular mechanism of potato tuber development regulated by jasmonic acid in vitro

BACKGROUND

Jasmonates (JAs) are one of important phytohormones regulating potato tuber development. It is a complex process and the underlying molecular mechanism regulating tuber development by JAs is still limited. This study attempted to illuminate it through the potential proteomic dynamics information about tuber development in vitro regulated by exogenous JA.

RESULTS

A combined analysis of physiological and iTRAQ (isobaric tags for relative and absolute quantification)-based proteomic approach was performed in tuber development in vitro under exogenous JA treatments (0, 0.5, 5 and 50 μΜ). Physiological results indicated that low JA concentration (especially 5 μM) promoted tuber development, whereas higher JA concentration (50 μM) showed inhibition effect. A total of 257 differentially expressed proteins (DEPs) were identified by iTRAQ, which provided a comprehensive overview on the functional protein profile changes of tuber development regulated by JA. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis indicated that low JA concentration (especially 5 μM) exhibited the promotion effects on tuber development in various cellular processes. Some cell wall polysaccharide synthesis and cytoskeleton formation-related proteins were up-regulated by JA to promote tuber cell expansion. Some primary carbon metabolism-related enzymes were up-regulated by JA to provide sufficient metabolism intermediates and energy for tuber development. And, a large number of protein biosynthesis, degradation and assembly-related were up-regulated by JA to promote tuber protein biosynthesis and maintain strict protein quality control during tuber development.

CONCLUSIONS

This study is the first to integrate physiological and proteomic data to provide useful information about the JA-signaling response mechanism of potato tuber development in vitro. The results revealed that the levels of a number of proteins involved in various cellular processes were regulated by JA during tuber development. The proposed hypothetical model would explain the interaction of these DEPs that associated with tuber development in vitro regulated by JA.

Identification of novel associations of candidate genes with resistance to Rhizoctonia solani AG-3PT in Solanum tuberosum stem canker

To develop an understanding mechanism to define responding of potatoes to R. solani, we analyzed the expression of ten novel candidate gene-markers using reverse-transcription-quantitative PCR (RT-qPCR) in resistant 'Savalan' and partially resistant 'Agria' in contrast to susceptible 'Sagita', and partially susceptible 'Pashandi'. In addition, oxidant-enzymatic-activity of catalase and superoxide-dismutase, as well as biomass-growth-parameters; shoot and root length, fresh and dry weight, and root volume were considered as complementary factors to the involving mechanism accordingly. Gene-markers up-regulated maximum up to 3.5-fold with the highest correlation, r = 0.939** following R. solani-inoculation, predominantly in resistant genotypes. Surprisingly, WRKY8-gene, basically resistant to late-blight-Phytophtora infestans was also up-regulated to 2.3-fold in resistant 'Savalan' followed by 'Agria'. Similar results with 3.1-fold were obtained on Osmotin-gene resistant to early-blight-Alternaria alternata. Enzymatic-activity of catalase with 1.6-fold and superoxide-dismutase, 6.8-fold also showed the highest level of activity in resistant genotypes, and had a high significant correlation, r = 773** and r = 0.881** with expression levels of related gene-markers respectively. Similarly, there were significant differences in biomass-growth-parameters, but with reductions in partially susceptible 'Sagita' and susceptible 'Pashandi'. Conclusively, S. tuberosum-R. solani interaction revealed that certain gene-markers can cover resistance to more than one disease simultaneously.

Transcriptome and Small RNA Sequencing Reveal the Mechanisms Regulating Harvest Index in Brassica napus

Harvest index (HI), the ratio of harvested seed weight to total aboveground biomass weight, is an economically critical value reflecting the convergence of complex agronomic traits. HI values in rapeseed (Brassica napus) remain much lower than in other major crops, and the underlying regulatory network is largely unknown. In this study, we performed mRNA and small RNA sequencing to reveal the mechanisms shaping HI in B. napus during the seed-filling stage. A total of 8,410 differentially expressed genes (DEGs) between high-HI and low-HI accessions in four tissues (silique pericarp, seed, leaves, and stem) were identified. Combining with co-expression network, 72 gene modules were identified, and a key gene BnaSTY46 was found to participate in retarded establishment of photosynthetic capacity to influence HI. Further research found that the genes involved in circadian rhythms and response to stimulus may play important roles in HI and that their transcript levels were modulated by differentially expressed microRNAs (DEMs), and we identified 903 microRNAs (miRNAs), including 46 known miRNAs and 857 novel miRNAs. Furthermore, transporter activity-related genes were critical to enhancing HI in good cultivation environments. Of 903 miRNAs, we found that the bna-miR396-Bna.A06SRp34a/Bna.A01EMB3119 pair may control the seed development and the accumulation of storage compounds, thus contributing to higher HI. Our findings uncovered the underlying complex regulatory network behind HI and offer potential approaches to rapeseed improvement.

DNA methylation affects photoperiodic tuberization in potato (Solanum tuberosum L.) by mediating the expression of genes related to the photoperiod and GA pathways

Overcoming short-day-dependent tuberization to adapt to long-day conditions is critical for the widespread geographical success of potato. The genetic pathways of photoperiodic tuberization are similar to those of photoperiodic flowering. DNA methylation plays an important role in photoperiodic flowering. However, little is known about how DNA methylation affects photoperiodic tuberization in potato. Here, we verified the effect of a DNA methylation inhibitor on photoperiodic tuberization and compared the DNA methylation levels and differentially methylated genes (DMGs) in the photoperiodic tuberization process between photoperiod-sensitive and photoperiod-insensitive genotypes, aiming to dissect the role of DNA methylation in the photoperiodic tuberization of potato. We found that a DNA methylation inhibitor could promote tuber initiation in strict short-day genotypes. Whole-genome DNA methylation sequencing showed that the photoperiod-sensitive and photoperiod-insensitive genotypes had distinct DNA methylation modes in which few differentially methylated genes were shared. Transcriptome analysis confirmed that the DNA methylation inhibitor regulated the expression of the key genes involved in the photoperiod and GA pathways to promote tuber initiation in the photoperiod-sensitive genotype. Comparison of the DNA methylation levels and transcriptome levels identified 52 candidate genes regulated by DNA methylation that were predicted to be involved in photoperiodic tuberization. Our findings provide a new perspective for understanding the relationship between photoperiod-dependent and GA-regulated tuberization. Uncovering the epigenomic signatures of these pathways will greatly enhance potato breeding for adaptation to a wide range of environments.

Genome-wide identification and expression profiling of basic leucine zipper transcription factors following abiotic stresses in potato (Solanum tuberosum L.)

Potato (Solanum tuberosum L.) is an important food crop that is grown and consumed worldwide. The growth and productivity of this crop are severely affected by various abiotic stresses. Basic leucine zipper (bZIP) transcription factors (TFs) in plants are well known for their function during growth and development. However, systematic and in-depth identification and functional characterization of the bZIP gene family of potato is lacking. In the current study, we identified a total of 90 bZIPs (StbZIP) distributed on 12 linkage groups of potato. Based on the previous functional annotation and classification of bZIPs in Arabidopsis, wheat, and rice, a phylogenetic tree of potato bZIPs was constructed and genes were categorized into various functional groups (A to I, S, and U) as previously annotated in Arabidopsis thaliana. Analyses of the transcript sequence (RNA-seq) data led to identifying a total of 18 candidate StbZIPs [four in roots, eight in the tuber, six in mesocarp and endocarp] that were expressed in a tissue-specific manner. Differential expression analysis under the various abiotic conditions (salt, mannitol, water, and heat stress) and treatment with phytohormones (ABA, GA, IAA, and BAP) led to the identification of forty-two [thirteen under salt stress, two under mannitol stress, ten under water stress, and eighteen under heat stress], and eleven [eight and three StbZIPs upon treatment with ABA, and IAA, respectively] candidate StbZIPs, respectively. Using sequence information of candidate StbZIPs, a total of 22 SSR markers were also identified in this study. In conclusion, the genome-wide identification analysis coupled with RNA-Seq expression data led to identifying candidate StbZIPs, which are dysregulated, and may play a pivotal role under various abiotic stress conditions. This study will pave the way for future functional studies using forward and reverse genetics to improve abiotic stress tolerance in potato.

Genome-wide identification and expression profiling of basic leucine zipper transcription factors following abiotic stresses in potato (Solanum tuberosum L.)

Potato (Solanum tuberosum L.) is an important food crop that is grown and consumed worldwide. The growth and productivity of this crop are severely affected by various abiotic stresses. Basic leucine zipper (bZIP) transcription factors (TFs) in plants are well known for their function during growth and development. However, systematic and in-depth identification and functional characterization of the bZIP gene family of potato is lacking. In the current study, we identified a total of 90 bZIPs (StbZIP) distributed on 12 linkage groups of potato. Based on the previous functional annotation and classification of bZIPs in Arabidopsis, wheat, and rice, a phylogenetic tree of potato bZIPs was constructed and genes were categorized into various functional groups (A to I, S, and U) as previously annotated in Arabidopsis thaliana. Analyses of the transcript sequence (RNA-seq) data led to identifying a total of 18 candidate StbZIPs [four in roots, eight in the tuber, six in mesocarp and endocarp] that were expressed in a tissue-specific manner. Differential expression analysis under the various abiotic conditions (salt, mannitol, water, and heat stress) and treatment with phytohormones (ABA, GA, IAA, and BAP) led to the identification of forty-two [thirteen under salt stress, two under mannitol stress, ten under water stress, and eighteen under heat stress], and eleven [eight and three StbZIPs upon treatment with ABA, and IAA, respectively] candidate StbZIPs, respectively. Using sequence information of candidate StbZIPs, a total of 22 SSR markers were also identified in this study. In conclusion, the genome-wide identification analysis coupled with RNA-Seq expression data led to identifying candidate StbZIPs, which are dysregulated, and may play a pivotal role under various abiotic stress conditions. This study will pave the way for future functional studies using forward and reverse genetics to improve abiotic stress tolerance in potato.

Transcriptome Analysis Reveals Key Pathways and Hormone Activities Involved in Early Microtuber Formation of Dioscorea opposita

Chinese yam (Dioscorea opposita) is an important tuberous crop used for both food and medicine. Despite a long history of cultivation, the understanding of D. opposita genetics and molecular biology remains scant, which has limited its genetic improvement. This work presents a de novo transcriptome sequencing analysis of microtuber formation in D. opposita. We assembled cDNA libraries from different stages during the process of microtuber formation, designated as initial explants (EXP), axillary bud proliferation after three weeks (BUD), and microtuber visible after four weeks (MTV). More differentially expressed genes (DEGs) and pathways were identified between BUD vs. EXP than in MTV vs. BUD, indicating that proliferation of the axillary bud is the key stage of microtuber induction. Gene classification and pathway enrichment analysis showed that microtuber formation is tightly coordinated with primary metabolism, such as amino acid biosynthesis, ribosomal component biosynthesis, and starch and sucrose metabolism. The formation of the microtuber is regulated by a variety of plant hormones, including ABA. Combined with analysis of physiological data, we suggest that ABA positively regulates tuberization in D. opposita. This study will serve as an empirical foundation for future molecular studies and for the propagation of D. opposita germplasm in field crops.

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.

The origins and adaptation of European potatoes reconstructed from historical genomes

Potato, one of the most important staple crops, originates from the highlands of the equatorial Andes. There, potatoes propagate vegetatively via tubers under short days, constant throughout the year. After their introduction to Europe in the sixteenth century, potatoes adapted to a shorter growing season and to tuber formation under long days. Here, we traced the demographic and adaptive history of potato introduction to Europe. To this end, we sequenced 88 individuals that comprise landraces, modern cultivars and historical herbarium samples, including specimens collected by Darwin during the voyage of the Beagle. Our findings show that European potatoes collected during the period 1650-1750 were closely related to Andean landraces. After their introduction to Europe, potatoes admixed with Chilean genotypes. We identified candidate genes putatively involved in long-day pre-adaptation, and showed that the 1650-1750 European individuals were not long-day adapted through previously described allelic variants of the CYCLING DOF FACTOR1 gene. Such allelic variants were detected in Europe during the nineteenth century. Our study highlights the power of combining contemporary and historical genomes to understand the complex evolutionary history of crop adaptation to new environments.

The RNA-Seq-based high resolution gene expression atlas of chickpea (Cicer arietinum L.) reveals dynamic spatio-temporal changes associated with growth and development

Chickpea is one of the world's largest cultivated food legumes and is an excellent source of high-quality protein to the human diet. Plant growth and development are controlled by programmed expression of a suite of genes at the given time, stage, and tissue. Understanding how the underlying genome sequence translates into specific plant phenotypes at key developmental stages, information on gene expression patterns is crucial. Here, we present a comprehensive Cicer arietinum Gene Expression Atlas (CaGEA) across different plant developmental stages and organs covering the entire life cycle of chickpea. One of the widely used drought tolerant cultivars, ICC 4958 has been used to generate RNA-Seq data from 27 samples at 5 major developmental stages of the plant. A total of 816 million raw reads were generated and of these, 794 million filtered reads after quality control (QC) were subjected to downstream analysis. A total of 15,947 unique number of differentially expressed genes across different pairwise tissue combinations were identified. Significant differences in gene expression patterns contributing in the process of flowering, nodulation, and seed and root development were inferred in this study. Furthermore, differentially expressed candidate genes from "QTL-hotspot" region associated with drought stress response in chickpea were validated.

Solanum tuberosum ZPR1 encodes a light-regulated nuclear DNA-binding protein adjusting the circadian expression of StBBX24 to light cycle

ZPR1 proteins belong to the C4-type of zinc finger coordinators known in animal cells to interact with other proteins and participate in cell growth and proliferation. In contrast, the current knowledge regarding plant ZPR1 proteins is very scarce. Here, we identify a novel potato nuclear factor belonging to this family and named StZPR1. StZPR1 is specifically expressed in photosynthetic organs during the light period, and the ZPR1 protein is located in the nuclear chromatin fraction. From modelling and experimental analyses, we reveal the StZPR1 ability to bind the circadian DNA cis motif 'CAACAGCATC', named CIRC and present in the promoter of the clock-controlled double B-box StBBX24 gene, the expression of which peaks in the middle of the day. We found that transgenic lines silenced for StZPR1 expression still display a 24 h period for the oscillation of StBBX24 expression but delayed by 4 h towards the night. Importantly, other BBX genes exhibit altered circadian regulation in these lines. Our data demonstrate that StZPR1 allows fitting of the StBBX24 circadian rhythm to the light period and provide evidence that ZPR1 is a novel clock-associated protein in plants necessary for the accurate rhythmic expression of specific circadian-regulated genes.

Transcriptome Analysis of Storage Roots and Fibrous Roots of the Traditional Medicinal Herb Callerya speciosa (Champ.) ScHot

Callerya speciosa (Champ.) ScHot is a woody perennial plant in Fabaceae, the roots of which are used medicinally. The storage roots of C. speciosa are derived from fibrous roots, but not all fibrous roots can develop into storage roots. To detect key genes involved in storage roots formation, we performed Illumina sequencing of the C. speciosa storage roots and fibrous roots. De novo assembly resulted in 161,926 unigenes, which were subsequently annotated by BLAST, GO and KEGG analyses. After expression profiling, 4538 differentially expressed genes were identified. The KEGG pathway enrichment analysis revealed changes in the biosynthesis of cytokinin, phenylpropanoid, starch, sucrose, flavone and other secondary metabolites. Transcription factor-related differentially expressed genes (DEGs) were also identified, including such gene families as GRAS, COL, MIKC, ERF, LBD, and NAC. The DEGs related to light signaling, starch, sugar, photohormones and cell wall-loosening might be involved in the formation of storage roots. This study provides the first transcriptome profiling of C. speciosa roots, data that will facilitate future research of root development and metabolites with medicinal value as well as the breeding of C. speciosa.

Mapping of quantitative trait loci for tuber starch and leaf sucrose contents in diploid potato

Most QTL for leaf sucrose content map to positions that are similar to positions of QTL for tuber starch content in diploid potato. In the present study, using a diploid potato mapping population and Diversity Array Technology (DArT) markers, we identified twelve quantitative trait loci (QTL) for tuber starch content on seven potato chromosomes: I, II, III, VIII, X, XI, and XII. The most important QTL spanned a wide region of chromosome I (42.0–104.6 cM) with peaks at 63 and 84 cM which explained 17.6 and 19.2% of the phenotypic variation, respectively. ADP-glucose pyrophosphorylase (AGPase) is the key enzyme for starch biosynthesis. The gene encoding the large subunit of this enzyme, AGPaseS-a, was localized to chromosome I at 102.3 cM and accounted for 15.2% of the variance in tuber starch content. A more than 100-fold higher expression of this gene was observed in RT-qPCR assay in plants with the marker allele AGPaseS-a1334. This study is the first to report QTL for sucrose content in potato leaves. QTL for sucrose content in leaves were located on eight potato chromosomes: I, II, III, V, VIII, IX, X and XII. In 5-week-old plants, only one QTL for leaf sucrose content was detected after 8 h of darkness; four QTL were detected after 8 h of illumination. In 11-week-old plants, 6 and 3 QTL were identified after dark and light phases, respectively. Of fourteen QTL for leaf sucrose content, eleven mapped to positions that were similar to QTL for tuber starch content. These results provide genetic information for further research examining the relationships between metabolic carbon molecule sources and sinks in potato plants.

Candidate Resistant Genes of Sand Pear (Pyrus pyrifolia Nakai) to Alternaria alternata Revealed by Transcriptome Sequencing

Pear black spot (PBS) disease, which is caused by Alternaria alternata (Aa), is one of the most serious diseases affecting sand pear (Pyrus pyrifolia Nakai) cultivation worldwide. To investigate the defense mechanisms of sand pear in response to Aa, the transcriptome of a sand pear germplasm with differential resistance to Aa was analyzed using Illumina paired-end sequencing. Four libraries derived from PBS-resistant and PBS-susceptible sand pear leaves were characterized through inoculation or mock-inoculation. In total, 20.5 Gbp of sequence data and 101,632,565 reads were generated, representing 44717 genes. Approximately 66% of the genes or sequenced reads could be aligned to the pear reference genome. A large number (5213) of differentially expressed genes related to PBS resistance were obtained; 34 microsatellites were detected in these genes, and 28 genes were found to be closely related to PBS resistance. Using a transcriptome analysis in response to PBS inoculation and comparison analysis to the PHI database, 4 genes (Pbr039001, Pbr001627, Pbr025080 and Pbr023112) were considered to be promising candidates for sand pear resistance to PBS. This study provides insight into changes in the transcriptome of sand pear in response to PBS infection, and the findings have improved our understanding of the resistance mechanism of sand pear to PBS and will facilitate future gene discovery and functional genome studies of sand pear.

Microarray analysis of gene expression patterns in the leaf during potato tuberization in the potato somatic hybrid Solanum tuberosum and Solanum etuberosum

Genes involved in photoassimilate partitioning and changes in hormonal balance are important for potato tuberization. In the present study, we investigated gene expression patterns in the tuber-bearing potato somatic hybrid (E1-3) and control non-tuberous wild species Solanum etuberosum (Etb) by microarray. Plants were grown under controlled conditions and leaves were collected at eight tuber developmental stages for microarray analysis. A t-test analysis identified a total of 468 genes (94 up-regulated and 374 down-regulated) that were statistically significant (p ≤ 0.05) and differentially expressed in E1-3 and Etb. Gene Ontology (GO) characterization of the 468 genes revealed that 145 were annotated and 323 were of unknown function. Further, these 145 genes were grouped based on GO biological processes followed by molecular function and (or) PGSC description into 15 gene sets, namely (1) transport, (2) metabolic process, (3) biological process, (4) photosynthesis, (5) oxidation-reduction, (6) transcription, (7) translation, (8) binding, (9) protein phosphorylation, (10) protein folding, (11) ubiquitin-dependent protein catabolic process, (12) RNA processing, (13) negative regulation of protein, (14) methylation, and (15) mitosis. RT-PCR analysis of 10 selected highly significant genes (p ≤ 0.01) confirmed the microarray results. Overall, we show that candidate genes induced in leaves of E1-3 were implicated in tuberization processes such as transport, carbohydrate metabolism, phytohormones, and transcription/translation/binding functions. Hence, our results provide an insight into the candidate genes induced in leaf tissues during tuberization in E1-3.

Thermodynamic and structural properties of tuber starches from transgenic potato plants grown in vitro and in vivo

Potato plants harboring Phytochrome B (PHYB) gene from Arabidopsis thaliana or rol genes from Agrobacterium rhizogenes were used to study the effect of transgene expression on structure and properties of starch in tubers. Thermodynamic characteristics of starch (melting temperature, enthalpy of melting, thickness of crystalline lamellae) were shown to be variable depending on the transgene expression and plant culturing mode: in vitro or in soil. The expression of rolB or rolC genes in in vitro cultured plants evoked opposite effects on starch melting temperature and crystalline lamellae thickness. AtPHYB or rolB expression in the soil-grown potato led to the formation of more defective or more ordered starch structures, respectively, in comparison with starches of the same lines grown in vitro. On the whole, our study revealed genotype-dependent differences between starches extracted from tubers of in vitro or in vivo grown plants.