Heritability and identification of QTLs and underlying candidate genes associated with the architecture of the grapevine cluster (Vitis vinifera L.)

Genetic analysis of grapevine root system architecture and loci associated gene networks

Own-rooted grapevines and grapevine rootstocks are vegetatively propagated from cuttings and have an adventitious root system. Unraveling the genetic underpinnings of the adventitious root system architecture (RSA) is important for improving own-rooted and grafted grapevine sustainability for a changing climate. Grapevine RSA genetic analysis was conducted in an Vitis sp. 'VRS-F2' population. Nine root morphology, three total root system morphology, and two biomass traits that contribute to root anchorage and water and nutrient uptake were phenotyped. Quantitative trait loci (QTL) analysis was performed using a high density integrated GBS and rhAmpSeq genetic map. Thirty-one QTL were detected for eleven of the RSA traits (surface area, root volume, total root length, fresh weight, number of tips, forks or links, longest root and average root diameter, link length, and link surface area) revealing many small effects. Several QTL were colocated on chromosomes 1, 9, 13, 18, and 19. QTL with identical peak positions on chromosomes 1 or 13 were enriched for AP2-EREBP, AS2, C2C2-CO, HMG, and MYB transcription factors, and QTL on chromosomes 9 or 13 were enriched for the ALFIN-LIKE transcription factor and regulation of autophagy pathways. QTL modeling for individual root traits identified eight models explaining 13.2 to 31.8% of the phenotypic variation. 'Seyval blanc' was the grandparent contributing to the allele models that included a greater surface area, total root length, and branching (number of forks and links) traits promoting a greater root density. In contrast, V. riparia 'Manitoba 37' contributed the allele for greater average branch length (link length) and diameter, promoting a less dense elongated root system with thicker roots. LATERAL ORGAN BOUNDARY DOMAIN (LBD or AS2/LOB) and the PROTODERMAL FACTOR (PFD2 and ANL2) were identified as important candidate genes in the enriched pathways underlying the hotspots for grapevine adventitious RSA. The combined QTL hotspot and trait modeling identified transcription factors, cell cycle and circadian rhythm genes with a known role in root cell and epidermal layer differentiation, lateral root development and cortex thickness. These genes are candidates for tailoring grapevine root system texture, density and length in breeding programs.

Review: Status and prospects of association mapping in grapevine

Thanks to current advances in sequencing technologies, novel bioinformatics tools, and efficient modeling solutions, association mapping has become a widely accepted approach to unravel the link between genotype and phenotype diversity in numerous crops. In grapevine, this strategy has been used in the last decades to understand the genetic basis of traits of agronomic interest (fruit quality, crop yield, biotic and abiotic resistance), of special relevance nowadays to improve crop resilience to cope with future climate scenarios. Genome-wide association studies have identified many putative causative loci for different traits, some of them overlapping well-known causal genes identified by conventional quantitative trait loci studies in biparental progenies, and/or validated by functional approaches. In addition, candidate-gene association studies have been useful to pinpoint the causal mutation underlying phenotypic variation for several traits of high interest in breeding programs (like berry color, seedlessness, and muscat flavor), information that has been used to develop highly informative and useful markers already in use in marker-assisted selection processes. Thus, association mapping has proved to represent a valuable step towards high quality and sustainable grape production. This review summarizes current applications of association mapping in grapevine research and discusses future prospects in view of current viticulture challenges.

Perspectives and recent progress of genome-wide association studies (GWAS) in fruits

BACKGROUND

Earlier next-generation sequencing technologies are being vastly used to explore, administer, and investigate the gene space with accurate profiling of nucleotide variations in the germplasm. OVERVIEW AND PROGRESS: Recently, novel advancements in high-throughput sequencing technologies allow a genotyping-by-sequencing approach that has opened up new horizons for extensive genotyping exploiting single-nucleotide-polymorphisms (SNPs). This method acts as a bridge to support and minimize a genotype to phenotype gap allowing genetic selection at the genome-wide level, named genomic selection that could facilitate the selection of traits also in the pomology sector. In addition to this, genome-wide genotyping is a prerequisite for genome-wide association studies that have been used successfully to discover the genes, which control polygenic traits including the genetic loci, associated with the trait of interest in fruit crops.

AIMS AND PROSPECTS

This review article emphasizes the role of genome-wide approaches to unlock and explore the genetic potential along with the detection of SNPs affecting the phenotype of fruit crops and highlights the prospects of genome-wide association studies in fruits.

Genetic variation and association analyses identify genes linked to fruit set-related traits in grapevine

Grapevine is one of the most valuable fruit crops in the world. Adverse environmental conditions reduce fruit quality and crop yield, so understanding the genetic and molecular mechanisms determining crop yield components is essential to optimize grape production. The analysis of a diverse collection of grapevine cultivars allowed us to evaluate the relationship between fruit set-related components of yield, including the incidence of reproductive disorders such as coulure and millerandage. The collection displayed a great phenotypic variation that we surveyed in a genetics association study using 15,309 single nucleotide polymorphisms (SNPs) detected in the sequence of 289 candidate genes scattered across the 19 grapevine linkage groups. After correcting statistical models for population structure and linkage disequilibrium effects, 164 SNPs from 34 of these genes were found to associate with fruit set-related traits, supporting a complex polygenic determinism. Many of them were found in the sequence of different putative MADS-box transcription factors, a gene family related with plant reproductive development control. In addition, we observed an additive effect of some of the associated SNPs on the phenotype, suggesting that advantageous alleles from different loci could be pyramided to generate superior cultivars with optimized fruit production.

High-Density Genetic Linkage Map Construction and White Rot Resistance Quantitative Trait Loci Mapping for Genus Vitis Based on Restriction Site-Associated DNA Sequencing

Grape white rot (Coniothyrium diplodiella) is a major fungal disease affecting grape yield and quality. Quantitative trait locus (QTL) analysis is an important method for studying important horticultural traits of grapevine. This study was conducted to construct a high-density map and conduct QTL mapping for grapevine white rot resistance. A mapping population with 177 genotypes was developed from interspecific hybridization of a white rot-resistant cultivar (Vitis vinifera × V. labrusca 'Zhuosexiang') and white rot-susceptible cultivar (V. vinifera 'Victoria'). Single-nucleotide polymorphism (SNP) markers were developed by restriction site-associated DNA sequencing. The female, male, and integrated maps contained 2,501, 4,110, and 6,249 SNP markers with average genetic distances of adjacent markers of 1.25, 0.77, and 0.50 cM, respectively. QTL mapping was conducted based on white rot resistance identification of 177 individuals in July and August of 2017 and 2018. Notably, one stable QTL related to white rot resistance was detected and located on linkage group LG14. The phenotypic variance ranged from 12.93 to 13.43%. An SNP marker (chr14_3929380), which cosegregated with white rot resistance, was discovered and shows potential for use in marker-assisted selection to generate new grapevine cultivars with resistance to white rot.

Somatic variants for seed and fruit set in grapevine

BACKGROUND

Grapevine reproductive development has direct implications on yield. It also impacts on berry and wine quality by affecting traits like seedlessness, berry and bunch size, cluster compactness and berry skin to pulp ratio. Seasonal fluctuations in yield, fruit composition and wine attributes, which are largely driven by climatic factors, are major challenges for worldwide table grape and wine industry. Accordingly, a better understanding of reproductive processes such as gamete development, fertilization, seed and fruit set is of paramount relevance for managing yield and quality. With the aim of providing new insights into this field, we searched for clones with contrasting seed content in two germplasm collections.

RESULTS

We identified eight variant pairs that seemingly differ only in seed-related characteristics while showing identical genotype when tested with the GrapeReSeq_Illumina_20K_SNP_chip and several microsatellites. We performed multi-year observations on seed and fruit set deriving from different pollination treatments, with special emphasis on the pair composed by Sangiovese and its seedless variant locally named Corinto Nero. The pollen of Corinto Nero failed to germinate in vitro and gave poor berry set when used to pollinate other varieties. Most berries from both open- and cross-pollinated Corinto Nero inflorescences did not contain seeds. The genetic analysis of seedlings derived from occasional Corinto Nero normal seeds revealed that the few Corinto Nero functional gametes are mostly unreduced. Moreover, three genotypes, including Sangiovese and Corinto Nero, were unexpectedly found to develop fruits without pollen contribution and occasionally showed normal-like seeds. Five missense single nucleotide polymorphisms were identified between Corinto Nero and Sangiovese from transcriptomic data.

CONCLUSIONS

Our observations allowed us to attribute a seedlessness type to some variants for which it was not documented in the literature. Interestingly, the VvAGL11 mutation responsible for Sultanina stenospermocarpy was also discovered in a seedless mutant of Gouais Blanc. We suggest that Corinto Nero parthenocarpy is driven by pollen and/or embryo sac defects, and both events likely arise from meiotic anomalies. The single nucleotide polymorphisms identified between Sangiovese and Corinto Nero are suitable for testing as traceability markers for propagated material and as functional candidates for the seedless phenotype.

Molecular Mapping of Water-Stress Responsive Genomic Loci in Lettuce (Lactuca spp.) Using Kinetics Chlorophyll Fluorescence, Hyperspectral Imaging and Machine Learning

Deep understanding of genetic architecture of water-stress tolerance is critical for efficient and optimal development of water-stress tolerant cultivars, which is the most economical and environmentally sound approach to maintain lettuce production with limited irrigation. Lettuce (Lactuca sativa L.) production in areas with limited precipitation relies heavily on the use of ground water for irrigation. Lettuce plants are highly susceptible to water-stress, which also affects their nutrient uptake efficiency. Water stressed plants show reduced growth, lower biomass, and early bolting and flowering resulting in bitter flavors. Traditional phenotyping methods to evaluate water-stress are labor intensive, time-consuming and prone to errors. High throughput phenotyping platforms using kinetic chlorophyll fluorescence and hyperspectral imaging can effectively attain physiological traits related to photosynthesis and secondary metabolites that can enhance breeding efficiency for water-stress tolerance. Kinetic chlorophyll fluorescence and hyperspectral imaging along with traditional horticultural traits identified genomic loci affected by water-stress. Supervised machine learning models were evaluated for their accuracy to distinguish water-stressed plants and to identify the most important water-stress related parameters in lettuce. Random Forest (RF) had classification accuracy of 89.7% using kinetic chlorophyll fluorescence parameters and Neural Network (NN) had classification accuracy of 89.8% using hyperspectral imaging derived vegetation indices. The top ten chlorophyll fluorescence parameters and vegetation indices selected by sequential forward selection by RF and NN were genetically mapped using a L. sativa × L. serriola interspecific recombinant inbred line (RIL) population. A total of 25 quantitative trait loci (QTL) segregating for water-stress related horticultural traits, 26 QTL for the chlorophyll fluorescence traits and 34 QTL for spectral vegetation indices (VI) were identified. The percent phenotypic variation (PV) explained by the horticultural QTL ranged from 6.41 to 19.5%, PV explained by chlorophyll fluorescence QTL ranged from 6.93 to 13.26% while the PV explained by the VI QTL ranged from 7.2 to 17.19%. Eight QTL clusters harboring co-localized QTL for horticultural traits, chlorophyll fluorescence parameters and VI were identified on six lettuce chromosomes. Molecular markers linked to the mapped QTL clusters can be targeted for marker-assisted selection to develop water-stress tolerant lettuce.

Towards Sensor-Based Phenotyping of Physical Barriers of Grapes to Improve Resilience to Botrytis Bunch Rot

Botrytis bunch rot is one of the economically most important fungal diseases in viticulture (aside from powdery mildew and downy mildew). So far, no active defense mechanisms and resistance loci against the necrotrophic pathogen are known. Since long, breeders are mostly selecting phenotypically for loose grape bunches, which is recently the most evident trait to decrease the infection risk of Botrytis bunch rot. This study focused on plant phenomics of multiple traits by applying fast sensor technologies to measure berry impedance (Z _REL ), berry texture, and 3D bunch architecture. As references, microscopic determined cuticle thickness (MS _CT ) and infestation of grapes with Botrytis bunch rot were used. Z _REL hereby is correlated to grape bunch density OIV204 (r = -0.6), cuticle thickness of berries (r = 0.61), mean berry diameter (r = -0.63), and Botrytis bunch rot (r = -0.7). However, no correlation between Z _REL and berry maturity or berry texture was observed. In comparison to the category of traditional varieties (mostly susceptible), elite breeding lines show an impressive increased Z _REL value (+317) and a 1-μm thicker berry cuticle. Quantitative trait loci (QTLs) on LGs 2, 6, 11, 15, and 16 were identified for Z _REL and berry texture explaining a phenotypic variance of between 3 and 10.9%. These QTLs providing a starting point for the development of molecular markers. Modeling of Z _REL and berry texture to predict Botrytis bunch rot resilience revealed McFadden R ² = 0.99. Taken together, this study shows that in addition to loose grape bunch architecture, berry diameter, Z _REL , and berry texture values are probably additional parameters that could be used to identify and select Botrytis-resilient wine grape varieties. Furthermore, grapevine breeding will benefit from these reliable methodologies permitting high-throughput screening for additional resilience traits of mechanical and physical barriers to Botrytis bunch rot. The findings might also be applicable to table grapes and other fruit crops like tomato or blueberry.

Differential expression of transcription factor- and further growth-related genes correlates with contrasting cluster architecture in Vitis vinifera 'Pinot Noir' and Vitis spp. genotypes

Grapevine (Vitis vinifera L.) is an economically important crop that needs to comply with high quality standards for fruit, juice and wine production. Intense plant protection is required to avoid fungal damage. Grapevine cultivars with loose cluster architecture enable reducing protective treatments due to their enhanced resilience against fungal infections, such as Botrytis cinerea-induced gray mold. A recent study identified transcription factor gene VvGRF4 as determinant of pedicel length, an important component of cluster architecture, in samples of two loose and two compact quasi-isogenic 'Pinot Noir' clones. Here, we extended the analysis to 12 differently clustered 'Pinot Noir' clones from five diverse clonal selection programs. Differential gene expression of these clones was studied in three different locations over three seasons. Two phenotypically opposite clones were grown at all three locations and served for standardization. Data were correlated with the phenotypic variation of cluster architecture sub-traits. A set of 14 genes with consistent expression differences between loosely and compactly clustered clones-independent from season and location-was newly identified. These genes have annotations related to cellular growth, cell division and auxin metabolism and include two more transcription factor genes, PRE6 and SEP1-like. The differential expression of VvGRF4 in relation to loose clusters was exclusively found in 'Pinot Noir' clones. Gene expression studies were further broadened to phenotypically contrasting F1 individuals of an interspecific cross and OIV reference varieties of loose cluster architecture. This investigation confirmed PRE6 and six growth-related genes to show differential expression related to cluster architecture over genetically divergent backgrounds.

High-density genetic linkage map construction and cane cold hardiness QTL mapping for Vitis based on restriction site-associated DNA sequencing

Background

Cold hardiness is an important agronomic trait and can significantly affect grape production and quality. Until now, there are no reports focusing on cold hardiness quantitative trait loci (QTL) mapping. In this study, grapevine interspecific hybridisation was carried out with the maternal parent ‘Cabernet sauvignon’ and paternal parent ‘Zuoyouhong’. A total of 181 hybrid offspring and their parents were used as samples for restriction-site associated DNA sequencing (RAD). Grapevine cane phloem and xylem cold hardiness of the experimental material was detected using the low-temperature exotherm method in 2016, 2017 and 2018. QTL mapping was then conducted based on the integrated map.

Results

We constructed a high-density genetic linkage map with 16,076, 11,643, and 25,917 single-nucleotide polymorphism (SNP) markers anchored in the maternal, paternal, and integrated maps, respectively. The average genetic distances of adjacent markers in the maps were 0.65 cM, 0.77 cM, and 0.41 cM, respectively. Colinearity analysis was conducted by comparison with the grape reference genome and showed good performance. Six QTLs were identified based on the phenotypic data of 3 years and they were mapped on linkage group (LG) 2, LG3, and LG15. Based on QTL results, candidate genes which may be involved in grapevine cold hardiness were selected.

Conclusions

High-density linkage maps can facilitate grapevine fine QTL mapping, genome comparison, and sequence assembly. The cold hardiness QTL mapping and candidate gene discovery performed in this study provide an important reference for molecular-assisted selection in grapevine cold hardiness breeding.

VviUCC1 Nucleotide Diversity, Linkage Disequilibrium and Association with Rachis Architecture Traits in Grapevine

Cluster compactness is a trait with high agronomic relevance, affecting crop yield and grape composition. Rachis architecture is a major component of cluster compactness determinism, and is a target trait toward the breeding of grapevine varieties less susceptible to pests and diseases. Although its genetic basis is scarcely understood, a preliminary result indicated a possible involvement of the VviUCC1 gene. The aim of this study was to characterize the VviUCC1 gene in grapevine and to test the association between the natural variation observed for a series of rachis architecture traits and the polymorphisms detected in the VviUCC1 sequence. This gene encodes an uclacyanin plant-specific cell-wall protein involved in fiber formation and/or lignification processes. A high nucleotide diversity in the VviUCC1 gene promoter and coding regions was observed, but no critical effects were predicted in the protein domains, indicating a high level of conservation of its function in the cultivated grapevine. After correcting statistical models for genetic stratification and linkage disequilibrium effects, marker-trait association results revealed a series of single nucleotide polymorphisms (SNPs) significantly associated with cluster compactness and rachis traits variation. Two of them (Y-984 and K-88) affected two common cis-transcriptional regulatory elements, suggesting an effect on phenotype via gene expression regulation. This work reinforces the interest of further studies aiming to reveal the functional effect of the detected VviUCC1 variants on grapevine rachis architecture.

Characterizing 3D inflorescence architecture in grapevine using X-ray imaging and advanced morphometrics: implications for understanding cluster density

Inflorescence architecture provides the scaffold on which flowers and fruits develop, and consequently is a primary trait under investigation in many crop systems. Yet the challenge remains to analyse these complex 3D branching structures with appropriate tools. High information content datasets are required to represent the actual structure and facilitate full analysis of both the geometric and the topological features relevant to phenotypic variation in order to clarify evolutionary and developmental inflorescence patterns. We combined advanced imaging (X-ray tomography) and computational approaches (topological and geometric data analysis and structural simulations) to comprehensively characterize grapevine inflorescence architecture (the rachis and all branches without berries) among 10 wild Vitis species. Clustering and correlation analyses revealed unexpected relationships, for example pedicel branch angles were largely independent of other traits. We identified multivariate traits that typified species, which allowed us to classify species with 78.3% accuracy, versus 10% by chance. Twelve traits had strong signals across phylogenetic clades, providing insight into the evolution of inflorescence architecture. We provide an advanced framework to quantify 3D inflorescence and other branched plant structures that can be used to tease apart subtle, heritable features for a better understanding of genetic and environmental effects on plant phenotypes.

Selection of candidate genes controlling veraison time in grapevine through integration of meta-QTL and transcriptomic data

BACKGROUND

High temperature during grape berry ripening impairs the quality of fruits and wines. Veraison time, which marks ripening onset, is a key factor for determining climatic conditions during berry ripening. Understanding its genetic control is crucial to successfully breed varieties more adapted to a changing climate. Quantitative trait loci (QTL) studies attempting to elucidate the genetic determinism of developmental stages in grapevine have identified wide genomic regions. Broad scale transcriptomic studies, by identifying sets of genes modulated during berry development and ripening, also highlighted a huge number of putative candidates.

RESULTS

With the final aim of providing an overview about available information on the genetic control of grapevine veraison time, and prioritizing candidates, we applied a meta-QTL analysis for grapevine phenology-related traits and checked for co-localization of transcriptomic candidates. A consensus genetic map including 3130 markers anchored to the grapevine genome assembly was compiled starting from 39 genetic maps. Two thousand ninety-three QTLs from 47 QTL studies were projected onto the consensus map, providing a comprehensive overview about distribution of available QTLs and revealing extensive co-localization especially across phenology related traits. From 141 phenology related QTLs we generated 4 veraison meta-QTLs located on linkage group (LG) 1 and 2, and 13 additional meta-QTLs connected to the veraison time genetic control, among which the most relevant were located on LG 14, 16 and 18. Functional candidates in these intervals were inspected. Lastly, taking advantage of available transcriptomic datasets, expression data along berry development were integrated, in order to pinpoint among positional candidates, those differentially expressed across the veraison transition.

CONCLUSION

Integration of meta-QTLs analysis on available phenology related QTLs and data from transcriptomic dataset allowed to strongly reduce the number of candidate genes for the genetic control of the veraison transition, prioritizing a list of 272 genes, among which 78 involved in regulation of gene expression, signal transduction or development.

A QTL associated with leaf trichome traits has a major influence on the abundance of the predatory mite Typhlodromus pyri in a hybrid grapevine population

The abundance of predatory phytoseiid mites, Typhlodromus pyri, important biological control agents of spider mite pests in numerous crops, is positively influenced by the density of leaf trichomes and tuft-form domatia in vein axils. Identification of the genetic regions controlling both trophic levels could facilitate the improvement of predatory mite habitat in breeding programs. The abundance of T. pyri and non-glandular trichomes was measured in a segregating F₁ family derived from the cross of the complex Vitis hybrid, 'Horizon', with Illinois 547-1 (V. rupestris B38 × V. cinerea B9), finding positive correlation among traits. High density genetic maps were used to localize one major quantitative trait locus (QTL) on chromosome 1 of Illinois 547-1 associated with both predatory mite abundance and leaf trichomes. This QTL explained 23% of the variation in phytoseiid abundance and similar amounts of variance in domatia rating (21%), domatia size (16%), leaf bristle density (37% in veins and 33% in blades), and leaf hair density (20% in veins and 15% in blades). Another nine QTL distributed among chromosomes 1, 2, 5, 8, and 15 were associated solely with trichome density, and explained 7-17% of the phenotypic variation. Combined, our results provide evidence of the genetic architecture of non-glandular trichomes in Vitis, with a major locus influencing trichome densities, domatia size and predatory mite abundance. This information is relevant for breeding grapevines with a more favorable habitat for biological control agents.

Identification of co-located QTLs and genomic regions affecting grapevine cluster architecture

Loose cluster architecture is an important aim in grapevine breeding since it has high impact on the phytosanitary status of grapes. This investigation analyzed the contributions of individual cluster sub-traits to the overall trait of cluster architecture. Six sub-traits showed large impact on cluster architecture as major determinants. They explained 57% of the OIV204 descriptor for cluster compactness rating in a highly diverse cross-population of 149 genotypes. Genetic analysis revealed several genomic regions involved in the expression of this trait. Based on the linkage of phenotypic features to molecular markers, QTL calculations shed new light on the genetic determinants of cluster architecture. Eight QTL clusters harbor overlapping confidence intervals of up to four co-located QTLs. A physical projection of the QTL clusters by confidence interval-flanking markers onto the PN40024 reference genome sequence revealed genes enriched in these regions.

Prioritization of Candidate Genes in QTL Regions for Physiological and Biochemical Traits Underlying Drought Response in Barley (Hordeum vulgare L.)

Drought is one of the most adverse abiotic factors limiting growth and productivity of crops. Among them is barley, ranked fourth cereal worldwide in terms of harvested acreage and production. Plants have evolved various mechanisms to cope with water deficit at different biological levels, but there is an enormous challenge to decipher genes responsible for particular complex phenotypic traits, in order to develop drought tolerant crops. This work presents a comprehensive approach for elucidation of molecular mechanisms of drought tolerance in barley at the seedling stage of development. The study includes mapping of QTLs for physiological and biochemical traits associated with drought tolerance on a high-density function map, projection of QTL confidence intervals on barley physical map, and the retrievement of positional candidate genes (CGs), followed by their prioritization based on Gene Ontology (GO) enrichment analysis. A total of 64 QTLs for 25 physiological and biochemical traits that describe plant water status, photosynthetic efficiency, osmoprotectant and hormone content, as well as antioxidant activity, were positioned on a consensus map, constructed using RIL populations developed from the crosses between European and Syrian genotypes. The map contained a total of 875 SNP, SSR and CGs, spanning 941.86 cM with resolution of 1.1 cM. For the first time, QTLs for ethylene, glucose, sucrose, maltose, raffinose, α-tocopherol, γ-tocotrienol content, and catalase activity, have been mapped in barley. Based on overlapping confidence intervals of QTLs, 11 hotspots were identified that enclosed more than 60% of mapped QTLs. Genetic and physical map integration allowed the identification of 1,101 positional CGs within the confidence intervals of drought response-specific QTLs. Prioritization resulted in the designation of 143 CGs, among them were genes encoding antioxidants, carboxylic acid biosynthesis enzymes, heat shock proteins, small auxin up-regulated RNAs, nitric oxide synthase, ATP sulfurylases, and proteins involved in regulation of flowering time. This global approach may be proposed for identification of new CGs that underlies QTLs responsible for complex traits.

GA3 application in grapes (Vitis vinifera L.) modulates different sets of genes at cluster emergence, full bloom, and berry stage as revealed by RNA sequence-based transcriptome analysis

In grapes (Vitis vinifera L.), exogenous gibberellic acid (GA₃) is applied at different stages of bunch development to achieve desirable bunch shape and berry size in seedless grapes used for table purpose. RNA sequence-based transcriptome analysis was used to understand the mechanism of GA₃ action at cluster emergence, full bloom, and berry stage in table grape variety Thompson Seedless. At cluster emergence, rachis samples were collected at 6 and 24 h after application of GA₃, whereas flower clusters and berry samples were collected at 6, 24, and 48 h after application at full bloom and 3-4 mm berry stages. Seven hundred thirty-three genes were differentially expressed in GA₃-treated samples. At rachis and flower cluster stage respectively, 126 and 264 genes were found to be significantly differentially expressed within 6 h of GA₃ application. The number of DEG reduced considerably at 24 h. However, at berry stage, major changes occurred even at 24 h and a number of DEGs at 6 and 24 h were 174 and 191, respectively. As compared to upregulated genes, larger numbers of genes were downregulated. Stage-specific response to the GA₃ application was observed as evident from the unique set of DEGs at each stage and only a few common genes among three stages. Among the DEGs, 67 were transcription factors. Functional categorization and enrichment analysis revealed that several transcripts involved in sucrose and hexose metabolism, hormone and secondary metabolism, and abiotic and biotic stimuli were enriched in response to application of GA₃. A high correlation was recorded for real-time PCR and transcriptome data for selected DEGs, thus indicating the robustness of transcriptome data obtained in this study for understanding the GA₃ response at different stages of berry development in grape. Chromosomal localization of DEGs and identification of polymorphic microsatellite markers in selected genes have potential for their use in breeding for varieties with improved bunch architecture.

Extended diversity analysis of cultivated grapevine Vitis vinifera with 10K genome-wide SNPs

Grapevine is a very important crop species that is mainly cultivated worldwide for fruits, wine and juice. Identification of the genetic bases of performance traits through association mapping studies requires a precise knowledge of the available diversity and how this diversity is structured and varies across the whole genome. An 18k SNP genotyping array was evaluated on a panel of Vitis vinifera cultivars and we obtained a data set with no missing values for a total of 10207 SNPs and 783 different genotypes. The average inter-SNP spacing was ~47 kbp, the mean minor allele frequency (MAF) was 0.23 and the genetic diversity in the sample was high (He = 0.32). Fourteen SNPs, chosen from those with the highest MAF values, were sufficient to identify each genotype in the sample. Parentage analysis revealed 118 full parentages and 490 parent-offspring duos, thus confirming the close pedigree relationships within the cultivated grapevine. Structure analyses also confirmed the main divisions due to an eastern-western gradient and human usage (table vs. wine). Using a multivariate approach, we refined the structure and identified a total of eight clusters. Both the genetic diversity (He, 0.26-0.32) and linkage disequilibrium (LD, 28.8-58.2 kbp) varied between clusters. Despite the short span LD, we also identified some non-recombining haplotype blocks that may complicate association mapping. Finally, we performed a genome-wide association study that confirmed previous works and also identified new regions for important performance traits such as acidity. Taken together, all the results contribute to a better knowledge of the genetics of the cultivated grapevine.

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

BACKGROUND

Berry size is considered as one of the main selection criteria in table grape breeding programs. However, this is a quantitative and polygenic trait, and its genetic determination is still poorly understood. Considering its economic importance, it is relevant to determine its genetic architecture and elucidate the mechanisms involved in its expression. To approach this issue, an RNA-Seq experiment based on Illumina platform was performed (14 libraries), including seedless segregants with contrasting phenotypes for berry weight at fruit setting (FST) and 6-8 mm berries (B68) phenological stages.

RESULTS

A group of 526 differentially expressed (DE) genes were identified, by comparing seedless segregants with contrasting phenotypes for berry weight: 101 genes from the FST stage and 463 from the B68 stage. Also, we integrated differential expression, principal components analysis (PCA), correlations and network co-expression analyses to characterize the transcriptome profiling observed in segregants with contrasting phenotypes for berry weight. After this, 68 DE genes were selected as candidate genes, and seven candidate genes were validated by real time-PCR, confirming their expression profiles.

CONCLUSIONS

We have carried out the first transcriptome analysis focused on table grape seedless segregants with contrasting phenotypes for berry weight. Our findings contributed to the understanding of the mechanisms involved in berry weight determination. Also, this comparative transcriptome profiling revealed candidate genes for berry weight which could be evaluated as selection tools in table grape breeding programs.

Association analysis of grapevine bunch traits using a comprehensive approach

KEY MESSAGE

A set of SNP markers associated to bunch compactness and related traits were identified in grapevine.

ABSTRACT

Bunch compactness plays an important role in the sanitary status and perceived quality of table and wine grapes, being influenced by cultural practices and by environmental and genetic factors, which are mostly unknown. In this work, we took advantage of genetic, genomic and bioinformatic advances to analyze part of its molecular basis through a combination of transcriptomic and association analyses. Results from different transcriptomic comparisons between loose and compact grapevine clones were analyzed to select a set of candidate genes likely involved in the observed variation for bunch compactness. Up to 183 genes were sequenced in a grapevine collection, and 7032 single nucleotide polymorphisms (SNPs) were detected in more than 100 varieties with a frequency of the minor allele over 5%. They were used to test their association in three consecutive seasons with bunch compactness and two of its most influencing factors: total berry number and length of the first ramification of the rachis. Only one SNP was associated with berry number in two seasons, suggesting the high sensitiveness of this trait to seasonal environmental changes. On the other hand, we found a set of SNPs associated with both the first ramification length and bunch compactness in various seasons, in several genes which had not previously related to bunch compactness or bunch compactness-related traits. They are proposed as interesting candidates for further functional analyses aimed to verify the results obtained in this work, as a previous step to their inclusion in marker-assisted selection strategies.

Genomics-assisted breeding in fruit trees

Recent advancements in genomic analysis technologies have opened up new avenues to promote the efficiency of plant breeding. Novel genomics-based approaches for plant breeding and genetics research, such as genome-wide association studies (GWAS) and genomic selection (GS), are useful, especially in fruit tree breeding. The breeding of fruit trees is hindered by their long generation time, large plant size, long juvenile phase, and the necessity to wait for the physiological maturity of the plant to assess the marketable product (fruit). In this article, we describe the potential of genomics-assisted breeding, which uses these novel genomics-based approaches, to break through these barriers in conventional fruit tree breeding. We first introduce the molecular marker systems and whole-genome sequence data that are available for fruit tree breeding. Next we introduce the statistical methods for biparental linkage and quantitative trait locus (QTL) mapping as well as GWAS and GS. We then review QTL mapping, GWAS, and GS studies conducted on fruit trees. We also review novel technologies for rapid generation advancement. Finally, we note the future prospects of genomics-assisted fruit tree breeding and problems that need to be overcome in the breeding.

Identification of stable QTLs for vegetative and reproductive traits in the microvine (Vitis vinifera L.) using the 18 K Infinium chip

BACKGROUND

The increasing temperature associated with climate change impacts grapevine phenology and development with critical effects on grape yield and composition. Plant breeding has the potential to deliver new cultivars with stable yield and quality under warmer climate conditions, but this requires the identification of stable genetic determinants. This study tested the potentialities of the microvine to boost genetics in grapevine. A mapping population of 129 microvines derived from Picovine x Ugni Blanc flb, was genotyped with the Illumina® 18 K SNP (Single Nucleotide Polymorphism) chip. Forty-three vegetative and reproductive traits were phenotyped outdoors over four cropping cycles, and a subset of 22 traits over two cropping cycles in growth rooms with two contrasted temperatures, in order to map stable QTLs (Quantitative Trait Loci).

RESULTS

Ten stable QTLs for berry development and quality or leaf area were identified on the parental maps. A new major QTL explaining up to 44 % of total variance of berry weight was identified on chromosome 7 in Ugni Blanc flb, and co-localized with QTLs for seed number (up to 76 % total variance), major berry acids at green lag phase (up to 35 %), and other yield components (up to 25 %). In addition, a minor QTL for leaf area was found on chromosome 4 of the same parent. In contrast, only minor QTLs for berry acidity and leaf area could be found as moderately stable in Picovine. None of the transporters recently identified as mutated in low acidity apples or Cucurbits were included in the several hundreds of candidate genes underlying the above berry QTLs, which could be reduced to a few dozen candidate genes when a priori pertinent biological functions and organ specific expression were considered.

CONCLUSIONS

This study combining the use of microvine and a high throughput genotyping technology was innovative for grapevine genetics. It allowed the identification of 10 stable QTLs, including the first berry acidity QTLs reported so far in a Vitis vinifera intra-specific cross. Robustness of a set of QTLs was assessed with respect to temperature variation.

Construction of a high-density genetic map and QTLs mapping for sugars and acids in grape berries

BACKGROUND

QTLs controlling individual sugars and acids (fructose, glucose, malic acid and tartaric acid) in grape berries have not yet been identified. The present study aimed to construct a high-density, high-quality genetic map of a winemaking grape cross with a complex parentage (V. vinifera × V. amurensis) × ((V. labrusca × V. riparia) × V. vinifera), using next-generation restriction site-associated DNA sequencing, and then to identify loci related to phenotypic variability over three years.

RESULTS

In total, 1 826 SNP-based markers were developed. Of these, 621 markers were assembled into 19 linkage groups (LGs) for the maternal map, 696 for the paternal map, and 1 254 for the integrated map. Markers showed good linear agreement on most chromosomes between our genetic maps and the previously published V. vinifera reference sequence. However marker order was different in some chromosome regions, indicating both conservation and variation within the genome. Despite the identification of a range of QTLs controlling the traits of interest, these QTLs explained a relatively small percentage of the observed phenotypic variance. Although they exhibited a large degree of instability from year to year, QTLs were identified for all traits but tartaric acid and titratable acidity in the three years of the study; however only the QTLs for malic acid and β ratio (tartaric acid-to-malic acid ratio) were stable in two years. QTLs related to sugars were located within ten LGs (01, 02, 03, 04, 07, 09, 11, 14, 17, 18), and those related to acids within three LGs (06, 13, 18). Overlapping QTLs in LG14 were observed for fructose, glucose and total sugar. Malic acid, total acid and β ratio each had several QTLs in LG18, and malic acid also had a QTL in LG06. A set of 10 genes underlying these QTLs may be involved in determining the malic acid content of berries.

CONCLUSION

The genetic map constructed in this study is potentially a high-density, high-quality map, which could be used for QTL detection, genome comparison, and sequence assembly. It may also serve to broaden our understanding of the grape genome.

Evolution of Catkins: Inflorescence Morphology of Selected Salicaceae in an Evolutionary and Developmental Context

Poplars (Populus sp.) and willows (Salix sp.) are well known woody plants common throughout the northern hemisphere, both with fully sequenced genomes. They bear compact unisexual inflorescences known as "catkins." Closely related genera of the "salicoid clade" within the family Salicaceae include the Asian genera Bennettiodendron, Idesia, Itoa, Poliothyrsis, and Carrierea and the Central American genera Olmediella and Macrohasseltia. Like willow and poplar, most of these genera are dioecious, but unlike willow and poplar they generally have loosely branched panicles rather than catkins, and less highly reduced flowers. However, the early developing inflorescences of Carrierea and Idesia show similarities to catkins which suggest possible pathways by which the salicoid catkin may have evolved.