Genetic diversity, linkage disequilibrium and power of a large grapevine (Vitis vinifera L) diversity panel newly designed for association studies

Single nucleotide polymorphism (SNP) analysis reveals ancestry and genetic diversity of cultivated and wild grapevines in Croatia

BACKGROUND

Croatia is a geographically small country with a remarkable diversity of cultivated and spontaneous grapevines. Local germplasm has been characterised by microsatellite markers, but a detailed analysis based on single nucleotide polymorphisms (SNPs) is still lacking. Here we characterize the genetic diversity of 149 accessions from three germplasm repositories and four natural sites using 516,101 SNPs to identify complete parent-offspring trios and their relations with spontaneous populations, offering a proof-of-concept for the use of reduced-representation genome sequencing in population genetics and genome-wide association studies (GWAS).

RESULTS

Principal component analysis revealed a clear discontinuity between cultivated (V. vinifera subsp. sativa) and spontaneous grapevines, supporting the notion that the latter represent local populations of the wild progenitor (V. vinifera subsp. sylvestris). ADMIXTURE identified three ancestry components. Two sativa components are alternatively predominant in cultivars grown either in northern Adriatic Croatia and Continental Croatia or in Dalmatia (i.e. central and southern Adriatic Croatia). A sylvestris component, which is predominant in accessions from spontaneous populations, is a minor ancestry component in cultivated accessions. TREEMIX provided evidence of unidirectional migration from the vineyards to natural sites, suggesting that gene flow has gone preferentially from the introduced domesticated germplasm into local wild populations rather than vice versa. Identity-by-descent analysis indicated an extensive kinship network, including 14 complete parent-offspring trios, involving only cultivated accessions, six full-sibling relationships and invalidated a presumed pedigree of one of the most important varieties in Croatia, 'Plavac Mali'. Despite this strong population structure, significant association was found between 143 SNPs and berry skin colour and between 2 SNPs and leaf hairiness, across two previously known genomic regions.

CONCLUSIONS

The clear genetic separation between Croatian cultivars and sylvestris ruled out the hypothesis that those cultivars originated from local domestication events. On the other hand, the evidence of a crop-to-wild gene flow signals the need for an urgent adoption of conservation strategies that preserve the residual genetic integrity of wild relatives. The use of this reduced-representation genome sequencing protocol in grapevine enables an accurate pedigree reconstruction and can be recommended for GWAS experiments.

The genetic characterization of grapevines prospected in old Serbian vineyards reveals multiple relationships between traditional varieties of the Balkans

Serbia preserves a high number of local grape varieties, which have been cultivated across the country for centuries. Now, these ancient varieties are in the spotlight, and there is a global trend towards their recovery and characterization because they can revitalize regional, national and international grape and wine sectors. In addition, their genetic study can be useful to find new pedigree relationships to reveal how local varietal assortment evolved over time. Here, the genetic characterization of 138 grapevines from old Serbian vineyards revealed 59 different genetic profiles, 49 of which were identified as grapevine varieties whose origin in the country could be linked to some major Serbian historical periods. Most of the genetic profiles found in this work arranged in a complex pedigree network that integrates numerous grapevine varieties from diverse Balkan countries, agreeing with an intense exchange of plant material among Balkan regions for centuries. This analysis identified some varieties as important founders of Balkan genetic resources, like 'Alba Imputotato', 'Braghina Rosie', 'Coarna Alba', and 'Vulpea'. After deepening into their genealogy, these major direct founders might have ultimately derived from 'Visparola', an ancient variety of likely Balkan origin with a major founding role in some European regions. Our results also indicated the genetic singularity of the grapevine resources from the Balkans when compared to those from other relevant winemaking regions, supporting the interest of their detailed study to evaluate their oenological potential and for the eventual identification of useful traits to counteract current viticulture challenges.

Integration of machine learning models with microsatellite markers: New avenue in world grapevine germplasm characterization

Development of efficient analytical techniques is required for effective interpretation of biological data to take novel hypotheses and finding the critical predictive patterns. Machine Learning algorithms provide a novel opportunity for development of low-cost and practical solutions in biology. In this study, we proposed a new integrated analytical approach using supervised machine learning algorithms and microsatellites data of worldwide vitis populations. A total of 1378 wild (V. vinifera spp. sylvestris) and cultivated (V. vinifera spp. sativa) accessions of grapevine were investigated using 20 microsatellite markers. Data cleaning, feature selection, and supervised machine learning classification models vis, Naive Bayes, Support Vector Machine (SVM) and Tree Induction methods were implied to find most indicative and diagnostic alleles to represent wild/cultivated and originated geography of each population. Our combined approaches showed microsatellite markers with the highest differentiating capacity and proved efficiency for our pipeline of classification and prediction of vitis accessions. Moreover, our study proposed the best combination of markers for better distinguishing of populations, which can be exploited in future germplasm conservation and breeding programs.

Enhancing grapevine breeding efficiency through genomic prediction and selection index

Grapevine (Vitis vinifera) breeding reaches a critical point. New cultivars are released every year with resistance to powdery and downy mildews. However, the traditional process remains time-consuming, taking 20-25 years, and demands the evaluation of new traits to enhance grapevine adaptation to climate change. Until now, the selection process has relied on phenotypic data and a limited number of molecular markers for simple genetic traits such as resistance to pathogens, without a clearly defined ideotype, and was carried out on a large scale. To accelerate the breeding process and address these challenges, we investigated the use of genomic prediction, a methodology using molecular markers to predict genotypic values. In our study, we focused on 2 existing grapevine breeding programs: Rosé wine and Cognac production. In these programs, several families were created through crosses of emblematic and interspecific resistant varieties to powdery and downy mildews. Thirty traits were evaluated for each program, using 2 genomic prediction methods: Genomic Best Linear Unbiased Predictor and Least Absolute Shrinkage Selection Operator. The results revealed substantial variability in predictive abilities across traits, ranging from 0 to 0.9. These discrepancies could be attributed to factors such as trait heritability and trait characteristics. Moreover, we explored the potential of across-population genomic prediction by leveraging other grapevine populations as training sets. Integrating genomic prediction allowed us to identify superior individuals for each program, using multivariate selection index method. The ideotype for each breeding program was defined collaboratively with representatives from the wine-growing sector.

A super-pangenome of the North American wild grape species

BACKGROUND

Capturing the genetic diversity of wild relatives is crucial for improving crops because wild species are valuable sources of agronomic traits that are essential to enhance the sustainability and adaptability of domesticated cultivars. Genetic diversity across a genus can be captured in super-pangenomes, which provide a framework for interpreting genomic variations.

RESULTS

Here we report the sequencing, assembly, and annotation of nine wild North American grape genomes, which are phased and scaffolded at chromosome scale. We generate a reference-unbiased super-pangenome using pairwise whole-genome alignment methods, revealing the extent of the genomic diversity among wild grape species from sequence to gene level. The pangenome graph captures genomic variation between haplotypes within a species and across the different species, and it accurately assesses the similarity of hybrids to their parents. The species selected to build the pangenome are a great representation of the genus, as illustrated by capturing known allelic variants in the sex-determining region and for Pierce's disease resistance loci. Using pangenome-wide association analysis, we demonstrate the utility of the super-pangenome by effectively mapping short reads from genus-wide samples and identifying loci associated with salt tolerance in natural populations of grapes.

CONCLUSIONS

This study highlights how a reference-unbiased super-pangenome can reveal the genetic basis of adaptive traits from wild relatives and accelerate crop breeding research.

Ripening dynamics revisited: an automated method to track the development of asynchronous berries on time-lapse images

BACKGROUND

Grapevine berries undergo asynchronous growth and ripening dynamics within the same bunch. Due to the lack of efficient methods to perform sequential non-destructive measurements on a representative number of individual berries, the genetic and environmental origins of this heterogeneity, remain nearly unknown. To address these limitations, we propose a method to track the growth and coloration kinetics of individual berries on time-lapse images of grapevine bunches.

RESULTS

First, a deep-learning approach is used to detect berries with at least 50 ± 10% of visible contours, and infer the shape they would have in the absence of occlusions. Second, a tracking algorithm was developed to assign a common label to shapes representing the same berry along the time-series. Training and validation of the methods were performed on challenging image datasets acquired in a robotised high-throughput phenotyping platform. Berries were detected on various genotypes with a F1-score of 91.8%, and segmented with a mean absolute error of 4.1% on their area. Tracking allowed to label and retrieve the temporal identity of more than half of the segmented berries, with an accuracy of 98.1%. This method was used to extract individual growth and colour kinetics of various berries from the same bunch, allowing us to propose the first statistically relevant analysis of berry ripening kinetics, with a time resolution lower than one day.

CONCLUSIONS

We successfully developed a fully-automated open-source method to detect, segment and track overlapping berries in time-series of grapevine bunch images acquired in laboratory conditions. This makes it possible to quantify fine aspects of individual berry development, and to characterise the asynchrony within the bunch. The interest of such analysis was illustrated here for one cultivar, but the method has the potential to be applied in a high throughput phenotyping context. This opens the way for revisiting the genetic and environmental variations of the ripening dynamics. Such variations could be considered both from the point of view of fruit development and the phenological structure of the population, which would constitute a paradigm shift.

Boosting grapevine breeding for climate-smart viticulture: from genetic resources to predictive genomics

The multifaceted nature of climate change is increasing the urgency to select resilient grapevine varieties, or generate new, fitter cultivars, to withstand a multitude of new challenging conditions. The attainment of this goal is hindered by the limiting pace of traditional breeding approaches, which require decades to result in new selections. On the other hand, marker-assisted breeding has proved useful when it comes to traits governed by one or few genes with great effects on the phenotype, but its efficacy is still restricted for complex traits controlled by many loci. On these premises, innovative strategies are emerging which could help guide selection, taking advantage of the genetic diversity within the Vitis genus in its entirety. Multiple germplasm collections are also available as a source of genetic material for the introgression of alleles of interest via adapted and pioneering transformation protocols, which present themselves as promising tools for future applications on a notably recalcitrant species such as grapevine. Genome editing intersects both these strategies, not only by being an alternative to obtain focused changes in a relatively rapid way, but also by supporting a fine-tuning of new genotypes developed with other methods. A review on the state of the art concerning the available genetic resources and the possibilities of use of innovative techniques in aid of selection is presented here to support the production of climate-smart grapevine genotypes.

Exploring genetic diversity and population structure of a large grapevine (Vitis vinifera L.) germplasm collection in Türkiye

Grapevine (Vitis Vinifera L.) has been one of the significant perennial crops in widespread temperate climate regions since its domestication around 6000 years ago. Grapevine and its products, particularly wine, table grapes, and raisins, have significant economic importance not only in grapevine-growing countries but also worldwide. Grapevine cultivation in Türkiye dates back to ancient times, and Anatolia is considered one of the main grapevine migration routes around the Mediterranean basin. Turkish germplasm collection, conserved at the Turkish Viticulture Research Institutes, includes cultivars and wild relatives mainly collected in Türkiye, breeding lines, rootstock varieties, and mutants, but also cultivars of international origin. Genotyping with high-throughput markers enables the investigation of genetic diversity, population structure, and linkage disequilibrium, which are crucial for applying genomic-assisted breeding. Here, we present the results of a high-throughput genotyping-by-sequencing (GBS) study of 341 genotypes from grapevine germplasm collection at Manisa Viticulture Research Institute. A total of 272,962 high-quality single nucleotide polymorphisms (SNP) markers on the nineteen chromosomes were identified using genotyping-by-sequencing (GBS) technology. The high-density coverage of SNPs resulted in an average of 14,366 markers per chromosome, an average polymorphism information content (PIC) value of 0.23 and an expected heterozygosity (He) value of 0.28 indicating the genetic diversity within 341 genotypes. LD decayed very fast when r² was between 0.45 and 0.2 and became flat when r² was 0.05. The average LD decay for the entire genome was 30 kb when r^2 = 0.2. The PCA and structure analysis did not distinguish the grapevine genotypes based on different origins, highlighting the occurrence of gene flow and a high amount of admixture. Analysis of molecular variance (AMOVA) results indicated a high level of genetic differentiation within populations, while variation among populations was extremely low. This study provides comprehensive information on the genetic diversity and population structure of Turkish grapevine genotypes.

Vitis vinifera genotyping toolbox to highlight diversity and germplasm identification

The contribution of vine cultivation to human welfare as well as the stimulation of basic social and cultural features of civilization has been great. The wide temporal and regional distribution created a wide array of genetic variants that have been used as propagating material to promote cultivation. Information on the origin and relationships among cultivars is of great interest from a phylogenetics and biotechnology perspective. Fingerprinting and exploration of the complicated genetic background of varieties may contribute to future breeding programs. In this review, we present the most frequently used molecular markers, which have been used on Vitis germplasm. We discuss the scientific progress that led to the new strategies being implemented utilizing state-of-the-art next generation sequencing technologies. Additionally, we attempted to delimit the discussion on the algorithms used in phylogenetic analyses and differentiation of grape varieties. Lastly, the contribution of epigenetics is highlighted to tackle future roadmaps for breeding and exploitation of Vitis germplasm. The latter will remain in the top of the edge for future breeding and cultivation and the molecular tools presented herein, will serve as a reference point in the challenging years to come.

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.

Interest of phenomic prediction as an alternative to genomic prediction in grapevine

BACKGROUND

Phenomic prediction has been defined as an alternative to genomic prediction by using spectra instead of molecular markers. A reflectance spectrum provides information on the biochemical composition within a tissue, itself being under genetic determinism. Thus, a relationship matrix built from spectra could potentially capture genetic signal. This new methodology has been mainly applied in several annual crop species but little is known so far about its interest in perennial species. Besides, phenomic prediction has only been tested for a restricted set of traits, mainly related to yield or phenology. This study aims at applying phenomic prediction for the first time in grapevine, using spectra collected on two tissues and over two consecutive years, on two populations and for 15 traits, related to berry composition, phenology, morphological and vigour. A major novelty of this study was to collect spectra and phenotypes several years apart from each other. First, we characterized the genetic signal in spectra and under which condition it could be maximized, then phenomic predictive ability was compared to genomic predictive ability.

RESULTS

For the first time, we showed that the similarity between spectra and genomic relationship matrices was stable across tissues or years, but variable across populations, with co-inertia around 0.3 and 0.6 for diversity panel and half-diallel populations, respectively. Applying a mixed model on spectra data increased phenomic predictive ability, while using spectra collected on wood or leaves from one year or another had less impact. Differences between populations were also observed for predictive ability of phenomic prediction, with an average of 0.27 for the diversity panel and 0.35 for the half-diallel. For both populations, a significant positive correlation was found across traits between predictive ability of genomic and phenomic predictions.

CONCLUSION

NIRS is a new low-cost alternative to genotyping for predicting complex traits in perennial species such as grapevine. Having spectra and phenotypes from different years allowed us to exclude genotype-by-environment interactions and confirms that phenomic prediction can rely only on genetics.

A genome-wide association and prediction study in grapevine deciphers the genetic architecture of multiple traits and identifies genes under many new QTLs

To cope with the challenges facing agriculture, speeding-up breeding programs is a worthy endeavor, especially for perennial species such as grapevine, but requires understanding the genetic architecture of target traits. To go beyond the mapping of quantitative trait loci in bi-parental crosses, we exploited a diversity panel of 279 Vitis vinifera L. cultivars planted in 5 blocks in the vineyard. This panel was phenotyped over several years for 127 traits including yield components, organic acids, aroma precursors, polyphenols, and a water stress indicator. The panel was genotyped for 63k single nucleotide polymorphisms by combining an 18K microarray and genotyping-by-sequencing. The experimental design allowed to reliably assess the genotypic values for most traits. Marker densification via genotyping-by-sequencing markedly increased the proportion of genetic variance explained by single nucleotide polymorphisms, and 2 multi-single nucleotide polymorphism models identified quantitative trait loci not found by a single nucleotide polymorphism-by-single nucleotide polymorphism model. Overall, 489 reliable quantitative trait loci were detected for 41% more response variables than by a single nucleotide polymorphism-by-single nucleotide polymorphism model with microarray-only single nucleotide polymorphisms, many new ones compared with the results from bi-parental crosses. A prediction accuracy higher than 0.42 was obtained for 50% of the response variables. Our overall approach as well as quantitative trait locus and prediction results provide insights into the genetic architecture of target traits. New candidate genes and the application into breeding are discussed.

Bioinformatic approach for the discovery of cis-eQTL signals during fruit ripening of a woody species as grape (Vitis vinifera L.)

Expression quantitative trait loci (eQTLs) are associations between genetic variants, such as Single Nucleotide Polymorphisms (SNPs), and gene expression. eQTLs are an important tool to understand the genetic variance of gene expression of complex phenotypes. eQTLs analyses are common in biomedical models but are scarce in woody crop species such as fruit trees or grapes. In this study, a comprehensive bioinformatic analysis was conducted leveraging with expression data from two different growth stages, around ripening onset, of 10 genotypes of grape (Vitis vinifera L.). A total of 2170 cis-eQTL were identified in 212 gene modulated at ripening onset. The 48% of these DEGs have a known function. Among the annotated protein-coding genes, terpene synthase, auxin-regulatory factors, GRFS, ANK_REP_REGION domain-containing protein, Kinesin motor domain-containing protein and flavonol synthase were noted. This new inventory of cis-eQTLs influencing gene expression during fruit ripening will be an important resource to examine variation for this trait and will help to elucidate the complex genetic architecture underlying this process in grape.

Across-population genomic prediction in grapevine opens up promising prospects for breeding

Crop breeding involves two selection steps: choosing progenitors and selecting individuals within progenies. Genomic prediction, based on genome-wide marker estimation of genetic values, could facilitate these steps. However, its potential usefulness in grapevine (Vitis vinifera L.) has only been evaluated in non-breeding contexts mainly through cross-validation within a single population. We tested across-population genomic prediction in a more realistic breeding configuration, from a diversity panel to ten bi-parental crosses connected within a half-diallel mating design. Prediction quality was evaluated over 15 traits of interest (related to yield, berry composition, phenology and vigour), for both the average genetic value of each cross (cross mean) and the genetic values of individuals within each cross (individual values). Genomic prediction in these conditions was found useful: for cross mean, average per-trait predictive ability was 0.6, while per-cross predictive ability was halved on average, but reached a maximum of 0.7. Mean predictive ability for individual values within crosses was 0.26, about half the within-half-diallel value taken as a reference. For some traits and/or crosses, these across-population predictive ability values are promising for implementing genomic selection in grapevine breeding. This study also provided key insights on variables affecting predictive ability. Per-cross predictive ability was well predicted by genetic distance between parents and when this predictive ability was below 0.6, it was improved by training set optimization. For individual values, predictive ability mostly depended on trait-related variables (magnitude of the cross effect and heritability). These results will greatly help designing grapevine breeding programs assisted by genomic prediction.

Harnessing multivariate, penalized regression methods for genomic prediction and QTL detection of drought-related traits in grapevine

Viticulture has to cope with climate change and to decrease pesticide inputs, while maintaining yield and wine quality. Breeding is a key lever to meet this challenge, and genomic prediction a promising tool to accelerate breeding programs. Multivariate methods are potentially more accurate than univariate ones. Moreover, some prediction methods also provide marker selection, thus allowing quantitative trait loci (QTLs) detection and the identification of positional candidate genes. To study both genomic prediction and QTL detection for drought-related traits in grapevine, we applied several methods, interval mapping (IM) as well as univariate and multivariate penalized regression, in a bi-parental progeny. With a dense genetic map, we simulated two traits under four QTL configurations. The penalized regression method Elastic Net (EN) for genomic prediction, and controlling the marginal False Discovery Rate on EN selected markers to prioritize the QTLs. Indeed, penalized methods were more powerful than IM for QTL detection across various genetic architectures. Multivariate prediction did not perform better than its univariate counterpart, despite strong genetic correlation between traits. Using 14 traits measured in semi-controlled conditions under different watering conditions, penalized regression methods proved very efficient for intra-population prediction whatever the genetic architecture of the trait, with predictive abilities reaching 0.68. Compared to a previous study on the same traits, these methods applied on a denser map found new QTLs controlling traits linked to drought tolerance and provided relevant candidate genes. Overall, these findings provide a strong evidence base for implementing genomic prediction in grapevine breeding.

Genome-Wide Identification of Sequence Variations and SSR Marker Development in the Munake Grape Cultivar

The Munake grape cultivar produces uniquely flavored and high-quality fruits. Despite the numerous beneficial agronomic traits of Munake, there are few genetic resources available for this cultivar. To address this knowledge gap, the entire genome was sequenced using whole-genome sequencing approaches and compared with a Vitis vinifera L. reference genome. This study describes more than 3 million single nucleotide polymorphism (SNP), 300,000 insertion and deletion (InDel), 14,000 structural variation (SV), and 80,000 simple sequence repeat (SSR) markers (one SSR per 4.23 kb), as well as their primers. Among the SSRs, 44 SSR primer pairs were randomly selected and validated by polymerase chain reaction (PCR), allowing discrimination between the different Munake cultivar genotypes. The genetic data provided allow a deeper understanding of Munake cultivar genomic sequence and contribute to better knowledge of the genetic basis behind its key agronomic traits.

Genetic diversity and population structure in Vitis species illustrate phylogeographic patterns in eastern North America

Geographical distribution and diversity of current plant species have been strongly shaped by climatic oscillations during the Quaternary. Analysing the resulting divergence among species and differentiation within species is crucial to understand the evolution of taxa like the Vitis genus, which provides very useful genetic resources for grapevine improvement and might reveal original recolonization patterns due to growth habit and dispersal mode. Here, we studied the genetic structure in natural populations of three species from eastern North America: Vitis aestivalis, V. cinerea and V. riparia using different marker types. Vitis aestivalis and V. cinerea showed higher diversity than V. riparia. The two former species are less differentiated, confirming an earlier divergence of V. riparia. V. aestivalis and V. riparia exhibited different genetic groups on both sides of the Appalachian Mountains that could mirror different recolonization routes from southern refugia. Genetic structure was stronger in V. cinerea, for which two varieties (var. berlandieri and var. cinerea) are morphologically recognized. Our results confirm this distinction and suggest the existence of three other lineages within var. cinerea. These discontinuities appear linked to adaptation of var. berlandieri to dry and limy areas of Texas and partially to the Mississippi River Valley. Rapid range expansions from refugia upon climate warming are also suggested by the low linkage disequilibrium values observed. Furthermore, large variation for downy mildew resistance was observed in the three species. Our findings appeared consistent with the vegetation history of eastern North America.

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.

Recent Large-Scale Genotyping and Phenotyping of Plant Genetic Resources of Vegetatively Propagated Crops

Several recent national and international projects have focused on large-scale genotyping of plant genetic resources in vegetatively propagated crops like fruit and berries, potatoes and woody ornamentals. The primary goal is usually to identify true-to-type plant material, detect possible synonyms, and investigate genetic diversity and relatedness among accessions. A secondary goal may be to create sustainable databases that can be utilized in research and breeding for several years ahead. Commonly applied DNA markers (like microsatellite DNA and SNPs) and next-generation sequencing each have their pros and cons for these purposes. Methods for large-scale phenotyping have lagged behind, which is unfortunate since many commercially important traits (yield, growth habit, storability, and disease resistance) are difficult to score. Nevertheless, the analysis of gene action and development of robust DNA markers depends on environmentally controlled screening of very large sets of plant material. Although more time-consuming, co-operative projects with broad-scale data collection are likely to produce more reliable results. In this review, we will describe some of the approaches taken in genotyping and/or phenotyping projects concerning a wide variety of vegetatively propagated crops.

Molecular Tools for Adapting Viticulture to Climate Change

Adaptation of viticulture to climate change includes exploration of new geographical areas, new training systems, new management practices, or new varieties, both for rootstocks and scions. Molecular tools can be defined as molecular approaches used to study DNAs, RNAs, and proteins in all living organisms. We present here the current knowledge about molecular tools and their potential usefulness in three aspects of grapevine adaptation to the ongoing climate change. (i) Molecular tools for understanding grapevine response to environmental stresses. A fine description of the regulation of gene expression is a powerful tool to understand the physiological mechanisms set up by the grapevine to respond to abiotic stress such as high temperatures or drought. The current knowledge on gene expression is continuously evolving with increasing evidence of the role of alternative splicing, small RNAs, long non-coding RNAs, DNA methylation, or chromatin activity. (ii) Genetics and genomics of grapevine stress tolerance. The description of the grapevine genome is more and more precise. The genetic variations among genotypes are now revealed with new technologies with the sequencing of very long DNA molecules. High throughput technologies for DNA sequencing also allow now the genetic characterization at the same time of hundreds of genotypes for thousands of points in the genome, which provides unprecedented datasets for genotype-phenotype associations studies. We review the current knowledge on the genetic determinism of traits for the adaptation to climate change. We focus on quantitative trait loci and molecular markers available for developmental stages, tolerance to water stress/water use efficiency, sugar content, acidity, and secondary metabolism of the berries. (iii) Controlling the genome and its expression to allow breeding of better-adapted genotypes. High-density DNA genotyping can be used to select genotypes with specific interesting alleles but genomic selection is also a powerful method able to take into account the genetic information along the whole genome to predict a phenotype. Modern technologies are also able to generate mutations that are possibly interesting for generating new phenotypes but the most promising one is the direct editing of the genome at a precise location.

Population structure, genetic diversity and genomic selection signatures among a Brazilian common bean germplasm

Brazil is the world's largest producer of common bean. Knowledge of the genetic diversity and relatedness of accessions adapted to Brazilian conditions is of great importance for the conservation of germplasm and for directing breeding programs aimed at the development of new cultivars. In this context, the objective of this study was to analyze the genetic diversity, population structure, and linkage disequilibrium (LD) of a diversity panel consisting of 219 common bean accessions, most of which belonging to the Mesoamerican gene pool. Genotyping by sequencing (GBS) of these accessions allowed the identification of 49,817 SNPs with minor allele frequency > 0.05. Of these, 17,149 and 12,876 were exclusive to the Mesoamerican and Andean pools, respectively, and 11,805 SNPs could differentiate the two gene pools. Further the separation according to the gene pool, bayesian analysis of the population structure showed a subdivision of the Mesoamerican accessions based on the origin and color of the seed tegument. LD analysis revealed the occurrence of long linkage blocks and low LD decay with physical distance between SNPs (LD half decay in 249 kb, corrected for population structure and relatedness). The GBS technique could effectively characterize the Brazilian common bean germplasms, and the diversity panel used in this study may be of great use in future genome-wide association studies.

Mining Grapevine Downy Mildew Susceptibility Genes: A Resource for Genomics-Based Breeding and Tailored Gene Editing

Several pathogens continuously threaten viticulture worldwide. Until now, the investigation on resistance loci has been the main trend to understand the interaction between grapevine and the mildew causal agents. Dominantly inherited gene-based resistance has shown to be race-specific in some cases, to confer partial immunity, and to be potentially overcome within a few years since its introgression. Recently, on the footprint of research conducted in Arabidopsis, putative genes associated with downy mildew susceptibility have been discovered also in the grapevine genome. In this work, we deep-sequenced four putative susceptibility genes-namely VvDMR6.1, VvDMR6.2, VvDLO1, VvDLO2-in 190 genetically diverse grapevine genotypes to discover new sources of broad-spectrum and recessively inherited resistance. Identified Single Nucleotide Polymorphisms were screened in a bottleneck analysis from the genetic sequence to their impact on protein structure. Fifty-five genotypes showed at least one impacting mutation in one or more of the scouted genes. Haplotypes were inferred for each gene and two of them at the VvDMR6.2 gene were found significantly more represented in downy mildew resistant genotypes. The current results provide a resource for grapevine and plant genetics and could corroborate genomic-assisted breeding programs as well as tailored gene editing approaches for resistance to biotic stresses.

Candidate genes and SNPs associated with stomatal conductance under drought stress in Vitis

BACKGROUND

Understanding the complexity of the vine plant's response to water deficit represents a major challenge for sustainable winegrowing. Regulation of water use requires a coordinated action between scions and rootstocks on which cultivars are generally grafted to cope with phylloxera infestations. In this regard, a genome-wide association study (GWAS) approach was applied on an 'ad hoc' association mapping panel including different Vitis species, in order to dissect the genetic basis of transpiration-related traits and to identify genomic regions of grape rootstocks associated with drought tolerance mechanisms. The panel was genotyped with the GrapeReSeq Illumina 20 K SNP array and SSR markers, and infrared thermography was applied to estimate stomatal conductance values during progressive water deficit.

RESULTS

In the association panel the level of genetic diversity was substantially lower for SNPs loci (0.32) than for SSR (0.87). GWAS detected 24 significant marker-trait associations along the various stages of drought-stress experiment and 13 candidate genes with a feasible role in drought response were identified. Gene expression analysis proved that three of these genes (VIT_13s0019g03040, VIT_17s0000g08960, VIT_18s0001g15390) were actually induced by drought stress. Genetic variation of VIT_17s0000g08960 coding for a raffinose synthase was further investigated by resequencing the gene of 85 individuals since a SNP located in the region (chr17_10,497,222_C_T) was significantly associated with stomatal conductance.

CONCLUSIONS

Our results represent a step forward towards the dissection of genetic basis that modulate the response to water deprivation in grape rootstocks. The knowledge derived from this study may be useful to exploit genotypic and phenotypic diversity in practical applications and to assist further investigations.

Population demographic history and population structure for Pakistani Nili-Ravi breeding bulls based on SNP genotyping to identify genomic regions associated with male effects for milk yield and body weight

The domestic Nili-Ravi water buffalo (Bubalus bubalis) is the best dairy animal contributing 68% to total milk production in Pakistan. In this study, we identified genome-wide single nucleotide polymorphisms (SNPs) to estimate various population genetic parameters such as diversity, pairwise population differentiation, linkage disequilibrium (LD) distribution and for genome-wide association study for milk yield and body weight traits in the Nili-Ravi dairy bulls that they may pass on to their daughters who are retained for milking purposes. The genotyping by sequencing approach revealed 13,039 reference genome-anchored SNPs with minor allele frequency of 0.05 among 167 buffalos. Population structure analysis revealed that the bulls were grouped into two clusters (K = 2), which indicates the presence of two different lineages in the Pakistani Nili-Ravi water buffalo population, and we showed the extent of admixture of these two lineages in our bull collection. LD analysis revealed 4169 significant SNP associations, with an average LD decay of 90 kb for these buffalo genome. Genome-wide association study involved a multi-locus mixed linear model for milk yield and body weight to identify genome-wide male effects. Our study further illustrates the utility of the genotyping by sequencing approach for identifying genomic regions to uncover additional demographic complexity and to improve the complex dairy traits of the Pakistani Nili-Ravi water buffalo population that would provide the lot of economic benefits to dairy industry.

Induction of cell cycle arrest and apoptosis by tomentosin in hepatocellular carcinoma HepG2 and Huh7 cells

Tomentosin, a sesquiterpene lactone, is known to possess various biological activities. However, its anticarcinogenic activity against human hepatocellular carcinoma (HCC) cells has not been investigated in detail. Thus, this study aimed to elucidate the cytotoxic mechanism of tomentosin in human HCC cell lines HepG2 and Huh7. WST-1, cell counting, and colony formation assay results showed that treatment with tomentosin decreased the viability and suppressed the proliferation rate of HepG2 and Huh7 cells in a dose- and time-dependent manner. Cell cycle analysis revealed increased population of cells at the SubG1 and G2/M stage, and decreased population of cells at the G0/1 stage in HepG2 and Huh7 cells treated with tomentosin. Annexin V/propidium iodide double staining and TUNEL assay results showed increased apoptotic cell population and DNA fragmentation in HepG2 and Huh7 cells treated with tomentosin. Western blotting analysis results showed that tomentosin treatment significantly increased the expression level of Bax, Bim (short form), cleaved PARP1, FOXO3, p53, pSer15p53, pSer20p53, pSer46p53, p21, and p27, but decreased the expression of Bcl2, caspase3, caspase7, caspase9, cyclin-dependent kinase 2 (CDK2), CDK4, CDK6, cyclinB1, cyclinD1, cyclinD2, cyclinD3, and cyclinE in a dose-dependent manner. Taken together, this study revealed that tomentosin, which acted through cell cycle arrest and apoptosis, may be a useful therapeutic option against HCC.

Genetic Diversity, Population Structure, and Parentage Analysis of Croatian Grapevine Germplasm

Croatian viticulture was most extensive at the beginning of the 20th century, when about 400 varieties were in use. Autochthonous varieties are the result of spontaneous hybridization from the pre-phylloxera era and are still cultivated today on about 35 % of vineyard area, while some exist only in repositories. We present what is the most comprehensive genetic analysis of all major Croatian national repositories, with a large number of microsatellite, or simple sequence repeat (SSR) markers, and it is also the first study to apply single nucleotide polymorphism (SNP) markers. After 212 accessions were fingerprinted, 95 were classified as unique to Croatian germplasm. Genetic diversity of Croatian germplasm is rather high considering its size. SNP markers proved useful for fingerprinting but less informative and practical than SSRs. Analysis of the genetic structure showed that Croatian germplasm is predominantly part of the Balkan grape gene pool. A high number of admixed varieties and synonyms is a consequence of complex pedigrees and migrations. Parentage analysis confirmed 24 full parentages, as well as 113 half-kinships. Unexpectedly, several key genitors could not be detected within the present Croatian germplasm. The low number of reconstructed parentages (19%) points to severe genetic erosion and stresses the importance of germplasm repositories.

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.

Recommendations for Choosing the Genotyping Method and Best Practices for Quality Control in Crop Genome-Wide Association Studies

High-throughput genotyping boosts genome-wide association studies (GWAS) in crop species, leading to the identification of single-nucleotide polymorphisms (SNPs) associated with economically important traits. Choosing a cost-effective genotyping method for crop GWAS requires careful examination of several aspects, namely, the purpose and the scale of the study, crop-specific genomic features, and technical and economic matters associated with each genotyping option. Once genotypic data have been obtained, quality control (QC) procedures must be applied to avoid bias and false signals in genotype–phenotype association tests. QC for human GWAS has been extensively reviewed; however, QC for crop GWAS may require different actions, depending on the GWAS population type. Here, we review most popular genotyping methods based on next-generation sequencing (NGS) and array hybridization, and report observations that should guide the investigator in the choice of the genotyping method for crop GWAS. We provide recommendations to perform QC in crop species, and deliver an overview of bioinformatics tools that can be used to accomplish all needed tasks. Overall, this work aims to provide guidelines to harmonize those procedures leading to SNP datasets ready for crop GWAS.

Identification of Selective Sweeps in the Domesticated Table and Wine Grape (Vitis vinifera L.)

Grapevine (Vitis vinifera) is one of the most important fruit species in the Classical Mediterranean world. It is thought to have been domesticated 6,000-8,000 years ago in the Near East. However, the domestication of its wild relative into wine grapes or table grapes remains largely unknown. In this study, we analyzed 30 table grapes, 30 wine grapes, 30 dual-purpose grape accessions, as well as 30 wild relatives (Vitis vinifera ssp. sylvestris). The phenotypic comparison showed striking differences in berry weight, acidity and the content of aroma. Based on a total of 7,522,958 single-nucleotide polymorphisms, we identified several significant selective sweep regions for table and wine grapes. Besides the well-known sex-determination locus on chromosome 2, the other four highest signals shared by table and wine grapes could not be linked to the known QTLs. The identification of these genomic regions under selection sweep may reveal agronomically important traits that have been selected during grape domestication. This information not only sheds light on the mechanisms of adaptions and diversification, but also guide the genetic improvement in breeding programs.

Genome-Wide Identification and Expression Analysis of GA2ox, GA3ox, and GA20ox Are Related to Gibberellin Oxidase Genes in Grape (Vitis Vinifera L.)

Gibberellin (GAs) plays the important role in the regulation of grape developmental and growth processes. The bioinformatics analysis confirmed the differential expression of GA2, GA3, and GA20 gibberellin oxidase genes (VvGA2oxs, VvGA3oxs, and VvGA20oxs) in the grape genome, and laid a theoretical basis for exploring its role in grape. Based on the Arabidopsis GA2oxs, GA3oxs, and GA20oxs genes already reported, the VvGA2oxs, VvGA3oxs, and VvGA20oxs genes in the grape genome were identified using the BLAST software in the grape genome database. Bioinformatics analysis was performed using software such as DNAMAN v.5.0, Clustalx, MapGene2Chrom, MEME, GSDS v.2.0, ExPASy, DNAsp v.5.0, and MEGA v.7.0. Chip expression profiles were generated using grape Affymetrix GeneChip 16K and Grape eFP Browser gene chip data in PLEXdb. The expression of VvGA2oxs, VvGA3oxs, and VvGA20oxs gene families in stress was examined by qRT-PCR (Quantitative real-time-PCR). There are 24 GAoxs genes identified with the grape genome that can be classified into seven subgroups based on a phylogenetic tree, gene structures, and conserved Motifs in our research. The gene family has higher codon preference, while selectivity is negative selection of codon bias and selective stress was analyzed. The expression profiles indicated that the most of VvGAox genes were highly expressed under different time lengths of ABA (Abscisic Acid) treatment, NaCl, PEG and 5 °C. Tissue expression analysis showed that the expression levels of VvGA2oxs and VvGA20oxs in different tissues at different developmental stages of grapes were relatively higher than that of VvGA3oxs. Last but not least, qRT-PCR (Real-time fluorescent quantitative PCR) was used to determine the relative expression of the GAoxs gene family under the treatment of GA3 (gibberellin 3) and uniconazole, which can find that some VvGA2oxs was upregulated under GA3 treatment. Simultaneously, some VvGA3oxs and VvGA20oxs were upregulated under uniconazole treatment. In a nutshell, the GA2ox gene mainly functions to inactivate biologically active GAs, while GA20ox mainly degrades C20 gibberellins, and GA3ox is mainly composed of biologically active GAs. The comprehensive analysis of the three classes of VvGAoxs would provide a basis for understanding the evolution and function of the VvGAox gene family in a grape plant.

Genotyping-by-sequencing and SNP-arrays are complementary for detecting quantitative trait loci by tagging different haplotypes in association studies

BACKGROUND

Single Nucleotide Polymorphism (SNP) array and re-sequencing technologies have different properties (e.g. calling rate, minor allele frequency profile) and drawbacks (e.g. ascertainment bias). This lead us to study their complementarity and the consequences of using them separately or combined in diversity analyses and Genome-Wide Association Studies (GWAS). We performed GWAS on three traits (grain yield, plant height and male flowering time) measured in 22 environments on a panel of 247 F1 hybrids obtained by crossing 247 diverse dent maize inbred lines with a same flint line. The 247 lines were genotyped using three genotyping technologies (Genotyping-By-Sequencing, Illumina Infinium 50 K and Affymetrix Axiom 600 K arrays).

RESULTS

The effects of ascertainment bias of the 50 K and 600 K arrays were negligible for deciphering global genetic trends of diversity and for estimating relatedness in this panel. We developed an original approach based on linkage disequilibrium (LD) extent in order to determine whether SNPs significantly associated with a trait and that are physically linked should be considered as a single Quantitative Trait Locus (QTL) or several independent QTLs. Using this approach, we showed that the combination of the three technologies, which have different SNP distributions and densities, allowed us to detect more QTLs (gain in power) and potentially refine the localization of the causal polymorphisms (gain in resolution).

CONCLUSIONS

Conceptually different technologies are complementary for detecting QTLs by tagging different haplotypes in association studies. Considering LD, marker density and the combination of different technologies (SNP-arrays and re-sequencing), the genotypic data available were most likely enough to well represent polymorphisms in the centromeric regions, whereas using more markers would be beneficial for telomeric regions.

Integration of expert knowledge in the definition of Swiss pear core collection

Core collections (CCs) constitute a key tool for the characterization and management of genetic resources (GR). When the institutions involved in GR preservation decide to define a CC, they frequently prefer to select accessions based not only on strictly objective criteria, but also to add others following expert knowledge considerations (popularity, prestige, role in breeding history, or presence of phenotypic features of interest). The aim of this study was to evaluate the implications of approaches that combine formal analytical procedures and expert knowledge on the efficiency of CC definition through a case study to establish a pear CC from the Swiss National Pear Inventory. The CC had to represent a maximum of the genetic diversity, not to exceed 150 accessions, and required to include a priority set (SPPS) with 86 genotypes selected based on expert knowledge. In total, nine strategies were evaluated, resulting of combining compositions of the dataset sampled, sampling sizes and methods. The CCs sampled by mixed approaches provided similar scores, irrespective of the approach considered, and obtained similar efficiency in optimizing the genetic diversity retained. Therefore, mixed approaches can be an appropriate choice for applications involving genetic conservation in tree germplasm collections.

Facing Climate Change: Biotechnology of Iconic Mediterranean Woody Crops

The Mediterranean basin is especially sensitive to the adverse outcomes of climate change and especially to variations in rainfall patterns and the incidence of extremely high temperatures. These two concurring adverse environmental conditions will surely have a detrimental effect on crop performance and productivity that will be particularly severe on woody crops such as citrus, olive and grapevine that define the backbone of traditional Mediterranean agriculture. These woody species have been traditionally selected for traits such as improved fruit yield and quality or alteration in harvesting periods, leaving out traits related to plant field performance. This is currently a crucial aspect due to the progressive and imminent effects of global climate change. Although complete genome sequence exists for sweet orange (Citrus sinensis) and clementine (Citrus clementina), olive tree (Olea europaea) and grapevine (Vitis vinifera), the development of biotechnological tools to improve stress tolerance still relies on the study of the available genetic resources including interspecific hybrids, naturally occurring (or induced) polyploids and wild relatives under field conditions. To this respect, post-genomic era studies including transcriptomics, metabolomics and proteomics provide a wide and unbiased view of plant physiology and biochemistry under adverse environmental conditions that, along with high-throughput phenotyping, could contribute to the characterization of plant genotypes exhibiting physiological and/or genetic traits that are correlated to abiotic stress tolerance. The ultimate goal of precision agriculture is to improve crop productivity, in terms of yield and quality, making a sustainable use of land and water resources under adverse environmental conditions using all available biotechnological tools and high-throughput phenotyping. This review focuses on the current state-of-the-art of biotechnological tools such as high throughput -omics and phenotyping on grapevine, citrus and olive and their contribution to plant breeding programs.

Characterization of molecular diversity and genome-wide association study of stripe rust resistance at the adult plant stage in Northern Chinese wheat landraces

BACKGROUND

Stripe rust is a serious fungal disease of wheat (Triticum aestivum L.) caused by Puccinia striiformis f. sp. tritici (Pst), which results in yield reduction and decreased grain quality. Breeding for genetic resistance to stripe rust is the most cost-effective method to control the disease. In the present study, a genome-wide association study (GWAS) was conducted to identify markers linked to stripe rust resistance genes (or loci) in 93 Northern Chinese wheat landraces, using Diversity Arrays Technology (DArT) and simple sequence repeat (SSR) molecular marker technology based on phenotypic data from two field locations over two growing seasons in China.

RESULTS

Seventeen accessions were verified to display stable and high levels of adult plant resistance (APR) to stripe rust via multi-environment field assessments. Significant correlations among environments and high heritability were observed for stripe rust infection type (IT) and disease severity (DS). Using mixed linear models (MLM) for the GWAS, a total of 32 significantly associated loci (P < 0.001) were detected. In combination with the linkage disequilibrium (LD) decay distance (6.4 cM), 25 quantitative trait loci (QTL) were identified. Based on the integrated map of previously reported genes and QTL, six QTL located on chromosomes 4A, 6A and 7D were mapped far from resistance regions identified previously, and represent potentially novel stripe rust resistance loci at the adult plant stage.

CONCLUSIONS

The present findings demonstrated that identification of genes or loci linked to significant markers in wheat by GWAS is feasible. Seventeen elite accessions conferred with stable and high resistance to stripe rust, and six putative newly detected APR loci were identified among the 93 Northern Chinese wheat landraces. The results illustrate the potential for acceleration of molecular breeding of wheat, and also provide novel sources of stripe rust resistance with potential utility in the breeding of improved wheat cultivars.

Whole-genome resequencing of 472 Vitis accessions for grapevine diversity and demographic history analyses

Understanding the Vitis species at the genomic level is important for cultivar improvement of grapevine. Here we report whole-genome genetic variation at single-base resolution of 472 Vitis accessions, which cover 48 out of 60 extant Vitis species from a wide geographic distribution. The variation helps to identify a recent dramatic expansion and contraction of effective population size in the domesticated grapevines and that cultivars from the pan-Black Sea region have a unique demographic history in comparison to the other domesticated cultivars. We also find selective sweeps for berry edibility and stress resistance improvement. Furthermore, we find associations between candidate genes and important agronomic traits, such as berry shape and aromatic compounds. These results demonstrate resource value of the resequencing data for illuminating the evolutionary biology of Vitis species and providing targets for grapevine genetic improvement.

Despite the importance of grapevine cultivation in human history and the economic values of cultivar improvement, large-scale genomic variation data are lacking. Here the authors resequence 472 Vitis accessions and use the identified genetic variations for domestication history, demography, and GWAS analyses.

Linking Terpene Synthases to Sesquiterpene Metabolism in Grapevine Flowers

Grapevine (Vitis vinifera L.) terpene synthases (VviTPS) are responsible for the biosynthesis of terpenic volatiles. Volatile profiling of nine commercial wine cultivars showed unique cultivar-specific variation in volatile terpenes emitted from grapevine flowers. The flower chemotypes of three divergent cultivars, Muscat of Alexandria, Sauvignon Blanc and Shiraz were subsequently investigated at two flower developmental stages (EL-18 and -26). The cultivars displayed unique flower sesquiterpene compositions that changed during flower organogenesis and the profiles were dominated by either (E)-β-farnesene, (E,E)-α-farnesene or (+)-valencene. In silico remapping of microarray probes to VviTPS gene models allowed for a meta-analysis of VviTPS expression patterns in the grape gene atlas to identify genes that could regulate terpene biosynthesis in flowers. Selected sesquiterpene synthase genes were isolated and functionally characterized in three cultivars. Genotypic differences that could be linked to the function of a targeted gene model resulted in the isolation of a novel and cultivar-specific single product sesquiterpene synthase from Muscat of Alexandria flowers (VvivMATPS10), synthesizing (E)-β-farnesene as its major volatile. Furthermore, we identified structural variations (SNPs, InDels and splice variations) in the characterized VviTPS genes that potentially impact enzyme function and/or volatile sesquiterpene production in a cultivar-specific manner.

Genetic Mapping of the Incompatibility Locus in Olive and Development of a Linked Sequence-Tagged Site Marker

The genetic control of self-incompatibility (SI) has been recently disclosed in olive. Inter-varietal crossing confirmed the presence of only two incompatibility groups (G1 and G2), suggesting a simple Mendelian inheritance of the trait. A double digest restriction associated DNA (ddRAD) sequencing of a biparental population segregating for incompatibility groups has been performed and high-density linkage maps were constructed in order to map the SI locus and identify gene candidates and linked markers. The progeny consisted of a full-sib family of 229 individuals derived from the cross 'Leccino' (G1) × 'Dolce Agogia' (G2) varieties, segregating 1:1 (G1:G2), in accordance with a diallelic self-incompatibility (DSI) model. A total of 16,743 single nucleotide polymorphisms was identified, 7,006 in the female parent 'Leccino' and 9,737 in the male parent 'Dolce Agogia.' Each parental map consisted of 23 linkage groups and showed an unusual large size (5,680 cM in 'Leccino' and 3,538 cM in 'Dolce Agogia'). Recombination was decreased across all linkage groups in pollen mother cells of 'Dolce Agogia,' the parent with higher heterozygosity, compared to megaspore mother cells of 'Leccino,' in a context of a species that showed exceptionally high recombination rates. A subset of 109 adult plants was assigned to either incompatibility group by a stigma test and the diallelic self-incompatibility (DSI) locus was mapped to an interval of 5.4 cM on linkage group 18. This region spanned a size of approximately 300 Kb in the olive genome assembly. We developed a sequence-tagged site marker in the DSI locus and identified five haplotypes in 57 cultivars with known incompatibility group assignment. A combination of two single-nucleotide polymorphisms (SNPs) was sufficient to predict G1 or G2 phenotypes in olive cultivars, enabling early marker-assisted selection of compatible genotypes and allowing for a rapid screening of inter-compatibility among cultivars in order to guarantee effective fertilization and increase olive production. The construction of high-density linkage maps has led to the development of the first functional marker in olive and provided positional candidate genes in the SI locus.

Genome-wide association study of berry-related traits in grape [Vitis vinifera L.] based on genotyping-by-sequencing markers

Deciphering the genetic control of grape berry traits is crucial for optimizing yield, fruit quality, and consumer acceptability. In this study, an association panel of 179 grape genotypes comprising a mixture of ancient cultivars, landraces, and modern varieties collected worldwide were genotyped with genotyping-by-sequencing using a genome-wide association approach based on 32,311 single-nucleotide polymorphism (SNP) markers. Genome-wide efficient mixed-model association was selected as the optimal statistical model based on the results of known control loci of grape berry color traits. Many of the associated SNPs identified in this study were in accordance with the previous QTL analyses using biparental mapping. The grape skin color locus was found to be associated with a mybA transcription factor on chromosome 2. Two strong and distinct association signals associated with berry development periods were found on chromosome 16. Most candidate genes of the interval were highlighted as receptor-like protein kinase. For berry weight, significant association loci were identified on chromosome 18, as previously known, and on chromosome 19 and chromosome 17, as newly mapped. Berry flesh texture was newly located on chromosome 16; candidate genes in the interval were related to calcium. Berry flavor was determined on chromosome 5. Genomic regions were further investigated to reveal candidate genes. In this work, we identified interesting genetic determinants of grape berry-related traits. The identification of the markers closely associated with these berry traits may be useful for grape molecular breeding.

Evidence for Strong Kinship Influence on the Extent of Linkage Disequilibrium in Cultivated Common Beans

Phaseolus vulgaris is an important grain legume for human consumption. Recently, association mapping studies have been performed for the species aiming to identify loci underlying quantitative variation of traits. It is now imperative to know whether the linkage disequilibrium (LD) reflects the true association between a marker and causative loci. The aim of this study was to estimate and analyze LD on a diversity panel of common beans using ordinary r² and r2 extensions which correct bias due to population structure (rS²), kinship (rV²), and both (rVS²). A total of 10,362 single nucleotide polymorphisms (SNPs) were identified by genotyping by sequencing (GBS), and polymorphisms were found to be widely distributed along the 11 chromosomes. In terms of r2, high values of LD (over 0.8) were identified between SNPs located at opposite chromosomal ends. Estimates for rV² were lower than those for rS². Results for rV² and rVS² were similar, suggesting that kinship may also include information on population structure. Over genetic distance, LD decayed to 0.1 at a distance of 1 Mb for rVS². Inter-chromosomal LD was also evidenced. This study showed that LD estimates decay dramatically according to the population structure, and especially the degree of kinship. Importantly, the LD estimates reported herein may influence our ability to perform association mapping studies on P. vulgaris.

Indole-3-Carbinol Induces Apoptosis in Human Osteosarcoma MG-63 and U2OS Cells

This study was focused on investigating the anticancer potential of indole-3-carbinol (I3C) against osteosarcoma MG-63 and U2OS cells. A wound healing assay indicated that IC3 inhibited migration of MG-63 and U2OS cells. MTT, WST-1, and colony formation assays revealed that treatment of MG-63 and U2OS cells with I3C decreased cell viability. Fluorescence-activated cell sorting (FACS) analysis showed that I3C induced apoptosis in a dose- and time-dependent manner in MG-63 and U2OS cells. Moreover, via terminal deoxynucleotidyl transferase- (TdT-) mediated dUTP-biotin nick-end labeling (TUNEL) assay, we detected that I3C induced DNA fragmentation. Western blotting demonstrated that activated forms of caspase-3, caspase-7, and caspase-9, as well as poly (ADP-ribose) polymerase (PARP) were increased in MG-63 and U2OS cells, following treatment with I3C. Furthermore, protein expression levels of FOXO3, Bax, and Bim extra-large form were increased while those of Akt, JNK, p38, phosphorylated ERK, and Bcl-xL were decreased by I3C treatment in MG-63 and U2OS cells. Thus, the study indicates that I3C may induce apoptosis in human osteosarcoma MG-63 and U2OS cells via the activation of apoptotic signaling pathways by FOXO3.

Genotype-Environment Interaction Shapes the Microbial Assemblage in Grapevine's Phyllosphere and Carposphere: An NGS Approach

Plant surface or phyllosphere is the habitat of hyperdiverse microbial communities and it is always exposed to the fluctuating environmental factors, which is thought to be one of the potential drivers of microbial community structuring. Impact of grapevine genotypes in variable environmental factors (i.e., at different geographic locations) on the phyllosphere has never been studied and is the main objective of this report. Using high throughput short amplicon sequencing of 16S rRNA genes and internal transcribed spacer (ITS), we analyzed the impacts of genotypes of Vitis Vinifera (coming from three genetic pool), on the microbial (bacterial and fungal) assemblage in the phyllosphere. First, we performed the analysis of the phyllosphere microbiome while using fifteen genotypes that were chosen to maximize intra-specific diversity and grown in two Mediterranean vineyards. Then, the same analysis was performed on five commercially important varieties of Vitis vinifera that were sampled from three different French agro-climatic zones (or terroir: a combination of climate, soils, and human practices). Our study revealed that, at a particular geographic location, genotypes have an impact on microbial assemblage in the phyllosphere and carposphere of leaf and fruit (or berries), respectively, which is more prominent on the carposphere but the effect of terroir was much stronger than the genotype when the leaf phyllosphere of five grapevine varieties grown in different agro-climatic zones was compared. Impacts of the season and exterior plant organs (leaf and berries) on microbial taxa structuring in the phyllosphere was also assessed and presented in this report.

Genomic signatures of different adaptations to environmental stimuli between wild and cultivated Vitis vinifera L

The application of population genetic methods in combination with gene mapping strategies can help to identify genes and mutations selected during the evolution from wild plants to crops and to explore the considerable genetic variation still maintained in natural populations. We genotyped a grapevine germplasm collection of 44 wild (Vitis vinifera subsp. sylvestris) and 48 cultivated (V. vinifera subsp. sativa) accessions at 54 K single-nucleotide polymorphisms (SNPs) to perform a whole-genome comparison of the main population genetic statistics. The analysis of Wright Fixation Index (FST) along the whole genome allowed us to identify several putative "signatures of selection" spanning over two thousand SNPs significantly differentiated between sativa and sylvestris. Many of these genomic regions included genes involved in the adaptation to environmental changes. An overall reduction of nucleotide diversity was observed across the whole genome within sylvestris, supporting a small effective population size of the wild grapevine. Tajima's D resulted positive in both wild and cultivated subgroups, which may indicate an ongoing balancing selection. Association mapping for six domestication-related traits was performed in combination with population genetics, providing further evidence of different perception and response to environmental stresses between sativa and sylvestris.

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.

Methods and Tools in Genome-wide Association Studies

Many traits, such as height, the response to a given drug, or the susceptibility to certain diseases are presumably co-determined by genetics. Especially in the field of medicine, it is of major interest to identify genetic aberrations that alter an individual's risk to develop a certain phenotypic trait. Addressing this question requires the availability of comprehensive, high-quality genetic datasets. The technological advancements and the decreasing cost of genotyping in the last decade led to an increase in such datasets. Parallel to and in line with this technological progress, an analysis framework under the name of genome-wide association studies was developed to properly collect and analyze these data. Genome-wide association studies aim at finding statistical dependencies-or associations-between a trait of interest and point-mutations in the DNA. The statistical models used to detect such associations are diverse, spanning the whole range from the frequentist to the Bayesian setting.Since genetic datasets are inherently high-dimensional, the search for associations poses not only a statistical but also a computational challenge. As a result, a variety of toolboxes and software packages have been developed, each implementing different statistical methods while using various optimizations and mathematical techniques to enhance the computations.This chapter is devoted to the discussion of widely used methods and tools in genome-wide association studies. We present the different statistical models and the assumptions on which they are based, explain peculiarities of the data that have to be accounted for and, most importantly, introduce commonly used tools and software packages for the different tasks in a genome-wide association study, complemented with examples for their application.

The long non-coding RNA MIAT regulates zinc finger E-box binding homeobox 1 expression by sponging miR-150 and promoteing cell invasion in non-small-cell lung cancer

The myocardial infarction associated transcript (MIAT), a long non-coding RNA (lncRNA), was originally identified as a candidate gene for myocardial infarction, and was recently shown to participate in the progression of cancer and the process of metastasis. However, the biological role of MIAT and the underlying mechanisms that mediate its role in non-small-cell lung cancer (NSCLC) remain unclear. Here, we have shown that the expression of MIAT in NSCLC tissues was upregulated. Knockdown of MIAT substantially inhibited the invasive ability of NSCLC cells. Moreover, the knockdown of MIAT significantly downregulated the expression of the zinc finger E-box binding homeobox 1 (ZEB1), that was upregulated in NSCLC and that promoted cell invasion. Rather than by direct interactions, we found that MIAT indirectly regulated ZEB1 expression through sponging and suppressing microRNA (miR)-150, which represses ZEB1 and interacts with MIAT in a sequence-specific manner. Thus, MIAT may inhibit ZEB1 expression and promote cell invasion of NSCLC cells via the miR-150/ZEB1 pathway. Taken together, our findings suggested that MIAT plays an oncogenic role in NSCLC through the ZEB1 signaling pathway by sponging miR-150, and MIAT may therefore serve as a valuable prognostic biomarker and therapeutic target for NSCLC.

Use of Natural Diversity and Biotechnology to Increase the Quality and Nutritional Content of Tomato and Grape

Improving fruit quality has become a major goal in plant breeding. Direct approaches to tackling fruit quality traits specifically linked to consumer preferences and environmental friendliness, such as improved flavor, nutraceutical compounds, and sustainability, have slowly been added to a breeder priority list that already includes traits like productivity, efficiency, and, especially, pest and disease control. Breeders already use molecular genetic tools to improve fruit quality although most advances have been made in producer and industrial quality standards. Furthermore, progress has largely been limited to simple agronomic traits easy-to-observe, whereas the vast majority of quality attributes, specifically those relating to flavor and nutrition, are complex and have mostly been neglected. Fortunately, wild germplasm, which is used for resistance against/tolerance of environmental stresses (including pathogens), is still available and harbors significant genetic variation for taste and health-promoting traits. Similarly, heirloom/traditional varieties could be used to identify which genes contribute to flavor and health quality and, at the same time, serve as a good source of the best alleles for organoleptic quality improvement. Grape (Vitis vinifera L.) and tomato (Solanum lycopersicum L.) produce fleshy, berry-type fruits, among the most consumed in the world. Both have undergone important domestication and selection processes, that have dramatically reduced their genetic variability, and strongly standardized fruit traits. Moreover, more and more consumers are asking for sustainable production, incompatible with the wide range of chemical inputs. In the present paper, we review the genetic resources available to tomato/grape breeders, and the recent technological progresses that facilitate the identification of genes/alleles of interest within the natural or generated variability gene pool. These technologies include omics, high-throughput phenotyping/phenomics, and biotech approaches. Our review also covers a range of technologies used to transfer to tomato and grape those alleles considered of interest for fruit quality. These include traditional breeding, TILLING (Targeting Induced Local Lesions in Genomes), genetic engineering, or NPBT (New Plant Breeding Technologies). Altogether, the combined exploitation of genetic variability and innovative biotechnological tools may facilitate breeders to improve fruit quality tacking more into account the consumer standards and the needs to move forward into more sustainable farming practices.

New Strategies and Tools in Quantitative Genetics: How to Go from the Phenotype to the Genotype

Quantitative genetics has a long history in plants: It has been used to study specific biological processes, identify the factors important for trait evolution, and breed new crop varieties. These classical approaches to quantitative trait locus mapping have naturally improved with technology. In this review, we show how quantitative genetics has evolved recently in plants and how new developments in phenotyping, population generation, sequencing, gene manipulation, and statistics are rejuvenating both the classical linkage mapping approaches (for example, through nested association mapping) as well as the more recently developed genome-wide association studies. These strategies are complementary in most instances, and indeed, one is often used to confirm the results of the other. Despite significant advances, an emerging trend is that the outcome and efficiency of the different approaches depend greatly on the genetic architecture of the trait in the genetic material under study.