Transcriptome profiling of grapevine seedless segregants during berry development reveals candidate genes associated with berry weight

Dissecting the complex genetic basis of pre- and post-harvest traits in Vitis vinifera L. using genome-wide association studies

Addressing the pressing challenges in agriculture necessitates swift advancements in breeding programs, particularly for perennial crops like grapevines. Moving beyond the traditional biparental quantitative trait loci (QTL) mapping, we conducted a genome-wide association study (GWAS) encompassing 588 Vitis vinifera L. cultivars from a Chilean breeding program, spanning three seasons and testing 13 key yield-related traits. A strong candidate gene, Vitvi11g000454, located on chromosome 11 and related to plant response to biotic and abiotic stresses through jasmonic acid signaling, was associated with berry width and holds potential for enhancing berry size in grape breeding. We also mapped novel QTL associated with post-harvest traits across chromosomes 2, 4, 9, 11, 15, 18, and 19, broadening our grasp on the genetic intricacies dictating fruit post-harvest behavior, including decay, shriveling, and weight loss. Leveraging gene ontology annotations, we drew parallels between traits and scrutinized candidate genes, laying a robust groundwork for future trait-feature identification endeavors in plant breeding. We also highlighted the importance of carefully considering the choice of the response variable in GWAS analyses, as the use of best linear unbiased estimators (BLUEs) corrections in our study may have led to the suppression of some common QTL in grapevine traits. Our results underscore the imperative of pioneering non-destructive evaluation techniques for long-term conservation traits, offering grape breeders and cultivators insights to improve post-harvest table grape quality and minimize waste.

Genetic changes in the genus Vitis and the domestication of vine

The genus Vitis belongs to the Vitaceae family and is divided into two subgenera: Muscadinia and Vitis, the main difference between these subgenera being the number of chromosomes. There are many hypotheses about the origin of the genus, which have been formed with archaeological studies and lately with molecular analyses. Even though there is no consensus on the place of origin, these studies have shown that grapes have been used by man since ancient times, starting later on its domestication. Most studies point to the Near East and Greece as the beginning of domestication, current research suggests it took place in parallel in different sites, but in all cases Vitis vinifera (L.) subsp. sylvestris [Vitis vinifera (L.) subsp. sylvestris (Gmelin) Hagi] seems to be the species chosen by our ancestors to give rise to the now known Vitis vinifera (L.) subsp. vinifera [=sativa (Hegi)= caucasica (Vavilov)]. Its evolution and expansion into other territories followed the formation of new empires and their expansion, and this is where the historical importance of this crop lies. In this process, plants with hermaphrodite flowers were preferentially selected, with firmer, sweeter, larger fruits of different colors, thus favoring the selection of genes associated with these traits, also resulting in a change in seed morphology. Currently, genetic improvement programs have made use of wild species for the introgression of disease resistance genes and tolerance to diverse soil and climate environments. In addition, the mapping of genes of interest, both linked to agronomic and fruit quality traits, has allowed the use of molecular markers for assisted selection. Information on the domestication process and genetic resources help to understand the gene pool available for the development of cultivars that respond to producer and consumer requirements.

Integrative genomic and transcriptomic analyses of a bud sport mutant 'Jinzao Wuhe' with the phenotype of large berries in grapevines

Background

Bud sport mutation occurs frequently in fruit plants and acts as an important approach for grapevine improvement and breeding. 'Jinzao Wuhe' is a bud sport of the elite cultivar 'Himord Seedless' with obviously enlarged organs and berries. To date, the molecular mechanisms underlying berry enlargement caused by bud sport in grapevines remain unclear.

Methods

Whole genome resequencing (WGRS) was performed for two pairs of bud sports and their maternal plants with similar phenotype to identify SNPs, InDels and structural variations (SVs) as well as related genes. Furthermore, transcriptomic sequencing at different developmental stages and weighted gene co-expression network analysis (WGCNA) for 'Jinzao Wuhe' and its maternal plant 'Himord Seedless' were carried out to identify the differentially expressed genes (DEGs), which were subsequently analyzed for Gene Ontology (GO) and function annotation.

Results

In two pairs of enlarged berry bud sports, a total of 1,334 SNPs, 272 InDels and 74 SVs, corresponding to 1,022 target genes related to symbiotic microorganisms, cell death and other processes were identified. Meanwhile, 1,149 DEGs associated with cell wall modification, stress-response and cell killing might be responsible for the phenotypic variation were also determined. As a result, 42 DEGs between 'Himord Seedless' and 'Jinzao Wuhe' harboring genetic variations were further investigated, including pectin esterase, cellulase A, cytochromes P450 (CYP), UDP-glycosyltransferase (UGT), zinc finger protein, auxin response factor (ARF), NAC transcription factor (TF), protein kinase, etc. These candidate genes offer important clues for a better understanding of developmental regulations of berry enlargement in grapevine.

Conclusion

Our results provide candidate genes and valuable information for dissecting the underlying mechanisms of berry development and contribute to future improvement of grapevine cultivars.

Integration of genomics, transcriptomics and metabolomics identifies candidate loci underlying fruit weight in loquat

Fruit weight is an integral part of fruit-quality traits and directly influences commodity values and economic returns of fruit crops. Despite its importance, the molecular mechanisms underlying fruit weight remain understudied, especially for perennial fruit tree crops such as cultivated loquat (Eriobotrya japonica Lindl.). Auxin is known to regulate fruit development, whereas its role and metabolism in fruit development remain obscure in loquat. In this study, we applied a multi-omics approach, integrating whole-genome resequencing-based quantitative trait locus (QTL) mapping with an F1 population, population genomics analysis using germplasm accessions, transcriptome analysis, and metabolic profiling to identify the genomic regions potentially associated with fruit weight in loquat. We identified three major loci associated with fruit weight, supported by both QTL mapping and comparative genomic analysis between small- and big-fruited loquat cultivars. Comparison between two genotypes with contrasting fruit weight performance through transcriptomic and metabolic profiling revealed an important role of auxin in regulating fruit development, especially at the fruit enlarging stage. The multi-omics approach identified two homologs of ETHYLENE INSENSITIVE 4 (EjEIN4) and TORNADO 1 (EjTRN1) as promising candidates controlling fruit weight. Moreover, three single nucleotide polymorphism (SNP) markers were closely associated with fruit weight. Results from this study provided insights from multiple perspectives into the genetic and metabolic controls of fruit weight in loquat. The candidate genomic regions, genes, and sequence variants will facilitate understanding the molecular basis of fruit weight and lay a foundation for future breeding and manipulation of fruit weight in loquat.

Pre-Anthesis Cytokinin Applications Increase Table Grape Berry Firmness by Modulating Cell Wall Polysaccharides

The use of plant growth regulators (PGRs) is widespread in commercial table grape vineyards. The synthetic cytokinin CPPU is a PGR that is extensively used to obtain higher quality grapes. However, the effect of CPPU on berry firmness is not clear. The current study investigated the effects of pre-anthesis applications (BBCH15 and BBCH55 stages) of CPPU on 'Thompson Seedless' berry firmness at harvest through a combination of cytological, morphological, and biochemical analyses. Ovaries in CPPU-treated plants presented morphological changes related to cell division and cell wall modification at the anthesis stage (BBCH65). Moreover, immunofluorescence analysis with monoclonal antibodies 2F4 and LM15 against pectin and xyloglucan demonstrated that CPPU treatment resulted in cell wall modifications at anthesis. These early changes have major repercussions regarding the hemicellulose and pectin cell wall composition of mature fruits, and are associated with increased calcium content and a higher berry firmness at harvest.

Identification of SNPs and InDels associated with berry size in table grapes integrating genetic and transcriptomic approaches

BACKGROUND

Berry size is considered as one of the main selection criteria in table grapes breeding programs, due to the consumer preferences. However, berry size is a complex quantitive trait under polygenic control, and its genetic determination of berry weight is not yet fully understood. The aim of this work was to perform marker discovery using a transcriptomic approach, in order to identify and characterize SNP and InDel markers associated with berry size in table grapes. We used an integrative analysis based on RNA-Seq, SNP/InDel search and validation on table grape segregants and varieties with different genetic backgrounds.

RESULTS

Thirty SNPs and eight InDels were identified using a transcriptomic approach (RNA-Seq). These markers were selected from SNP/InDel found among segregants from a Ruby x Sultanina population with contrasting phenotypes for berry size. The set of 38 SNP and InDel markers was distributed in eight chromosomes. Genotype-phenotype association analyses were performed using a set of 13 RxS segregants and 41 table grapes varieties with different genetic backgrounds during three seasons. The results showed several degrees of association of these markers with berry size (10.2 to 30.7%) as other berry-related traits such as length and width. The co-localization of SNP and /or InDel markers and previously reported QTLs and candidate genes associated with berry size were analysed.

CONCLUSIONS

We identified a set of informative and transferable SNP and InDel markers associated with berry size. Our results suggest the suitability of SNPs and InDels as candidate markers for berry weight in seedless table grape breeding. The identification of genomic regions associated with berry weight in chromosomes 8, 15 and 17 was achieved with supporting evidence derived from a transcriptome experiment focused on SNP/InDel search, as well as from a QTL-linkage mapping approach. New regions possibly associated with berry weight in chromosomes 3, 6, 9 and 14 were identified.

Transcriptome Analysis of Elm (Ulmus pumila) Fruit to Identify Phytonutrients Associated Genes and Pathways

Plant fruit is an important source of natural active phytonutrients that are profitable for human health. Elm (Ulmus pumila) fruit is considered as natural plant food in China that is rich in nutrients. In the present study, high-throughput RNA sequencing was performed in U. pumila edible fruits and leaves and 11,386 unigenes were filtered as dysregulated genes in fruit samples, including 5231 up- and 6155 downregulated genes. Hundreds of pathways were predicted to participate in seed development and phytonutrient biosynthesis in U. pumila by GO, MapMan, and KEGG enrichment analysis, including “seed maturation”, “glycine, serine, and threonine metabolism” and “phenylpropanoid biosynthesis”. ABA-mediated glucose response-related ethylene-activated signaling pathway (e.g., ABI4) were supposed to associate with elm fruit development; unsaturated fatty acids pathway (e.g., ACX2 and SAD) were predicted to participate in determination of fatty acid composition in elm fruit; flavonoid and coumarins biosynthesis (e.g., CYP98A3 and CCoAOMT1) were demonstrated to correlate with the bioactivity of elm fruits in human cancer and inflammation resistance. To provide more information about fruit developmental status, the qRT-PCR analysis for key genes of “phenylpropanoid biosynthesis” and “alpha-Linolenic acid metabolism” were conducted in samples of young fruits, ripe fruit, old fruit, and leaves. Two biosynthetic pathways for unsaturated fatty acid and Jasmonic acid (JA) were deduced to be involved in fruit development in U. pumila and the phenylpropanoid glycoside, syringin, was speculated to accumulate in the early development stages of elm fruit. Our transcriptome data supports molecular clues for seed development and biologically active substances in elm fruits.

Identification of transcription factor genes involved in anthocyanin biosynthesis in carrot (Daucus carota L.) using RNA-Seq

BACKGROUND

Anthocyanins are water-soluble colored flavonoids present in multiple organs of various plant species including flowers, fruits, leaves, stems and roots. DNA-binding R2R3-MYB transcription factors, basic helix-loop-helix (bHLH) transcription factors, and WD40 repeat proteins are known to form MYB-bHLH-WD repeat (MBW) complexes, which activates the transcription of structural genes in the anthocyanin pathway. Although black cultivars of carrots (Daucus carota L.) can accumulate large quantities of anthocyanin in their storage roots, the regulatory genes responsible for their biosynthesis are not well characterized. The current study aimed to analyze global transcription profiles based on RNA sequencing (RNA-Seq), and mine MYB, bHLH and WD40 genes that may function as positive or negative regulators in the carrot anthocyanin biosynthesis pathways.

RESULTS

RNA was isolated from differently colored calli, as well as tissue samples from taproots of various black carrot cultivars across the course of development, and gene expression levels of colored and non-colored tissue and callus samples were compared. The expression of 32 MYB, bHLH and WD40 genes were significantly correlated with anthocyanin content in black carrot taproot. Of those, 11 genes were consistently up- or downregulated in a purple color-specific manner across various calli and cultivar comparisons. The expression of 10 out of these 11 genes was validated using real-time quantitative reverse transcriptase polymerase chain reaction (qRT-PCR).

CONCLUSIONS

The results of this study provide insights into regulatory genes that may be responsible for carrot anthocyanin biosynthesis, and suggest that future focus on them may help improve our overall understanding of the anthocyanin synthesis pathway.

Influence of natural variation in berry size on the volatile profiles of Vitis vinifera L. cv. Merlot and Cabernet Gernischt grapes

This study was conducted during the 2014 and 2015 vintages on Vitis vinifera L. cv. Merlot and Cabernet Gernischt to investigate whether natural variation in berry size could affect grape aromatic compounds. Grape berries were separated into three size categories based on their diameter: small, middle and large. The results showed that berry size exerted a significant influence on the volatile profiles of both winegrape varieties. Hierarchical clustering analysis demonstrated that the volatile profiles of middle berries were different from those of large and small berries. Middle berries had the greatest abundance of aroma compounds, followed by small and large berries. Especially, C6/C9 compounds, norisoprenoids, terpenoids showed markedly different concentrations among differently sized Merlot berries and C6/C9 compounds, terpenoids among differently sized Cabernet Gernischt berries. Middle berries of both grape varieties may possess the greatest intensity of fresh-green, fruity and floral aromas due to the high odour activity values (OAVs) of decanal, hexanal, (E)-2-hexenal, (E)-β-damascenone and β-ionone in middle sizes of Merlot berries and the high OAVs of (E)-2-hexenal and (E)-β-damascenone in middle sizes of Cabernet Gernischt berries. This knowledge could be important for winemakers to conduct targeted berry sorting, thereby improving the aromatic quality of grapes.

Transcriptome analysis of grapevine under salinity and identification of key genes responsible for salt tolerance

The negative effects of soil salinity towards grape yield depend upon salt concentration, cultivar type, developmental stage, and rootstock. Thompson Seedless variety of grape plant is considered moderately sensitive to salinity when grown upon its own root stock. In recent epoch, identification of key genes responsive to salinity offers hope to generate salinity-tolerant crop plants by their overexpression through genetic manipulation. In the present report, salt responsive transcriptome analysis of Thompson Seedless grape variety was done to identify vital genes involved in salinity tolerance which could be used further to generate salt liberal grape plant or other crop plants. Transcriptome libraries for control and 150-mM-NaCl-treated grape leaves were sequenced on Illumina platform where 714 genes were found to be differentially expressed. Gene ontology analysis indicated that under salinity conditions, the genes involved in metabolic process were highly enriched. Keto Encyclopedia of Genes and Genomes analysis revealed that, among the top 22 enriched pathways for the salt stress upregulated genes, the carbohydrate metabolism, signal transduction, energy metabolism, amino acid metabolism, biosynthesis of secondary metabolite, and lipid metabolism pathways possessed the largest number of transcripts. Key salinity-induced genes were selected and validated through qRT-PCR analysis which was comparable to RNA-seq results. Real-time PCR analysis also revealed that after 24 days of salinity, the expression of most of the selected key genes was highest. These salinity-induced genes will be characterized further in a model plant and also in Vitis vinifera through transgenic approach to disclose their role towards salt tolerance.

Inter-varietal structural variation in grapevine genomes

Grapevine (Vitis vinifera L.) is one of the world's most important crop plants, which is of large economic value for fruit and wine production. There is much interest in identifying genomic variations and their functional effects on inter-varietal, phenotypic differences. Using an approach developed for the analysis of human and mammalian genomes, which combines high-throughput sequencing, array comparative genomic hybridization, fluorescent in situ hybridization and quantitative PCR, we created an inter-varietal atlas of structural variations and single nucleotide variants (SNVs) for the grapevine genome analyzing four economically and genetically relevant table grapevine varieties. We found 4.8 million SNVs and detected 8% of the grapevine genome to be affected by genomic variations. We identified more than 700 copy number variation (CNV) regions and more than 2000 genes subjected to CNV as potential candidates for phenotypic differences between varieties.