Comparative transcriptomics reveals shared gene expression changes during independent evolutionary origins of stem and hypocotyl/root tubers in Brassica (Brassicaceae)

RNA-Seq Reveals That Multiple Pathways Are Involved in Tuber Expansion in Tiger Nuts (Cyperus esculentus L.)

The tiger nut (Cyperus esculentus L.) is a usable tuber and edible oil plant. The size of the tubers is a key trait that determines the yield and the mechanical harvesting of tiger nut tubers. However, little is known about the anatomical and molecular mechanisms of tuber expansion in tiger nut plants. This study conducted anatomical and comprehensive transcriptomics analyses of tiger nut tubers at the following days after sowing: 40 d (S1); 50 d (S2); 60 d (S3); 70 d (S4); 90 d (S5); and 110 d (S6). The results showed that, at the initiation stage of a tiger nut tuber (S1), the primary thickening meristem (PTM) surrounded the periphery of the stele and was initially responsible for the proliferation of parenchyma cells of the cortex (before S1) and then the stele (S2-S3). The increase in cell size of the parenchyma cells occurred mainly from S1 to S3 in the cortex and from S3 to S4 in the stele. A total of 12,472 differentially expressed genes (DEGs) were expressed to a greater extent in the S1-S3 phase than in S4-S6 phase. DEGs related to tuber expansion were involved in cell wall modification, vesicle transport, cell membrane components, cell division, the regulation of plant hormone levels, signal transduction, and metabolism. DEGs involved in the biosynthesis and the signaling of indole-3-acetic acid (IAA) and jasmonic acid (JA) were expressed highly in S1-S3. The endogenous changes in IAA and JAs during tuber development showed that the highest concentrations were found at S1 and S1-S3, respectively. In addition, several DEGs were related to brassinosteroid (BR) signaling and the G-protein, MAPK, and ubiquitin-proteasome pathways, suggesting that these signaling pathways have roles in the tuber expansion of tiger nut. Finally, we come to the conclusion that the cortex development preceding stele development in tiger nut tubers. The auxin signaling pathway promotes the division of cortical cells, while the jasmonic acid pathway, brassinosteroid signaling, G-protein pathway, MAPK pathway, and ubiquitin protein pathway regulate cell division and the expansion of the tuber cortex and stele. This finding will facilitate searches for genes that influence tuber expansion and the regulatory networks in developing tubers.

Functional and ecological diversification of underground organs in Solanum

The evolution of geophytes in response to different environmental stressors is poorly understood largely due to the great morphological variation in underground plant organs, which includes species with rhizomatous structures or underground storage organs (USOs). Here we compare the evolution and ecological niche patterns of different geophytic organs in Solanum L., classified based on a functional definition and using a clade-based approach with an expert-verified specimen occurrence dataset. Results from PERMANOVA and Phylogenetic ANOVAs indicate that geophytic species occupy drier areas, with rhizomatous species found in the hottest areas whereas species with USOs are restricted to cooler areas in the montane tropics. In addition, rhizomatous species appear to be adapted to fire-driven disturbance, in contrast to species with USOs that appear to be adapted to prolonged climatic disturbance such as unfavorable growing conditions due to drought and cold. We also show that the evolution of rhizome-like structures leads to changes in the relationship between range size and niche breadth. Ancestral state reconstruction shows that in Solanum rhizomatous species are evolutionarily more labile compared to species with USOs. Our results suggest that underground organs enable plants to shift their niches towards distinct extreme environmental conditions and have different evolutionary constraints.

Comparative transcriptome analysis during tuberous stem formation in Kohlrabi (B. oleracea var. gongylodes) at early growth periods (seedling stages)

Tuberous stem of kohlrabi is an important agronomic trait, however, the molecular basis of tuberization is poorly understood. To elucidate the tuberization mechanism, we conducted a comparative transcriptomic analysis between kohlrabi and broccoli at 10 and 20 days after germination (DAG) as tuberous stem initiated between these time points. A total of 5580 and 2866 differentially expressed transcripts (DETs) were identified between genotypes (kohlrabi vs. broccoli) and growth stages (10 DAG vs. 20 DAG), respectively, and most of the DETs were down-regulated in kohlrabi. Gene ontology (GO) and KEGG pathway enrichment analyses showed that the DETs between genotypes are involved in cell wall loosening and expansion, cell cycle and division, carbohydrate metabolism, hormone transport, hormone signal transduction and in several transcription factors. The DETs identified in those categories may directly/indirectly relate to the initiation and development of tuberous stem in kohlrabi. In addition, the expression pattern of the hormone synthesis related DETs coincided with the endogenous IAA, IAAsp, GA, ABA, and tZ profiles in kohlrabi and broccoli seedlings, that were revealed in our phytohormone analysis. This is the first report on comparative transcriptome analysis for tuberous stem formation in kohlrabi at early growth periods. The resulting data could provide significant insights into the molecular mechanism underlying tuberous stem development in kohlrabi as well as in other tuberous organ forming crops.

Transcriptomic Analysis of Radish (Raphanus sativus L.) Spontaneous Tumor

Spontaneous tumors can develop in different organs of various plant species without any pathogen infection and, as a rule, appear in plants with a certain genotype: Mutants, interspecific hybrids, etc. In particular, among the inbred lines of radish (Raphanus sativus L.), lines that form spontaneous tumors on the taproot during the flowering period were obtained many years ago. In this work, we analyzed the differential gene expression in the spontaneous tumors of radish versus the lateral roots using the RNA-seq method. Data were obtained indicating the increased expression of genes associated with cell division and growth (especially genes that regulate G2-M transition and cytokinesis) in the spontaneous tumor. Among genes downregulated in the tumor tissue, genes participating in the response to stress and wounding, mainly involved in the biosynthesis of jasmonic acid and glucosinolates, were enriched. Our data will help elucidate the mechanisms of spontaneous tumor development in higher plants.

Comparative transcriptomics of a monocotyledonous geophyte reveals shared molecular mechanisms of underground storage organ formation

Many species from across the vascular plant tree-of-life have modified standard plant tissues into tubers, bulbs, corms, and other underground storage organs (USOs), unique innovations which allow these plants to retreat underground. Our ability to understand the developmental and evolutionary forces that shape these morphologies is limited by a lack of studies on certain USOs and plant clades. We take a comparative transcriptomics approach to characterizing the molecular mechanisms of tuberous root formation in Bomarea multiflora (Alstroemeriaceae) and compare these mechanisms to those identified in other USOs across diverse plant lineages; B. multiflora fills a key gap in our understanding of USO molecular development as the first monocot with tuberous roots to be the focus of this kind of research. We sequenced transcriptomes from the growing tip of four tissue types (aerial shoot, rhizome, fibrous root, and root tuber) of three individuals of B. multiflora. We identified differentially expressed isoforms between tuberous and non-tuberous roots and tested the expression of a priori candidate genes implicated in underground storage in other taxa. We identify 271 genes that are differentially expressed in root tubers versus non-tuberous roots, including genes implicated in cell wall modification, defense response, and starch biosynthesis. We also identify a phosphatidylethanolamine-binding protein, which has been implicated in tuberization signalling in other taxa and, through gene-tree analysis, place this copy in a phylogenetic context. These findings suggest that some similar molecular processes underlie the formation of USOs across flowering plants despite the long evolutionary distances among taxa and non-homologous morphologies (e.g., bulbs vs. tubers). (Plant development, tuberous roots, comparative transcriptomics, geophytes).

QTL Mapping and Candidate Gene Identification of Swollen Root Formation in Turnip

The swollen root is an important agronomic trait and is a determinant of yield for turnips, which are cultivated as both vegetables and fodder. However, the genetic mechanism of swollen root formation is poorly understood. In this study, we analyzed the F₂ and BC₁P₂ populations derived from a cross between “10601” (European turnip with swollen root, Brassica rapa ssp. rapifera, AA, 2n = 2× = 20) and “10603” (Chinese cabbage with normal root, Brassica rapa ssp. pekinensis, AA, 2n = 2× = 20), and suggested that the swollen root is a quantitative trait. Two major quantitative trait loci (QTLs), FR1.1 (Fleshy root 1.1) and FR7.1 (Fleshy root 7.1), were identified by QTL-seq analysis and further confirmed by QTL mapping in F₂ and BC₁P₂ populations. The QTL FR1.1 with a likelihood of odd (LOD) of 7.01 explained 17.2% of the total phenotypic variations for root diameter and the QTL FR7.1 explained 23.0% (LOD = 9.38) and 31.0% (LOD = 13.27) of the total phenotypic variations in root diameter and root weight, respectively. After a recombinant screening, the major QTL FR7.1 was further narrowed down to a 220 kb region containing 47 putative genes. A candidate gene, Bra003652, which is a homolog of AT1G78240 that plays an essential role in cell adhesion and disorganized tumor-like formation in Arabidopsis thaliana, was identified in this region. In addition, expression and parental allele analysis supported that Bra003652 was a possible candidate gene of QTL FR7.1 for swollen root formation in turnip. Our research may provide new insight into the molecular mechanism of swollen root formation in root crops.

Turing-like mechanism in a stochastic reaction-diffusion model recreates three dimensional vascular patterning of plant stems

Vascular tissue in plants provides a resource distribution network for water and nutrients that exhibits remarkable diversity in patterning among different species. In many succulent plants, the vascular network includes longitudinally-oriented supplemental vascular bundles (SVBs) in the central core of the succulent stems and roots in addition to the more typical zone of vascular tissue development (vascular cambium) in a cylinder at the periphery of the succulent organ. Plant SVBs evolved in over 38 plant families often in tandem with evolutionary increases in stem and root parenchyma storage tissue, so it is of interest to understand the evolutionary-developmental processes responsible for their recurrent evolution and patterning. Previous mathematical models have successfully recreated the two-dimensional vascular patterns in stem and root cross sections, but such models have yet to recreate three-dimensional vascular patterning. Here, a stochastic reaction-diffusion model of plant vascular bundle patterning is developed in an effort to highlight a potential mechanism of three dimensional patterning-Turing pattern formation coupled with longitudinal efflux of a regulatory molecule. A relatively simple model of four or five molecules recreated empirical SVB patterns and many other common vascular arrangements. SVBs failed to develop below a threshold width of parenchymatous tissues, suggesting a mechanism of evolutionary character loss due to changes in the spatial context in which development takes place. Altered diffusion rates of the modeled activator and substrate molecules affected the number and size of the simulated SVBs. This work provides a first mathematical model employing a stochastic Turing-type mechanism that recreates three dimensional vascular patterns seen in plant stems. The model offers predictions that can be tested using molecular-genetic approaches. Evolutionary-developmental ramifications concerning evolution of diffusion rates, organ size and geometry are discussed.

Expression profiles of organogenesis-related genes over the time course of one-step de novo shoot organogenesis from intact seedlings of kohlrabi

Kohlrabi (Brassica oleracea var. gongylodes) is an important vegetable crop that is able to undergo shoot regeneration in culture from intact seedlings in a single-step regeneration process, using cytokinin as the only plant growth regulator. In this work, we present the expression profiles of seven organogenesis-related genes over the time course of shoot regeneration from intact seedlings of kohlrabi cv. Vienna Purple on shoot regeneration media containing trans-zeatin, cis-zeatin, benzyl adenine or thidiazuron. Two auxin transporter genes - PIN3 and PIN4, a cytokinin response regulator - ARR5, two shoot apical meristem-related transcription factors - CUC1 and RGD3, and two cell cycle-related genes - CDKB2;1 and CYCB2;4 - displayed bimodal expression patterns on most cytokinin-containing media when their expression levels were normalized against control plants grown on hormone-free media. The first expression peak corresponded to direct upregulation by cytokinin from the growth media, and the second one reflected transcriptional events related to callus formation and/or acquisition of organogenic competence, corresponding to the shoot regeneration phases that have already been characterized in Arabidopsis thaliana. We demonstrate that the genes involved in the two-step shoot regeneration of Arabidopsis display their expected expression profiles during the single-step shoot regeneration of its close phylogenetic relative kohlrabi confirming the universality of their roles in the distinct phases of the regeneration process in Brassicaceae. The results presented here represent a first step towards genetic characterization of the morphogenetic processes in this important crop species.