Next generation sequencing and omics in cucumber (Cucumis sativus L.) breeding directed research

Analysis of the Complete Mitochondrial Genome of the Bitter Gourd (Momordica charantia)

Bitter gourd (Momordica charantia L.) is a significant vegetable. Although it has a special bitter taste, it is still popular with the public. The industrialization of bitter gourd could be hampered by a lack of genetic resources. The bitter gourd's mitochondrial and chloroplast genomes have not been extensively studied. In the present study, the mitochondrial genome of bitter gourd was sequenced and assembled, and its substructure was investigated. The mitochondrial genome of bitter gourd is 331,440 bp with 24 unique core genes, 16 variable genes, 3 rRNAs, and 23 tRNAs. We identified 134 SSRs and 15 tandem repeats in the entire mitochondrial genome of bitter gourd. Moreover, 402 pairs of repeats with a length greater than or equal to 30 were observed in total. The longest palindromic repeat was 523 bp, and the longest forward repeat was 342 bp. We found 20 homologous DNA fragments in bitter gourd, and the summary insert length was 19,427 bp, accounting for 5.86% of the mitochondrial genome. We predicted a total of 447 potential RNA editing sites in 39 unique PCGs and also discovered that the ccmFN gene has been edited the most often, at 38 times. This study provides a basis for a better understanding and analysis of differences in the evolution and inheritance patterns of cucurbit mitochondrial genomes.

Multiple Stressors in Vegetable Production: Insights for Trait-Based Crop Improvement in Cucurbits

Vegetable production is a key determinant of contribution from the agricultural sector toward national Gross Domestic Product in a country like India, the second largest producer of fresh vegetables in the world. This calls for a careful scrutiny of the threats to vegetable farming in the event of climate extremes, environmental degradation and incidence of plant pests/diseases. Cucurbits are a vast group of vegetables grown almost throughout the world, which contribute to the daily diet on a global scale. Increasing food supply to cater to the ever-increasing world population, calls for intensive, off-season and year-round cultivation of cucurbits. Current situation predisposes these crops to a multitude of stressors, often simultaneously, under field conditions. This scenario warrants a systematic understanding of the different stress specific traits/mechanisms/pathways and their crosstalk that have been examined in cucurbits and identification of gaps and formulation of perspectives on prospective research directions. The careful dissection of plant responses under specific production environments will help in trait identification for genotype selection, germplasm screens to identify superior donors or for direct genetic manipulation by modern tools for crop improvement. Cucurbits exhibit a wide range of acclimatory responses to both biotic and abiotic stresses, among which a few like morphological characters like waxiness of cuticle; primary and secondary metabolic adjustments; membrane thermostability, osmoregulation and, protein and reactive oxygen species homeostasis and turnover contributing to cellular tolerance, appear to be common and involved in cross talk under combinatorial stress exposures. This is assumed to have profound influence in triggering system level acclimation responses that safeguard growth and metabolism. The possible strategies attempted such as grafting initiatives, molecular breeding, novel genetic manipulation avenues like gene editing and ameliorative stress mitigation approaches, have paved way to unravel the prospects for combined stress tolerance. The advent of next generation sequencing technologies and big data management of the omics output generated have added to the mettle of such emanated concepts and ideas. In this review, we attempt to compile the progress made in deciphering the biotic and abiotic stress responses of cucurbits and their associated traits, both individually and in combination.

A novel insight into transcriptional and epigenetic regulation underlying sex expression and flower development in melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important cucurbit and has been considered as a model plant for studying sex determination. The four most common sexual morphotypes in melon are monoecious (A-G-M), gynoecious (--ggM-), andromonoecious (A-G-mm), and hermaphrodite (--ggmm). Sex expression in melons is complex, as the genes and associated networks that govern the sex expression are not fully explored. Recently, RNA-seq transcriptomic profiling, ChIP-qPCR analysis integrated with gene ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathways predicted the differentially expressed genes including sex-specific ACS and ACO genes, in regulating the sex-expression, phytohormonal cross-talk, signal transduction, and secondary metabolism in melons. Integration of transcriptional control through genetic interaction in between the ACS7, ACS11, and WIP1 in epistatic or hypostatic manner, along with the recruitment of H3K9ac and H3K27me3, epigenetically, overall determine sex expression. Alignment of protein sequences for establishing phylogenetic evolution, motif comparison, and protein-protein interaction supported the structural conservation while presence of the conserved hydrophilic and charged residues across the diverged evolutionary group predicted the functional conservation of the ACS protein. Presence of the putative cis-binding elements or DNA motifs, and its further comparison with DAP-seq-based cistrome and epicistrome of Arabidopsis, unraveled strong ancestry of melons with Arabidopsis. Motif comparison analysis also characterized putative genes and transcription factors involved in ethylene biosynthesis, signal transduction, and hormonal cross-talk related to sex expression. Overall, we have comprehensively reviewed research findings for a deeper insight into transcriptional and epigenetic regulation of sex expression and flower development in melons.

Agarose-resolvable InDel markers based on whole genome re-sequencing in cucumber

Insertion and Deletion (InDel) are common features in genomes and are associated with genetic variation. The whole-genome re-sequencing data from two parents (X1 and X2) of the elite cucumber (Cucumis sativus) hybrid variety Lvmei No.1 was used for genome-wide InDel polymorphisms analysis. Obtained sequence reads were mapped to the genome reference sequence of Chinese fresh market type inbred line ‘9930’ and gaps conforming to InDel were pinpointed. Further, the level of cross-parents polymorphism among five pairs of cucumber breeding parents and their corresponding hybrid varieties were used for evaluating hybrid seeds purity test efficiency of InDel markers. A panel of 48 cucumber breeding lines was utilized for PCR amplification versatility and phylogenetic analysis of these markers. In total, 10,470 candidate InDel markers were identified for X1 and X2. Among these, 385 markers with more than 30 nucleotide difference were arbitrary chosen. These markers were selected for experimental resolvability through electrophoresis on an Agarose gel. Two hundred and eleven (211) accounting for 54.81% of markers could be validated as single and clear polymorphic pattern while 174 (45.19%) showed unclear or monomorphic genetic bands between X1 and X2. Cross-parents polymorphism evaluation recorded 68 (32.23%) of these markers, which were designated as cross-parents transferable (CPT) InDel markers. Interestingly, the marker InDel114 presented experimental transferability between cucumber and melon. A panel of 48 cucumber breeding lines including parents of Lvmei No. 1 subjected to PCR amplification versatility using CPT InDel markers successfully clustered them into fruit and common cucumber varieties based on phylogenetic analysis. It is worth noting that 16 of these markers were predominately associated to enzymatic activities in cucumber. These agarose-based InDel markers could constitute a valuable resource for hybrid seeds purity testing, germplasm classification and marker-assisted breeding in cucumber.

A high-quality cucumber genome assembly enhances computational comparative genomics

Genetic variation is expressed by the presence of polymorphisms in compared genomes of individuals that can be transferred to next generations. The aim of this work was to reveal genome dynamics by predicting polymorphisms among the genomes of three individuals of the highly inbred B10 cucumber (Cucumis sativus L.) line. In this study, bioinformatic comparative genomics was used to uncover cucumber genome dynamics (also called real-time evolution). We obtained a new genome draft assembly from long single molecule real-time (SMRT) sequencing reads and used short paired-end read data from three individuals to analyse the polymorphisms. Using this approach, we uncovered differentiation aspects in the genomes of the inbred B10 line. The newly assembled genome sequence (B10v3) has the highest contiguity and quality characteristics among the currently available cucumber genome draft sequences. Standard and newly designed approaches were used to predict single nucleotide and structural variants that were unique among the three individual genomes. Some of the variant predictions spanned protein-coding genes and their promoters, and some were in the neighbourhood of annotated interspersed repetitive elements, indicating that the highly inbred homozygous plants remained genetically dynamic. This is the first bioinformatic comparative genomics study of a single highly inbred plant line. For this project, we developed a polymorphism prediction method with optimized precision parameters, which allowed the effective detection of small nucleotide variants (SNVs). This methodology could significantly improve bioinformatic pipelines for comparative genomics and thus has great practical potential in genomic metadata handling.

Comparative transcriptome analysis reveals new molecular pathways for cucumber genes related to sex determination

Transcriptome data and qPCR analysis revealed new insight into genes regulatory mechanism related to cucumber sex determination. Cucumber (Cucumis sativus L.) is an economically important crop cultivated worldwide. Enhancing the genomic resources for cucumber may enable the regulation of traits relevant to crop productivity and quality. Sequencing technologies and bioinformatics tools provide opportunities for the development of such resources. The aims of this study were to identify and characterize the genes involved in sex determination and flower morphogenesis in cucumber isogenic lines that differed regarding flower sex type. We obtained transcripts for 933 genes related to shoot apex development, among which 310 were differentially expressed genes (DEGs) among the male, female, and hermaphroditic lines. We performed gene ontology and molecular network analyses and explored the DEGs related to already known processes like: hormone synthesis and signaling, lipid and sugar metabolism; and also newly discovered processes related to cell wall, membrane, and cytoskeleton modifications; ion homeostasis which appears to be important for ethylene perception and signaling, and genes expression mediated by transcription factors related to floral organ identities. We proposed a new model of regulatory mechanism network of sex development in cucumber. Our results may be useful for clarifying the molecular genetics and the functional mechanisms underlying the sex determination processes.

Transcriptome Analyses Provide Novel Insights into Heat Stress Responses in Chieh-Qua (Benincasa hispida Cogn. var. Chieh-Qua How)

Temperature rising caused by global warming has imposed significant negative effects on crop qualities and yields. To get the well-known molecular mechanism upon the higher temperature, we carefully analyzed the RNA sequencing-based transcriptomic responses of two contrasting chieh-qua genotypes: A39 (heat-tolerant) and H5 (heat-sensitive). In this study, twelve cDNA libraries generated from A39 and H5 were performed with a transcriptome assay under normal and heat stress conditions, respectively. A total of 8705 differentially expressed genes (DEGs) were detected under normal conditions (3676 up-regulated and 5029 down-regulated) and 1505 genes under heat stress (914 up-regulated and 591 down-regulated), respectively. A significant positive correlation between RNA-Seq data and qRT-PCR results was identified. DEGs related to heat shock proteins (HSPs), ubiquitin-protein ligase, transcriptional factors, and pentatricopeptide repeat-containing proteins were significantly changed after heat stress. Several genes, which encoded HSPs (CL2311.Contig3 and CL6612.Contig2), cytochrome P450 (CL4517.Contig4 and CL683.Contig7), and bHLH TFs (CL914.Contig2 and CL8321.Contig1) were specifically induced after four days of heat stress. DEGs detected in our study between these two contrasting cultivars would provide a novel basis for isolating useful candidate genes of heat stress responses in chieh-qua.

Genetic mapping of psl locus and quantitative trait loci for angular leaf spot resistance in cucumber (Cucumis sativus L.)

One of the most important cucumber diseases is bacterial angular leaf spot (ALS), whose increased occurrence in open-field production has been observed over the last years. To map ALS resistance genes, a recombinant inbred line (RIL) mapping population was developed from a narrow cross of cucumber line Gy14 carrying psl resistance gene and susceptible B10 line. Parental lines and RILs were tested under growth chamber conditions as well as in the field for angular leaf spot symptoms. Based on simple sequence repeat and DArTseq, genotyping a genetic map was constructed, which contained 717 loci in seven linkage groups, spanning 599.7 cM with 0.84 cM on average between markers. Monogenic inheritance of the lack of chlorotic halo around the lesions, which is typical for ALS resistance and related with the presence of recessive psl resistance gene, was confirmed. The psl locus was mapped on cucumber chromosome 5. Two major quantitative trait loci (QTL) psl5.1 and psl5.2 related to disease severity were found and located next to each other on chromosome 5; moreover, psl5.1 was co-located with psl locus. Identified QTL were validated in the field experiment. Constructed genetic map and markers linked to ALS resistance loci are novel resources that can contribute to cucumber breeding programs.

Whole-Transcriptome Sequence Analysis of Verbena bonariensis in Response to Drought Stress

Drought is an important abiotic factor that threatens the growth and development of plants. Verbena bonariensis is a widely used landscape plant with a very high ornamental value. We found that Verbena has drought tolerance in production practice, so in order to delve into its mechanism of drought resistance and screen out its drought-resistance genes, we used the RNA-Seq platform to perform a de novo transcriptome assembly to analyze Verbena transcription response to drought stress. By high-throughput sequencing with Illumina Hiseq Xten, a total of 44.59 Gb clean data was obtained from T01 (control group) and T02 (drought experiment group). After assembly, 111,313 unigenes were obtained, and 53,757 of them were annotated by compared databases. In this study, 4829 differentially expressed genes were obtained, of which 4165 were annotated. We performed GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analyses, and explored a lot of differently expressed genes related to plant energy production, hormone synthesis, cell signal transduction, and metabolism to understand the stress response of Verbena in drought stress. In addition, we also found that a series of TFs related to drought-resistance of Verbena and provide excellent genetic resources for improving the drought tolerance of crops.

Transcriptomic analysis of differentially expressed genes in flower-buds of genetic male sterile and wild type cucumber by RNA sequencing

Cucumber (Cucumis sativus L.) pollen development involves a diverse range of gene interactions between sporophytic and gametophytic tissues. Previous studies in our laboratory showed that male sterility was controlled by a single recessive nuclear gene, and occurred in pollen mother cell meiophase. To fully explore the global gene expression and identify genes related to male sterility, a RNA-seq analysis was adopted in this study. Young male flower-buds (1-2 mm in length) from genetic male sterility (GMS) mutant and homozygous fertile cucumber (WT) were collected for two sequencing libraries. Total 545 differentially expressed genes (DEGs), including 142 up-regulated DEGs and 403 down-regulated DEGs, were detected in two libraries (Fold Change ≥ 2, FDR < 0.01). These genes were involved in a variety of metabolic pathways, like ethylene-activated signaling pathway, sporopollenin biosynthetic pathway, cell cycle and DNA damage repair pathway. qRT-PCR analysis was performed and showed that the correlation between RNA-Seq and qRT-PCR was 0.876. These findings contribute to a better understanding of the mechanism that leads to GMS in cucumber.

Comparative analysis of miRNA and mRNA abundance in determinate cucumber by high-throughput sequencing

Determinate cucumber is a type of short vines, fewer nodes, and terminal flowers, it is suitable for high-density planting and available harvesting in field cultivation, whereas the indeterminate cucumber is preferred to cultivate in greenhouses. However, many biotic or abiotic stresses could lead indeterminate cucumber to be determinate in greenhouse cultivation, which may decrease yield and fruit quality. Therefore, it is urgent and essential to investigate the key factors forming determinate and terminal flowering in cucumber. In this study, two close background inbred lines were selected and conducted the miRNA and mRNA high throughput sequencing. Interestingly, ethylene-associated miRNAs and mRNAs were intensively obtained, indicating that the plant hormone ethylene is a key factor impacting determinate and terminal flowering in cucumber. The ethylene metabolites analysis showed that significant higher ethylene was observed in determinate line than that in the indeterminate line. The RT-qPCR validation of ethylene related miRNAs Cas-miR172, Cas-miR396, and Cas-miR414 and their target mRNAs showed a significant negative correlation. These data suggested that ethylene-associated miRNAs might affect determinate and terminal flower phenotypes by regulating their target genes expression. This study not only provides a potential molecular mechanism for determinate formation in cucumber but also establishes a method to demonstrate important physiological processes through the comprehensive association of miRNA and mRNA high-throughput sequencing.

Post genomics era for orchid research

Among 300,000 species in angiosperms, Orchidaceae containing 30,000 species is one of the largest families. Almost every habitats on earth have orchid plants successfully colonized, and it indicates that orchids are among the plants with significant ecological and evolutionary importance. So far, four orchid genomes have been sequenced, including Phalaenopsis equestris, Dendrobium catenatum, Dendrobium officinale, and Apostaceae shengen. Here, we review the current progress and the direction of orchid research in the post genomics era. These include the orchid genome evolution, genome mapping (genome-wide association analysis, genetic map, physical map), comparative genomics (especially receptor-like kinase and terpene synthase), secondary metabolomics, and genome editing.