Analysis of the Citrullus colocynthis transcriptome during water deficit stress

Chromosome-level genome assembly and functional annotation of Citrullus colocynthis: unlocking genetic resources for drought-resilient crop development

MAIN CONCLUSION

The chromosome-level genome assembly of Citrullus colocynthis reveals its genetic potential for enhancing drought tolerance, paving the way for innovative crop improvement strategies. This study presents the first comprehensive genome assembly and annotation of Citrullus colocynthis, a drought-tolerant wild close relative of cultivated watermelon, highlighting its potential for enhancing agricultural resilience to climate change. The study achieved a chromosome-level assembly using advanced sequencing technologies, including PacBio HiFi and Hi-C, revealing a genome size of approximately 366 Mb with low heterozygosity and substantial repetitive content. Our analysis identified 23,327 gene models, that could encode stress response mechanisms for species' adaptation to arid environments. Comparative genomics with closely related species illuminated the evolutionary dynamics within the Cucurbitaceae family. In addition, resequencing of 27 accessions from the United Arab Emirates (UAE) identified genetic diversity, suggesting a foundation for future breeding programs. This genomic resource opens new avenues for the de novo domestication of C. colocynthis, offering a blueprint for developing crops with enhanced drought tolerance, disease resistance, and nutritional profiles, crucial for sustaining future food security in the face of escalating climate challenges.

Metatranscriptome analysis of symptomatic bitter apple plants revealed mixed viral infections with a putative novel polerovirus

BACKGROUND

Next-generation Sequencing (NGS) combined with bioinformatic analyses constitutes a powerful approach for identifying and characterizing previously unknown viral genomes. In this study, leaf samples from bitter apple plants (Citrullus colocynthis (L.) Schrad) exhibiting symptoms such as dwarfing, leaf crinkling, and chlorosis were collected from the southern part of Kerman province, Iran.

RESULTS

Putative infecting viruses were identified through de novo assembly of sequencing reads using various tools, followed by BLAST analysis. Complete genomes for Squash vein yellowing virus (SqVYV), Citrus-associated rhabdovirus (CiaRV), and a novel polerovirus-related strain termed Bitter apple aphid-borne yellows virus (BaABYV) were assembled and characterized. Additionally, a partial genome for Watermelon mosaic virus (WMV) was assembled. The genomic organization of the BaABYV was determined to be 5'-ORF0-ORF1-ORF1,2-ORF3a-ORF3-ORF3,5-ORF4-3'. Amino acid sequence identities for inferred proteins (P0 and P1, P1,2) with known poleroviruses were found to be the 90% species delineation limit, implying that BaABYV should be considered a new member of the genus Polerovirus. Recombination events were observed in the BaABYV and WMV strains; such events were not found in the CiaRV strain.

CONCLUSIONS

Molecular evidence from this study suggests that C. colocynthis is a reservoir host of several plant viruses. Among them, BaABYV is proposed as a new member of the genus Polerovirus. Furthermore, the CiaRV strain has been reported for the first time from Iran.

Integrated Bulk Segregant Analysis, Fine Mapping, and Transcriptome Revealed QTLs and Candidate Genes Associated with Drought Adaptation in Wild Watermelon

Drought stress has detrimental effects on crop productivity worldwide. A strong root system is crucial for maintaining water and nutrients uptake under drought stress. Wild watermelons possess resilient roots with excellent drought adaptability. However, the genetic factors controlling this trait remain uninvestigated. In this study, we conducted a bulk segregant analysis (BSA) on an F2 population consisting of two watermelon genotypes, wild and domesticated, which differ in their lateral root development under drought conditions. We identified two quantitative trait loci (qNLR_Dr. Chr01 and qNLR_Dr. Chr02) associated with the lateral root response to drought. Furthermore, we determined that a small region (0.93 Mb in qNLR_Dr. Chr01) is closely linked to drought adaptation through quantitative trait loci (QTL) validation and fine mapping. Transcriptome analysis of the parent roots under drought stress revealed unique effects on numerous genes in the sensitive genotype but not in the tolerant genotype. By integrating BSA, fine mapping, and the transcriptome, we identified six genes, namely L-Ascorbate Oxidase (AO), Cellulose Synthase-Interactive Protein 1 (CSI1), Late Embryogenesis Abundant Protein (LEA), Zinc-Finger Homeodomain Protein 2 (ZHD2), Pericycle Factor Type-A 5 (PFA5), and bZIP transcription factor 53-like (bZIP53-like), that might be involved in the drought adaptation. Our findings provide valuable QTLs and genes for marker-assisted selection in improving water-use efficiency and drought tolerance in watermelon. They also lay the groundwork for the genetic manipulation of drought-adapting genes in watermelon and other Cucurbitacea species.

Plant hormone signals regulate trehalose accumulation against osmotic stress in watermelon cells

Trehalose, one of the most chemically stable sugars, can effectively improve the tolerance of various plants against abiotic stress by protecting and stabilizing protein and cell membranes. However, the signaling pathway in trehalose biosynthesis triggered by abiotic stresses is still unclear. In the study, it can be shown that exogenous trehalose can alleviate the inhibitory effect of osmotic stress on cell growth, suppress extracellular alkalization, ROS burst, and maintain the integrity of the microtubular cytoskeleton. Trehalose-6-phosphate synthase (TPS) is the key limiting enzyme for trehalose synthesis and is encoded by 7 ClTPS genes, located in 7 different chromosomes of the watermelon genome. Expression analysis by qRT-PCR indicated that osmotic stress could upregulate the expression of all the family members of ClTPS and promote the accumulation of trehalose in watermelon cells accordingly. Exogenous methyl jasmonate (MeJA), ethephon (ETH), abscisic acid (ABA), or salicylic acid (SA) induced trehalose accumulation, with MeJA being the most effective treatment. When fluridone (FL), an ABA biosynthesis inhibitor, was pre-perfused into the cells before osmotic stress, trehalose accumulation and packed cell volume were suppressed significantly, whereas inhibition of ethylene biosynthesis could even restore cell growth. Moreover, inhibition of trehalose hydrolysis could also increase the tolerance against osmotic stress. This study shows that trehalose biosynthesis is phytohormone-dependent and the hydrolysis of trehalose is involved in osmotic tolerance regulation.

MicroRNAs expression profiles in early responses to different levels of water deficit in Setaria viridis

Water deficit is a major constraint for crops of economic importance in almost all agricultural regions. However, plants have an active defense system to adapt to these adverse conditions, acting in the reprogramming of gene expression responsible for encoding microRNAs (miRNAs). These miRNAs promote the regulation to the target gene expression by the post-transcriptional (PTGS) and transcriptional gene silencing (TGS), modulating several pathways including defense response to water deficit. The broader knowledge of the miRNA expression profile and its regulatory networks in response to water deficit can provide evidence for the development of new biotechnological tools for genetic improvement of several important crops. In this study, we used Setaria viridis accession A10.1 as a C4 model plant to widely investigate the miRNA expression profile in early responses to different levels of water deficit. Ecophysiological studies in Setaria viridis under water deficit and after rewatering demonstrated a drought tolerant accession, capable of a rapid recovery from the stress. Deep small RNA sequencing and degradome studies were performed in plants submitted to drought to identify differentially expressed miRNA genes and their predicted targets, using in silico analysis. Our findings showed that several miRNAs were differentially modulated in response to distinctive levels of water deficit and after rewatering. The predicted mRNA targets mainly corresponded to genes related to cell wall remodeling, antioxidant system and drought-related transcription factors, indicating that these genes are rapidly regulated in early responses to drought stress. The implications of these modulations are extensively discussed, and higher-effect miRNAs are suggested as major players for potential use in genetic engineering to improve drought tolerance in economically important crops, such as sugarcane, maize, and sorghum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01226-z.

Pangenomes as a Resource to Accelerate Breeding of Under-Utilised Crop Species

Pangenomes are a rich resource to examine the genomic variation observed within a species or genera, supporting population genetics studies, with applications for the improvement of crop traits. Major crop species such as maize (Zea mays), rice (Oryza sativa), Brassica (Brassica spp.), and soybean (Glycine max) have had pangenomes constructed and released, and this has led to the discovery of valuable genes associated with disease resistance and yield components. However, pangenome data are not available for many less prominent crop species that are currently under-utilised. Despite many under-utilised species being important food sources in regional populations, the scarcity of genomic data for these species hinders their improvement. Here, we assess several under-utilised crops and review the pangenome approaches that could be used to build resources for their improvement. Many of these under-utilised crops are cultivated in arid or semi-arid environments, suggesting that novel genes related to drought tolerance may be identified and used for introgression into related major crop species. In addition, we discuss how previously collected data could be used to enrich pangenome functional analysis in genome-wide association studies (GWAS) based on studies in major crops. Considering the technological advances in genome sequencing, pangenome references for under-utilised species are becoming more obtainable, offering the opportunity to identify novel genes related to agro-morphological traits in these species.

Integrated Transcriptome Analysis and Single-Base Resolution Methylomes of Watermelon (Citrullus lanatus) Reveal Epigenome Modifications in Response to Osmotic Stress

DNA methylation plays an important role against adverse environment by reshaping transcriptional profile in plants. To better understand the molecular mechanisms of watermelon response to osmotic stress, the suspension cultured watermelon cells were treated with 100mM mannitol, and then a methylated cytosines map was generated by whole genome bisulfite sequencing (WGBS). Combined with transcriptome sequencing, the effects of osmotic stress on differentially methylated expressed genes (DMEGs) were assessed. It was found that genes related to plant hormone synthesis, signal transduction, osmoregulatory substance-related and reactive oxygen species scavenging-related enzyme could rapidly respond to osmotic stress. The overall methylation level of watermelon decreased after osmotic stress treatment, and demethylation occurred in CG, CHG, and CHH contexts. Moreover, differentially methylated expressed genes (DMEGs) were significantly enriched in RNA transport, starch and sucrose metabolism, plant hormone signal transduction and biosynthesis of secondary metabolites, especially in biosynthesis of osmolytes synthase genes. Interestingly, demethylation of a key enzyme gene Cla014489 in biosynthesis of inositol upregulated its expression and promoted accumulation of inositol, which could alleviate the inhibition of cell growth caused by osmotic stress. Meanwhile, a recombinant plasmid pET28a-Cla014489 was constructed and transferred into Escherichia coli BL21 for prokaryotic expression and the expression of ClMIPS protein could improve the tolerance of E. coli to osmotic stress. The effect of methylation level on the expression properties of inositol and its related genes was further confirmed by application of DNA methylation inhibitor 5-azacytidine. These results provide a preliminary insight into the altered methylation levels of watermelon cells in response to osmotic stress and suggest a new mechanism that how watermelon cells adapt to osmotic stress.

DNA methylation and histone modifications induced by abiotic stressors in plants

BACKGROUND

A review of research shows that methylation in plants is more complex and sophisticated than in microorganisms and animals. Overall, studies on the effects of abiotic stress on epigenetic modifications in plants are still scarce and limited to few species. Epigenetic regulation of plant responses to environmental stresses has not been elucidated. This study summarizes key effects of abiotic stressors on DNA methylation and histone modifications in plants.

DISCUSSION

Plant DNA methylation and histone modifications in responses to abiotic stressors varied and depended on the type and level of stress, plant tissues, age, and species. A critical analysis of the literature available revealed that 44% of the epigenetic modifications induced by abiotic stressors in plants involved DNA hypomethylation, 40% DNA hypermethylation, and 16% histone modification. The epigenetic changes in plants might be underestimated since most authors used methods such as methylation-sensitive amplification polymorphism (MSAP), High performance liquid chromatography (HPLC), and immunolabeling that are less sensitive compared to bisulfite sequencing and single-base resolution methylome analyses. More over, mechanisms underlying epigenetic changes in plants have not yet been determined since most reports showed only the level or/and distribution of DNA methylation and histone modifications.

CONCLUSIONS

Various epigenetic mechanisms are involved in response to abiotic stressors, and several of them are still unknown. Integrated analysis of the changes in the genome by omic approaches should help to identify novel components underlying mechanisms involved in DNA methylation and histone modifications associated with plant response to environmental stressors.

Transcriptomic analysis reveals key transcription factors associated to drought tolerance in a wild papaya (Carica papaya) genotype

Most of the commercial papaya genotypes show susceptibility to water deficit stress and require high volumes of irrigation water to yield properly. To tackle this problem, we have collected wild native genotypes of Carica papaya that have proved to show better physiological performance under water deficit stress than the commercial cultivar grown in Mexico. In the present study, plants from a wild Carica papaya genotype and a commercial genotype were subjected to water deficit stress (WDS), and their response was characterized in physiological and molecular terms. The physiological parameters measured (water potential, photosynthesis, Fv/Fm and electrolyte leakage) confirmed that the papaya wild genotype showed better physiological responses than the commercial one when exposed to WDS. Subsequently, RNA-Seq was performed for 4 cDNA libraries in both genotypes (susceptible and tolerant) under well-watered conditions, and when they were subjected to WDS for 14 days. Consistently, differential expression analysis revealed that after 14 days of WDS, the wild tolerant genotype had a higher number of up-regulated genes, and a higher number of transcription factors (TF) that were differentially expressed in response to WDS, than the commercial genotype. Thus, six TF genes (CpHSF, CpMYB, CpNAC, CpNFY-A, CpERF and CpWRKY) were selected for further qRT-PCR analysis as they were highly expressed in response to WDS in the wild papaya genotype. qRT-PCR results confirmed that the wild genotype had higher expression levels (REL) in all 6 TF genes than the commercial genotype. Our transcriptomic analysis should help to unravel candidate genes that may be useful in the development of new drought-tolerant cultivars of this important tropical crop.

Functional Relevance of Citrulline in the Vegetative Tissues of Watermelon During Abiotic Stresses

A non-protein amino acid, citrulline, is a compatible solute involved in the maintenance of cellular osmolarity during abiotic stresses. Despite its significance, a coherent model indicating the role of citrulline during stress conditions has not yet emerged. We have used watermelon, naturally rich in citrulline, as a model to understand its accumulation during drought stress and nitrogen perturbation using transcriptomic and metabolomic analysis. Experiments were performed in the semi-controlled environment, and open field to study the accumulation of drought-induced citrulline in the vegetative tissues of watermelon by monitoring the stress treatments using physiological measurements. The amino acid profiling of leaves and stems in response to drought stress showed up to a 38 and 16-fold increase in citrulline content, respectively. Correlation between amino acids indicated a concomitant activation of a metabolic pathway that included citrulline, its precursor (ornithine), and catabolic product (arginine). Consistent with its accumulation, the gene expression analysis and RNA-Sequencing confirmed activation of citrulline biosynthesis-related genes - Ornithine carbamoyl-transferase (OTC), N-acetylornithine deacetylase (AOD) and Carbamoyl phosphate synthases (CPS), and down-regulation of catabolic genes; Arginosuccinate lyase (ASL) and Arginosuccinate synthases (ASS) in drought-stressed leaf tissues. Based on the relative abundance in the nitrogen-depleted vegetative tissues and down-regulation of genes involved in citrulline biosynthesis, we also demonstrated that the nitrogen status of the plant regulates citrulline. Taken together, these data provide further insights into the metabolic and molecular mechanisms underlying the amino acid metabolism under environmental stress and the significance of non-protein amino acid citrulline in plants.

Citrulline and Arginine Content of Taxa of Cucurbitaceae

Watermelon is the most significant, natural plant source of L-citrulline, a non-proteinaceous amino acid that benefits cardiovascular health and increases vasodilation in many tissues of the body. Watermelon is a member of the Cucurbitaceae, which includes squash, melon, pumpkin, and cucumber. It is possible that other cucurbits could be good sources of citrulline or of arginine, its direct precursor. Twenty-one cultigens were evaluated in triplicate at two locations in North Carolina to estimate citrulline and arginine amounts and variation due to cultigen, replication, and environment. Cultigens containing the highest amount of citrulline (based on LS means) in g/kg fresh weight were ’Crimson Sweet’ watermelon (2.85), ’Dixielee’ watermelon (2.43), casaba-type melon (0.86), mouse melon (0.64), and horned melon rind (0.45). Additionally, mouse melon, horned melon, and bitter gourd (arils) may be interesting sources of arginine-family amino acids, perhaps because of their large seed and aril content relative to mesocarp.

Transcriptomic and physiological analyses reveal drought adaptation strategies in drought-tolerant and -susceptible watermelon genotypes

Drought stress has become one of the most urgent environmental hazards for horticultural crops. In this research, we analyzed watermelon adaptation strategies to drought stress in drought-tolerant (M20) and -susceptible (Y34) genotypes via transcriptomic and physiological analyses. After drought stress, a total of 6228 and 4311 differentially expressed genes (DEGs) were observed in Y34 and M20, respectively. Numerous DEGs were involved in various defense responses such as antioxidation, protein protection, osmotic adjustment, wax accumulation, hormone signaling, and melatonin biosynthesis. Accordingly, the contents of ABA, melatonin, wax, and some osmoprotectants were increased by drought stress in both Y34 and especially M20. Exogenous application of melatonin or ABA induced wax accumulation and drought tolerance and melatonin may function upstream of ABA. In comparison to Y34, M20 was more able to activate ABA signaling, melatonin biosynthesis, osmotic adjustment, antioxidation, and wax accumulation under drought stress. These stronger responses confer M20 tolerance to drought. Photosynthesis and most DEGs involved in photosynthesis and cell growth were decreased by drought stress in both M20 and especially Y34. For drought-susceptible genotypes, growth retardation may be an important mechanism for saving and redistributing resources in order to reprogram stress signaling networks.

Integrated mRNA and microRNA analysis identifies genes and small miRNA molecules associated with transcriptional and post-transcriptional-level responses to both drought stress and re-watering treatment in tobacco

BACKGROUND

Drought stress is one of the most severe problem limited agricultural productivity worldwide. It has been reported that plants response to drought-stress by sophisticated mechanisms at both transcriptional and post-transcriptional levels. However, the precise molecular mechanisms governing the responses of tobacco leaves to drought stress and water status are not well understood. To identify genes and miRNAs involved in drought-stress responses in tobacco, we performed both mRNA and small RNA sequencing on tobacco leaf samples from the following three treatments: untreated-control (CL), drought stress (DL), and re-watering (WL).

RESULTS

In total, we identified 798 differentially expressed genes (DEGs) between the DL and CL (DL vs. CL) treatments and identified 571 DEGs between the WL and DL (WL vs. DL) treatments. Further analysis revealed 443 overlapping DEGs between the DL vs. CL and WL vs. DL comparisons, and, strikingly, all of these genes exhibited opposing expression trends between these two comparisons, strongly suggesting that these overlapping DEGs are somehow involved in the responses of tobacco leaves to drought stress. Functional annotation analysis showed significant up-regulation of genes annotated to be involved in responses to stimulus and stress, (e.g., late embryogenesis abundant proteins and heat-shock proteins) antioxidant defense (e.g., peroxidases and glutathione S-transferases), down regulation of genes related to the cell cycle pathway, and photosynthesis processes. We also found 69 and 56 transcription factors (TFs) among the DEGs in, respectively, the DL vs. CL and the WL vs. DL comparisons. In addition, small RNA sequencing revealed 63 known microRNAs (miRNA) from 32 families and 368 novel miRNA candidates in tobacco. We also found that five known miRNA families (miR398, miR390, miR162, miR166, and miR168) showed differential regulation under drought conditions. Analysis to identify negative correlations between the differentially expressed miRNAs (DEMs) and DEGs revealed 92 mRNA-miRNA interactions between CL and DL plants, and 32 mRNA-miRNA interactions between DL and WL plants.

CONCLUSIONS

This study provides a global view of the transcriptional and the post-transcriptional responses of tobacco under drought stress and re-watering conditions. Our results establish an empirical foundation that should prove valuable for further investigations into the molecular mechanisms through which tobacco, and plants more generally, respond to drought stress at multiple molecular genetic levels.

Development and Application of Transcriptome-Derived Microsatellites in Actinidia eriantha (Actinidiaceae)

Actinidia eriantha Benth. is a diploid perennial woody vine native to China and is recognized as a valuable species for commercial kiwifruit improvement with high levels of ascorbic acid as well as having been used in traditional Chinese medicine. Due to the lack of genomic resources for the species, microsatellite markers for population genetics studies are scarce. In this study, RNASeq was conducted on fruit tissue of A. eriantha, yielding 5,678,129 reads with a total output of 3.41 Gb. De novo assembly yielded 69,783 non-redundant unigenes (41.3 Mb), of which 21,730 were annotated using protein databases. A total of 8,658 EST-SSR loci were identified in 7,495 unigene sequences, for which primer pairs were successfully designed for 3,842 loci (44.4%). Among these, 183 primer pairs were assayed for PCR amplification, yielding 69 with detectable polymorphism in A. eriantha. Additionally, 61 of the 69 polymorphic loci could be successfully amplified in at least one other Actinidia species. Of these, 14 polymorphic loci (mean N_A = 6.07 ± 2.30) were randomly selected for assessing levels of genetic diversity and population structure within A. eriantha. Finally, a neighbor-joining tree and Bayesian clustering analysis showed distinct clustering into two groups (K = 2), agreeing with the geographical distributions of these populations. Overall, our results will facilitate further studies of genetic diversity within A. eriantha and will aid in discriminating outlier loci involved in local adaptation.

Transcriptome Profiling of Watermelon Root in Response to Short-Term Osmotic Stress

Osmotic stress adversely affects the growth, fruit quality and yield of watermelon (Citrullus lanatus (Thunb.) Matsum. & Nakai). Increasing the tolerance of watermelon to osmotic stress caused by factors such as high salt and water deficit is an effective way to improve crop survival in osmotic stress environments. Roots are important organs in water absorption and are involved in the initial response to osmosis stress; however, few studies have examined the underlying mechanism of tolerance to osmotic stress in watermelon roots. For better understanding of this mechanism, the inbred watermelon accession M08, which exhibits relatively high tolerance to water deficits, was treated with 20% polyethylene glycol (PEG) 6000. The root samples were harvested at 6 h after PEG treatment and untreated samples were used as controls. Transcriptome analyses were carried out by Illumina RNA sequencing. A total of 5246 differentially expressed genes were identified. Gene ontology enrichment and biochemical pathway analyses of these 5246 genes showed that short-term osmotic stress affected osmotic adjustment, signal transduction, hormone responses, cell division, cell cycle and ribosome, and M08 may repress root growth to adapt osmotic stress. The results of this study describe the watermelon root transcriptome under osmotic stress and propose new insight into watermelon root responses to osmotic stress at the transcriptome level. Accordingly, these results allow us to better understand the molecular mechanisms of watermelon in response to drought stress and will facilitate watermelon breeding projects to improve drought tolerance.

De novo assembly and functional annotation of Myrciaria dubia fruit transcriptome reveals multiple metabolic pathways for L-ascorbic acid biosynthesis

Background

Myrciaria dubia is an Amazonian fruit shrub that produces numerous bioactive phytochemicals, but is best known by its high L-ascorbic acid (AsA) content in fruits. Pronounced variation in AsA content has been observed both within and among individuals, but the genetic factors responsible for this variation are largely unknown. The goals of this research, therefore, were to assemble, characterize, and annotate the fruit transcriptome of M. dubia in order to reconstruct metabolic pathways and determine if multiple pathways contribute to AsA biosynthesis.

Results

In total 24,551,882 high-quality sequence reads were de novo assembled into 70,048 unigenes (mean length = 1150 bp, N50 = 1775 bp). Assembled sequences were annotated using BLASTX against public databases such as TAIR, GR-protein, FB, MGI, RGD, ZFIN, SGN, WB, TIGR_CMR, and JCVI-CMR with 75.2 % of unigenes having annotations. Of the three core GO annotation categories, biological processes comprised 53.6 % of the total assigned annotations, whereas cellular components and molecular functions comprised 23.3 and 23.1 %, respectively. Based on the KEGG pathway assignment of the functionally annotated transcripts, five metabolic pathways for AsA biosynthesis were identified: animal-like pathway, myo-inositol pathway, L-gulose pathway, D-mannose/L-galactose pathway, and uronic acid pathway. All transcripts coding enzymes involved in the ascorbate-glutathione cycle were also identified. Finally, we used the assembly to identified 6314 genic microsatellites and 23,481 high quality SNPs.

Conclusions

This study describes the first next-generation sequencing effort and transcriptome annotation of a non-model Amazonian plant that is relevant for AsA production and other bioactive phytochemicals. Genes encoding key enzymes were successfully identified and metabolic pathways involved in biosynthesis of AsA, anthocyanins, and other metabolic pathways have been reconstructed. The identification of these genes and pathways is in agreement with the empirically observed capability of M. dubia to synthesize and accumulate AsA and other important molecules, and adds to our current knowledge of the molecular biology and biochemistry of their production in plants. By providing insights into the mechanisms underpinning these metabolic processes, these results can be used to direct efforts to genetically manipulate this organism in order to enhance the production of these bioactive phytochemicals.

The accumulation of AsA precursor and discovery of genes associated with their biosynthesis and metabolism in M. dubia is intriguing and worthy of further investigation. The sequences and pathways produced here present the genetic framework required for further studies. Quantitative transcriptomics in concert with studies of the genome, proteome, and metabolome under conditions that stimulate production and accumulation of AsA and their precursors are needed to provide a more comprehensive view of how these pathways for AsA metabolism are regulated and linked in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2225-6) contains supplementary material, which is available to authorized users.

Transcriptome profiling of differentially expressed genes in floral buds and flowers of male sterile and fertile lines in watermelon

Background

Male sterility is an important mechanism for the production of hybrid seeds in watermelon. Although fruit development has been studied extensively in watermelon, there are no reports on gene expression in floral organs. In this study, RNA-sequencing (RNA-seq) was performed in two near-isogenic watermelon lines (genic male sterile [GMS] line, DAH3615-MS and male fertile line, DAH3615) to identify the differentially expressed genes (DEGs) related to male sterility.

Results

DEG analysis showed that 1259 genes were significantly associated with male sterility at a FDR P-value of < 0.01. Most of these genes were only expressed in the male fertile line. In addition, 11 functional clusters were identified using DAVID functional classification analysis. Of detected genes in RNA-seq analysis, 19 were successfully validated by qRT-PCR.

Conclusions

In this study, we carried out a comprehensive floral transcriptome sequence comparison of a male fertile line and its near-isogenic male sterile line in watermelon. This analysis revealed essential genes responsible for stamen development, including pollen development and pollen tube elongation, and allowed their functional classification. These results provided new information on global mechanisms related to male sterility in watermelon.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2186-9) contains supplementary material, which is available to authorized users.