A bulk segregant transcriptome analysis reveals metabolic and cellular processes associated with Orange allelic variation and fruit β-carotene accumulation in melon fruit

Carotenoid sequestration protein FIBRILLIN participates in CmOR-regulated β-carotene accumulation in melon

Chromoplasts are plant organelles with a unique ability to sequester and store massive carotenoids. Chromoplasts have been hypothesized to enable high levels of carotenoid accumulation due to enhanced sequestration ability or sequestration substructure formation. However, the regulators that control the substructure component accumulation and substructure formation in chromoplasts remain unknown. In melon (Cucumis melo) fruit, β-carotene accumulation in chromoplasts is governed by ORANGE (OR), a key regulator for carotenoid accumulation in chromoplasts. By using comparative proteomic analysis of a high β-carotene melon variety and its isogenic line low-β mutant that is defective in CmOr with impaired chromoplast formation, we identified carotenoid sequestration protein FIBRILLIN1 (CmFBN1) as differentially expressed. CmFBN1 expresses highly in melon fruit tissue. Overexpression of CmFBN1 in transgenic Arabidopsis (Arabidopsis thaliana) containing ORHis that genetically mimics CmOr significantly enhances carotenoid accumulation, demonstrating its involvement in CmOR-induced carotenoid accumulation. Both in vitro and in vivo evidence showed that CmOR physically interacts with CmFBN1. Such an interaction occurs in plastoglobules and results in promoting CmFBN1 accumulation. CmOR greatly stabilizes CmFBN1, which stimulates plastoglobule proliferation and subsequently carotenoid accumulation in chromoplasts. Our findings show that CmOR directly regulates CmFBN1 protein levels and suggest a fundamental role of CmFBN1 in facilitating plastoglobule proliferation for carotenoid sequestration. This study also reveals an important genetic tool to further enhance OR-induced carotenoid accumulation in chromoplasts in crops.

Pan-genome of Citrullus genus highlights the extent of presence/absence variation during domestication and selection

The rich genetic diversity in Citrullus lanatus and the other six species in the Citrullus genus provides important sources in watermelon breeding. Here, we present the Citrullus genus pan-genome based on the 400 Citrullus genus resequencing data, showing that 477 Mb contigs and 6249 protein-coding genes were absent in the Citrullus lanatus reference genome. In the Citrullus genus pan-genome, there are a total of 8795 (30.5%) genes that exhibit presence/absence variations (PAVs). Presence/absence variation (PAV) analysis showed that a lot of gene PAV were selected during the domestication and improvement, such as 53 favorable genes and 40 unfavorable genes were identified during the C. mucosospermus to C. lanatus landrace domestication. We also identified 661 resistance gene analogs (RGAs) in the Citrullus genus pan-genome, which contains 90 RGAs (89 variable and 1 core gene) located on the pangenome additional contigs. By gene PAV-based GWAS, 8 gene presence/absence variations were found associated with flesh color. Finally, based on the results of gene PAV selection analysis between watermelon populations with different fruit colors, we identified four non-reference candidate genes associated with carotenoid accumulation, which had a significantly higher frequency in the white flesh. These results will provide an important source for watermelon breeding.

The win-win effects of nitrification inhibitors on soil-crop systems: Decreasing carbendazim residues but promoting soil bacterial community diversities and stabilities and crop yields

Applying nitrogen (N)-cycling inhibitors is an effective measure to improve N fertilizer utilization efficiency, but the effects of N-cycling inhibitors on fungicide residues in soil-crop systems are unclear. In this study, nitrification inhibitors dicyandiamide (DCD) and 3, 4-dimethylpyrazole phosphate (DMPP) and urease inhibitor N-(n-butyl) thiophosphoric triamide (NBPT) were applied into agricultural soils with fungicide carbendazim applications. The soil abiotic properties, carrot yields, carbendazim residues, bacterial communities and their comprehensive relationships were also quantified. Compared to the control treatment, the DCD and DMPP significantly decreased soil carbendazim residues by 96.2% and 96.0%, and the DMPP and NBPT significantly reduced carrot carbendazim residues by 74.3% and 60.3%, respectively. The nitrification inhibitor applications also generated significant and positive effects on carrot yields and soil bacterial community diversities. The DCD application significantly stimulated soil Bacteroidota and endophytic Myxococcota and modified soil and endophytic bacterial communities. Meanwhile, the DCD and DMPP applications also positively stimulated the co-occurrence network edges of soil bacterial communities by 32.6% and 35.2%, respectively. The linear correlation coefficients between soil carbendazim residues and pH, ETSA and NH₄⁺-N contents were - 0.84, - 0.57 and - 0.80, respectively. The nitrification inhibitor applications generated win-win effects on the soil-crop systems by decreasing carbendazim residues but promoting soil bacterial community diversities and stabilities and crop yields.

Genome-Wide Identification and Chilling Stress Analysis of the NF-Y Gene Family in Melon

The nuclear factor Y (NF-Y) transcription factor contains three subfamilies: NF-YA, NF-YB, and NF-YC. The NF-Y family have been reported to be key regulators in plant growth and stress responses. However, little attention has been given to these genes in melon (Cucumis melo L.). In this study, twenty-five NF-Ys were identified in the melon genome, including six CmNF-YAs, eleven CmNF-YBs, and eight CmNF-YCs. Their basic information (gene location, protein characteristics, and subcellular localization), conserved domains and motifs, and phylogeny and gene structure were subsequently analyzed. Results showed highly conserved motifs exist in each subfamily, which are distinct between subfamilies. Most CmNF-Ys were expressed in five tissues and exhibited distinct expression patterns. However, CmNF-YA6, CmNF-YB1/B2/B3/B8, and CmNF-YC6 were not expressed and might be pseudogenes. Twelve CmNF-Ys were induced by cold stress, indicating the NF-Y family plays a key role in melon cold tolerance. Taken together, our findings provide a comprehensive understanding of CmNF-Y genes in the development and stress response of melon and provide genetic resources for solving the practical problems of melon production.

Transcriptome analysis reveals association of carotenoid metabolism pathway with fruit color in melon

Flesh color is an important quality of melon (Cucumis melo L.) and is determined mainly by carotenoid content, awarding them with colors, aromas, and nutrients. enhancing the nutritional and health benefits of fruits and vegetables for humans. In this study, we performed transcriptomic analysis of two melon inbred line “B-14” (orange-flesh) and “B-6” (white-flesh) at three developmental stages. We observed that the β-carotene content of inbred line “B-6” (14.232 μg/g) was significantly lower than that of inbred line “B-14” (0.534 μg/g). RNA-sequencing and quantitative reverse transcription PCR analyses were performed to identify differentially expressed genes (DEGs) between the two inbred lines at different stages; the DEGs were analyzed using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes databases (KEGG). We identified 33 structural DEGs in different developmental periods of the two lines that were related to carotenoid metabolism. Among them, PSY, Z-ISO, ZDS, CRTISO, CCD4, VDE1, and NCED2 were highly correlated with carotenoid content. Thus, this study provides a basis for molecular mechanism of carotenoid biosynthesis and flesh color in melon fruit.

R3-MYB repressor Mybr97 is a candidate gene associated with the Anthocyanin3 locus and enhanced anthocyanin accumulation in maize

Anthocyanin3 inhibits the anthocyanin and monolignol pathways in maize. Transposon-tagging, RNA-sequencing, and GST-pulldown assays determine Anthocyanin3 may be R3-MYB repressor gene Mybr97. Anthocyanins are colorful molecules receiving recent attention due to their numerous health benefits and applications as natural colorants and nutraceuticals. Purple corn is being investigated as a more economical source of anthocyanins. Anthocyanin3 (A3) is a known recessive intensifier of anthocyanin pigmentation in maize. In this study, anthocyanin content was elevated 100-fold in recessive a3 plants. Two approaches were used to discover candidates involved with the a3 intense purple plant phenotype. First, a large-scale transposon-tagging population was created with a Dissociation (Ds) insertion in the nearby Anthocyanin1 gene. A de novo a3-m1::Ds mutant was generated, and the transposon insertion was found to be located in the promoter of Mybr97, which has homology to R3-MYB repressor CAPRICE in Arabidopsis. Second, a bulked segregant RNA-sequencing population found expression differences between pools of green A3 plants and purple a3 plants. All characterized anthocyanin biosynthetic genes were upregulated in a3 plants along with several genes of the monolignol pathway. Mybr97 was highly downregulated in a3 plants, suggesting its role as a negative regulator of the anthocyanin pathway. Photosynthesis-related gene expression was reduced in a3 plants through an unknown mechanism. Numerous transcription factors and biosynthetic genes were also upregulated and need further investigation. Mybr97 may inhibit anthocyanin synthesis by associating with basic helix-loop helix transcription factors like Booster1. Overall, Mybr97 is the most likely candidate gene for the A3 locus. A3 has a profound effect on the maize plant and has many favorable implications for crop protection, human health, and natural colorant production.

New Insights into MdSPS4-Mediated Sucrose Accumulation under Different Nitrogen Levels Revealed by Physiological and Transcriptomic Analysis

Nitrogen nutrition participates in many physiological processes and understanding the physiological and molecular mechanisms of apple responses to nitrogen is very significant for improving apple quality. This study excavated crucial genes that regulates sugar metabolism in response to nitrogen in apples through physiology and transcriptome analysis, so as to lay a theoretical foundation for improving fruit quality. In this paper, the content of sugar and organic acid in apple fruit at different developmental periods under different nitrogen levels (0, 150, 300, and 600 kg·hm^-2) were determined. Then, the transcriptomic analysis was performed in 120 days after bloom (DAB) and 150 DAB. The results showed that the fructose and glucose content were the highest at 120 DAB under 600 kg·hm^-2 nitrogen level. Meanwhile, different nitrogen treatments decreased malate content in 30 and 60 DAB. RNA-seq analysis revealed a total of 4537 UniGenes were identified as differentially expressed genes (DEGs) under nitrogen treatments. Among these DEGs, 2362 (52.06%) were up-regulated and 2175 (47.94%) were down-regulated. The gene co-expression clusters revealed that most DEGs were significantly annotated in the photosynthesis, glycolysis/gluconeogenesis, pyruvate metabolism, carbon metabolism, carbon fixation in photosynthetic organisms and plant hormone signal transduction pathways. The key transcription factor genes (ERF, NAC, WRKY, and C2H2 genes) were differentially expressed in apple fruit. Sugar and acid metabolism-related genes (e.g., HXK1, SPS4, SS2, PPC16-2, and MDH2 genes) exhibited significantly up-regulated expression at 120 DAB, whereas they were down-regulated at 150 DAB. Furthermore, the MdSPS4 gene overexpression positively promoted sucrose accumulation in apple callus and fruit. In conclusion, the combinational analysis of transcriptome and the functional validation of the MdSPS4 gene provides new insights into apple responses to different nitrogen levels.

Multispecies transcriptomes reveal core fruit development genes

During angiosperm evolution there have been repeated transitions from an ancestral dry fruit to a derived fleshy fruit, often with dramatic ecological and economic consequences. Following the transition to fleshy fruits, domestication may also dramatically alter the fruit phenotype via artificial selection. Although the morphologies of these fruits are well documented, relatively less is known about the molecular basis of these developmental and evolutionary shifts. We generated RNA-seq libraries from pericarp tissue of desert tobacco and both cultivated and wild tomato species at common developmental time points and combined this with corresponding, publicly available data from Arabidopsis and melon. With this broadly sampled dataset consisting of dry/fleshy fruits and wild/domesticated species, we applied novel bioinformatic methods to investigate conserved and divergent patterns of gene expression during fruit development and evolution. A small set of 121 orthologous "core" fruit development genes show a common pattern of expression across all five species. These include key players in developmental patterning such as orthologs of KNOLLE, PERIANTHIA, and ARGONAUTE7. GO term enrichment suggests that these genes function in basic cell division processes, cell wall biosynthesis, and developmental patterning. We furthermore uncovered a number of "accessory" genes with conserved expression patterns within but not among fruit types, and whose functional enrichment highlights the conspicuous differences between these phenotypic classes. We observe striking conservation of gene expression patterns despite large evolutionary distances, and dramatic phenotypic shifts, suggesting a conserved function for a small subset of core fruit development genes.

Role of reactive oxygen species in lesion mimic formation and conferred basal resistance to Fusarium graminearum in barley lesion mimic mutant 5386

This study investigated the barley lesion mimic mutant (LMM) 5386, evidenced by a leaf brown spot phenotype localized on the chromosome 3H, and its conferred basal resistance to Fusarium graminearum. RNA-seq analysis identified 1453 genes that were differentially expressed in LMM 5386 compared to those in the wild type. GO and KEGG functional annotations suggested that lesion mimic formation was mediated by pathways involving oxidation reduction and glutathione metabolism. Additionally, reactive oxygen species (ROS) accumulation in brown spots was substantially higher in LMM 5386 than in the wild-type plant; therefore, antioxidant competence, which is indicated by ROS accumulation, was significantly lower in LMM 5386. Furthermore, the reduction of glycine in LMM 5386 inhibited glutathione biosynthesis. These results suggest that the decrease in antioxidant competence and glutathione biosynthesis caused considerable ROS accumulation, leading to programmed cell death, which eventually reduced the yield components in LMM 5386.

Transcription Factor CmNAC34 Regulated CmLCYB-Mediated β-Carotene Accumulation during Oriental Melon Fruit Ripening

Ripened oriental melon (Cucumis melo) with orange-colored flesh is rich in β-carotene. Lycopene β-cyclase (LCYB) is the synthetic enzyme that directly controls the massive accumulation of β-carotene. However, the regulatory mechanism underlying the CmLCYB-mediated β-carotene accumulation in oriental melon is fairly unknown. Here, we screened and identified a transcription factor, CmNAC34, by combining bioinformatics analysis and yeast one-hybrid screen with CmLCYB promoter. CmNAC34 was located in the nucleus and acted as a transcriptional activator. The expression profile of CmNAC34 was consistent with that of CmLCYB during the fruit ripening. Additionally, the transient overexpression of CmNAC34 in oriental melon fruit promoted the expression of CmLCYB and enhanced β-carotene concentration, while transient silence of CmNAC34 in fruit was an opposite trend, which indicated CmNAC34 could modulate CmLCYB-mediated β-carotene biosynthesis in oriental melon. Finally, the yeast one-hybrid (Y1H), electrophoretic mobility shift assay (EMSA), β-glucuronidase (GUS) analysis assay, and luciferase reporter (LUC) assay indicated that CmNAC34 could bind to the promoter of CmLCYB and positively regulated the CmLCYB transcription level. These findings suggested that CmNAC34 acted as an activator to regulate β-carotene accumulation by directly binding the promoter of CmLCYB, which provides new insight into the regulatory mechanism of carotenoid metabolism during the development and ripening of oriental melon.

Harnessing the potential of bulk segregant analysis sequencing and its related approaches in crop breeding

Most plant traits are governed by polygenes including both major and minor genes. Linkage mapping and positional cloning have contributed greatly to mapping genomic loci controlling important traits in crop species. However, they are low-throughput, time-consuming, and have low resolution due to which their efficiency in crop breeding is reduced. In this regard, the bulk segregant analysis sequencing (BSA-seq) and its related approaches, viz., quantitative trait locus (QTL)-seq, bulk segregant RNA-Seq (BSR)-seq, and MutMap, have emerged as efficient methods to identify the genomic loci/QTLs controlling specific traits at high resolution, accuracy, reduced time span, and in a high-throughput manner. These approaches combine BSA with next-generation sequencing (NGS) and enable the rapid identification of genetic loci for qualitative and quantitative assessments. Many previous studies have shown the successful identification of the genetic loci for different plant traits using BSA-seq and its related approaches, as discussed in the text with details. However, the efficiency and accuracy of the BSA-seq depend upon factors like sequencing depth and coverage, which enhance the sequencing cost. Recently, the rapid reduction in the cost of NGS together with the expected cost reduction of third-generation sequencing in the future has further increased the accuracy and commercial applicability of these approaches in crop improvement programs. This review article provides an overview of BSA-seq and its related approaches in crop breeding together with their merits and challenges in trait mapping.

Identification of Key Gene Networks Controlling Soluble Sugar and Organic Acid Metabolism During Oriental Melon Fruit Development by Integrated Analysis of Metabolic and Transcriptomic Analyses

Oriental melon (Cucumis melo var. acidulus) is one of the most economically important fruit crops worldwide. To elucidate the molecular basis related to soluble sugar and organic acid metabolism in the fruits of two oriental melon cultivars with different sweetness, we performed integrated metabolomic and transcriptomic analyses of the fruits of 'Tianbao' (A) with high sweetness and 'Xiaocuigua' (B) with low sweetness at different ripening stages. The high accumulation of sucrose, D-glucose, D-(+)-raffinose, and the relatively lower citric acid and malic acid might contribute to the sweet taste of A. By screening the differentially expressed genes (DEGs) and correlation analysis of the DEGs and differentially accumulated metabolites, we deduced that the B cultivar might promote the conversion of glucose and fructose into intermediate compounds for downstream processes such as glycolysis. The tricarboxylic acid (TCA) cycle might also be enhanced compared to A, thus resulting in the differential accumulation of soluble sugars and organic acids, ultimately causing the taste difference between the two oriental melon cultivars. Our finding provides important information for further exploring the metabolic mechanisms of soluble sugars and organic acids in oriental melon.

Plant carotenoids: recent advances and future perspectives

Carotenoids are isoprenoid metabolites synthesized de novo in all photosynthetic organisms. Carotenoids are essential for plants with diverse functions in photosynthesis, photoprotection, pigmentation, phytohormone synthesis, and signaling. They are also critically important for humans as precursors of vitamin A synthesis and as dietary antioxidants. The vital roles of carotenoids to plants and humans have prompted significant progress toward our understanding of carotenoid metabolism and regulation. New regulators and novel roles of carotenoid metabolites are continuously revealed. This review focuses on current status of carotenoid metabolism and highlights recent advances in comprehension of the intrinsic and multi-dimensional regulation of carotenoid accumulation. We also discuss the functional evolution of carotenoids, the agricultural and horticultural application, and some key areas for future research.

Multi-strategy engineering greatly enhances provitamin A carotenoid accumulation and stability in Arabidopsis seeds

Staple grains with low levels of provitamin A carotenoids contribute to the global prevalence of vitamin A deficiency and therefore are the main targets for provitamin A biofortification. However, carotenoid stability during both seed maturation and postharvest storage is a serious concern for the full benefits of carotenoid biofortified grains. In this study, we utilized Arabidopsis as a model to establish carotenoid biofortification strategies in seeds. We discovered that manipulation of carotenoid biosynthetic activity by seed-specific expression of Phytoene synthase (PSY) increases both provitamin A and total carotenoid levels but the increased carotenoids are prone to degradation during seed maturation and storage, consistent with previous studies of provitamin A biofortified grains. In contrast, stacking with Orange (OR ^His ), a gene that initiates chromoplast biogenesis, dramatically enhances provitamin A and total carotenoid content and stability. Up to 65- and 10-fold increases of β-carotene and total carotenoids, respectively, with provitamin A carotenoids composing over 63% were observed in the seeds containing OR ^His and PSY. Co-expression of Homogentisate geranylgeranyl transferase (HGGT) with OR ^His and PSY further increases carotenoid accumulation and stability during seed maturation and storage. Moreover, knocking-out of β-carotene hydroxylase 2 (BCH2) by CRISPR/Cas9 not only potentially facilitates β-carotene accumulation but also minimizes the negative effect of carotenoid over production on seed germination. Our findings provide new insights into various processes on carotenoid accumulation and stability in seeds and establish a multiplexed strategy to simultaneously target carotenoid biosynthesis, turnover, and stable storage for carotenoid biofortification in crop seeds.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-021-00046-1.

Comparative transcriptome analyses shed light on carotenoid production and plastid development in melon fruit

Carotenoids, such as β-carotene, accumulate in chromoplasts of various fleshy fruits, awarding them with colors, aromas, and nutrients. The Orange (CmOr) gene controls β-carotene accumulation in melon fruit by posttranslationally enhancing carotenogenesis and repressing β-carotene turnover in chromoplasts. Carotenoid isomerase (CRTISO) isomerizes yellow prolycopene into red lycopene, a prerequisite for further metabolism into β-carotene. We comparatively analyzed the developing fruit transcriptomes of orange-colored melon and its two isogenic EMS-induced mutants, low-β (Cmor) and yofi (Cmcrtiso). The Cmor mutation in low-β caused a major transcriptomic change in the mature fruit. In contrast, the Cmcrtiso mutation in yofi significantly changed the transcriptome only in early fruit developmental stages. These findings indicate that melon fruit transcriptome is primarily altered by changes in carotenoid metabolic flux and plastid conversion, but minimally by carotenoid composition in the ripe fruit. Clustering of the differentially expressed genes into functional groups revealed an association between fruit carotenoid metabolic flux with the maintenance of the photosynthetic apparatus in fruit chloroplasts. Moreover, large numbers of thylakoid localized photosynthetic genes were differentially expressed in low-β. CmOR family proteins were found to physically interact with light-harvesting chlorophyll a-b binding proteins, suggesting a new role of CmOR for chloroplast maintenance in melon fruit. This study brings more insights into the cellular and metabolic processes associated with fruit carotenoid accumulation in melon fruit and reveals a new maintenance mechanism of the photosynthetic apparatus for plastid development.

Encapsulation of pea protein in an alginate matrix by cold set gelation method and use of the capsules in fruit juices

Plant-based proteins gained importance in recent years due to the increase in the awareness of healthy diet and in the consumption of plant-based foods. However, some features of plant-based proteins like the undesirable odor and flavor affect the sensorial properties of protein containing foods. Therefore, encapsulation of these proteins could be a good strategy to tackle with this problem. The objective of this study was to design microcapsules (beads) consisting of pea protein by using sodium alginate and to investigate the effect of different alginate concentrations (1.0, 1.5, and 2.0%) on the protein content, encapsulation efficiency, particle size, bead stability, and the morphology of the capsules and then add them to different fruit juices (pomegranate and melon) and examine the release behavior from the capsules. Rheological behavior of the juices including pectin were also investigated. TD- nuclear magnetic resonance relaxometry analysis through T₂ relaxation times was conducted on the capsules to observe the changes in the beads. In conclusion, alginate was found to be a suitable encapsulation coating for pea protein. Beads containing 1% alginate concentration was found to be the most effective with respect to protein content and bead stability. PRACTICAL APPLICATION: This study aims to design and characterize pea protein containing microcapsules capsules and their utilization in fruit juices. The study itself focused on a specific application on the fruit juices.

Genome-wide identification of microRNAs involved in the regulation of fruit ripening and climacteric stages in melon (Cucumis melo)

Fruit ripening is influenced by multiple plant hormones and the regulation of genes. However, studies on posttranscriptional regulators (e.g., miRNAs) of fruit growth and ripening are limited. We used miRNA sequencing and degradome methods to identify miRNAs and their target genes in melon (Cucumis melo cv. Hetao melon). A total of 61 conserved miRNAs and 36 novel miRNAs were identified from fruit growth, ripening, climacteric, and postclimacteric developmental stage samples, of which 32 conserved miRNAs were differentially expressed between developmental stage samples. Sixty-two target genes of 43 conserved miRNAs and 1 novel miRNA were identified from degradome sequencing. To further investigate miRNA influencing fruit ripening, transgenic melon plants overexpressing pre-cme-miR393 (cme-miR393-OE) were generated and characterized. The results showed that fruit ripening was delayed in cme-miR393-OE transgenic lines compared to nontransgenic fruits. The target of cme-miR393 was also identified, and the expression of CmAFB2 was repressed in transgenic plants. These results provide evidence that miRNA regulates melon fruit ripening and provide potential targets to improve the horticultural traits of melon fruit.

Transcriptome profiling of non-climacteric 'yellow' melon during ripening: insights on sugar metabolism

BACKGROUND

The non-climacteric 'Yellow' melon (Cucumis melo, inodorus group) is an economically important crop and its quality is mainly determined by the sugar content. Thus, knowledge of sugar metabolism and its related pathways can contribute to the development of new field management and post-harvest practices, making it possible to deliver better quality fruits to consumers.

RESULTS

The RNA-seq associated with RT-qPCR analyses of four maturation stages were performed to identify important enzymes and pathways that are involved in the ripening profile of non-climacteric 'Yellow' melon fruit focusing on sugar metabolism. We identified 895 genes 10 days after pollination (DAP)-biased and 909 genes 40 DAP-biased. The KEGG pathway enrichment analysis of these differentially expressed (DE) genes revealed that 'hormone signal transduction', 'carbon metabolism', 'sucrose metabolism', 'protein processing in endoplasmic reticulum' and 'spliceosome' were the most differentially regulated processes occurring during melon development. In the sucrose metabolism, five DE genes are up-regulated and 12 are down-regulated during fruit ripening.

CONCLUSIONS

The results demonstrated important enzymes in the sugar pathway that are responsible for the sucrose content and maturation profile in non-climacteric 'Yellow' melon. New DE genes were first detected for melon in this study such as invertase inhibitor LIKE 3 (CmINH3), trehalose phosphate phosphatase (CmTPP1) and trehalose phosphate synthases (CmTPS5, CmTPS7, CmTPS9). Furthermore, the results of the protein-protein network interaction demonstrated general characteristics of the transcriptome of young and full-ripe melon and provide new perspectives for the understanding of ripening.

Quantitative Trait Loci (QTL) Analysis of Fruit and Agronomic Traits of Tropical Pumpkin (Cucurbita moschata) in an Organic Production System

Interest in the development of organically grown vegetable crops has risen over the past decades due to consumer preferences. However, most crops that have desirable consumer traits have been bred in conventional growing conditions, and their transfer to an organic setting is challenging. Here, the organically grown Hawaiian pumpkin (Cucurbita moschata) accession ‘Shima’ was crossed with the conventionally grown Puerto Rican variety ‘Taina Dorada’ to develop a backcross (BC1) population, where ‘Shima’ was the recurrent parent. A total of 202 BC1 (‘Shima’ X F1) progenies were planted in a certified organic field, and twelve traits were evaluated. We used genotype-by-sequencing (GBS) to identify the Quantitative Trait Loci (QTL) associated with insect tolerance along with commercially desirable traits. A total of 1582 single nucleotide polymorphisms (SNPs) were identified, from which 711 SNPs were used to develop a genetic map and perform QTL mapping. Reads associated with significant QTLs were aligned to the publicly available Cucurbita moschata genome and identified several markers linked to genes that have been previously reported to be associated with that trait in other crop systems, such as melon (Cucumis melo L.). This research provides a resource for marker-assisted selection (MAS) efforts in Cucurbita moschata, as well as serving as a model study to improve cultivars that are transitioning from a conventional to an organic setting.

Toward the 'golden' era: The status in uncovering the regulatory control of carotenoid accumulation in plants

Carotenoids are essential pigments to plants and important natural products to humans. Carotenoids as both primary and specialized metabolites fulfill multifaceted functions in plants. As such, carotenoid accumulation (a net process of biosynthesis, degradation and sequestration) is subjected to complicated regulation throughout plant life cycle in response to developmental and environmental signals. Investigation of transcriptional regulation of carotenoid metabolic genes remains the focus in understanding the regulatory control of carotenoid accumulation. While discovery of bona fide carotenoid metabolic regulators is still challenging, the recent progress of identification of various transcription factors and regulators helps us to construct hierarchical regulatory network of carotenoid accumulation. The elucidation of carotenoid regulatory mechanisms at protein level and in chromoplast provides some insights into post-translational regulation of carotenogenic enzymes and carotenoid sequestration in plastid sink. This review briefly describes the pathways and main flux-controlling steps for carotenoid accumulation in plants. It highlights our recent understanding of the regulatory mechanisms underlying carotenoid accumulation at both transcriptional and post-translational levels. It also discusses the opportunities to expand toolbox for further shedding light upon the intrinsic regulation of carotenoid accumulation in plants.

The anther-specific CYP704B is potentially responsible for MSG26 male sterility in barley

Plant male sterility is a valuable trait in breeding and hybrid seed production. The barley male-sterility gene msg26 was mapped to a 0.02-cM region that anchors to a 506-kb low-quality assembly between two cleaved amplified polymorphic sequence (CAPS) markers, SP1M14 and SP1M49. The barley gene HORVU4Hr1G074840, which encodes a putative cytochrome P450 CYP704B protein, appears to be a strong candidate for the MSG26 trait. Barley (Hordeum vulgare L.) is an important cereal crop worldwide. Traditional breeding in barley is time-consuming and labor-intensive. The use of male-sterile genotypes may significantly improve the efficacy of hybrid breeding and seed production. The barley accession 'GSHO745' is a spontaneous male-sterile mutant from the barley variety, 'Unitan'. The male sterility in 'GSHO745' is controlled by the recessive gene, msg26 (originally named as ms-u). We revealed that the barley plants homozygous for msg26 proceeded normally through Meiosis II until the tetrad stage, but became fully defective in the late uninucleate microspores and developed pollen-less anthers. Using seven barley F₂ populations, we mapped MSG26 to a 0.02-cM region that anchored to a 506-kb low-quality assembly between two cleaved amplified polymorphic sequence markers, SP1M14 and SP1M49. The HORVU4Hr1G074840 gene that encodes a putative cytochrome P450 protein (CYP704B) was identified as the most plausible candidate for MSG26. First, HORVU4Hr1G074840 is located in a collinear region of the rice CYP704B2 and the maize CYP704B1. Both of these genes are essential for male gamete production. Second, the male-sterile allele of HORVU4Hr1G074840 in GSHO745 contained a 4-bp deletion in the last exon. The resulting frame shift causes a Gly436Gln substitution, scrambles the sequence of the remainder of the protein, and forms a new termination site at the 70th triplet of the shifted reading frame. We thus called the variant protein CYP704B:p.G436Qfs*70. Third, the barley HORVU4Hr1G074840 gene was specifically expressed in anthers. Altogether, HORVU4Hr1G074840 represents a strong candidate for MSG26 in barley.

Quantitative trait loci analysis of hormone levels in Arabidopsis roots

Quantitative trait loci (QTL) analyses for five groups of hormones, including cytokinins in Arabidopsis roots were performed using recombinant inbred lines (Ler×Cvi). Significant QTLs were detected for cytokinins, jasmonic acid and salicylic acid. Separate analysis of two sub-populations, viz., vegetative and flowering plants revealed that many of the QTLs were development-specific. Using near-isogenic lines, several significant QTLs were confirmed; three co-localized QTL regions were responsible for determining several cytokinin metabolites. Using a knock-out plant, a functional role of zeatin N-glucosyltransferase gene (UGT76C2) underlying a large-effect QTL for levels of tZ-N-glucosides and tZRMP was evaluated in the metabolism of cytokinins. Pleotropic effects of this gene were found for cytokinin levels in both roots and leaves, but significant changes of morphological traits were observed only in roots. Hormone QTL analysis reveals development-specific and organ-dependent aspects of the regulation of plant hormone content and metabolism.

Rapid identification of a candidate nicosulfuron sensitivity gene (Nss) in maize (Zea mays L.) via combining bulked segregant analysis and RNA-seq

A candidate nicosulfuron sensitivity gene Nss was identified by combining bulked segregant analysis and RNA-seq. Multiple mutations of this gene were discovered in nicosulfuron-sensitive maize compared with the tolerant. It has been demonstrated that variabilities exist in maize response to nicosulfuron. Two nicosulfuron-sensitive inbred lines (HB39, HB41) and two tolerant inbred lines (HB05, HB09) were identified via greenhouse and field trials. Genetic analysis indicated that the sensitivity to nicosulfuron in maize was controlled by a single, recessive gene. To precisely and rapidly map the nicosulfuron sensitivity gene (Nss), two independent F₂ segregating populations, Population A (HB41 × HB09) and Population B (HB39 × HB05), were constructed. By applying bulked segregant RNA-Seq (BSR-Seq), the Nss gene was, respectively, mapped on the short arm of chromosome 5 (chr5: 1.1-15.3 Mb) and (chr5: 0.5-18.2 Mb) using two populations, with 14.2 Mb region in common. Further analysis revealed that there were 43 and 119 differentially expressed genes in the mapping intervals, with 18 genes in common. Gene annotation results showed that a cytochrome P450 gene (CYP81A9) appeared to be the candidate gene of Nss associated with nicosulfuron sensitivity in maize. Sequence analysis demonstrated that two common deletion mutations existed in the sensitive maize, which might lead to the nicosulfuron sensitivity in maize. The results will make valuable contributions to the understanding of molecular mechanism of herbicide sensitivity in maize.

Fine Mapping of Lycopene Content and Flesh Color Related Gene and Development of Molecular Marker-Assisted Selection for Flesh Color in Watermelon (Citrullus lanatus)

Lycopene content and flesh color are important traits determined by a network of carotenoid metabolic pathways in watermelon. Based on our previous study of genetic inheritance and initial mapping using F₂ populations of LSW-177 (red flesh) × cream of Saskatchewan (pale yellow flesh), red flesh color was controlled by one recessive gene regulating red and pale yellow pigmentation, and a candidate region related to lycopene content was detected spanning a 392,077-bp region on chromosome 4. To obtain a more precise result for further study, three genetic populations and a natural panel of 81 watermelon accessions with different flesh colors were used in this research. Herein, we narrowed the preliminary mapping region to 41,233 bp with the linkage map generated from F₂ populations of LSW-177 (red flesh) × cream of Saskatchewan (pale yellow flesh) with 1,202 individuals. Two candidate genes, Cla005011 and Cla005012, were found in the fine mapping region; therein Cla005011 was a key locus annotated as a lycopene β-cyclase gene. Phylogenetic tree analysis showed that Cla005011 was the closest relative gene in gourd. LSW-177 × PI 186490 (white flesh) and another BC₁ population derived from garden female (red flesh) × PI 186490 were generated to verify the accuracy of the red flesh candidate gene region. By analyzing the expression levels of candidate genes in different developmental stages of different color watermelon varieties, Cla005011 for the expression differences was not the main reason for the flesh color variation between COS and LSW-177. This indicated that the LCYB gene might regulate fruit color changes at the protein level. A new marker-assisted selection system to identify red and yellow flesh colors in watermelon was developed with flesh color-specific CAPS markers and tested in 81 watermelon accessions.

Ectopic expression of ORANGE promotes carotenoid accumulation and fruit development in tomato

Carotenoids are critically important to plants and humans. The ORANGE (OR) gene is a key regulator for carotenoid accumulation, but its physiological roles in crops remain elusive. In this study, we generated transgenic tomato ectopically overexpressing the Arabidopsis wild-type OR (AtOR^WT ) and a 'golden SNP'-containing OR (AtOR^His ). We found that AtOR^His initiated chromoplast formation in very young fruit and stimulated carotenoid accumulation at all fruit developmental stages, uncoupled from other ripening activities. The elevated levels of carotenoids in the AtOR lines were distributed in the same subplastidial fractions as in wild-type tomato, indicating an adaptive response of plastids to sequester the increased carotenoids. Microscopic analysis revealed that the plastid sizes were increased in both AtOR^WT and AtOR^His lines at early fruit developmental stages. Moreover, AtOR overexpression promoted early flowering, fruit set and seed production. Ethylene production and the expression of ripening-associated genes were also significantly increased in the AtOR transgenic fruit at ripening stages. RNA-Seq transcriptomic profiling highlighted the primary effects of OR overexpression on the genes in the processes related to RNA, protein and signalling in tomato fruit. Taken together, these results expand our understanding of OR in mediating carotenoid accumulation in plants and suggest additional roles of OR in affecting plastid size as well as flower and fruit development, thus making OR a target gene not only for nutritional biofortification of agricultural products but also for alteration of horticultural traits.

The Role of Carotenogenic Metabolic Flux in Carotenoid Accumulation and Chromoplast Differentiation: Lessons From the Melon Fruit

Carotenoids have various roles in plant physiology. Plant carotenoids are synthesized in plastids and are highly abundant in the chromoplasts of ripening fleshy fruits. Considerable research efforts have been devoted to elucidating mechanisms that regulate carotenoid biosynthesis, yet, little is known about the mechanism that triggers storage capacity, mainly through chromoplast differentiation. The Orange gene (OR) product stabilizes phytoene synthase protein (PSY) and triggers chromoplast differentiation. OR underlies carotenoid accumulation in orange cauliflower and melon. The OR's 'golden SNP', found in melon, alters the highly evolutionary conserved Arginine¹⁰⁸ to Histidine and controls β-carotene accumulation in melon fruit, in a mechanism yet to be elucidated. We have recently shown that similar carotenogenic metabolic flux is active in non-orange and orange melon fruit. This flux probably leads to carotenoid turnover but known carotenoid turnover products are not detected in non-orange fruit. Arrest of this metabolic flux, using chemical inhibitors or mutations, induces carotenoid accumulation and biogenesis of chromoplasts, regardless of the allelic state of OR. We suggest that the 'golden SNP' induces β-carotene accumulation probably by negatively affecting the capacity to synthesize downstream compounds. The accumulation of carotenoids induces chromoplast biogenesis through a metabolite-induced mechanism. Carotenogenic turnover flux can occur in non-photosynthetic tissues, which do not accumulate carotenoids. Arrest of this flux by the 'golden SNP' or other flux-arrest mutations is a potential tool for the biofortification of agricultural products with carotenoids.

The Role of Orange Gene in Carotenoid Accumulation: Manipulating Chromoplasts Toward a Colored Future

Carotenoids are isoprenoid pigments synthesized in plants, algae, and photosynthetic bacteria and fungus. Their role is essential in light capture, photoprotection, pollinator attraction, and phytohormone production. Furthermore, they can regulate plant development when they are processed as small signaling molecules. Due to their importance for human health, as promoters of the immune system and antioxidant activity, carotenoids have been used in the pharmaceutical, food, and nutraceutical industries. Regulation of carotenoid synthesis and accumulation has been extensively studied. Excellent work has been done unraveling the mode of action of phytoene synthase (PSY), a rate-limiting enzyme of carotenoid biosynthesis pathway, in model species and staple crops. Lately, interest has been turned to Orange protein and its interaction with PSY during carotenoid biosynthesis. Discovered as a dominant mutation in Brassica oleracea, Orange protein regulates carotenoid accumulation by posttranscriptionally regulating PSY, promoting the formation of carotenoid-sequestering structures, and also preventing carotenoid degradation. Furthermore, Orange protein contributes to homeostasis regulation, improving plant tolerance to abiotic stress. In this mini review, the focus is made on recent evidence that elucidates Orange protein mode of action and expression in different plant species. Additionally, strategies are proposed to modify Orange gene by utilization of genome editing techniques. A better understanding of carotenoid biosynthesis and accumulation will lead to a positive impact on the development of healthy food for a growing population.

Fruit flesh volatile and carotenoid profile analysis within the Cucumis melo L. species reveals unexploited variability for future genetic breeding

BACKGROUND

Aroma profile and carotenoids content of melon flesh are two important aspects influencing the quality of this fruit that have been characterized using only selected genotypes. However, the extant variability of the whole species remains unknown.

RESULTS

A complete view of the volatile/carotenoid profiles of melon flesh was obtained analyzing 71 accessions, representing the whole diversity of the species. Gas chromatography-mass spectrometry and high-performance liquid chromatography were used to analyze 200 volatile compounds and five carotenoids. Genotypes were classified into two main clusters (high/low aroma), but with a large diversity of differential profiles within each cluster, consistent with the ripening behavior, flesh color and proposed evolutionary and breeding history of the different horticultural groups.

CONCLUSION

Our results highlight the huge amount of untapped aroma diversity of melon germplasm, especially of non-commercial types. Also, landraces with high nutritional value with regard to carotenoids have been identified. All this knowledge will encourage melon breeding, facilitating the selection of the genetic resources more appropriate to develop cultivars with new aromatic profiles or to minimize the impact of breeding on melon quality. The newly characterized sources provide the basis for further investigations into specific genes/alleles contributing to melon flesh quality. © 2018 Society of Chemical Industry.

Orange: a target gene for regulating carotenoid homeostasis and increasing plant tolerance to environmental stress in marginal lands

Carotenoids play essential roles in various light-harvesting processes in plants and help protect the photosynthetic machinery from photo-oxidative damage. Orange genes, which play a role in carotenoid accumulation, have recently been isolated from several plant species, and their functions have been intensively investigated. The Orange gene (IbOr) of sweet potato [Ipomoea batatas (L.) Lam] helps maintain carotenoid homeostasis to improve plant tolerance to environmental stress. IbOr, a protein with strong holdase chaperone activity, directly interacts with phytoene synthase, a key enzyme involved in carotenoid biosynthesis, in plants under stress conditions, resulting in increased carotenoid accumulation and abiotic stress tolerance. In addition, IbOr interacts with the oxygen-evolving enhancer protein 2-1, a member of a protein complex in photosystem II that is denatured under heat stress. Transgenic sweet potato plants overexpressing IbOr showed enhanced tolerance to high temperatures (47 °C). These findings indicate that IbOr protects plants from environmental stress not only by controlling carotenoid biosynthesis, but also by directly stabilizing photosystem II. In this review, we discuss the functions of IbOr and Or proteins in other plant species and their possible biotechnological applications for molecular breeding for sustainable development on marginal lands.

Deciphering genetic factors that determine melon fruit-quality traits using RNA-Seq-based high-resolution QTL and eQTL mapping

Combined quantitative trait loci (QTL) and expression-QTL (eQTL) mapping analysis was performed to identify genetic factors affecting melon (Cucumis melo) fruit quality, by linking genotypic, metabolic and transcriptomic data from a melon recombinant inbred line (RIL) population. RNA sequencing (RNA-Seq) of fruit from 96 RILs yielded a highly saturated collection of > 58 000 single-nucleotide polymorphisms, identifying 6636 recombination events that separated the genome into 3663 genomic bins. Bin-based QTL analysis of 79 RILs and 129 fruit-quality traits affecting taste, aroma and color resulted in the mapping of 241 QTL. Thiol acyltransferase (CmThAT1) gene was identified within the QTL interval of its product, S-methyl-thioacetate, a key component of melon fruit aroma. Metabolic activity of CmThAT1-encoded protein was validated in bacteria and in vitro. QTL analysis of flesh color intensity identified a candidate white-flesh gene (CmPPR1), one of two major loci determining fruit flesh color in melon. CmPPR1 encodes a member of the pentatricopeptide protein family, involved in processing of RNA in plastids, where carotenoid and chlorophyll pigments accumulate. Network analysis of > 12 000 eQTL mapped for > 8000 differentially expressed fruit genes supported the role of CmPPR1 in determining the expression level of plastid targeted genes. We highlight the potential of RNA-Seq-based QTL analysis of small to moderate size, advanced RIL populations for precise marker-assisted breeding and gene discovery. We provide the following resources: a RIL population genotyped with a unique set of SNP markers, confined genomic segments that harbor QTL governing 129 traits and a saturated set of melon eQTLs.

Carotenoid Metabolism in Plants: The Role of Plastids

Carotenoids are indispensable to plants and critical in human diets. Plastids are the organelles for carotenoid biosynthesis and storage in plant cells. They exist in various types, which include proplastids, etioplasts, chloroplasts, amyloplasts, and chromoplasts. These plastids have dramatic differences in their capacity to synthesize and sequester carotenoids. Clearly, plastids play a central role in governing carotenogenic activity, carotenoid stability, and pigment diversity. Understanding of carotenoid metabolism and accumulation in various plastids expands our view on the multifaceted regulation of carotenogenesis and facilitates our efforts toward developing nutrient-enriched food crops. In this review, we provide a comprehensive overview of the impact of various types of plastids on carotenoid biosynthesis and accumulation, and discuss recent advances in our understanding of the regulatory control of carotenogenesis and metabolic engineering of carotenoids in light of plastid types in plants.

A high-density linkage map and QTL mapping of fruit-related traits in pumpkin (Cucurbita moschata Duch.)

Pumpkin (Cucurbita moschata) is an economically worldwide crop. Few quantitative trait loci (QTLs) were reported previously due to the lack of genomic and genetic resources. In this study, a high-density linkage map of C. moschata was structured by double-digest restriction site-associated DNA sequencing, using 200 F2 individuals of CMO-1 × CMO-97. By filtering 74,899 SNPs, a total of 3,470 high quality SNP markers were assigned to the map spanning a total genetic distance of 3087.03 cM on 20 linkage groups (LGs) with an average genetic distance of 0.89 cM. Based on this map, both pericarp color and strip were fined mapped to a novel single locus on LG8 in the same region of 0.31 cM with phenotypic variance explained (PVE) of 93.6% and 90.2%, respectively. QTL analysis was also performed on carotenoids, sugars, tuberculate fruit, fruit diameter, thickness and chamber width with a total of 12 traits. 29 QTLs distributed in 9 LGs were detected with PVE from 9.6% to 28.6%. It was the first high-density linkage SNP map for C. moschata which was proved to be a valuable tool for gene or QTL mapping. This information will serve as significant basis for map-based gene cloning, draft genome assembling and molecular breeding.

Bulked Segregant RNA-seq Reveals Differential Expression and SNPs of Candidate Genes Associated with Waterlogging Tolerance in Maize

Waterlogging has increasingly become one of the major constraints to maize productivity in some maize production zones because it causes serious yield loss. Bulked segregant RNA-seq (BSR-seq) has been widely applied to profile candidate genes and map associated Single Nucleotide Polymorphism (SNP) markers in many species. In this study, 10 waterlogging sensitive and eight tolerant inbred lines were selected from 60 maize inbred lines with waterlogging response determined and preselected by the International Maize and Wheat Improvement Center (CIMMYT) from over 400 tropical maize inbred lines. BSR-seq was performed to identify differentially expressed genes and SNPs associated with waterlogging tolerance. Upon waterlogging stress, 354 and 1094 genes were differentially expressed in the tolerant and sensitive pools, respectively, compared to untreated controls. When tolerant and sensitive pools were compared, 593 genes were differentially expressed under untreated and 431 genes under waterlogged conditions, of which 122 genes overlapped. To validate the BSR-seq results, the expression levels of six genes were determined by qRT-PCR. The qRT-PCR results were consistent with BSR-seq results. Comparison of allelic polymorphism in mRNA sequences between tolerant and sensitive pools revealed 165 (normal condition) and 128 (waterlogged condition) high-probability SNPs. We found 18 overlapping SNPs with genomic positions mapped. Eighteen SNPs were contained in 18 genes, and eight and nine of 18 genes were responsive to waterlogging stress in tolerant and sensitive lines, respectively. Six alleles of the 18 originated from tolerant pool were significantly up-regulated under waterlogging, but not those from sensitive pool. Importantly, one allele (GRMZM2G055704) of the six genes was mapped between umc1619 and umc1948 on chromosome 1 where a QTL associated with waterlogging tolerance was identified in a previous research, strongly indicating that GRMZM2G055704 is a candidate gene responsive to waterlogging. Our research contributes to the knowledge of the molecular mechanism for waterlogging tolerance in maize.

Distinct Mechanisms of the ORANGE Protein in Controlling Carotenoid Flux

β-Carotene adds nutritious value and determines the color of many fruits, including melon (Cucumis melo). In melon mesocarp, β-carotene accumulation is governed by the Orange gene (CmOr) golden single-nucleotide polymorphism (SNP) through a yet to be discovered mechanism. In Arabidopsis (Arabidopsis thaliana), OR increases carotenoid levels by posttranscriptionally regulating phytoene synthase (PSY). Here, we identified a CmOr nonsense mutation (Cmor-lowβ) that lowered fruit β-carotene levels with impaired chromoplast biogenesis. Cmor-lowβ exerted a minimal effect on PSY transcripts but dramatically decreased PSY protein levels and enzymatic activity, leading to reduced carotenoid metabolic flux and accumulation. However, the golden SNP was discovered to not affect PSY protein levels and carotenoid metabolic flux in melon fruit, as shown by carotenoid and immunoblot analyses of selected melon genotypes and by using chemical pathway inhibitors. The high β-carotene accumulation in golden SNP melons was found to be due to a reduced further metabolism of β-carotene. This was revealed by genetic studies with double mutants including carotenoid isomerase (yofi), a carotenoid-isomerase nonsense mutant, which arrests the turnover of prolycopene. The yofi F2 segregants accumulated prolycopene independently of the golden SNP Moreover, Cmor-lowβ was found to inhibit chromoplast formation and chloroplast disintegration in fruits from 30 d after anthesis until ripening, suggesting that CmOr regulates the chloroplast-to-chromoplast transition. Taken together, our results demonstrate that CmOr is required to achieve PSY protein levels to maintain carotenoid biosynthesis metabolic flux but that the mechanism of the CmOr golden SNP involves an inhibited metabolism downstream of β-carotene to dramatically affect both carotenoid content and plastid fate.

A Combined Comparative Transcriptomic, Metabolomic, and Anatomical Analyses of Two Key Domestication Traits: Pod Dehiscence and Seed Dormancy in Pea (Pisum sp.)

The origin of the agriculture was one of the turning points in human history, and a central part of this was the evolution of new plant forms, domesticated crops. Seed dispersal and germination are two key traits which have been selected to facilitate cultivation and harvesting of crops. The objective of this study was to analyze anatomical structure of seed coat and pod, identify metabolic compounds associated with water-impermeable seed coat and differentially expressed genes involved in pea seed dormancy and pod dehiscence. Comparative anatomical, metabolomics, and transcriptomic analyses were carried out on wild dormant, dehiscent Pisum elatius (JI64, VIR320) and cultivated, indehiscent Pisum sativum non-dormant (JI92, Cameor) and recombinant inbred lines (RILs). Considerable differences were found in texture of testa surface, length of macrosclereids, and seed coat thickness. Histochemical and biochemical analyses indicated genotype related variation in composition and heterogeneity of seed coat cell walls within macrosclereids. Liquid chromatography-electrospray ionization/mass spectrometry and Laser desorption/ionization-mass spectrometry of separated seed coats revealed significantly higher contents of proanthocyanidins (dimer and trimer of gallocatechin), quercetin, and myricetin rhamnosides and hydroxylated fatty acids in dormant compared to non-dormant genotypes. Bulk Segregant Analysis coupled to high throughput RNA sequencing resulted in identification of 770 and 148 differentially expressed genes between dormant and non-dormant seeds or dehiscent and indehiscent pods, respectively. The expression of 14 selected dormancy-related genes was studied by qRT-PCR. Of these, expression pattern of four genes: porin (MACE-S082), peroxisomal membrane PEX14-like protein (MACE-S108), 4-coumarate CoA ligase (MACE-S131), and UDP-glucosyl transferase (MACE-S139) was in agreement in all four genotypes with Massive analysis of cDNA Ends (MACE) data. In case of pod dehiscence, the analysis of two candidate genes (SHATTERING and SHATTERPROOF) and three out of 20 MACE identified genes (MACE-P004, MACE-P013, MACE-P015) showed down-expression in dorsal and ventral pod suture of indehiscent genotypes. Moreover, MACE-P015, the homolog of peptidoglycan-binding domain or proline-rich extensin-like protein mapped correctly to predicted Dpo1 locus on PsLGIII. This integrated analysis of the seed coat in wild and cultivated pea provides new insight as well as raises new questions associated with domestication and seed dormancy and pod dehiscence.