Genotype and phenotype data standardization, utilization and integration in the big data era for agricultural sciences.

Large-scale genotype and phenotype data have been increasingly generated to identify genetic markers, understand gene function and evolution and facilitate genomic selection. These datasets hold immense value for both current and future studies, as they are vital for crop breeding, yield improvement and overall agricultural sustainability. However, integrating these datasets from heterogeneous sources presents significant challenges and hinders their effective utilization. We established the Genotype-Phenotype Working Group in November 2021 as a part of the AgBioData Consortium (https://www.agbiodata.org) to review current data types and resources that support archiving, analysis and visualization of genotype and phenotype data to understand the needs and challenges of the plant genomic research community. For 2021-22, we identified different types of datasets and examined metadata annotations related to experimental design/methods/sample collection, etc. Furthermore, we thoroughly reviewed publicly funded repositories for raw and processed data as well as secondary databases and knowledgebases that enable the integration of heterogeneous data in the context of the genome browser, pathway networks and tissue-specific gene expression. Based on our survey, we recommend a need for (i) additional infrastructural support for archiving many new data types, (ii) development of community standards for data annotation and formatting, (iii) resources for biocuration and (iv) analysis and visualization tools to connect genotype data with phenotype data to enhance knowledge synthesis and to foster translational research. Although this paper only covers the data and resources relevant to the plant research community, we expect that similar issues and needs are shared by researchers working on animals. Database URL: https://www.agbiodata.org.

Genome-wide association analysis and pathway enrichment provide insights into the genetic basis of photosynthetic responses to drought stress in Persian walnut.

Uncovering the genetic basis of photosynthetic trait variation under drought stress is essential for breeding climate-resilient walnut cultivars. To this end, we examined photosynthetic capacity in a diverse panel of 150 walnut families (1500 seedlings) from various agro-climatic zones in their habitats and grown in a common garden experiment. Photosynthetic traits were measured under well-watered (WW), water-stressed (WS) and recovery (WR) conditions. We performed genome-wide association studies (GWAS) using three genomic datasets: genotyping by sequencing data (∼43 K SNPs) on both mother trees (MGBS) and progeny (PGBS) and the Axiom™ Juglans regia 700 K SNP array data (∼295 K SNPs) on mother trees (MArray). We identified 578 unique genomic regions linked with at least one trait in a specific treatment, 874 predicted genes that fell within 20 kb of a significant or suggestive SNP in at least two of the three GWAS datasets (MArray, MGBS, and PGBS), and 67 genes that fell within 20 kb of a significant SNP in all three GWAS datasets. Functional annotation identified several candidate pathways and genes that play crucial roles in photosynthesis, amino acid and carbohydrate metabolism, and signal transduction. Further network analysis identified 15 hub genes under WW, WS and WR conditions including GAPB, PSAN, CRR1, NTRC, DGD1, CYP38, and PETC which are involved in the photosynthetic responses. These findings shed light on possible strategies for improving walnut productivity under drought stress.

Valorization of Traditional Italian Walnut (Juglans regia L.) Production: Genetic, Nutritional and Sensory Characterization of Locally Grown Varieties in the Trentino Region.

Juglans regia (L.) is cultivated worldwide for its nutrient-rich nuts. In Italy, despite the growing demand, walnut cultivation has gone through a strong decline in recent decades, which led to Italy being among the top five net importing countries. To promote the development of local high-quality Italian walnut production, we devised a multidisciplinary project to highlight the distinctive traits of three varieties grown in the mountainous region Trentino (northeast of Italy): the heirloom 'Bleggiana', a second local accession called local Franquette and the French cultivar 'Lara', recently introduced in the local production to increase yield. The genetic characterization confirmed the uniqueness of 'Bleggiana' and revealed local Franquette as a newly described autochthonous variety, thus named 'Blegette'. The metabolic profiles highlighted a valuable nutritional composition of the local varieties, richer in polyphenols and with a lower ω-6/ω-3 ratio than the commercial 'Lara'. 'Blegette' obtained the highest preference scores from consumers for both the visual aspect and tasting; however, the volatile organic compound profiles did not discriminate among the characterized cultivars. The described local varieties represent an interesting reservoir of walnut genetic diversity and quality properties, which deserve future investigation on agronomically useful traits (e.g., local adaptation and water usage) for a high-quality and sustainable production.

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.

Genetic Analysis of Walnut (Juglans regia L.) Pellicle Pigment Variation Through a Novel, High-Throughput Phenotyping Platform.

Walnut pellicle color is a key quality attribute that drives consumer preference and walnut sales. For the first time a high-throughput, computer vision-based phenotyping platform using a custom algorithm to quantitatively score each walnut pellicle in L* a* b* color space was deployed at large-scale. This was compared to traditional qualitative scoring by eye and was used to dissect the genetics of pellicle pigmentation. Progeny from both a bi-parental population of 168 trees ('Chandler' × 'Idaho') and a genome-wide association (GWAS) with 528 trees of the UC Davis Walnut Improvement Program were analyzed. Color phenotypes were found to have overlapping regions in the 'Chandler' genetic map on Chr01 suggesting complex genetic control. In the GWAS population, multiple, small effect QTL across Chr01, Chr07, Chr08, Chr09, Chr10, Chr12 and Chr13 were discovered. Marker trait associations were co-localized with QTL mapping on Chr01, Chr10, Chr14, and Chr16. Putative candidate genes controlling walnut pellicle pigmentation were postulated.

High-quality chromosome-scale assembly of the walnut (Juglans regia L.) reference genome.

BACKGROUND: The release of the first reference genome of walnut (Juglans regia L.) enabled many achievements in the characterization of walnut genetic and functional variation. However, it is highly fragmented, preventing the integration of genetic, transcriptomic, and proteomic information to fully elucidate walnut biological processes. FINDINGS: Here, we report the new chromosome-scale assembly of the walnut reference genome (Chandler v2.0) obtained by combining Oxford Nanopore long-read sequencing with chromosome conformation capture (Hi-C) technology. Relative to the previous reference genome, the new assembly features an 84.4-fold increase in N50 size, with the 16 chromosomal pseudomolecules assembled and representing 95% of its total length. Using full-length transcripts from single-molecule real-time sequencing, we predicted 37,554 gene models, with a mean gene length higher than the previous gene annotations. Most of the new protein-coding genes (90%) present both start and stop codons, which represents a significant improvement compared with Chandler v1.0 (only 48%). We then tested the potential impact of the new chromosome-level genome on different areas of walnut research. By studying the proteome changes occurring during male flower development, we observed that the virtual proteome obtained from Chandler v2.0 presents fewer artifacts than the previous reference genome, enabling the identification of a new potential pollen allergen in walnut. Also, the new chromosome-scale genome facilitates in-depth studies of intraspecies genetic diversity by revealing previously undetected autozygous regions in Chandler, likely resulting from inbreeding, and 195 genomic regions highly differentiated between Western and Eastern walnut cultivars. CONCLUSION: Overall, Chandler v2.0 will serve as a valuable resource to better understand and explore walnut biology.

Quantitative phenotyping of shell suture strength in walnut (Juglans regia L.) enhances precision for detection of QTL and genome-wide association mapping.

Walnut shell suture strength directly impacts the ability to maintain shell integrity during harvest and processing, susceptibility to insect damage and other contamination, and the proportion of kernel halves recovered during cracking. Suture strength is therefore an important breeding objective. Here, two methods of phenotyping this trait were investigated: 1) traditional, qualitative and rather subjective scoring on an interval scale by human observers, and; 2) quantitative and continuous measurements captured by a texturometer. The aim of this work was to increase the accuracy of suture strength phenotyping and to then apply two mapping approaches, quantitative trait loci (QTL) mapping and genome wide association (GWAS) models, in order to dissect the genetic basis of the walnut suture trait. Using data collected on trees within the UC Davis Walnut Improvement Program (n = 464), the genetic correlation between the texturometer method and qualitatively scored method was high (0.826). Narrow sense heritability calculated using quantitative measurements was 0.82. A major QTL for suture strength was detected on LG05, explaining 34% of the phenotypic variation; additionally, two minor QTLs were identified on LG01 and LG11. All three QTLs were confirmed with GWAS on corresponding chromosomes. The findings reported in this study are relevant for application towards a molecular breeding program in walnut.

Association and linkage mapping to unravel genetic architecture of phenological traits and lateral bearing in Persian walnut (Juglans regia L.).

BACKGROUND: Unravelling the genetic architecture of agronomic traits in walnut such as budbreak date and bearing habit, is crucial for climate change adaptation and yield improvement. A Genome-Wide Association Study (GWAS) using multi-locus models was conducted in a panel of 170 walnut accessions genotyped using the Axiom™ J. regia 700 K SNP array, with phenological data from 2018, 2019 and legacy data. These accessions come from the INRAE walnut germplasm collection which is the result of important prospecting work performed in many countries around the world. In parallel, an F1 progeny of 78 individuals segregating for phenology-related traits, was genotyped with the same array and phenotyped for the same traits, to construct linkage maps and perform Quantitative Trait Loci (QTLs) detection. RESULTS: Using GWAS, we found strong associations of SNPs located at the beginning of chromosome 1 with both budbreak and female flowering dates. These findings were supported by QTLs detected in the same genomic region. Highly significant associated SNPs were also detected using GWAS for heterodichogamy and lateral bearing habit, both on chromosome 11. We developed a Kompetitive Allele Specific PCR (KASP) marker for budbreak date in walnut, and validated it using plant material from the Walnut Improvement Program of the University of California, Davis, demonstrating its effectiveness for marker-assisted selection in Persian walnut. We found several candidate genes involved in flowering events in walnut, including a gene related to heterodichogamy encoding a sugar catabolism enzyme and a cell division related gene linked to female flowering date. CONCLUSIONS: This study enhances knowledge of the genetic architecture of important agronomic traits related to male and female flowering processes and lateral bearing in walnut. The new marker available for budbreak date, one of the most important traits for good fruiting, will facilitate the selection and development of new walnut cultivars suitable for specific climates.

Stacking up RADSeq assembly programs: From complete hit to completely abysmal.

Decreasing sequencing costs have driven a rapid expansion of novel genotyping methods. One of these methods is the exploitation of restriction enzyme cut sites to generate genome-wide but reduced representation sequencing libraries (RRLs), alternatively termed genotyping by sequencing or restriction-site associated DNA sequencing. Without a reference genome, the resulting short sequence reads must be assembled de novo. There are many possible assembly programs, most not explicitly developed for RRL data, and we know little of their effectiveness. In this issue of Molecular Ecology Resources, LaCava et al. (2020) systematically evaluate six commonly used programs and two commonly varied parameters for complete and accurate assembly of RRLs, using simulated double digests of Homo sapiens and Arabidopsis thaliana genomes with varied mutation rates and types. The authors find substantial variation in performance across assembly programs. The most consistently high-performing assembler is infrequently used in their literature survey (CD-HIT; Li and Godzik, 2006), while several others fail to produce complete, accurate assemblies under many conditions. LaCava et al. additionally recommend best practices in parameter choice and evaluation of future assembly programs-advice that molecular ecologists working to assemble sequences of all kinds should take to heart.

Genome-Wide Association Study Reveals Candidate Genes Involved in Fruit Trait Variation in Persian Walnut (Juglans regia L.).

Elucidating the genetic determinants of fruit quality traits in walnut is essential to breed new cultivars meeting the producers and consumers' needs. We conducted a genome-wide association study (GWAS) using multi-locus models in a panel of 170 accessions of Juglans regia from the INRAE walnut germplasm collection, previously genotyped using the AxiomTM J. regia 700K SNP array. We phenotyped the panel for 25 fruit traits related to morphometrics, shape, volume, weight, ease of cracking, and nutritional composition. We found more than 60 marker-trait associations (MTAs), including a highly significant SNP associated with nut face diameter, nut volume and kernel volume on chromosome 14, and 5 additional associations were detected for walnut weight. We proposed several candidate genes involved in nut characteristics, such as a gene coding for a beta-galactosidase linked to several size-related traits and known to be involved in fruit development in other species. We also confirmed associations on chromosomes 5 and 11 with nut suture strength, recently reported by the University of California, Davis. Our results enhance knowledge of the genetic control of important agronomic traits related to fruit quality in walnut, and pave the way for the development of molecular markers for future assisted selection.

Combining phenotype, genotype, and environment to uncover genetic components underlying water use efficiency in Persian walnut.

Walnut production is challenged by climate change and abiotic stresses. Elucidating the genomic basis of adaptation to climate is essential to breeding drought-tolerant cultivars for enhanced productivity in arid and semi-arid regions. Here, we aimed to identify loci potentially involved in water use efficiency (WUE) and adaptation to drought in Persian walnut using a diverse panel of 95 walnut families (950 seedlings) from Iran, which show contrasting levels of water availability in their native habitats. We analyzed associations between phenotypic, genotypic, and environmental variables from data sets of 609 000 high-quality single nucleotide polymorphisms (SNPs), three categories of phenotypic traits [WUE-related traits under drought, their drought stress index, and principal components (PCs)], and 21 climate variables and their combination (first three PCs). Our genotype-phenotype analysis identified 22 significant and 266 suggestive associations, some of which were for multiple traits, suggesting their correlation and a possible common genetic control. Also, genotype-environment association analysis found 115 significant and 265 suggestive SNP loci that displayed potential signals of local adaptation. Several sets of stress-responsive genes were found in the genomic regions significantly associated with the aforementioned traits. Most of the candidate genes identified are involved in abscisic acid signaling, stomatal regulation, transduction of environmental signals, antioxidant defense system, osmotic adjustment, and leaf growth and development. Upon validation, the marker-trait associations identified for drought tolerance-related traits would allow the selection and development of new walnut rootstocks or scion cultivars with superior WUE.

Deciphering of the Genetic Control of Phenology, Yield, and Pellicle Color in Persian Walnut (Juglans regia L.).

Yield, nut quality, and ability to adapt to specific climate conditions, are all important factors to consider in the development and selection of new Persian walnut (Juglans regia L.) varieties. The genetic control of these traits is still unknown in walnut, limiting the accuracy and rapidity of releasing new cultivars for commercial use. We studied the genetic architecture of five traits crucial for either marketing (i.e., yield, lateral fruit-bearing, and pellicle color) or selection of individuals with specific phenology (i.e., leafing and harvest date). By combining over 30 years of phenotypic data with genetic profiles generated using the latest Axiom™ J. regia 700K SNP array, we were able to identify and confirm major loci for all these traits. In particular, we revealed that a genomic region at the beginning of Chr1 controls both leafing and harvest date in walnut, consistent with the observed strong phenotypical correlation between these traits, and including candidate genes involved in plant development, leaf formation, and cell division. In addition, a large genomic region on Chr11 that includes genes with a central role in flowering control and shoot meristem growth underlies both lateral fruit-bearing and yield in walnut. We observed a more complex genetic architecture for pellicle color, strongly influenced by the environment (h 2 = 0.43). We identified two marker-trait associations on Chr6 and 7 for pellicle color, where genes encoding a UDP-glycosyltransferase or involved in the response to oxidation were found. In conclusion, by combining classical quantitative trait loci (QTL) mapping and genome-wide association mapping, we deciphered, for the first time, the molecular pathways controlling walnut phenology, yield, lateral fruitfulness, and pellicle color. Our findings represent a further milestone in the transition from conventional to genome-assisted breeding in Persian walnut.

Genome-wide patterns of population structure and association mapping of nut-related traits in Persian walnut populations from Iran using the Axiom J. regia 700K SNP array.

Persian plateau (including Iran) is considered as one of the primary centers of origin of walnut. Sampling walnut trees originating from this arena and exploiting the capabilities of next-generation sequencing (NGS) can provide new insights into the degree of genetic variation across the walnut genome. The present study aimed to explore the population structure and genomic variation of an Iranian collection of Persian walnut (Juglans regia L.) and identify loci underlying the variation in nut and kernel related traits using the new Axiom J. regia 700K SNP genotyping array. We genotyped a diversity panel including 95 walnut genotypes from eight Iranian provinces with a variety of climate zones. A majority of the SNPs (323,273, 53.03%) fell into the "Poly High Resolution" class of polymorphisms, which includes the highest quality variants. Genetic structure assessment, using several approaches, divided the Iranian walnut panel into four principal clusters, reflecting their geographic partitioning. We observed high genetic variation across all of the populations (HO = 0.34 and HE = 0.38). The overall level of genetic differentiation among populations was moderate (FST = 0.07). However, the Semnan population showed high divergence from the other Iranian populations (on average FST = 0.12), most likely due to its geographical isolation. Based on parentage analysis, the level of relatedness was very low among the Iranian walnuts examined, reflecting the geographical distance between the Iranian provinces considered in our study. Finally, we performed a genome-wide association study (GWAS), identifying 55 SNPs significantly associated with nut and kernel-related traits. In conclusion, by applying the novel Axiom J. regia 700K SNP array we uncovered new unexplored genetic diversity and identified significant marker-trait associations for nut-related traits in Persian walnut that will be useful for future breeding programs in Iran and other countries.

A new genomic tool for walnut (Juglans regia L.): development and validation of the high-density Axiom™ J. regia 700K SNP genotyping array.

Over the last 20 years, global production of Persian walnut (Juglans regia L.) has grown enormously, likely reflecting increased consumption due to its numerous benefits to human health. However, advances in genome-wide association (GWA) studies and genomic selection (GS) for agronomically important traits in walnut remain limited due to the lack of powerful genomic tools. Here, we present the development and validation of a high-density 700K single nucleotide polymorphism (SNP) array in Persian walnut. Over 609K high-quality SNPs have been thoroughly selected from a set of 9.6 m genome-wide variants, previously identified from the high-depth re-sequencing of 27 founders of the Walnut Improvement Program (WIP) of University of California, Davis. To validate the effectiveness of the array, we genotyped a collection of 1284 walnut trees, including 1167 progeny of 48 WIP families and 26 walnut cultivars. More than half of the SNPs (55.7%) fell in the highest quality class of 'Poly High Resolution' (PHR) polymorphisms, which were used to assess the WIP pedigree integrity. We identified 151 new parent-offspring relationships, all confirmed with the Mendelian inheritance test. In addition, we explored the genetic variability among cultivars of different origin, revealing how the varieties from Europe and California were differentiated from Asian accessions. Both the reconstruction of the WIP pedigree and population structure analysis confirmed the effectiveness of the Applied Biosystems™ Axiom™ J. regia 700K SNP array, which initiates a novel genomic and advanced phase in walnut genetics and breeding.

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.

Collection and characterization of grapevine genetic resources (Vitis vinifera) in the Holy Land, towards the renewal of ancient winemaking practices.

The importance and extent of wine consumption in all life aspects at the Holy Land is well documented. The Muslim influence in this region led to the abandonment of winemaking practices, and possible loss of indigenous wine varieties. Here we present a country wide collection of the local grapevine population including wild and cultivated forms, and its characterization by genetic, ampelographic and enological methods. The ampelographic analysis shows clear differences between Sativa and Sylvestris groups in flower, leaf and cluster parameters, and that most Sativa belong to proles orientalis. Genetic population analysis was conducted by analyzing 22 common SSR markers, determining first the unique genotypes, and internally assessing the population's structure, showing the existence of two distinct Sativa and Sylvestris populations, and a third mixed one. Likewise, the relationship between the Israeli grapevine population and grapevine populations in Europe and parts of Asia was investigated, showing that the Israeli Sativa and Sylvestris populations cluster closely together, suggesting a common genetic source. Lastly, the enological characteristics of selected Sativa and Sylvestris genotypes are presented, demonstrating their potential for quality wine production. This research significantly contributes toward the re-establishment of indigenous and traditional local grapevine varieties into the modern international wine industry.

SNP-Discovery by RAD-Sequencing in a Germplasm Collection of Wild and Cultivated Grapevines (V. vinifera L.).

Whole-genome comparisons of Vitis vinifera subsp. sativa and V. vinifera subsp. sylvestris are expected to provide a better estimate of the valuable genetic diversity still present in grapevine, and help to reconstruct the evolutionary history of a major crop worldwide. To this aim, the increase of molecular marker density across the grapevine genome is fundamental. Here we describe the SNP discovery in a grapevine germplasm collection of 51 cultivars and 44 wild accessions through a novel protocol of restriction-site associated DNA (RAD) sequencing. By resequencing 1.1% of the grapevine genome at a high coverage, we recovered 34K BamHI unique restriction sites, of which 6.8% were absent in the 'PN40024' reference genome. Moreover, we identified 37,748 single nucleotide polymorphisms (SNPs), 93% of which belonged to the 19 assembled chromosomes with an average of 1.8K SNPs per chromosome. Nearly half of the SNPs fell in genic regions mostly assigned to the functional categories of metabolism and regulation, whereas some nonsynonymous variants were identified in genes related with the detection and response to environmental stimuli. SNP validation was carried-out, showing the ability of RAD-seq to accurately determine genotypes in a highly heterozygous species. To test the usefulness of our SNP panel, the main diversity statistics were evaluated, highlighting how the wild grapevine retained less genetic variability than the cultivated form. Furthermore, the analysis of Linkage Disequilibrium (LD) in the two subspecies separately revealed how the LD decays faster within the domesticated grapevine compared to its wild relative. Being the first application of RAD-seq in a diverse grapevine germplasm collection, our approach holds great promise for exploiting the genetic resources available in one of the most economically important fruit crops.