Prunus transcription factors: breeding perspectives

Monitoring Fruit Growth and Development in Apricot (Prunus armeniaca L.) through Gene Expression Analysis

The main objective of this study was to monitor apricot development and ripening through gene expression analysis of key candidate genes using the RT-qPCR technique. Eight apricot cultivars were selected to analyze phenological and genetic patterns from pre-ripening stages through to postharvest. In addition, 19 selected genes were analyzed in the contrasting cultivars 'Cebas Red' and 'Rojo Pasión' in different stages (two preharvest stages S1 and S2, one harvest stage S3, and two postharvest stages S4 and S5). This pool of genes included genes related to fruit growth and ripening, genes associated with fruit color, and genes linked to the fruit's nutraceutical aspects. Among the studied genes, Polygalacturonase (PG), Pectin methylesterase (PME), Aminocyclopropane-1-carboxylate synthase (ACS), and Myo-inositol-1-phosphate synthase (INO1) were directly related to fruit maturation and quality. Significant differential expression was observed between the cultivars, which correlated with variations in firmness, shelf life, and sensory characteristics of the apricots. 'Rojo Pasión' displayed high levels of PG, associated with rapid maturation and shorter postharvest shelf life, whereas 'Cebas Red' exhibited lower levels of this gene, resulting in greater firmness and extended shelf life. Genes CCD4, CRTZ, and ZDS, related to carotenoids, showed varied expression patterns during growth and postharvest stages, with higher levels in 'Rojo Pasión'. On the other hand, Sucrose synthase (SUSY) and Lipoxygenase (LOX2) were prominent during the postharvest and growth stages, respectively. Additionally, GDP-L-galactose phosphorylase (VTC2_5) was linked to better postharvest performance. This research provides valuable insights for future breeding initiatives aimed at enhancing the quality and sustainability of apricot cultivation.

Transcriptomic analysis reveals candidate genes associated with salinity stress tolerance during the early vegetative stage in fababean genotype, Hassawi-2

Abiotic stresses are a significant constraint to plant production globally. Identifying stress-related genes can aid in the development of stress-tolerant elite genotypes and facilitate trait and crop manipulation. The primary aim of this study was to conduct whole transcriptome analyses of the salt-tolerant faba bean genotype, Hassawi-2, under different durations of salt stress (6 h, 12 h, 24 h, 48 h, and 72 h) at the early vegetative stage, to better understand the molecular basis of salt tolerance. After de novo assembly, a total of 140,308 unigenes were obtained. The up-regulated differentially expressed genes (DEGs) were 2380, 2863, 3057, 3484, and 4820 at 6 h, 12 h, 24 h, 48 h, and 72 h of salt stress, respectively. Meanwhile, 1974, 3436, 2371, 3502, and 5958 genes were downregulated at 6 h, 12 h, 24 h, 48 h, and 72 h of salt stress, respectively. These DEGs encoded various regulatory and functional proteins, including kinases, plant hormone proteins, transcriptional factors (TFs) basic helix-loop-helix (bHLH), Myeloblastosis (MYB), and (WRKY), heat shock proteins (HSPs), late embryogenesis abundant (LEA) proteins, dehydrin, antioxidant enzymes, and aquaporin proteins. This suggests that the faba bean genome possesses an abundance of salinity resistance genes, which trigger different adaptive mechanisms under salt stress. Some selected DEGs validated the RNA sequencing results, thus confirming similar gene expression levels. This study represents the first transcriptome analysis of faba bean leaves subjected to salinity stress offering valuable insights into the mechanisms governing salt tolerance in faba bean during the vegetative stage. This comprehensive investigation enhances our understanding of precise gene regulatory mechanisms and holds promise for the development of novel salt-tolerant faba bean salt-tolerant cultivars.

Molecular mechanisms involved in fruit cracking: A review

Several fleshy fruits are highly affected by cracking, a severe physiological disorder that compromises their quality and causes high economical losses to the producers. Cracking can occur due to physiological, genetic or environmental factors and may happen during fruit growth, development and ripening. Moreover, in fleshy fruits, exocarp plays an important role, acting as a mechanical protective barrier, defending against biotic or abiotic factors. Thus, when biochemical properties of the cuticle + epidermis + hypodermis are affected, cracks appear in the fruit skin. The identification of genes involved in development such as cell wall modifications, biosynthesis and transport of cuticular waxes, cuticular membrane deposition and associated transcription factors provides new insights to better understand how fruit cracking is affected by genetic factors. Amongst the major environmental stresses causing cracking are excessive water during fruit development, leading to imbalances in cations such as Ca. This review focus on expression of key genes in these pathways, in their influence in affected fruits and the potential for molecular breeding programs, aiming to develop cultivars more resistant to cracking under adverse environmental conditions.

Prunus Knotted-like Genes: Genome-Wide Analysis, Transcriptional Response to Cytokinin in Micropropagation, and Rootstock Transformation

Knotted1-like homeobox (KNOX) transcription factors are involved in plant development, playing complex roles in aerial organs. As Prunus species include important fruit tree crops of Italy, an exhaustive investigation of KNOX genes was performed using genomic and RNA-seq meta-analyses. Micropropagation is an essential technology for rootstock multiplication; hence, we investigated KNOX transcriptional behavior upon increasing 6-benzylaminopurine (BA) doses and the effects on GF677 propagules. Moreover, gene function in Prunus spp. was assessed by Gisela 6 rootstock transformation using fluorescence and peach KNOX transgenes. Based on ten Prunus spp., KNOX proteins fit into I-II-M classes named after Arabidopsis. Gene number, class member distribution, and chromosome positions were maintained, and exceptions supported the diversification of Prunus from Cerasus subgenera, and that of Armeniaca from the other sections within Prunus. Cytokinin (CK) cis-elements occurred in peach and almond KNOX promoters, suggesting a BA regulatory role in GF677 shoot multiplication as confirmed by KNOX expression variation dependent on dose, time, and interaction. The tripled BA concentration exacerbated stress, altered CK perception genes, and modified KNOX transcriptions, which are proposed to concur in in vitro anomalies. Finally, Gisela 6 transformation efficiency varied (2.6-0.6%) with the genetic construct, with 35S:GFP being more stable than 35S:KNOPE1 lines, which showed leaf modification typical of KNOX overexpression.

Meta-QTL analysis and candidate genes identification for various abiotic stresses in maize (Zea mays L.) and their implications in breeding programs

Global climate change leads to the concurrence of a number of abiotic stresses including moisture stress (drought, waterlogging), temperature stress (heat, cold), and salinity stress, which are the major factors affecting maize production. To develop abiotic stress tolerance in maize, many quantitative trait loci (QTL) have been identified, but very few of them have been utilized successfully in breeding programs. In this context, the meta-QTL analysis of the reported QTL will enable the identification of stable/real QTL which will pave a reliable way to introgress these QTL into elite cultivars through marker-assisted selection. In this study, a total of 542 QTL were summarized from 33 published studies for tolerance to different abiotic stresses in maize to conduct meta-QTL analysis using BiomercatorV4.2.3. Among those, only 244 major QTL with more than 10% phenotypic variance were preferably utilised to carry out meta-QTL analysis. In total, 32 meta-QTL possessing 1907 candidate genes were detected for different abiotic stresses over diverse genetic and environmental backgrounds. The MQTL2.1, 5.1, 5.2, 5.6, 7.1, 9.1, and 9.2 control different stress-related traits for combined abiotic stress tolerance. The candidate genes for important transcription factor families such as ERF, MYB, bZIP, bHLH, NAC, LRR, ZF, MAPK, HSP, peroxidase, and WRKY have been detected for different stress tolerances. The identified meta-QTL are valuable for future climate-resilient maize breeding programs and functional validation of candidate genes studies, which will help to deepen our understanding of the complexity of these abiotic stresses.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01294-9.

Effects of Auxin (Indole-3-butyric Acid) on Adventitious Root Formation in Peach-Based Prunus Rootstocks

Several Prunus species are among the most important cultivated stone fruits in the Mediterranean region, and there is an urgent need to obtain rootstocks with specific adaptations to challenging environmental conditions. The development of adventitious roots (ARs) is an evolutionary mechanism of high relevance for stress tolerance, which has led to the development of environmentally resilient plants. As a first step towards understanding the genetic determinants involved in AR formation in Prunus sp., we evaluated the rooting of hardwood cuttings from five Prunus rootstocks (Adafuel, Adarcias, Cadaman, Garnem, and GF 677) grown in hydroponics. We found that auxin-induced callus and rooting responses were strongly genotype-dependent. To investigate the molecular mechanisms involved in these differential responses, we performed a time-series study of AR formation in two rootstocks with contrasting rooting performance, Garnem and GF 677, by culturing in vitro microcuttings with and without auxin treatment (0.9 mg/L of indole-3-butyric acid [IBA]). Despite showing a similar histological structure, Garnem and GF677 rootstocks displayed dynamic changes in endogenous hormone homeostasis involving metabolites such as indole-3-acetic acid (IAA) conjugated to aspartic acid (IAA-Asp), and these changes could explain the differences observed during rooting.

Genome-Wide Prediction of Transcription Start Sites in Conifers

The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.

Long-read transcriptome sequencing provides insight into lignan biosynthesis during fruit development in Schisandra chinensis

BACKGROUND

Schisandra chinensis, an ancient member of the most basal angiosperm lineage which is known as the ANITA, is a fruit-bearing vine with the pharmacological effects of a multidrug system, such as antioxidant, anti-inflammatory, cardioprotective, neuroprotective, anti-osteoporosis effects. Its major bioactive compound is represented by lignans such as schisandrin. Molecular characterization of lignan biosynthesis in S. chinensis is of great importance for improving the production of this class of active compound. However, the biosynthetic mechanism of schisandrin remains largely unknown.

RESULTS

To understand the potential key catalytic steps and their regulation of schisandrin biosynthesis, we generated genome-wide transcriptome data from three different tissues of S. chinensis cultivar Cheongsoon, including leaf, root, and fruit, via long- and short-read sequencing technologies. A total of 132,856 assembled transcripts were generated with an average length of 1.9 kb and high assembly completeness. Overall, our data presented effective, accurate gene annotation in the prediction of functional pathways. In particular, the annotation revealed the abundance of transcripts related to phenylpropanoid biosynthesis. Remarkably, transcriptome profiling during fruit development of S. chinensis cultivar Cheongsoon revealed that the phenylpropanoid biosynthetic pathway, specific to coniferyl alcohol biosynthesis, showed a tendency to be upregulated at the postfruit development stage. Further the analysis also revealed that the pathway forms a transcriptional network with fruit ripening-related genes, especially the ABA signaling-related pathway. Finally, candidate unigenes homologous to isoeugenol synthase 1 (IGS1) and dirigent-like protein (DIR), which are subsequently activated by phenylpropanoid biosynthesis and thus catalyze key upstream steps in schisandrin biosynthesis, were identified. Their expression was increased at the postfruit development stage, suggesting that they may be involved in the regulation of schisandrin biosynthesis in S. chinensis.

CONCLUSIONS

Our results provide new insights into the production and accumulation of schisandrin in S. chinensis berries and will be utilized as a valuable transcriptomic resource for improving the schisandrin content.

Whole transcriptome analysis resulted in the identification of Chinese sprangletop (Leptochloa chinensis) genes involved in cyhalofop-butyl tolerance

Background

Chinese sprangletop [Leptochloa chinensis (L.) Nees] is an annual malignant weed, which can often be found in paddy fields. Cyhalofop-butyl is a specialized herbicide which is utilized to control L. chinensis. However, in many areas, L. chinensis has become tolerant to this key herbicide due to its continuous long-term use.

Results

In this study, we utilized a tolerant (LC18002) and a sensitive (LC17041) L. chinensis populations previously identified in our laboratory, which were divided into four different groups. We then employed whole transcriptome analysis to identify candidate genes which may be involved in cyhalofop-butyl tolerance. This analysis resulted in the identification of six possible candidate genes, including three cytochrome P450 genes and three ATP-binding cassette transporter genes. We then carried out a phylogenetic analysis to identify homologs of the differentially expressed cytochrome P450 genes. This phylogenetic analysis indicated that all genes have close homologs in other species, some of which have been implicated in non-target site resistance (NTSR).

Conclusions

This study is the first to use whole transcriptome analysis to identify herbicide non-target resistance genes in L. chinensis. The differentially expressed genes represent promising targets for better understanding herbicide tolerance in L. chinensis. The six genes belonging to classes already associated in herbicide tolerance may play important roles in the metabolic resistance of L. chinensis to cyhalofop-butyl, although the exact mechanisms require further study.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07856-z.

Tuning promoter boundaries improves regulatory motif discovery in nonmodel plants: the peach example

The identification of functional elements encoded in plant genomes is necessary to understand gene regulation. Although much attention has been paid to model species like Arabidopsis (Arabidopsis thaliana), little is known about regulatory motifs in other plants. Here, we describe a bottom-up approach for de novo motif discovery using peach (Prunus persica) as an example. These predictions require pre-computed gene clusters grouped by their expression similarity. After optimizing the boundaries of proximal promoter regions, two motif discovery algorithms from RSAT::Plants (http://plants.rsat.eu) were tested (oligo and dyad analysis). Overall, 18 out of 45 co-expressed modules were enriched in motifs typical of well-known transcription factor (TF) families (bHLH, bZip, BZR, CAMTA, DOF, E2FE, AP2-ERF, Myb-like, NAC, TCP, and WRKY) and a few uncharacterized motifs. Our results indicate that small modules and promoter window of [-500 bp, +200 bp] relative to the transcription start site (TSS) maximize the number of motifs found and reduce low-complexity signals in peach. The distribution of discovered regulatory sites was unbalanced, as they accumulated around the TSS. This approach was benchmarked by testing two different expression-based clustering algorithms (network-based and hierarchical) and, as control, genes grouped for harboring ChIPseq peaks of the same Arabidopsis TF. The method was also verified on maize (Zea mays), a species with a large genome. In summary, this article presents a glimpse of the peach regulatory components at genome scale and provides a general protocol that can be applied to other species. A Docker software container is released to facilitate the reproduction of these analyses.

Shotgun proteomics of peach fruit reveals major metabolic pathways associated to ripening

BACKGROUND

Fruit ripening in Prunus persica melting varieties involves several physiological changes that have a direct impact on the fruit organoleptic quality and storage potential. By studying the proteomic differences between the mesocarp of mature and ripe fruit, it would be possible to highlight critical molecular processes involved in the fruit ripening.

RESULTS

To accomplish this goal, the proteome from mature and ripe fruit was assessed from the variety O'Henry through shotgun proteomics using 1D-gel (PAGE-SDS) as fractionation method followed by LC/MS-MS analysis. Data from the 131,435 spectra could be matched to 2740 proteins, using the peach genome reference v1. After data pre-treatment, 1663 proteins could be used for comparison with datasets assessed using transcriptomic approaches and for quantitative protein accumulation analysis. Close to 26% of the genes that code for the proteins assessed displayed higher expression at ripe fruit compared to other fruit developmental stages, based on published transcriptomic data. Differential accumulation analysis between mature and ripe fruit revealed that 15% of the proteins identified were modulated by the ripening process, with glycogen and isocitrate metabolism, and protein localization overrepresented in mature fruit, as well as cell wall modification in ripe fruit. Potential biomarkers for the ripening process, due to their differential accumulation and gene expression pattern, included a pectin methylesterase inhibitor, a gibbellerin 2-beta-dioxygenase, an omega-6 fatty acid desaturase, a homeobox-leucine zipper protein and an ACC oxidase. Transcription factors enriched in NAC and Myb protein domains would target preferentially the genes encoding proteins more abundant in mature and ripe fruit, respectively.

CONCLUSIONS

Shotgun proteomics is an unbiased approach to get deeper into the proteome allowing to detect differences in protein abundance between samples. This technique provided a resolution so that individual gene products could be identified. Many proteins likely involved in cell wall and sugar metabolism, aroma and color, change their abundance during the transition from mature to ripe fruit.

Molecular Bases of Fruit Quality in Prunus Species: An Integrated Genomic, Transcriptomic, and Metabolic Review with a Breeding Perspective

In plants, fruit ripening is a coordinated developmental process that requires the change in expression of hundreds to thousands of genes to modify many biochemical and physiological signal cascades such as carbohydrate and organic acid metabolism, cell wall restructuring, ethylene production, stress response, and organoleptic compound formation. In Prunus species (including peaches, apricots, plums, and cherries), fruit ripening leads to the breakdown of complex carbohydrates into sugars, fruit firmness reductions (softening by cell wall degradation and cuticle properties alteration), color changes (loss of green color by chlorophylls degradation and increase in non-photosynthetic pigments like anthocyanins and carotenoids), acidity decreases, and aroma increases (the production and release of organic volatile compounds). Actually, the level of information of molecular events at the transcriptional, biochemical, hormonal, and metabolite levels underlying ripening in Prunus fruits has increased considerably. However, we still poorly understand the molecular switch that occurs during the transition from unripe to ripe fruits. The objective of this review was to analyze of the molecular bases of fruit quality in Prunus species through an integrated metabolic, genomic, transcriptomic, and epigenetic approach to better understand the molecular switch involved in the ripening process with important consequences from a breeding point of view.

Comparative transcriptomic and physiological analyses of contrasting hybrid cultivars ND476 and ZX978 identify important differentially expressed genes and pathways regulating drought stress tolerance in maize

BACKGROUND

Drought is the major abiotic stress factor that negatively influences growth and yield in cereal grain crops such as maize (Zea mays L.). A multitude of genes and pathways tightly modulate plant growth, development and responses to environmental stresses including drought. Therefore, crop breeding efforts for enhanced drought resistance require improved knowledge of plant drought responses.

OBJECTIVE

Here, we sought to elucidate the molecular and physiological mechanisms underpinning maize drought stress tolerance.

METHODS

We therefore applied a 12-day water-deficit stress treatment to maize plants of two contrasting (drought tolerant ND476 and drought sensitive ZX978) hybrid cultivars at the late vegetative (V12) growth stage and performed a large-scale RNA sequencing (RNA-seq) transcriptome analysis of the leaf tissues.

RESULTS

A comparative analysis of the two genotypes leaf transcriptomes and physiological parameters revealed the key differentially expressed genes (DEGs) and metabolic pathways that respond to drought in a genotype-specific manner. A total of 3114 DEGs were identified, with 21 DEGs being specifically expressed in tolerant genotype ND476 in response to drought stress. Of these, genes involved in secondary metabolites biosynthesis, transcription factor regulation, detoxification and stress defense were highly expressed in ND476. Physiological analysis results substantiated our RNA-seq data, with ND476 exhibiting better cell water retention, higher soluble protein content and guaiacol peroxidase activity, along with low lipid peroxidation extent than the sensitive cultivar ZX978 under drought conditions.

CONCLUSION

Our findings enrich the maize genetic resources and enhance our further understanding of the molecular mechanisms regulating drought stress tolerance in maize. Additionally, the DEGs screened in this study may provide a foundational basis for our future targeted cloning studies.

Comparative transcriptomic analysis reveals novel roles of transcription factors and hormones during the flowering induction and floral bud differentiation in sweet cherry trees (Prunus avium L. cv. Bing)

In sweet cherry trees, flowering is commercially important because the flowers, after fertilization, will generate the fruits. In P. avium, the flowering induction and flower organogensis are the first developmental steps towards flower formation and they occur within specialized organs known as floral buds during the summer, nine months before blooming. During this period the number of floral buds per tree and the bud fruitfulness (number of flowers per bud) are stablished affecting the potential yield of orchards and the plant architecture. The floral bud development is sensitive to any type of stress and the hotter and drier summers will interfere with this process and are calling for new adapted cultivars. A better understanding of the underlying molecular and hormonal mechanisms would be of help, but unlike the model plant Arabidopsis, very little is known about floral induction in sweet cherry. To explore the molecular mechanism of floral bud differentiation, high-throughput RNA sequencing was used to detect differences in the gene expression of P. avium floral buds at five differentiation stages. We found 2,982 differentially expressed genes during floral bud development. We identified genes associated with floral initiation or floral organ identity that appear to be useful biomarkers of floral development and several transcription factor families (ERF, MYB, bHLH, MADS-box and NAC gene family) with novel potential roles during floral transition in this species. We analyzed in deep the MADS-box gene family and we shed light about their key role during floral bud and organs development in P. avium. Furthermore, the hormonal-related signatures in the gene regulatory networks and the dynamic changes of absicic acid, zeatin and indolacetic acid contents in buds suggest an important role for these hormones during floral bud differentiation in sweet cherry. These data provide a rich source of novel informacion for functional and evolutionary studies about floral bud development in sweet cherry and new tools for biotechnology and breeding.

A Prunus persica genome-wide RNA-seq approach uncovers major differences in the transcriptome among chilling injury sensitive and non-sensitive varieties

Chilling injury represents a major constrain for crops productivity. Prunus persica, one of the most relevant rosacea crops, have early season varieties that are resistant to chilling injury, in contrast to late season varieties, which display chilling symptoms such as mealiness (dry, sandy fruit mesocarp) after prolonged storage at chilling temperatures. To uncover the molecular processes related to the ability of early varieties to withstand mealiness, postharvest and genome-wide RNA-seq assessments were performed in two early and two late varieties. Differences in juice content and ethylene biosynthesis were detected among early and late season fruits that became mealy after exposed to prolonged chilling. Principal component and data distribution analysis revealed that cold-stored late variety fruit displayed an exacerbated and unique transcriptome profile when compared to any other postharvest condition. A differential expression analysis performed using an empirical Bayes mixture modeling approach followed by co-expression and functional enrichment analysis uncover processes related to ethylene, lipids, cell wall, carotenoids and DNA metabolism, light response, and plastid homeostasis associated to the susceptibility or resistance of P. persica varieties to chilling stress. Several of the genes related to these processes are in quantitative trait loci (QTL) associated to mealiness in P. persica. Together, these analyses exemplify how P. persica can be used as a model for studying chilling stress in plants.

Key Maize Drought-Responsive Genes and Pathways Revealed by Comparative Transcriptome and Physiological Analyses of Contrasting Inbred Lines

To unravel the molecular mechanisms underpinning maize (Zea mays L.) drought stress tolerance, we conducted comprehensive comparative transcriptome and physiological analyses of drought-tolerant YE8112 and drought-sensitive MO17 inbred line seedlings that had been exposed to drought treatment for seven days. Resultantly, YE8112 seedlings maintained comparatively higher leaf relative water and proline contents, greatly increased peroxidase activity, but decreased malondialdehyde content, than MO17 seedlings. Using an RNA sequencing (RNA-seq)-based approach, we identified a total of 10,612 differentially expressed genes (DEGs). From these, we mined out four critical sets of drought responsive DEGs, including 80 specific to YE8112, 5140 shared between the two lines after drought treatment (SD_TD), five DEGs of YE8112 also regulated in SD_TD, and four overlapping DEGs between the two lines. Drought-stressed YE8112 DEGs were primarily associated with nitrogen metabolism and amino-acid biosynthesis pathways, whereas MO17 DEGs were enriched in the ribosome pathway. Additionally, our physiological analyses results were consistent with the predicted RNA-seq-based findings. Furthermore, quantitative real-time polymerase chain reaction (qRT-PCR) analysis and the RNA-seq results of twenty representative DEGs were highly correlated (R² = 98.86%). Crucially, tolerant line YE8112 drought-responsive genes were predominantly implicated in stress signal transduction; cellular redox homeostasis maintenance; MYB, NAC, WRKY, and PLATZ transcriptional factor modulated; carbohydrate synthesis and cell-wall remodeling; amino acid biosynthesis; and protein ubiquitination processes. Our findings offer insights into the molecular networks mediating maize drought stress tolerance.

Characterization of drought stress-responsive root transcriptome of faba bean (Vicia faba L.) using RNA sequencing

Drought and salinity are the major factors that limit the faba bean (Vicia faba L.) production worldwide. The aim of this study is to identify the water stress differentially expressed genes (DEGs) through the root transcriptome analyses of the drought-tolerant Hassawi 2 genotype at vegetative and flowering stages. A total of 624.8 M high-quality Illumina reads were generated and assembled into 198,155 all-unigenes with a mean length of 738 bp and an N50 length of 1347 bp. Among all-unigenes, 78,262 were assigned to non-redundant (Nr), 66,254 to nucleotide (Nt), 54,034 to KEGG, and 43,913 to gene ontology (GO) annotations. A total of 36,834 and 35,510 unigenes were differentially expressed at the vegetative and flowering stages of Hassawi 2 under drought stress, respectively. The majority of unigenes were down-regulated at both developmental stages. However, the number of genes up-regulated (15,366) at the flowering stage exceeded the number of those up-regulated (14,097) at the vegetative stage, and the number of genes down-regulated (20,144) at the flowering stage was smaller than the number of those down-regulated (22,737) at the vegetative stage. The drought stress-responsive differentially expressed unigenes coded for various regulatory proteins, including protein kinases and phosphatases, transcription factors and plant hormones and functional proteins including enzymes for osmoprotectant, detoxification and transporters were differentially expressed, most of which were largely up-regulated. Moreover, a substantial proportion of the DEGs identified in this study were novel, most exhibited a significant change in their expression levels under water stress, making them an unexploited resource that might control specific responses to drought stress in the faba bean. Finally, qRT-PCR results were found almost consistent with the results of next-generation sequencing. Our data will help in understanding the drought tolerance mechanisms in plants and will provide resources for functional genomics.

Monitoring Dormancy Transition in Almond [Prunus Dulcis (Miller) Webb] during Cold and Warm Mediterranean Seasons through the Analysis of a DAM (Dormancy-Associated MADS-Box) Gene

For fruit tree (Prunus) species, flower bud dormancy completion determines the quality of bud break and the flowering time. In the present climate change and global warming context, the relationship between dormancy and flowering processes is a fundamental goal in molecular biology of these species. In almond [P. dulcis (Miller) Webb], flowering time is a trait of great interest in the development of new cultivars adapted to different climatic areas. Late flowering is related to a long dormancy period due to high chilling requirements of the cultivar. It is considered a quantitative and highly heritable character but a dominant gene (Late bloom, Lb) was also described. A major QTL (quantitative trait loci) in the linkage group (LG) 4 was associated with Lb, together with other three QTLs in LG1 and LG7. In addition, DAM (Dormancy-Associated MADS-Box) genes located in LG1 have been largely described as a gene family involved in bud dormancy in different Prunus species including peach [P. persica (L.) Batsch] and Japanese apricot (P. mume Sieb. et Zucc.). In this work, a DAM transcript was cloned and its expression was analysed by qPCR (quantitative Polymerase Chain Reaction) in almond flower buds during the dormancy release. For this purpose two almond cultivars (‘Desmayo Largueta’ and ‘Penta’) with different chilling requirements and flowering time were used, and the study was performed along two years. The complete coding sequence, designated PdDAM6 (Prunus dulcis DAM6), was subjected to a phylogenetic analysis with homologous sequences from other Prunus species. Finally, expression dynamics analysed by using qPCR showed a continuous decrease in transcript levels for both cultivars and years during the period analysed. Monitoring almond flower bud dormancy through DAM expression should be used to improve almond production in different climate conditions.

Transcriptomic signature of drought response in pearl millet (Pennisetum glaucum (L.) and development of web-genomic resources

Pearl millet, (Pennisetum glaucum L.), an efficient (C4) crop of arid/semi-arid regions is known for hardiness. Crop is valuable for bio-fortification combating malnutrition and diabetes, higher caloric value and wider climatic resilience. Limited studies are done in pot-based experiments for drought response at gene-expression level, but field-based experiment mimicking drought by withdrawal of irrigation is still warranted. We report de novo assembly-based transcriptomic signature of drought response induced by irrigation withdrawal in pearl millet. We found 19983 differentially expressed genes, 7595 transcription factors, gene regulatory network having 45 hub genes controlling drought response. We report 34652 putative markers (4192 simple sequence repeats, 12111 SNPs and 6249 InDels). Study reveals role of purine and tryptophan metabolism in ABA accumulation mediating abiotic response in which MAPK acts as major intracellular signal sensing drought. Results were validated by qPCR of 13 randomly selected genes. We report the first web-based genomic resource ( http://webtom.cabgrid.res.in/pmdtdb/ ) which can be used for candidate genes-based SNP discovery programs and trait-based association studies. Looking at climatic change, nutritional and pharmaceutical importance of this crop, present investigation has immense value in understanding drought response in field condition. This is important in germplasm management and improvement in endeavour of pearl millet productivity.

Discovery of Putative Herbicide Resistance Genes and Its Regulatory Network in Chickpea Using Transcriptome Sequencing

Background: Chickpea (Cicer arietinum L.) contributes 75% of total pulse production. Being cheaper than animal protein, makes it important in dietary requirement of developing countries. Weed not only competes with chickpea resulting into drastic yield reduction but also creates problem of harboring fungi, bacterial diseases and insect pests. Chemical approach having new herbicide discovery has constraint of limited lead molecule options, statutory regulations and environmental clearance. Through genetic approach, transgenic herbicide tolerant crop has given successful result but led to serious concern over ecological safety thus non-transgenic approach like marker assisted selection is desirable. Since large variability in tolerance limit of herbicide already exists in chickpea varieties, thus the genes offering herbicide tolerance can be introgressed in variety improvement programme. Transcriptome studies can discover such associated key genes with herbicide tolerance in chickpea. Results: This is first transcriptomic studies of chickpea or even any legume crop using two herbicide susceptible and tolerant genotypes exposed to imidazoline (Imazethapyr). Approximately 90 million paired-end reads generated from four samples were processed and assembled into 30,803 contigs using reference based assembly. We report 6,310 differentially expressed genes (DEGs), of which 3,037 were regulated by 980 miRNAs, 1,528 transcription factors associated with 897 DEGs, 47 Hub proteins, 3,540 putative Simple Sequence Repeat-Functional Domain Marker (SSR-FDM), 13,778 genic Single Nucleotide Polymorphism (SNP) putative markers and 1,174 Indels. Randomly selected 20 DEGs were validated using qPCR. Pathway analysis suggested that xenobiotic degradation related gene, glutathione S-transferase (GST) were only up-regulated in presence of herbicide. Down-regulation of DNA replication genes and up-regulation of abscisic acid pathway genes were observed. Study further reveals the role of cytochrome P450, xyloglucan endotransglucosylase/hydrolase, glutamate dehydrogenase, methyl crotonoyl carboxylase and of thaumatin-like genes in herbicide resistance. Conclusion: Reported DEGs can be used as genomic resource for future discovery of candidate genes associated with herbicide tolerance. Reported markers can be used for future association studies in order to develop marker assisted selection (MAS) for refinement. In endeavor of chickpea variety development programme, these findings can be of immense use in improving productivity of chickpea germplasm.

Transcriptional Responses in Root and Leaf of Prunus persica under Drought Stress Using RNA Sequencing

Prunus persica L. Batsch, or peach, is one of the most important crops and it is widely established in irrigated arid and semi-arid regions. However, due to variations in the climate and the increased aridity, drought has become a major constraint, causing crop losses worldwide. The use of drought-tolerant rootstocks in modern fruit production appears to be a useful method of alleviating water deficit problems. However, the transcriptomic variation and the major molecular mechanisms that underlie the adaptation of drought-tolerant rootstocks to water shortage remain unclear. Hence, in this study, high-throughput sequencing (RNA-seq) was performed to assess the transcriptomic changes and the key genes involved in the response to drought in root tissues (GF677 rootstock) and leaf tissues (graft, var. Catherina) subjected to 16 days of drought stress. In total, 12 RNA libraries were constructed and sequenced. This generated a total of 315 M raw reads from both tissues, which allowed the assembly of 22,079 and 17,854 genes associated with the root and leaf tissues, respectively. Subsets of 500 differentially expressed genes (DEGs) in roots and 236 in leaves were identified and functionally annotated with 56 gene ontology (GO) terms and 99 metabolic pathways, which were mostly associated with aminobenzoate degradation and phenylpropanoid biosynthesis. The GO analysis highlighted the biological functions that were exclusive to the root tissue, such as "locomotion," "hormone metabolic process," and "detection of stimulus," indicating the stress-buffering role of the GF677 rootstock. Furthermore, the complex regulatory network involved in the drought response was revealed, involving proteins that are associated with signaling transduction, transcription and hormone regulation, redox homeostasis, and frontline barriers. We identified two poorly characterized genes in P. persica: growth-regulating factor 5 (GRF5), which may be involved in cellular expansion, and AtHB12, which may be involved in root elongation. The reliability of the RNA-seq experiment was validated by analyzing the expression patterns of 34 DEGs potentially involved in drought tolerance using quantitative reverse transcription polymerase chain reaction. The transcriptomic resources generated in this study provide a broad characterization of the acclimation of P. persica to drought, shedding light on the major molecular responses to the most important environmental stressor.