The grapevine homeobox gene VvHB58 influences seed and fruit development through multiple hormonal signaling pathways

Integrative genomics reveals the polygenic basis of seedlessness in grapevine

Seedlessness is a crucial quality trait in table grape (Vitis vinifera L.) breeding. However, the development of seeds involved intricate regulations, and the polygenic basis of seed abortion remains unclear. Here, we combine comparative genomics, population genetics, quantitative genetics, and integrative genomics to unravel the evolution and polygenic basis of seedlessness in grapes. We generated the haplotype-resolved genomes for two seedless grape cultivars, "Thompson Seedless" (TS, syn. "Sultania") and "Black Monukka" (BM). Comparative genomics identified a ∼4.25 Mb hemizygous inversion on Chr10 specific in seedless cultivars, with seedless-associated genes VvTT16 and VvSUS2 located at breakpoints. Population genomic analyses of 548 grapevine accessions revealed two distinct clusters of seedless cultivars, and the identity-by-descent (IBD) results indicated that the origin of the seedlessness trait could be traced back to "Sultania." Introgression, rather than convergent selection, shaped the evolutionary history of seedlessness in grape improvement. Genome-wide association study (GWAS) analysis identified 110 quantitative trait loci (QTLs) associated with 634 candidate genes, including previously unidentified candidate genes, such as three 11S GLOBULIN SEED STORAGE PROTEIN and two CYTOCHROME P450 genes, and well-known genes like VviAGL11. Integrative genomic analyses resulted in 339 core candidate genes categorized into 13 functional categories related to seed development. Machine learning-based genomic selection achieved a remarkable prediction accuracy of 97% for seedlessness in grapevines. Our findings highlight the polygenic nature of seedlessness and provide candidate genes for molecular genetics and an effective prediction for seedlessness in grape genomic breeding.

In silico genome-wide analysis of homeodomain-leucine zipper transcription factors in Cannabis sativa L

HD-Zip (Homeodomain-Leucine Zipper) is a family of transcription factors unique to higher plants and plays a vital role in plant growth and development. Increasing research results show that HD-Zip transcription factors are widely involved in many life processes in plants. However, the HD-Zip transcription factor for cannabis, a valuable crop, has not yet been identified. The sequence characteristics, chromosome localization, system evolution, conservative motif, gene structure, and gene expression of the HD-Zip transcription factor in the cannabis genome were systematically studied. Real-time quantitative polymerase chain reaction (qRT-PCR) was used to verify its function. The results showed that cannabis contained 33 HD-Zip gene members. The number of amino acids is 136-849aa, the isoelectric point is 4.54-9.04, and the molecular weight is 23264.32-93147.87Da. Many cis-acting elements are corresponding to hormone and abiotic stress in the HD-Zip family promoter area of cannabis. Sequencing of the transcriptome at 5 tissue sites of hemp, stems, leaves, bracts, and seeds showed similar levels of expression of 33 members of the HD-Zip gene family at 5 tissue sites. Bioinformatics results show that HD-Zip expression is tissue-specific and may be influenced by hormones and environmental factors. This lays a foundation for further research on the gene function of HD-Zip.

Genome-wide identification of B-box zinc finger (BBX) gene family in Medicago sativa and their roles in abiotic stress responses

BACKGROUND

B-box (BBX) family is a class of zinc finger transcription factors (TFs) that play essential roles in regulating plant growth, development, as well as abiotic stress. However, no systematic analysis of BBX genes has yet been conducted in alfalfa (Medica go sativa L.), and their functions have not been elucidated up to now.

RESULTS

In this study, 28 MsBBX genes were identified from the alfalfa genome, which were clustered into 4 subfamilies according to an evolutionary tree of BBX proteins. Exon-intron structure and conserved motif analysis reflected the evolutionary conservation of MsBBXs in alfalfa. Collinearity analysis showed that segmental duplication promoted the expansion of the MsBBX family. Analysis of cis-regulatory elements suggested that the MsBBX genes possessed many growth/development-, light-, phytohormone-, and abiotic stress-related elements. MsBBX genes were differentially expressed in leaves, flowers, pre-elongated stems, elongated stems, roots and nodules, and most MsBBXs were remarkably induced by drought, salt and various plant growth regulators (ABA, JA, and SA). Further functional verification demonstrated that overexpressing of the MsBBX11 gene clearly promoted salt tolerance in transgenic Arabidopsis by regulating growth and physiological processes of seedlings.

CONCLUSIONS

This research provides insights into further functional research and regulatory mechanisms of MsBBX family genes under abiotic stress of alfalfa.

Evolutionary Consequences of Functional and Regulatory Divergence of HD-Zip I Transcription Factors as a Source of Diversity in Protein Interaction Networks in Plants

The HD superfamily has been studied in detail for several decades. The plant-specific HD-Zip I subfamily attracts the most attention because of its involvement in plant development and stress responses. In this review, we provide a comprehensive insight into the evolutionary events responsible for the functional redundancy and diversification of the HD-Zip I genes in regulating various biological processes. We summarized the evolutionary history of the HD-Zip family, highlighting the important role of WGDs in its expansion and divergence of retained duplicates in the genome. To determine the relationship between the evolutionary origin and functional conservation of HD-Zip I in different species, we performed a phylogenetic analysis, compared their expression profiles in different tissues and under stress and traced the role of orthologs and paralogs in regulating developmental processes. We found that HD-Zip I from different species have similar gene structures with a highly conserved HD and Zip, bind to the same DNA sequences and are involved in similar biological processes. However, they exhibit a functional diversity, which is manifested in altered expression patterns. Some of them are involved in the regulation of species-specific leaf morphology and phenotypes. Here, we discuss the role of changes in functional domains involved in DNA binding and protein interaction of HD-Zip I and in cis-regulated regions of its target genes in promoting adaptive innovations through the formation of de novo regulatory systems. Understanding the role of the HD-Zip I subfamily in organism-environment interactions remains a challenge for evolutionary developmental biology (evo-devo).

Fine mapping of the major gene BhHLS1 controlling seed size in wax gourd (Benincasa hispida)

Introduction/Background

The seed size of wax gourds is an important agronomic trait; however, the associated genes have not yet been reported.

Methods

In this study, we used a high-density genetic map constructed based on F8 recombinant inbred line populations derived from a cross between MY-1 (large seed) and GX-71 (small seed) strains to detect quantitative trait locis (QTLs) for seed-size-related traits in wax gourd over a two-year period.

Results

Two stable QTLs (qSL10 and qSW10) for seed length (SL) and seed width (SW) on chromosome 10 were repeatedly detected over two years (2021-2022). qSL10 had a phenotypic variation rate of 75.30% and 80.80% in 2021 and 2022, respectively. Whereas, qSW10 had a phenotypic variation rate of 66.60% and 73.80% in 2021 and 2022, respectively. Further, a single nucleotide polymorphism mutation was found to cause early termination of Bch10G006400 (BhHLS1) translation in GX-71 through sequencing analysis of candidate genes. Based on gene functional annotation and quantitative real-time PCR analyses, BhHLS1 encoded a probable N-acetyltransferase HLS1-like protein and its expression level was significantly different between parents. Therefore, BhHLS1 is a major candidate gene associated with a one-factor polymorphism regulating the SL and SW of wax gourds. Finally, based on variation in the BhHLS1 sequence, a cleaved amplified polymorphic sequence marker was developed for the molecular marker-assisted breeding of wax gourds.

Discussion

Overall, this study is of great significance for the genetic improvement of seed size, verification of gene functions, and cultivation of specific germplasm resources for wax gourds.

Construction of a high density genetic linkage map to define the locus conferring seedlessness from Mukaku Kishu mandarin

Mukaku Kishu ('MK'), a small sized mandarin, is an important source of seedlessness in citrus breeding. Identification and mapping the gene(s) governing 'MK' seedlessness will expedite seedless cultivar development. In this study, two 'MK'-derived mapping populations- LB8-9 Sugar Belle^® ('SB') × 'MK' (N=97) and Daisy ('D') × 'MK' (N=68) were genotyped using an Axiom_Citrus56 Array encompassing 58,433 SNP probe sets, and population specific male and female parent linkage maps were constructed. The parental maps of each population were integrated to produce sub-composite maps, which were further merged to develop a consensus linkage map. All the parental maps (except 'MK_D') had nine major linkage groups, and contained 930 ('SB'), 810 ('MK_SB'), 776 ('D') and 707 ('MK_D') SNPs. The linkage maps displayed 96.9 ('MK_D') to 98.5% ('SB') chromosomal synteny with the reference Clementine genome. The consensus map was comprised of 2588 markers including a phenotypic seedless (Fs)-locus and spanned a genetic distance of 1406.84 cM, with an average marker distance of 0.54 cM, which is substantially lower than the reference Clementine map. For the phenotypic Fs-locus, the distribution of seedy and seedless progenies in both 'SB' × 'MK' (55:42, χ² = 1.74) and 'D' × 'MK' populations (33:35, χ² = 0.06) followed a test cross pattern. The Fs-locus mapped on chromosome 5 with SNP marker 'AX-160417325' at 7.4 cM in 'MK_SB' map and between two SNP markers 'AX-160536283' and 'AX-160906995' at a distance of 2.4 and 4.9 cM, respectively in 'MK_D' map. The SNPs 'AX-160417325' and 'AX-160536283' correctly predicted seedlessness of 25-91.9% progenies in this study. Based on the alignment of flanking SNP markers to the Clementine reference genome, the candidate gene for seedlessness hovered in a ~ 6.0 Mb region between 3.97 Mb (AX-160906995) to 10.00 Mb (AX-160536283). This region has 131 genes of which 13 genes (belonging to seven gene families) reportedly express in seed coat or developing embryo. The findings of the study will prove helpful in directing future research for fine mapping this region and eventually underpinning the exact causative gene governing seedlessness in 'MK'.

Integrated microRNA and transcriptome profiling reveals the regulatory network of embryo abortion in jujube

Hybridization is an important approach to the production of new varieties with exceptional traits. Although the kernel rate of wild jujube (Ziziphus jujuba Mill. var. spinosa Hu.) is generally high, that of cultivated jujube (Z. jujuba Mill.) is low, greatly hampering the jujube breeding process. However, the mechanism by which this trait changed during jujube domestication remains unclear. Here, we explored the potential regulatory network that governs jujube embryo abortion using correlation analysis of population traits, artificial pollination, sugar content measurements and multi-omics analysis. The results showed that embryo abortion was an important reason for the low kernel rate of cultivated jujube, and kernel rate was negatively correlated with edible rate. Twenty-one days after pollination was a critical period for embryo abortion. At this time, the sugar content of cultivated 'Junzao' kernels decreased significantly compared with that of the pulp, but sugar content remained relatively stable in kernels of wild 'Suanzao'. A total of 1142 differentially expressed genes targeted by 93 microRNAs (miRNAs) were identified by transcriptome, miRNA and degradome sequencing, and may be involved in the regulation of embryo abortion during kernel development. Among them, DELLA protein, TCP14 and bHLH93 transcription factors have been shown to participate in the regulation of embryonic development. Our findings suggest that carbohydrate flow between different tissues of cultivated jujube exhibits a bias toward the pulp at 21 days after pollination, thereby restricting the process of kernel development. This information enhances our understanding of the embryo abortion process and reveals miRNA-target gene pairs that may be useful for molecular-assisted breeding.

RsCLE22a regulates taproot growth through an auxin signaling-related pathway in radish (Raphanus sativus L.)

CLAVATA3/EMBRYO SURROUNDING REGION-related (CLE) peptides are a class of small molecules involved in plant growth and development. Although radish (Raphanus sativus) is an important root vegetable crop worldwide, the functions of CLE peptides in its taproot formation remain elusive. Here, a total of 48 RsCLE genes were identified from the radish genome. RNA in situ hybridization showed that RsCLE22a gene was highly expressed in the vascular cambium. Overexpression of RsCLE22a inhibited root growth by impairing stem cell proliferation in Arabidopsis, and radish plants with exogenous supplementation of RsCLE22 peptide (CLE22p) showed a similar phenotype. The vascular cambial activity was increased in RsCLE22a-silenced plants. Transcriptome analysis revealed that CLE22p altered the expression of several genes involved in meristem development and hormone signal transduction in radish. Immunolocalization results showed that CLE22p increased auxin accumulation in vascular cambium. Yeast one-hybrid and dual-luciferase assays showed that the WUSCHEL-RELATED HOMEOBOX 4 (RsWOX4) binds to RsCLE22a promoter and activates its transcription. The expression level of RsWOX4 was related to vascular cambial activity and was regulated by auxin. Furthermore, a RsCLE22a-RsWOX4 module is proposed to regulate taproot vascular cambium activity through an auxin signaling-related pathway in radish. These findings provide novel insights into the regulation of root growth in a horticultural crop.

Genome-Wide Identification and Expression Profile of the HD-Zip Transcription Factor Family Associated with Seed Germination and Abiotic Stress Response in Miscanthus sinensis

Miscanthus sinensis is an ornamental grass, non-food bioenergy crop, and forage with high feeding value. It can adapt to many kinds of soil conditions due to its high level of resistance to various abiotic stresses. However, a low level of seed germination restricts the utilization and application of M. sinensis. It is reported that the Homeodomain-leucine zipper (HD-Zip) gene family participates in plant growth and development and ability to cope with outside environment stresses, which may potentially regulate seed germination and stress resistance in M. sinensis. In this study, a complete overview of M. sinensis HD-Zip genes was conducted, including gene structure, conserved motifs, chromosomal distribution, and gene duplication patterns. A total of 169 members were identified, and the HD-Zip proteins were divided into four subgroups. Inter-chromosomal evolutionary analysis revealed that four pairs of tandem duplicate genes and 72 segmental duplications were detected, suggesting the possible role of gene replication events in the amplification of the M. sinensis HD-Zip gene family. There was an uneven distribution of HD-Zip genes on 19 chromosomes of M. sinensis. Also, evolutionary analysis showed that M. sinensis HD-Zip gene family members had more collinearity with monocotyledons and less with dicotyledons. The gene structure analysis showed that there were 93.5% of proteins with motif 1 and motif 4, while motif 10 was only found in group IV. Based on the cis-elements analysis, it appeared that most of the genes were related to plant growth and development, various hormones, and abiotic stress. Furthermore, qRT-PCR analysis showed that Misin06G303300.1 was significantly expressed in seed germination and Misin05G030000.1 and Misin06G303300.1 were highly expressed under chromium, salt, and drought stress. Results in this study will provide a basis for further exploring the potential role of HD-Zip genes in stress responses and genetic improvement of M. sinensis seed germination.

Genome-Wide Analysis of KNOX Transcription Factors and Expression Pattern of Dwarf-Related KNOX Genes in Pear

KNOTTED1-like homeobox (KNOX) transcription factors (TFs) belonging to the homeobox TF family play important roles in plant growth, development, and responses to abiotic and biotic stress. However, little information is available on KNOX TF in pear (Pyrus). In this study, 19 PbKNOXs TFs were re-identified in pear (Pyrus bretschneideri Rehd.). Phylogenetic analysis revealed that the TFs were clustered into three groups with 10 conserved motifs, some of which were group- or subgroup-specific, implying that they are important for the functions of the KNOX in these clades. PbKNM1 and PbKNM2 are KNM (encodes a MEINOX domain but not a homeodomain) genes identified in pear for the first time. KNOX genes in Pyrus and Malus were closely related, and a collinear relationship among PbKNOX genes in Pyrus and Malus was observed. Analysis of the expression patterns of PbKNOX genes in different tissues, at various growth stages, and in response to abiotic and biotic stress revealed that PbKNOXs are involved in plant growth and development. Our comparative transcriptional analysis of dwarf mutant varieties revealed that genes belonging to class I are highly expressed compared with genes in other classes. Analysis of the expression of PbKNOX genes in the hybrid offspring of vigorous and dwarf varieties revealed that PbKNOX genes were highly expressed in the vigorous offspring and weakly expressed in the dwarf offspring. These findings provide new insight into the function of KNOX TFs in pear and will aid future studies of dwarf fruit trees.

Genome-wide analysis of poplar HD-Zip family and over-expression of PsnHDZ63 confers salt tolerance in transgenic Populus simonii × P.nigra

HD-Zip is a plant-specific HB transcription factor, which participates in plant development and stress response. In this study, we identified 63 poplar HD-Zip transcription factors, which were randomly distributed on 19 chromosomes of poplar. Based on the gene structure and phylogenetic relationship, these members are divided into four groups, which have a variety of collinear evolutionary relationships. They also have rich segmental replication events and experienced strong purification selection. Based on RNA-seq analysis, we profiled the expression pattern of the 63 HD-Zip members under salt stress. Subsequently, we carried out in-depth study on the significantly up-regulated PsnHDZ63 in the stems and leaves. The transgenic Populus simonii × P.nigra plants over-expressing PsnHDZ63 displayed better morphological and physiological indexes than WT under salt stress. In addition, PsnHDZ63 enhanced salt stress tolerance of transgenic lines by combining effective stress-resistant elements to improve reactive oxygen species scavenging ability. These studies laid a foundation for a comprehensive understanding of poplar HD-Zip family members, and revealed the important role of PsnHDZ63 in plant salt tolerance.

Elevated transcription of transposable elements is accompanied by het-siRNA-driven de novo DNA methylation in grapevine embryogenic callus

Background

Somatic variation is a valuable source of trait diversity in clonally propagated crops. In grapevine, which has been clonally propagated worldwide for centuries, important phenotypes such as white berry colour are the result of genetic changes caused by transposable elements. Additionally, epiallele formation may play a role in determining geo-specific (‘terroir’) differences in grapes and thus ultimately in wine. This genomic plasticity might be co-opted for crop improvement via somatic embryogenesis, but that depends on a species-specific understanding of the epigenetic regulation of transposable element (TE) expression and silencing in these cultures. For this reason, we used whole-genome bisulphite sequencing, mRNA sequencing and small RNA sequencing to study the epigenetic status and expression of TEs in embryogenic callus, in comparison with leaf tissue.

Results

We found that compared with leaf tissue, grapevine embryogenic callus cultures accumulate relatively high genome-wide CHH methylation, particularly across heterochromatic regions. This de novo methylation is associated with an abundance of transcripts from highly replicated TE families, as well as corresponding 24 nt heterochromatic siRNAs. Methylation in the TE-specific CHG context was relatively low over TEs located within genes, and the expression of TE loci within genes was highly correlated with the expression of those genes.

Conclusions

This multi-‘omics analysis of grapevine embryogenic callus in comparison with leaf tissues reveals a high level of genome-wide transcription of TEs accompanied by RNA-dependent DNA methylation of these sequences in trans. This provides insight into the genomic conditions underlying somaclonal variation and epiallele formation in plants regenerated from embryogenic cultures, which is an important consideration when using these tissues for plant propagation and genetic improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07973-9.

VvHDZ28 positively regulate salicylic acid biosynthesis during seed abortion in Thompson Seedless

Seedlessness in grapes is one of the features most appreciated by consumers. However, the mechanisms underlying seedlessness in grapes remain obscure. Here, we observe small globular embryos and globular embryos in Pinot Noir and Thompson Seedless from 20 to 30 days after flowering (DAF). From 40 to 50 DAF, we observe torpedo embryos and cotyledon embryos in Pinot Noir but aborted embryos and endosperm in Thompson Seedless. Thus, RNA-Seq analyses of seeds at these stages from Thompson Seedless and Pinot Noir were performed. A total of 6442 differentially expressed genes were identified. Among these, genes involved in SA biosynthesis, VvEDS1 and VvSARD1, were more highly expressed in Thompson Seedless than in Pinot Noir. Moreover, the content of endogenous SA is at least five times higher in Thompson Seedless than in Pinot Noir. Increased trimethylation of H3K27 of VvEDS1 and VvSARD1 may be correlated with lower SA content in Pinot Noir. We also demonstrate that VvHDZ28 positively regulates the expression of VvEDS1. Moreover, over-expression of VvHDZ28 results in seedless fruit and increased SA contents in Solanum lycopersicum. Our results reveal the potential role of SA and feedback regulation of VvHDZ28 in seedless grapes.

HD-ZIP Gene Family: Potential Roles in Improving Plant Growth and Regulating Stress-Responsive Mechanisms in Plants

Exploring the molecular foundation of the gene-regulatory systems underlying agronomic parameters or/and plant responses to both abiotic and biotic stresses is crucial for crop improvement. Thus, transcription factors, which alone or in combination directly regulated the targeted gene expression levels, are appropriate players for enlightening agronomic parameters through genetic engineering. In this regard, homeodomain leucine zipper (HD-ZIP) genes family concerned with enlightening plant growth and tolerance to environmental stresses are considered key players for crop improvement. This gene family containing HD and LZ domain belongs to the homeobox superfamily. It is further classified into four subfamilies, namely HD-ZIP I, HD-ZIP II, HD-ZIP III, and HD-ZIP IV. The first HD domain-containing gene was discovered in maize cells almost three decades ago. Since then, with advanced technologies, these genes were functionally characterized for their distinct roles in overall plant growth and development under adverse environmental conditions. This review summarized the different functions of HD-ZIP genes in plant growth and physiological-related activities from germination to fruit development. Additionally, the HD-ZIP genes also respond to various abiotic and biotic environmental stimuli by regulating defense response of plants. This review, therefore, highlighted the various significant aspects of this important gene family based on the recent findings. The practical application of HD-ZIP biomolecules in developing bioengineered plants will not only mitigate the negative effects of environmental stresses but also increase the overall production of crop plants.

Genome-wide identification and expression analysis of the B-box transcription factor gene family in grapevine (Vitis vinifera L.)

Background

B-box (BBX) zinc-finger transcription factors play important roles in plant growth, development, and stress response. Although these proteins have been studied in model plants such as Arabidopsis thaliana or Oryza sativa, little is known about the evolutionary history or expression patterns of BBX proteins in grapevine (Vitis vinifera L.).

Results

We identified a total of 25 VviBBX genes in the grapevine genome and named them according to the homology with Arabidopsis. These proteins were classified into five groups on the basis of their phylogenetic relationships, number of B-box domains, and presence or absence of a CCT domain or VP motif. BBX proteins within the same group showed similar exon-intron structures and were unevenly distributed in grapevine chromosomes. Synteny analyses suggested that only segmental duplication events contributed to the expansion of the VviBBX gene family in grapevine. The observed syntenic relationships between some BBX genes from grapevine and Arabidopsis suggest that they evolved from a common ancestor. Transcriptional analyses showed that the grapevine BBX genes were regulated distinctly in response to powdery mildew infection and various phytohormones. Moreover, the expression levels of a subset of BBX genes in ovules were much higher in seedless grapevine cultivars compared with seeded cultivars during ovule development, implying a potential role in seed abortion. Additionally, VviBBX8, VquBBX15a and VquBBX29b were all located in the nucleus and had transcriptional activity except for VquBBX29b.

Conclusions

The results of this study establish the genome-wide analysis of the grapevine BBX family and provide a framework for understanding the biological roles of BBX genes in grapevine.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07479-4.

A moso bamboo transcription factor, Phehdz1, positively regulates the drought stress response of transgenic rice

78 HD-Zip family genes in Phyllostachys edulis were analyzed. Overexpression of Phehdz1 can improve the drought tolerance of transgenic rice and affect its secondary metabolism. Many studies suggested homeodomain-leucine zipper (HD-Zip) transcription factors are important regulators of plant growth and development, signal transduction, and responses to environmental stresses. In this study, 78 moso bamboo (Phyllostachys edulis) HD-Zip genes were investigated and classified into four subfamilies (HD-Zip I-IV). Additionally, Phehdz1 (HD-Zip I gene) was isolated and confirmed to be highly expressed in the roots. A quantitative real-time PCR analysis indicated Phehdz1 expression was significantly induced by drought, high salinity, and abscisic acid (ABA). A transient expression assay proved that Phehdz1 was localized in the nucleus of tobacco cells. Moreover, it could bind to the core region encoded by the H-box sequence (CAATAATTG) in yeast. In response to mannitol treatments, the Phehdz1-overexpressing transgenic rice had a higher germination rate and longer shoots than the wild-type controls. Moreover, Phehdz1-overexpressing rice plants had a higher survival rate as well as higher relative water and proline contents, but a lower malondialdehyde content, than the WT plants after a 30% polyethylene glycol 6000 treatment. Accordingly, the overexpression of Phehdz1 enhances the drought tolerance of transgenic rice. Many of the differentially expressed genes identified by a transcriptome analysis are involved in MAPK signal transduction and the biosynthesis of secondary metabolites. Thus, the overexpression of Phehdz1 enhances the drought stress tolerance of transgenic rice, while also potentially modulating the expression of metabolism-related genes.

Quantitative Trait Loci for Seed Size Variation in Cucurbits - A Review

Cucurbits (Cucurbitaceae family) include many economically important fruit vegetable crops such as watermelon, pumpkin/squash, cucumber, and melon. Seed size (SS) is an important trait in cucurbits breeding, which is controlled by quantitative trait loci (QTL). Recent advances have deciphered several signaling pathways underlying seed size variation in model plants such as Arabidopsis and rice, but little is known on the genetic basis of SS variation in cucurbits. Here we conducted literature review on seed size QTL identified in watermelon, pumpkin/squash, cucumber and melon, and inferred 14, 9 and 13 consensus SS QTL based on their physical positions in respective draft genomes. Among them, four from watermelon (ClSS2.2, ClSS6.1, ClSS6.2, and ClSS8.2), two from cucumber (CsSS4.1 and CsSS5.1), and one from melon (CmSS11.1) were major-effect, stable QTL for seed size and weight. Whole genome sequence alignment revealed that these major-effect QTL were located in syntenic regions across different genomes suggesting possible structural and functional conservation of some important genes for seed size control in cucurbit crops. Annotation of genes in the four watermelon consensus SS QTL regions identified genes that are known to play important roles in seed size control including members of the zinc finger protein and the E3 ubiquitin-protein ligase families. The present work highlights the utility of comparative analysis in understanding the genetic basis of seed size variation, which may help future mapping and cloning of seed size QTL in cucurbits.

Genome-Wide Characterization and Expression Profiling of GASA Genes during Different Stages of Seed Development in Grapevine (Vitis vinifera L.) Predict Their Involvement in Seed Development

Members of the plant-specific GASA (gibberellic acid-stimulated Arabidopsis) gene family have multiple potential roles in plant growth and development, particularly in flower induction and seed development. However, limited information is available about the functions of these genes in fruit plants, particularly in grapes. We identified 14 GASA genes in grapevine (Vitis vinifera L.) and performed comprehensive bioinformatics and expression analyses. In the bioinformatics analysis, the locations of genes on chromosomes, physiochemical properties of proteins, protein structure, and subcellular positions were described. We evaluated GASA proteins in terms of domain structure, exon-intron distribution, motif arrangements, promoter analysis, phylogenetic, and evolutionary history. According to the results, the GASA domain is conserved in all proteins and the proteins are divided into three well-conserved subgroups. Synteny analysis proposed that segmental and tandem duplication have played a role in the expansion of the GASA gene family in grapes, and duplicated gene pairs have negative selection pressure. Most of the proteins were predicted to be in the extracellular region, chloroplasts, and the vacuole. In silico promoter analysis suggested that the GASA genes may influence different hormone signaling pathways and stress-related mechanisms. Additionally, we performed a comparison of the expression between seedless (Thompson seedless) and seeded (Red globe) cultivars in different plant parts, including the ovule during different stages of development. Furthermore, some genes were differentially expressed in different tissues, signifying their role in grapevine growth and development. Several genes (VvGASA2 and 7) showed different expression levels in later phases of seed development in Red globe and Thompson seedless, suggesting their involvement in seed development. Our study presents the first genome-wide identification and expression profiling of grapevine GASA genes and provides the basis for functional characterization of GASA genes in grapes. We surmise that this information may provide new potential resources for the molecular breeding of grapes.