A comprehensive RNA-Seq-based gene expression atlas of the summer squash (Cucurbita pepo) provides insights into fruit morphology and ripening mechanisms.

BACKGROUND: Summer squash (Cucurbita pepo: Cucurbitaceae) are a popular horticultural crop for which there is insufficient genomic and transcriptomic information. Gene expression atlases are crucial for the identification of genes expressed in different tissues at various plant developmental stages. Here, we present the first comprehensive gene expression atlas for a summer squash cultivar, including transcripts obtained from seeds, shoots, leaf stem, young and developed leaves, male and female flowers, fruits of seven developmental stages, as well as primary and lateral roots. RESULTS: In total, 27,868 genes and 2352 novel transcripts were annotated from these 16 tissues, with over 18,000 genes common to all tissue groups. Of these, 3812 were identified as housekeeping genes, half of which assigned to known gene ontologies. Flowers, seeds, and young fruits had the largest number of specific genes, whilst intermediate-age fruits the fewest. There also were genes that were differentially expressed in the various tissues, the male flower being the tissue with the most differentially expressed genes in pair-wise comparisons with the remaining tissues, and the leaf stem the least. The largest expression change during fruit development was early on, from female flower to fruit two days after pollination. A weighted correlation network analysis performed on the global gene expression dataset assigned 25,413 genes to 24 coexpression groups, and some of these groups exhibited strong tissue specificity. CONCLUSIONS: These findings enrich our understanding about the transcriptomic events associated with summer squash development and ripening. This comprehensive gene expression atlas is expected not only to provide a global view of gene expression patterns in all major tissues in C. pepo but to also serve as a valuable resource for functional genomics and gene discovery in Cucurbitaceae.

Expanding Phaseolus coccineus Genomic Resources: De Novo Transcriptome Assembly and Analysis of Landraces 'Gigantes' and 'Elephantes' Reveals Rich Functional Variation.

Beans are one of the most important staple crops in the world. Runner bean (Phaseolus coccineus L.) is a small-scale agriculture crop compared to common bean (Phaseolusvulgaris). Beans have been introduced to Europe from the Central America to Europe and since then they have been scattered to different geographical regions. This has resulted in the generation of numerous local cultivars and landraces with distinguished characters and adaptive potential. To identify and characterize the underlying genomic variation of two very closely related runner bean cultivars, we performed RNA-Seq with de novo transcriptome assembly in two landraces of P. coccineus, 'Gigantes' and 'Elephantes' phenotypically distinct, differing in seed size and shape. The cleaned reads generated 37,379 and 37,774 transcripts for 'Gigantes' and 'Elephantes,' respectively. A total of 1896 DEGs were identified between the two cultivars, 1248 upregulated in 'Elephantes' and 648 upregulated in 'Gigantes.' A significant upregulation of defense-related genes was observed in 'Elephantes,' among those, numerous members of the AP2-EREBP, WRKY, NAC, and bHLH transcription factor families. In total, 3956 and 4322 SSRs were identified in 'Gigantes' and 'Elephantes,' respectively. Trinucleotide repeats were the most dominant repeat motif, accounting for 41.9% in 'Gigantes' and 40.1% in 'Elephantes' of the SSRs identified, followed by dinucleotide repeats (29.1% in both cultivars). Additionally, 19,281 putative SNPs were identified, among those 3161 were non-synonymous, thus having potential functional implications. High-confidence non-synonymous SNPs were successfully validated with an HRM assay, which can be directly adopted for P. coccineus molecular breeding. These results significantly expand the number of polymorphic markers within P. coccineus genus, enabling the robust identification of runner bean cultivars, the construction of high-resolution genetic maps, potentiating genome-wide association studies. They finally contribute to the genetic reservoir for the improvement of the closely related and intercrossable Phaseolus vulgaris.

Identification and evidence of positive selection upon resistance gene analogs in cotton (Gossypium hirsutum L.).

Upland cotton (Gossypium hirsutum L.) is an important fiber crop species, which is intensively plagued by a plethora of phytopathogenic fungi such as Fusarium oxysporum f. sp. vasinfectum (Fov) causing severe wilt disease. Resistance gene analogs (RGAs) are the largest class of potential resistance (R) genes depicting highly conserved domains and structures in plants. Additionally, RGAs are pivotal components of breeding projects towards host disease resistance, serving as useful functional markers linked to R genes. In this study, a cloning approach based on conserved RGAs motifs was used in order to amplify 38 RGAs from two upland cotton cultivars differing in their Fov susceptibility. Besides, we assessed the phylogenetic expansion and the evolutionary pressures acting upon 127 RGA homologues, which were previously deposited in GenBank along with the 38 RGAs from this study. A total of 165 RGAs sequences were clustered according to their BLAST(P) similarities in ten paralogous genes groups (PGGs). These RGAs exhibited intensive signs of positive selection as it was revealed by inferring various maximum likelihood analyses. The results showed robust signs of positive selection, acting in almost all PGGs across the phylogeny. The evolutionary analysis revealed the existence of 42 positively selected residue sites across the PGG lineages, putatively affecting their ligand-binding specificities. As RGAs derived markers are in close linkage to R genes, these results could be used in ongoing breeding programs of upland cotton.

DNA barcode ITS2 coupled with high resolution melting (HRM) analysis for taxonomic identification of Sideritis species growing in Greece.

Identification of genotypes in Sideritis is complicated owing to the morphological similarity and common occurrence of natural hybridisation within Sideritis species. Species- and genotype-specific DNA markers are very useful for plant identification, breeding and preservation programs. Herein, a real-time polymerase chain reaction (PCR) of ITS2 barcode region coupled with high resolution melting-curve (HRM) analysis was evaluated for an accurate, rapid and sensitive tool for species identification focusing on seven Sideritis species growing in Greece. The HRM assay developed in this study is a rapid and straightforward method for the identification and discrimination of the investigated Sideritis species. This assay is simple compared to other genotyping methods as it does not require DNA sequencing or post-PCR processing. Therefore, this method offers a new alternative for rapid detection of Sideritis species.

Evaluation of a dill (Anethum graveolens L.) gene bank germplasm collection using multivariate analysis of morphological traits, molecular genotyping and chemical composition to identify novel genotypes for plant breeding.

Dill (Anethum graveolens L.) is an aromatic herb widely used in the food industry, with several commercial cultivars available with different qualitative characteristics. Commercial cultivars are usually preferred over landraces due to their higher yield and also the lack of improved landraces than can be commercialized. In Greece, however, traditional dill landraces are cultivated by local communities. Many are conserved in the Greek Gene Bank and the aim here was to investigate and compare the morphological, genetic, and chemical biodiversity of twenty-two Greek landraces and nine modern/commercial cultivars. Multivariate analysis of the morphological descriptors, molecular markers, and essential oil and polyphenol composition revealed that the Greek landraces were clearly distinguished compared with modern cultivars at the level of phenological, molecular and chemical traits. Landraces were typically taller, with larger umbels, denser foliage, and larger leaves. Plant height, density of foliage, density of feathering as well as aroma characteristics were desirable traits observed for some landraces, such as T538/06 and GRC-1348/04, which were similar or superior to those of some commercial cultivars. Polymorphic loci for inter-simple sequence repeat (ISSR) and start codon targeted (SCoT) molecular markers were 76.47% and 72.41% for landraces, and 68.24% and 43.10% for the modern cultivars, respectively. Genetic divergence was shown, but not complete isolation, indicating that some gene flow may have occurred between landraces and cultivars. The major constituent in all dill leaf essential oils was α-phellandrene (54.42-70.25%). Landraces had a higher α-phellandrene and dill ether content than cultivars. Two dill landraces were rich in chlorogenic acid, the main polyphenolic compound determined. The study highlighted for the first-time Greek landraces with desirable characteristics regarding quality, yield, and harvest time suitable for breeding programs to develop new dill cultivars with superior features.

Characterization of PROFILIN genes from allotetraploid (Gossypium hirsutum) cotton and its diploid progenitors and expression analysis in cotton genotypes differing in fiber characteristics.

The actin-binding protein profilin (PRF) plays an important role in cell growth and expansion by regulating the organization of the actin filaments. Recent studies have reported association between fiber elongation in cultivated cotton (Gossypium hirsutum) and PRF expression. In the present study, we cloned four genomic clones from allotetraploid cotton (G. hirsutum) and its putative diploid progenitors (G. arboreum and G. raimondii) designated GhPRF1_A, GhPRF1_D, GaPRF1, and GrPRF1 encoding cotton PRF and characterized their genomic structure, phylogenetic relationships and promoter structure. Sequence analysis of the coding regions of all clones resulted in a single protein product which revealed more than 80% similarity to most plant PRFs and a typical organization with an actin-binding and a polybasic phospholipid binding motif at the carboxy terminus. DNA blot hybridization suggested that PRF gene is present with more than one copy in the allotetraploid species G. hirsutum. Expression analysis performed in various organs of cultivated cotton revealed that the PRF gene was preferentially expressed in cotton fibers. Very low levels of expression were observed in whole flowers, while PRF transcripts were not detected in other organs examined. Furthermore, higher levels of expression were observed at the early stages of cotton fiber development (at 10 days post anthesis), indicative that this gene may play a major role in the early stages of cotton fiber development. Quantitation of the expression by real-time PCR revealed higher expression levels in a G. hirsutum variety with higher fiber percentage compared to a variety with lower percentage. In addition, higher levels of expression were found in cultivated allotetraploid G. barbadense cotton species with higher fiber length in comparison to cultivated allotetraploid G. hirsutum.

Isolation of a CENTRORADIALIS/TERMINAL FLOWER1 homolog in saffron (Crocus sativus L.): characterization and expression analysis.

Genes in the phosphatidyl-ethanolamine-binding protein (PEBP) family are instrumental in regulating the fate of meristems and flowering time. To investigate the role of these genes in the monocotyledonous plant Crocus (Crocus sativus L), an industrially important crop cultivated for its nutritional and medicinal properties, we have cloned and characterized a CENTRORADIALIS/TERMINAL FLOWER1 (CEN/TFL1) like gene, named CsatCEN/TFL1-like, the first reported CEN/TFL1 gene characterized from such a perennial geophyte. Sequence analysis revealed that CsatCEN/TFL1 shows high similarity to its homologous PEBP family genes CEN/TFL1, FT and MFT from a variety of plant species and maintains the same exon/intron organization. Phylogenetic analysis of the CsatCEN/TFL1 amino acid sequence confirmed that the isolated sequences belong to the CEN/TFL1 clade of the PEBP family. CsatCEN/TFL1 transcripts could be detected in corms, flower and flower organs but not in leaves. An alternative spliced transcript was also detected in the flower. Comparison of expression levels of CsatCEN/TFL1 and its alternative spliced transcript in wild type flower and a double flower mutant showed no significant differences. Overexpression of CsatCEN/TFL1 transcript in Arabidopsis tfl1 plants reversed the phenotype of early flowering and terminal flowering of the tfl1 plants to a normal one. Computational analysis of the obtained promoter sequences revealed, next to common binding motifs in CEN/TFL1-like genes as well as other flowering gene promoters, the presence of two CArG binding sites indicative of control of CEN/TFL1 by MADS-box transcription factors involved in crocus flowering and flower organ formation.

Characterization of the Genetic Diversity Present in a Diverse Sesame Landrace Collection Based on Phenotypic Traits and EST-SSR Markers Coupled With an HRM Analysis.

A selection of sesame (Sesamum indicum L.) landraces of different eco-geographical origin and breeding history have been characterized using 28 qualitative morpho-physiological descriptors and seven expressed sequence tag-simple sequence repeat (EST-SSR) markers coupled with a high-resolution melting (HRM) analysis. The most variable qualitative traits that could efficiently discriminate landraces, as revealed by the correlation analyses, were the plant growth type and position of the branches, leaf blade width, stem pubescence, flowering initiation, capsule traits and seed coat texture. The agglomerative hierarchical clustering analysis based on a dissimilarity matrix highlighted three main groups among the sesame landraces. An EST-SSR marker analysis revealed an average polymorphism information content (PIC) value of 0.82, which indicated that the selected markers were highly polymorphic. A principal coordinate analysis and dendrogram reconstruction based on the molecular data classified the sesame genotypes into four major clades. Both the morpho-physiological and molecular analyses showed that landraces from the same geographical origin were not always grouped in the same cluster, forming heterotic groups; however, clustering patterns were observed for the Greek landraces. The selective breeding of such traits could be employed to unlock the bottleneck of local phenotypic diversity and create new cultivars with desirable traits.

famRCA-RACE: a rolling circle amplification race for isolating a family of homologous cDNAs in one reaction and its application to obtain NAC genes transcription factors from crocus (Crocus sativus) flower.

We describe an improvement of the RCA-RACE (rolling circle amplification-rapid amplification of cDNA ends) method, called family RCA-RACE (famRCA-RACE). The method is based on the generation of circular cDNA fragments, followed by rolling circle amplification of the circular cDNA using phi29 DNA polymerase and the application of PCR using degenerate outworking primers, designed for a conserved region of homologous genes, that allows the isolation of homologous cDNA sequences expressed in the mRNA preparation in a single polymerase chain reaction (PCR). As an example we present the isolation of seven NAC-like transcription factors cDNA sequences expressed in Crocus sativus flower, used for saffron production. Sequence alignment revealed that CsatNAC proteins contain the typical domain structure of plant NAC proteins, consisting of the conserved N-terminal NAC domain used to design the primers and the five subdomains. Phylogenetic analysis revealed that CsatNAC proteins fall in subgroup I of the NAC family of proteins.

Effect οf Genotype and Growing Year on the Nutritional, Phytochemical, and Antioxidant Properties of Industrial Hemp (Cannabis sativa L.) Seeds.

Cannabis sativa L. seeds have been an important source of protein, oil, and dietary fiber for human and animals. Currently, there is a growing interest in the commercial products of these seeds, which are recognized as a legitimate source of medicaments, cosmeceuticals, and nutraceuticals. The objective of this study was to investigate the nutritional, phytochemical composition, and antioxidant properties of seeds from seven hemp cultivars grown in Greece for three consecutive years. All the measured parameters strongly varied under the influence of growing year and genotype. In particular, protein, oil, and carbohydrates' content of hemp seeds as well as fatty acids' composition were mainly affected by genotype, whereas the growing year had a major effect on phytochemical components and antioxidant activity, which was determined by the 2,2'-azino-bis (3-ethylbenzthiazoline sulfonate) (ABTS) and ferric-reducing antioxidant power (FRAP) assays. Moreover, a predominant effect of the year was observed for phenolic profiles as determined by high-performance liquid chromatography and total carotenoids' content. This study suggests that hemp seeds could be a promising food crop as a result of their high nutritive traits and antioxidant potential. A comparison of the studied cultivars, showed that Finola seeds had the highest oil and protein contents and, thus, appeared to be the most promising cultivar for cultivation in Greece.

Whole-genome resequencing of Cucurbita pepo morphotypes to discover genomic variants associated with morphology and horticulturally valuable traits.

Cucurbita pepo contains two cultivated subspecies, each of which encompasses four fruit-shape morphotypes (cultivar groups). The Pumpkin, Vegetable Marrow, Cocozelle, and Zucchini Groups are of subsp. pepo and the Acorn, Crookneck, Scallop, and Straightneck Groups are of subsp. ovifera. Recently, a de novo assembly of the C. pepo subsp. pepo Zucchini genome was published, providing insights into its evolution. To expand our knowledge of evolutionary processes within C. pepo and to identify variants associated with particular morphotypes, we performed whole-genome resequencing of seven of these eight C. pepo morphotypes. We report for the first time whole-genome resequencing of the four subsp. pepo (Pumpkin, Vegetable Marrow, Cocozelle, green Zucchini, and yellow Zucchini) morphotypes and three of the subsp. ovifera (Acorn, Crookneck, and Scallop) morphotypes. A high-depth resequencing approach was followed, using the BGISEQ-500 platform that enables the identification of rare variants, with an average of 33.5X. Approximately 94.5% of the clean reads were mapped against the reference Zucchini genome. In total, 3,823,977 high confidence single-nucleotide polymorphisms (SNPs) were identified. Within each accession, SNPs varied from 636,918 in green Zucchini to 2,656,513 in Crookneck, and were distributed homogeneously along the chromosomes. Clear differences between subspecies pepo and ovifera in genetic variation and linkage disequilibrium are highlighted. In fact, comparison between subspecies pepo and ovifera indicated 5710 genes (22.5%) with Fst > 0.80 and 1059 genes (4.1%) with Fst = 1.00 as potential candidate genes that were fixed during the independent evolution and domestication of the two subspecies. Linkage disequilibrium was greater in subsp. ovifera than in subsp. pepo, perhaps reflective of the earlier differentiation of morphotypes within subsp. ovifera. Some morphotype-specific genes have been localized. Our results offer new clues that may provide an improved understanding of the underlying genomic regions involved in the independent evolution and domestication of the two subspecies. Comparisons among SNPs unique to particular subspecies or morphotypes may provide candidate genes responsible for traits of high economic importance.

Heterotopic expression of B-class floral homeotic genes PISTILLATA/GLOBOSA supports a modified model for crocus (Crocus sativus L.) flower formation.

For uncovering and understanding the molecular mechanisms controlling flower development in cultivated Crocus sativus and particularly the transformation of sepals in outer whorl (whorl 1) tepals, we have cloned and characterized the expression of a family of five PISTILLATA/GLOBOSA-like (PI/GLO-like) MADS-box genes expressed in the C. sativus flower. The deduced amino acid sequences of the coded proteins indicated high homology with members of the MADS-box family of transcription factors, and particularly with other members of the PI/GLO family of MADS-box proteins that control floral organ identity. PI/GLO expression studies in cultivated C. sativus uncover the presence of PI/GLO transcripts not only in the second and third whorls of flower organs as expected, but also in the outer whorl tepals that are the sepals in most typical flowers. This heterotopic expression of both B-class genes: PI/GLO and AP3/DEF, known to form heterodimers for stamens and petals (petaloid inner whor l-whorl 2-tepals in C. sativus), explains the homeotic transformation of sepals into outer whorl tepals in this species. Analysis of PI/GLO sequences from C. sativus for putative targets to known micro-RNAs (miRNAs) showed that the target site for ath-miRNA167 found in Arabidopsis thaliana PI is not present in C. sativus, however, the PI/GLO sequences may be regulated by an ath-miRNA163.

Cloning, structural characterization, and phylogenetic analysis of flower MADS-box genes from crocus (Crocus sativus L.).

Crocus (Crocus sativus L.) is a crop species cultivated for its flowers and, more specifically, for its red stigmas. The flower of crocus is bisexual and sterile, since crocus is a triploid species. Its perianth consists of six petaloid tepals: three tepals in whorl 1 (outer tepals) and three tepals in whorl 2 (inner tepals). The androecium consists of three distinct stamens and the gynoecium consists of a single compound pistil with three carpels, a single three-branched style, and an inferior ovary. The dry form of the stigmas constitutes the commercial saffron used as a food additive, in the coloring industry, and in medicine. In order to uncover and understand the molecular mechanisms controlling flower development in cultivated crocus and its relative wild progenitor species, and characterize a number of crocus flower mutants, we have cloned and characterized different, full-length, cDNA sequences encoding MADS-box transcription factor proteins involved in flower formation. Here we review the different methods followed or developed for obtaining these sequences involving conventional 5 inverted exclamation markä 3 inverted exclamation markä RACE, as well as newly developed methods from our group, named Rolling Circle Amplification C RACE (RCA-RACE) and its modification named familyRCA-RACE (famRCA-RACE). Furthermore, the characteristics of the protein structure and their common and specific domains for each type of MADS-box transcription factors in this lower nongrass monocot belonging to the Iridaceae family are described. Finally, a phylogenetic tree of all the MADS-box sequences available in our lab is presented and discussed in relation to other data from studies of species of the Iridaceae group and closely related families from an evolutionary perspective. The structural and phylogenetic analyses are based on both published and unpublished data.

Comparative metagenomics reveals alterations in the soil bacterial community driven by N-fertilizer and Amino 16® application in lettuce.

Nutrients in the form of fertilizers and/or other additives such as amino acids, dramatically influence plant development and growth, plant nutrient composition and the level of soil pollution. Moreover, the treatment of soil microbiota is emerging as a new strategy in plant breeding to achieve desirable traits. Thus, integrated study of fertilizer application and soil microbiota might lead to a better understanding of soil-plant interactions and inform the design of novel ways to fertilize plants. Herein we report metagenomics data for soil microbiota in lettuce (Lactuca sativa) treated with fertilizer, amino acids or their combinations as follows: N-fertilizer + Amino16®, Amino16®, N-fertilizer and no treatment control. Data have been deposited in the NCBI Sequence Read Archive (SRA) (accession number: PRJNA388765).

De novo comparative transcriptome analysis of genes involved in fruit morphology of pumpkin cultivars with extreme size difference and development of EST-SSR markers.

The genetic basis of fruit size and shape was investigated for the first time in Cucurbita species and genetic loci associated with fruit morphology have been identified. Although extensive genomic resources are available at present for tomato (Solanum lycopersicum), cucumber (Cucumis sativus), melon (Cucumis melo) and watermelon (Citrullus lanatus), genomic databases for Cucurbita species are limited. Recently, our group reported the generation of pumpkin (Cucurbita pepo) transcriptome databases from two contrasting cultivars with extreme fruit sizes. In the current study we used these databases to perform comparative transcriptome analysis in order to identify genes with potential roles in fruit morphology and fruit size. Differential Gene Expression (DGE) analysis between cv. 'Munchkin' (small-fruit) and cv. 'Big Moose' (large-fruit) revealed a variety of candidate genes associated with fruit morphology with significant differences in gene expression between the two cultivars. In addition, we have set the framework for generating EST-SSR markers, which discriminate different C. pepo cultivars and show transferability to related Cucurbitaceae species. The results of the present study will contribute to both further understanding the molecular mechanisms regulating fruit morphology and furthermore identifying the factors that determine fruit size. Moreover, they may lead to the development of molecular marker tools for selecting genotypes with desired morphological traits.

The study of the E-class SEPALLATA3-like MADS-box genes in wild-type and mutant flowers of cultivated saffron crocus (Crocus sativus L.) and its putative progenitors.

To further understand flowering and flower organ formation in the monocot crop saffron crocus (Crocus sativus L.), we cloned four MIKC(c) type II MADS-box cDNA sequences of the E-class SEPALLATA3 (SEP3) subfamily designated CsatSEP3a/b/c/c_as as well as the three respective genomic sequences. Sequence analysis showed that cDNA sequences of CsatSEP3 c and c_as are the products of alternative splicing of the CsatSEP3c gene. Bioinformatics analysis with putative orthologous sequences from various plant species suggested that all four cDNA sequences encode for SEP3-like proteins with characteristic motifs and amino acids, and highlighted intriguing sequence features. Phylogenetically, the isolated sequences were closest to the SEP3-like genes from monocots such as Asparagus virgatus, Oryza sativa, Zea mays, and the dicot Arabidopsis SEP3 gene. All four isolated C. sativus sequences were strongly expressed in flowers and in all flower organs: whorl1 tepals, whorl2 tepals, stamens and carpels, but not in leaves. Expression of CsatSEP3a/b/c/c_as cDNAs was compared in wild-type and mutant flowers. Expression of the isolatedCsatSEP3-like genes in whorl1 tepals together with E-class CsatAP1/FUL subfamily and B-class CsatAP3 and CsatPI subfamilies of genes, fits the ABCE "quartet model," an extended form of the original ABC model proposed to explain the homeotic transformation of whorl1 sepals into whorl1 tepals in Liliales and Asparagales plants such as C. sativus. This conclusion was also supported by the interaction of the CsatSEP3b protein with CsatAP1/FUL and CsatAP3 proteins. In contrast, expression of both B-class CsatAP3 and CsatPI genes and the C-class CsatAGAMOUS genes together with E-class CsatSEP3-like genes in carpels, without any phenotypic effects on carpels, raises questions about the role of these gene classes in carpel formation in this non-grass monocot and requires further experimentation. Finally, taking advantage of the size and sequence differences in amplified genomic sequences of the triploid C. sativus and comparing them with the respective sequences from C. tomasii, C. hadriaticus and C. cartwrightianus, three putative wild-type diploid progenitor species, we examined the origin of CsatSEP3a sequence.