Comparative Transcriptomics of Root Development in Wild and Cultivated Carrots

Deciphering Carotenoid and Flowering Pathway Gene Variations in Eastern and Western Carrots (Daucus carota L.)

BACKGROUND/OBJECTIVES

Carrot is a major root vegetable in the Apiaceae owing to its abundant carotenoids, antioxidants, vitamins, and minerals. The modern dark orange western carrot was derived from sequential domestication events from the white-rooted wild form to the pale orange-, purple-, or yellow-rooted eastern carrot. Genetic and molecular studies between eastern and western carrots are meager despite their evolutionary relatedness.

METHODS

Twelve RNA seq libraries obtained from distinct eastern and western cultivars at vegetative and reproductive developmental stages were utilized to identify differentially expressed genes (DEGs) to decode the key molecular genetic changes in carotenoid and flowering pathways.

RESULTS

In the carotenoid pathway, an upregulation of the PSY, CRTISO, and LCYE genes was observed in the western cultivar, while the eastern cultivar exhibited a higher abundance of downstream enzymes, particularly CCD and NCED1. These later enzymes are crucial in linking apocarotenoids and xanthin-mediated ABA signaling. In the flowering pathway, we noted a greater expression of DEGs associated with the photoperiod and vernalization pathways in the western cultivar. In contrast, the eastern cultivar displayed a dominance of genes from the autonomous pathway (FLD, LD, FLK, and PEBP) that function to repress FLC. The experimental validation of 12 key genes through quantitative real-time PCR further confirms their functional role in carrots.

CONCLUSIONS

The identified key regulatory genes in these major pathways are valuable for designing breeding strategies for manipulating carotenoid content and flowering time while developing climate-specific carrots. The knowledge of carotenoid and flowering pathways is advantageous in producing nutritionally improved roots and seeds in carrots across diverse climates.

Comparative Analysis of Hulless Barley Transcriptomes to Regulatory Effects of Phosphorous Deficiency

Hulless barley is a cold-resistant crop widely planted in the northwest plateau of China. It is also the main food crop in this region. Phosphorus (P), as one of the important essential nutrient elements, regulates plant growth and defense. This study aimed to analyze the development and related molecular mechanisms of hulless barley under P deficiency and explore the regulatory genes so as to provide a basis for subsequent molecular breeding research. Transcriptome analysis was performed on the root and leaf samples of hulless barley cultured with different concentrations of KH2PO4 (1 mM and 10 μM) Hoagland solution. A total of 46,439 genes were finally obtained by the combined analysis of leaf and root samples. Among them, 325 and 453 genes had more than twofold differences in expression. These differentially expressed genes (DEGs) mainly participated in the abiotic stress biosynthetic process through Gene Ontology prediction. Moreover, the Kyoto Encyclopedia of Genes and Genomes showed that DEGs were mainly involved in photosynthesis, plant hormone signal transduction, glycolysis, phenylpropanoid biosynthesis, and synthesis of metabolites. These pathways also appeared in other abiotic stresses. Plants initiated multiple hormone synergistic regulatory mechanisms to maintain growth under P-deficient conditions. Transcription factors (TFs) also proved these predictions. The enrichment of ARR-B TFs, which positively regulated the phosphorelay-mediated cytokinin signal transduction, and some other TFs (AP2, GRAS, and ARF) was related to plant hormone regulation. Some DEGs showed different values in their FPKM (fragment per kilobase of transcript per million mapped reads), but the expression trends of genes responding to stress and phosphorylation remained highly consistent. Therefore, in the case of P deficiency, the first response of plants was the expression of stress-related genes. The effects of this stress on plant metabolites need to be further studied to improve the relevant regulatory mechanisms so as to further understand the importance of P in the development and stress resistance of hulless barley.

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.

Comparative Root Transcriptome Profiling and Gene Regulatory Network Analysis between Eastern and Western Carrot (Daucus carota L.) Cultivars Reveals Candidate Genes for Vascular Tissue Patterning

Carrot (Daucus carota L.) is a highly consumed vegetable rich in carotenoids, known for their potent antioxidant, anti-inflammatory, and immune-protecting properties. While genetic and molecular studies have largely focused on wild and Western carrot cultivars (cvs), little is known about the evolutionary interactions between closely related Eastern and Western cvs. In this study, we conducted comparative transcriptome profiling of root tissues from Eastern (UHSBC-23-1) and Western (UHSBC-100) carrot cv. to better understand differentially expressed genes (DEGs) associated with storage root development and vascular cambium (VC) tissue patterning. Through reference-guided TopHat mapping, we achieved an average mapping rate of 73.87% and identified a total of 3544 DEGs (p < 0.05). Functional annotation and gene ontology classification revealed 97 functional categories, including 33 biological processes, 19 cellular components, 45 metabolic processes, and 26 KEGG pathways. Notably, Eastern cv. exhibited enrichment in cell wall, plant-pathogen interaction, and signal transduction terms, while Western cv. showed dominance in photosynthesis, metabolic process, and carbon metabolism terms. Moreover, constructed gene regulatory network (GRN) for both cvs. obtained orthologs with 1222 VC-responsive genes of Arabidopsis thaliana. In Western cv, GRN revealed VC-responsive gene clusters primarily associated with photosynthetic processes and carbon metabolism. In contrast, Eastern cv. exhibited a higher number of stress-responsive genes, and transcription factors (e.g., MYB15, WRKY46, AP2/ERF TF connected via signaling pathways with NAC036) were identified as master regulators of xylem vessel differentiation and secondary cell wall thickening. By elucidating the comparative transcriptome profiles of Eastern and Western cvs. for the first time, our study provides valuable insights into the differentially expressed genes involved in root development and VC tissue patterning. The identification of key regulatory genes and their roles in these processes represents a significant advancement in our understanding of the evolutionary relations and molecular mechanisms underlying secondary growth of carrot and regulation by vascular cambium.

A Global Landscape of Miniature Inverted-Repeat Transposable Elements in the Carrot Genome

Miniature inverted-repeat transposable elements (MITEs) are the most abundant group of Class II mobile elements in plant genomes. Their presence in genic regions may alter gene structure and expression, providing a new source of functional diversity. Owing to their small size and lack of coding capacity, the identification of MITEs has been demanding. However, the increasing availability of reference genomes and bioinformatic tools provides better means for the genome-wide identification and analysis of MITEs and for the elucidation of their contribution to the evolution of plant genomes. We mined MITEs in the carrot reference genome DH1 using MITE-hunter and developed a curated carrot MITE repository comprising 428 families. Of the 31,025 MITE copies spanning 10.34 Mbp of the carrot genome, 54% were positioned in genic regions. Stowaways and Tourists were frequently present in the vicinity of genes, while Mutator-like MITEs were relatively more enriched in introns. hAT-like MITEs were relatively more frequently associated with transcribed regions, including untranslated regions (UTRs). Some carrot MITE families were shared with other Apiaceae species. We showed that hAT-like MITEs were involved in the formation of new splice variants of insertion-harboring genes. Thus, carrot MITEs contributed to the accretion of new diversity by altering transcripts and possibly affecting the regulation of many genes.

Characteristics of the AT-Hook Motif Containing Nuclear Localized (AHL) Genes in Carrot Provides Insight into Their Role in Plant Growth and Storage Root Development

The AT-hook motif containing nuclear localized (AHL) gene family, controlling various developmental processes, is conserved in land plants. They comprise Plant and Prokaryote Conserved (PPC) domain and one or two AT-hook motifs. DcAHLc1 has been proposed as a candidate gene governing the formation of the carrot storage root. We identified and in-silico characterized carrot AHL proteins, performed phylogenetic analyses, investigated their expression profiles and constructed gene coexpression networks. We found 47 AHL genes in carrot and grouped them into two clades, A and B, comprising 29 and 18 genes, respectively. Within Clade-A, we distinguished three subclades, one of them grouping noncanonical AHLs differing in their structure (two PPC domains) and/or cellular localization (not nucleus). Coexpression network analysis attributed AHLs expressed in carrot roots into four of the 72 clusters, some of them showing a large number of interactions. Determination of expression profiles of AHL genes in various tissues and samples provided basis to hypothesize on their possible roles in the development of the carrot storage root. We identified a group of rapidly evolving noncanonical AHLs, possibly differing functionally from typical AHLs, as suggested by their expression profiles and their predicted cellular localization. We pointed at several AHLs likely involved in the development of the carrot storage root.

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).

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

A de novo transcriptome analysis revealed that photomorphogenic genes are required for carotenoid synthesis in the dark-grown carrot taproot

Carotenoids are terpenoid pigments synthesized by all photosynthetic and some non-photosynthetic organisms. In plants, these lipophilic compounds are involved in photosynthesis, photoprotection, and phytohormone synthesis. In plants, carotenoid biosynthesis is induced by several environmental factors such as light including photoreceptors, such as phytochromes (PHYs) and negatively regulated by phytochrome interacting factors (PIFs). Daucus carota (carrot) is one of the few plant species that synthesize and accumulate carotenoids in the storage root that grows in darkness. Contrary to other plants, light inhibits secondary root growth and carotenoid accumulation suggesting the existence of new mechanisms repressed by light that regulate both processes. To identify genes induced by dark and repressed by light that regulate carotenoid synthesis and carrot root development, in this work an RNA-Seq analysis was performed from dark- and light-grown carrot roots. Using this high-throughput sequencing methodology, a de novo transcriptome model with 63,164 contigs was obtained, from which 18,488 were differentially expressed (DEG) between the two experimental conditions. Interestingly, light-regulated genes are preferably expressed in dark-grown roots. Enrichment analysis of GO terms with DEGs genes, validation of the transcriptome model and DEG analysis through qPCR allow us to hypothesize that genes involved in photomorphogenesis and light perception such as PHYA, PHYB, PIF3, PAR1, CRY2, FYH3, FAR1 and COP1 participate in the synthesis of carotenoids and carrot storage root development.

Gene Regulatory Network Guided Investigations and Engineering of Storage Root Development in Root Crops

The plasticity of plant development relies on its ability to balance growth and stress resistance. To do this, plants have established highly coordinated gene regulatory networks (GRNs) of the transcription factors and signaling components involved in developmental processes and stress responses. In root crops, yields of storage roots are mainly determined by secondary growth driven by the vascular cambium. In relation to this, a dynamic yet intricate GRN should operate in the vascular cambium, in coordination with environmental changes. Despite the significance of root crops as food sources, GRNs wired to mediate secondary growth in the storage root have just begun to emerge, specifically with the study of the radish. Gene expression data available with regard to other important root crops are not detailed enough for us directly to infer underlying molecular mechanisms. Thus, in this review, we provide a general overview of the regulatory programs governing the development and functions of the vascular cambium in model systems, and the role of the vascular cambium on the growth and yield potential of the storage roots in root crops. We then undertake a reanalysis of recent gene expression data generated for major root crops and discuss common GRNs involved in the vascular cambium-driven secondary growth in storage roots using the wealth of information available in Arabidopsis. Finally, we propose future engineering schemes for improving root crop yields by modifying potential key nodes in GRNs.

Mining for Candidate Genes Controlling Secondary Growth of the Carrot Storage Root

Background: Diverse groups of carrot cultivars have been developed to meet consumer demands and industry needs. Varietal groups of the cultivated carrot are defined based on the shape of roots. However, little is known about the genetic basis of root shape determination. Methods: Here, we used 307 carrot plants from 103 open-pollinated cultivars for a genome wide association study to identify genomic regions associated with the storage root morphology. Results: A 180 kb-long region on carrot chromosome 1 explained 10% of the total observed phenotypic variance in the shoulder diameter. Within that region, DcDCAF1 and DcBTAF1 genes were proposed as candidates controlling secondary growth of the carrot storage root. Their expression profiles differed between the cultivated and the wild carrots, likely indicating that their elevated expression was required for the development of edible roots. They also showed higher expression at the secondary root growth stage in cultivars producing thick roots, as compared to those developing thin roots. Conclusions: We provided evidence for a likely involvement of DcDCAF1 and/or DcBTAF1 in the development of the carrot storage root and developed a genotyping assay facilitating the identification of variants in the region on carrot chromosome 1 associated with secondary growth of the carrot root.

Comparative Analysis between Wild and Cultivated Cucumbers Reveals Transcriptional Changes during Domestication Process

The cultivated cucumber (Cucumis sativus L.) was reported to have been developed from a wild cucumber (Cucumis hystrix Chakrav.), nevertheless, these two organisms exhibit noteworthy differences. For example, the wild cucumber is known for its high resistance to different biotic and abiotic stresses. Moreover, the leaves and fruits of the wild cucumber have a bitter taste compared to the cultivated cucumber. These differences could be attributed mainly to the differences in gene expression levels. In the present investigation, we analyzed the RNA-sequencing data to show the differentially expressed genes (DEGs) between the wild and cultivated cucumbers. The identified DEGs were further utilized for Gene Ontology (GO) and pathway enrichment analysis and for identification of transcription factors and regulators. In the results, several enriched GO terms in the biological process, cellular component, and molecular functions categories were identified and various enriched pathways, especially the biosynthesis pathways of secondary products were recognized. Plant-specific transcription factor families were differentially expressed between the wild and cultivated cucumbers. The results obtained provide preliminary evidence for the transcriptional differences between the wild and cultivated cucumbers which developed during the domestication process as a result of natural and/or artificial selection, and they formulate the basis for future genetic research and improvement of the cultivated cucumber.

Transcriptome Profiling Provides Insight into the Genes in Carotenoid Biosynthesis during the Mesocarp and Seed Developmental Stages of Avocado (Persea americana)

Avocado (Persea americana Mill.) is an economically important crop because of its high nutritional value. However, the absence of a sequenced avocado reference genome has hindered investigations of secondary metabolism. For next-generation high-throughput transcriptome sequencing, we obtained 365,615,152 and 348,623,402 clean reads as well as 109.13 and 104.10 Gb of sequencing data for avocado mesocarp and seed, respectively, during five developmental stages. High-quality reads were assembled into 100,837 unigenes with an average length of 847.40 bp (N50 = 1725 bp). Additionally, 16,903 differentially expressed genes (DEGs) were detected, 17 of which were related to carotenoid biosynthesis. The expression levels of most of these 17 DEGs were higher in the mesocarp than in the seed during five developmental stages. In this study, the avocado mesocarp and seed transcriptome were also sequenced using single-molecule long-read sequencing to acquired 25.79 and 17.67 Gb clean data, respectively. We identified 233,014 and 238,219 consensus isoforms in avocado mesocarp and seed, respectively. Furthermore, 104 and 59 isoforms were found to correspond to the putative 11 carotenoid biosynthetic-related genes in the avocado mesocarp and seed, respectively. The isoform numbers of 10 out of the putative 11 genes involved in the carotenoid biosynthetic pathway were higher in the mesocarp than those in the seed. Besides, alpha- and beta-carotene contents in the avocado mesocarp and seed during five developmental stages were also measured, and they were higher in the mesocarp than in the seed, which validated the results of transcriptome profiling. Gene expression changes and the associated variations in gene dosage could influence carotenoid biosynthesis. These results will help to further elucidate carotenoid biosynthesis in avocado.