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

Complete mitochondrial genome assembly of Juglans regia unveiled its molecular characteristics, genome evolution, and phylogenetic implications

BACKGROUND

The Persian walnut (Juglans regia), an economically vital species within the Juglandaceae family, has seen its mitochondrial genome sequenced and assembled in the current study using advanced Illumina and Nanopore sequencing technology.

RESULTS

The 1,007,576 bp mitogenome of J. regia consisted of three circular chromosomes with a 44.52% GC content encoding 39 PCGs, 47 tRNA, and five rRNA genes. Extensive repetitive sequences, including 320 SSRs, 512 interspersed, and 83 tandem repeats, were identified, contributing to genomic complexity. The protein-coding sequences (PCGs) favored A/T-ending codons, and the codon usage bias was primarily shaped by selective pressure. Intracellular gene transfer occurred among the mitogenome, chloroplast, and nuclear genomes. Comparative genomic analysis unveiled abundant structure and sequence variation among J. regia and related species. The results of selective pressure analysis indicated that most PCGs underwent purifying selection, whereas the atp4 and ccmB genes had experienced positive selection between many species pairs. In addition, the phylogenetic examination, grounded in mitochondrial genome data, precisely delineated the evolutionary and taxonomic relationships of J. regia and its relatives. We identified a total of 539 RNA editing sites, among which 288 were corroborated by transcriptome sequencing data. Furthermore, expression profiling under temperature stress highlighted the complex regulation pattern of 28 differently expressed PCGs, wherein NADH dehydrogenase and ATP synthase genes might be critical in the mitochondria response to cold stress.

CONCLUSIONS

Our results provided valuable molecular resources for understanding the genetic characteristics of J. regia and offered novel perspectives for population genetics and evolutionary studies in Juglans and related woody species.

A multiplexed plant-animal SNP array for selective breeding and species conservation applications

Reliable and high-throughput genotyping platforms are of immense importance for identifying and dissecting genomic regions controlling important phenotypes, supporting selection processes in breeding programs, and managing wild populations and germplasm collections. Amongst available genotyping tools, single nucleotide polymorphism arrays have been shown to be comparatively easy to use and generate highly accurate genotypic data. Single-species arrays are the most commonly used type so far; however, some multi-species arrays have been developed for closely related species that share single nucleotide polymorphism markers, exploiting inter-species cross-amplification. In this study, the suitability of a multiplexed plant-animal single nucleotide polymorphism array, including both closely and distantly related species, was explored. The performance of the single nucleotide polymorphism array across species for diverse applications, ranging from intra-species diversity assessments to parentage analysis, was assessed. Moreover, the value of genotyping pooled DNA of distantly related species on the single nucleotide polymorphism array as a technique to further reduce costs was evaluated. Single nucleotide polymorphism performance was generally high, and species-specific single nucleotide polymorphisms proved suitable for diverse applications. The multi-species single nucleotide polymorphism array approach reported here could be transferred to other species to achieve cost savings resulting from the increased throughput when several projects use the same array, and the pooling technique adds another highly promising advancement to additionally decrease genotyping costs by half.

Genome-wide characterization and development of SSR markers for genetic diversity analysis in northwestern Himalayas Walnut (Juglans regia L.)

In the present study, we designed and validated genome-wide polymorphic SSR markers (110 SSRs) by mining the walnut genome. A total of 198,924 SSR loci were identified. Among these, successful primers were designed for 162,594 (81.73%) SSR loci. Dinucleotides were the most predominant accounting for 88.40% (175,075) of total SSRs. The SSR frequency was 377.312 SSR/Mb and it showed a decreasing trend from dinucleotide to octanucleotide motifs. We identified 20 highly polymorphic SSR markers and used them to genotype 72 walnut accessions. Over all, we obtained 118 alleles that ranged from 2 to 12 with an average value of 5.9. The higher SSR PIC values indicate their robustness in discriminating walnut genotypes. Heat map, PCA, and population structure categorized 72 walnut genotypes into 2 distinct clusters. The genetic variation within population was higher than among population as inferred by analysis of molecular variance (AMOVA). For walnut improvement, it is necessary to have a large repository of SSRs with high discriminative power. The present study reports 150,000 SSRs, which is the largest SSR repository for this important nut crop. Scientific communities may use this repository for walnut improvement such as QTL mapping, genetic studies, linkage map construction, and marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03563-6.

Large-scale population structure and genetic architecture of agronomic traits of garlic

Garlic, an asexually propagated crop, is the second important bulb crop after the onion and is used as a vegetable and medicinal plant. Abundant and diverse garlic resources have been formed over thousands of years of cultivation. However, genome variation, population structure and genetic architecture of garlic agronomic traits were still not well elucidated. Here, 1 100 258 single nucleotide polymorphisms (SNPs) were identified using genotyping-by-sequencing in 606 garlic accessions collected from 43 countries. Population structure, principal component and phylogenetic analysis showed that these accessions were divided into five subpopulations. Twenty agronomic traits, including above-ground growth traits, bulb-related and bolt-related traits in two consecutive years were implemented in a genome-wide association study. In total, 542 SNPs were associated with these agronomic traits, among which 188 SNPs were repeatedly associated with more than two traits. One SNP (chr6: 1896135972) was repeatedly associated with ten traits. These associated SNPs were located within or near 858 genes, 56 of which were transcription factors. Interestingly, one non-synonymous SNP (Chr4: 166524085) in ribosomal protein S5 was repeatedly associated with above-ground growth and bulb-related traits. Additionally, gene ontology enrichment analysis of candidate genes for genomic selection regions between complete-bolting and non-bolting accessions showed that these genes were significantly enriched in 'vegetative to reproductive phase transition of meristem', 'shoot system development', 'reproductive process', etc. These results provide valuable information for the reliable and efficient selection of candidate genes to achieve garlic genetic improvement and superior varieties.

Development and Evaluation of an AxiomTM 60K SNP Array for Almond (Prunus dulcis)

A high-density single nucleotide polymorphism (SNP) array is essential to enable faster progress in plant breeding for new cultivar development. In this regard, we have developed an Axiom 60K almond SNP array by resequencing 81 almond accessions. For the validation of the array, a set of 210 accessions were genotyped and 82.8% of the SNPs were classified in the best recommended SNPs. The rate of missing data was between 0.4% and 2.7% for the almond accessions and less than 15.5% for the few peach and wild accessions, suggesting that this array can be used for peach and interspecific peach × almond genetic studies. The values of the two SNPs linked to the RMja (nematode resistance) and SK (bitterness) genes were consistent. We also genotyped 49 hybrids from an almond F2 progeny and could build a genetic map with a set of 1159 SNPs. Error rates, less than 1%, were evaluated by comparing replicates and by detection of departures from Mendelian inheritance in the F2 progeny. This almond array is commercially available and should be a cost-effective genotyping tool useful in the search for new genes and quantitative traits loci (QTL) involved in the control of agronomic traits.

THE WALNUT - CONSTRAINTS AND ADVANTAGES FOR A SUSTAINABLE DEVELOPMENT

The Persian walnut or English walnut (Juglans regia L.) is a commonly grown species for nut production and noble wood. The nut is one of the oldest food and traditional medicine sources. The native and commercial walnut genotypes present a large diversity that differ widely in nut productivity and quality. However, genetic erosion poses a serious threat to this tree. Several researches of walnut genetic diversity are being carried out utilizing morphological, biochemical, and molecular approaches in order to select superior walnut cultivars of different agroclimatic areas to increase nut production and quality. Genetic resource evaluation and agrodiversity conservation have a major role in ensuring food security for future generations through a continuous supply of new rootstocks and improved cultivars.

Genome-wide association analysis of 101 accessions dissects the genetic basis of shell thickness for genetic improvement in Persian walnut (Juglans regia L.)

BACKGROUND

Understanding the underlying genetic mechanisms that drive phenotypic variations is essential for enhancing the efficacy of crop improvement. Persian walnut (Juglans regia L.), which is grown extensively worldwide, is an important economic tree fruit due to its horticultural, medicinal, and material value. The quality of the walnut fruit is related to the selection of traits such as thinner shells, larger filling rates, and better taste, which is very important for breeding in China. The complex quantitative fruit-related traits are influenced by a variety of physiological and environmental factors, which can vary widely between walnut genotypes.

RESULTS

For this study, a set of 101 Persian walnut accessions were re-sequenced, which generated a total of 906.2 Gb of Illumina sequence data with an average read depth of 13.8× for each accession. We performed the genome-wide association study (GWAS) using 10.9 Mb of high-quality single-nucleotide polymorphisms (SNPs) and 10 agronomic traits to explore the underlying genetic basis of the walnut fruit. Several candidate genes are proposed to be involved in walnut characteristics, including JrPXC1, JrWAKL8, JrGAMYB, and JrFRK1. Specifically, the JrPXC1 gene was confirmed to participate in the regulation of secondary wall cellulose thickening in the walnut shell.

CONCLUSION

In addition to providing considerable available genetic resources for walnut trees, this study revealed the underlying genetic basis involved in important walnut agronomic traits, particularly shell thickness, as well as providing clues for the improvement of genetic breeding and domestication in other perennial economic crops.

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.

Population-genomic analyses reveal bottlenecks and asymmetric introgression from Persian into iron walnut during domestication

Background

Persian walnut, Juglans regia, occurs naturally from Greece to western China, while its closest relative, the iron walnut, Juglans sigillata, is endemic in southwest China; both species are cultivated for their nuts and wood. Here, we infer their demographic histories and the time and direction of possible hybridization and introgression between them.

Results

We use whole-genome resequencing data, different population-genetic approaches (PSMC and GONE), and isolation-with-migration models (IMa3) on individuals from Europe, Iran, Kazakhstan, Pakistan, and China. IMa3 analyses indicate that the two species diverged from each other by 0.85 million years ago, with unidirectional gene flow from eastern J. regia and its ancestor into J. sigillata, including the shell-thickness gene. Within J. regia, a western group, located from Europe to Iran, and an eastern group with individuals from northern China, experienced dramatically declining population sizes about 80 generations ago (roughly 2400 to 4000 years), followed by an expansion at about 40 generations, while J. sigillata had a constant population size from about 100 to 20 generations ago, followed by a rapid decline.

Conclusions

Both J. regia and J. sigillata appear to have suffered sudden population declines during their domestication, suggesting that the bottleneck scenario of plant domestication may well apply in at least some perennial crop species. Introgression from introduced J. regia appears to have played a role in the domestication of J. sigillata.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02720-z.

Dissection of the practical soybean breeding pipeline by developing ZDX1, a high-throughput functional array

KEY MESSAGE

We developed the ZDX1 high-throughput functional soybean array for high accuracy evaluation and selection of both parents and progeny, which can greatly accelerate soybean breeding. Microarray technology facilitates rapid, accurate, and economical genotyping. Here, using resequencing data from 2214 representative soybean accessions, we developed the high-throughput functional array ZDX1, containing 158,959 SNPs, covering 90.92% of soybean genes and sites related to important traits. By application of the array, a total of 817 accessions were genotyped, including three subpopulations of candidate parental lines, parental lines and their progeny from practical breeding. The fixed SNPs were identified in progeny, indicating artificial selection during the breeding process. By identifying functional sites of target traits, novel soybean cyst nematode-resistant progeny and maturity-related novel sources were identified by allele combinations, demonstrating that functional sites provide an efficient method for the rapid screening of desirable traits or gene sources. Notably, we found that the breeding index (BI) was a good indicator for progeny selection. Superior progeny were derived from the combination of distantly related parents, with at least one parent having a higher BI. Furthermore, new combinations based on good performance were proposed for further breeding after excluding redundant and closely related parents. Genomic best linear unbiased prediction (GBLUP) analysis was the best analysis method and achieved the highest accuracy in predicting four traits when comparing SNPs in genic regions rather than whole genomic or intergenic SNPs. The prediction accuracy was improved by 32.1% by using progeny to expand the training population. Collectively, a versatile assay demonstrated that the functional ZDX1 array provided efficient information for the design and optimization of a breeding pipeline for accelerated soybean breeding.

Taming the massive genome of Scots pine with PiSy50k, a new genotyping array for conifer research

Pinus sylvestris (Scots pine) is the most widespread coniferous tree in the boreal forests of Eurasia, with major economic and ecological importance. However, its large and repetitive genome presents a challenge for conducting genome-wide analyses such as association studies, genetic mapping and genomic selection. We present a new 50K single-nucleotide polymorphism (SNP) genotyping array for Scots pine research, breeding and other applications. To select the SNP set, we first genotyped 480 Scots pine samples on a 407 540 SNP screening array and identified 47 712 high-quality SNPs for the final array (called 'PiSy50k'). Here, we provide details of the design and testing, as well as allele frequency estimates from the discovery panel, functional annotation, tissue-specific expression patterns and expression level information for the SNPs or corresponding genes, when available. We validated the performance of the PiSy50k array using samples from Finland and Scotland. Overall, 39 678 (83.2%) SNPs showed low error rates (mean = 0.9%). Relatedness estimates based on array genotypes were consistent with the expected pedigrees, and the level of Mendelian error was negligible. In addition, array genotypes successfully discriminate between Scots pine populations of Finnish and Scottish origins. The PiSy50k SNP array will be a valuable tool for a wide variety of future genetic studies and forestry applications.

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.

A genome variation map provides insights into the genetics of walnut adaptation and agronomic traits

BACKGROUND

Common walnut (Juglans regia L.) is one of the top four most consumed nuts in the world due to its health benefits and pleasant taste. Despite its economic importance, the evolutionary history and genetic control of its adaptation and agronomic traits remain largely unexplored.

RESULTS

We report a comprehensive walnut genomic variation map based on whole-genome resequencing of 815 walnut accessions. Evolutionary analyses suggest that Chinese J. regia diverged from J. sigillata with extensive hybridizations after the split of the two species. In contrast to annual crops, the genetic diversity and heterozygous deleterious mutations of Chinese common walnut trees have continued to increase during the improvement process. Selective sweep analyses identify 902 genes uniquely selected in the improved common walnut compared to its progenitor population. Five major-effect loci are identified to be involved in walnut adaptations to temperature, precipitation, and altitude. Genome-wide association studies reveal 27 genomic loci responsible for 18 important agronomic traits, among which JrFAD2 and JrANR are the potentially major-effect causative genes controlling linoleic acid content and color of the endopleura of the nut, respectively.

CONCLUSIONS

The largest genomic resource for walnuts to date has been generated and explored in this study, unveiling their evolutionary history and cracking the genetic code for agronomic traits and environmental adaptation of this economically crucial crop tree.

Transcriptome and proteome analysis of walnut (Juglans regia L.) fruit in response to infection by Colletotrichum gloeosporioides

BACKGROUND

Walnut anthracnose induced by Colletotrichum gloeosporioides is a disastrous disease affecting walnut production. The resistance of walnut fruit to C. gloeosporioides is a highly complicated and genetically programmed process. However, the underlying mechanisms have not yet been elucidated.

RESULTS

To understand the molecular mechanism underlying the defense of walnut to C. gloeosporioides, we used RNA sequencing and label-free quantitation technologies to generate transcriptomic and proteomic profiles of tissues at various lifestyle transitions of C. gloeosporioides, including 0 hpi, pathological tissues at 24 hpi, 48 hpi, and 72 hpi, and distal uninoculated tissues at 120 hpi, in anthracnose-resistant F26 fruit bracts and anthracnose-susceptible F423 fruit bracts, which were defined through scanning electron microscopy. A total of 21,798 differentially expressed genes (DEGs) and 1929 differentially expressed proteins (DEPs) were identified in F26 vs. F423 at five time points, and the numbers of DEGs and DEPs were significantly higher in the early infection stage. Using pairwise comparisons and weighted gene co-expression network analysis of the transcriptome, we identified two modules significantly related to disease resistance and nine hub genes in the transcription expression gene networks. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis of the DEGs and DEPs revealed that many genes were mainly related to immune response, plant hormone signal transduction, and secondary metabolites, and many DEPs were involved in carbon metabolism and photosynthesis. Correlation analysis between the transcriptome data and proteome data also showed that the consistency of the differential expression of the mRNA and corresponding proteins was relatively higher in the early stage of infection.

CONCLUSIONS

Collectively, these results help elucidate the molecular response of walnut fruit to C. gloeosporioides and provide a basis for the genetic improvement of walnut disease resistance.

Integration of Infinium and Axiom SNP array data in the outcrossing species Malus × domestica and causes for seemingly incompatible calls

Background

Single nucleotide polymorphism (SNP) array technology has been increasingly used to generate large quantities of SNP data for use in genetic studies. As new arrays are developed to take advantage of new technology and of improved probe design using new genome sequence and panel data, a need to integrate data from different arrays and array platforms has arisen. This study was undertaken in view of our need for an integrated high-quality dataset of Illumina Infinium® 20 K and Affymetrix Axiom® 480 K SNP array data in apple (Malus × domestica). In this study, we qualify and quantify the compatibility of SNP calling, defined as SNP calls that are both accurate and concordant, across both arrays by two approaches. First, the concordance of SNP calls was evaluated using a set of 417 duplicate individuals genotyped on both arrays starting from a set of 10,295 robust SNPs on the Infinium array. Next, the accuracy of the SNP calls was evaluated on additional germplasm (n = 3141) from both arrays using Mendelian inconsistent and consistent errors across thousands of pedigree links. While performing this work, we took the opportunity to evaluate reasons for probe failure and observed discordant SNP calls.

Results

Concordance among the duplicate individuals was on average of 97.1% across 10,295 SNPs. Of these SNPs, 35% had discordant call(s) that were further curated, leading to a final set of 8412 (81.7%) SNPs that were deemed compatible. Compatibility was highly influenced by the presence of alternate probe binding locations and secondary polymorphisms. The impact of the latter was highly influenced by their number and proximity to the 3′ end of the probe.

Conclusions

The Infinium and Axiom SNP array data were mostly compatible. However, data integration required intense data filtering and curation. This work resulted in a workflow and information that may be of use in other data integration efforts. Such an in-depth analysis of array concordance and accuracy as ours has not been previously described in the literature and will be useful in future work on SNP array data integration and interpretation, and in probe/platform development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07565-7.

Diversity of Chemical Structures and Biosynthesis of Polyphenols in Nut-Bearing Species

Nuts, such as peanut, almond, and chestnut, are valuable food crops for humans being important sources of fatty acids, vitamins, minerals, and polyphenols. Polyphenols, such as flavonoids, stilbenoids, and hydroxycinnamates, represent a group of plant-specialized (secondary) metabolites which are characterized as health-beneficial antioxidants within the human diet as well as physiological stress protectants within the plant. In food chemistry research, a multitude of polyphenols contained in culinary nuts have been studied leading to the identification of their chemical properties and bioactivities. Although functional elucidation of the biosynthetic genes of polyphenols in nut species is crucially important for crop improvement in the creation of higher-quality nuts and stress-tolerant cultivars, the chemical diversity of nut polyphenols and the key biosynthetic genes responsible for their production are still largely uncharacterized. However, current technical advances in whole-genome sequencing have facilitated that nut plant species became model plants for omics-based approaches. Here, we review the chemical diversity of seed polyphenols in majorly consumed nut species coupled to insights into their biological activities. Furthermore, we present an example of the annotation of key genes involved in polyphenolic biosynthesis in peanut using comparative genomics as a case study outlining how we are approaching omics-based approaches of the nut plant species.

Genetic Mapping by Integration of 55K SNP Array and KASP Markers Reveals Candidate Genes for Important Agronomic Traits in Hexaploid Wheat

Agronomic traits such as heading date (HD), plant height (PH), thousand grain weight (TGW), and spike length (SL) are important factors affecting wheat yield. In this study, we constructed a high-density genetic linkage map using the Wheat55K SNP Array to map quantitative trait loci (QTLs) for these traits in 207 recombinant inbred lines (RILs). A total of 37 QTLs were identified, including 9 QTLs for HD, 7 QTLs for PH, 12 QTLs for TGW, and 9 QTLs for SL, which explained 3.0-48.8% of the phenotypic variation. Kompetitive Allele Specific PCR (KASP) markers were developed based on sequencing data and used for validation of the stably detected QTLs on chromosomes 3A, 4B and 6A using 400 RILs. A QTL cluster on chromosome 4B for PH and TGW was delimited to a 0.8 Mb physical interval explaining 12.2-22.8% of the phenotypic variation. Gene annotations and analyses of SNP effects suggested that a gene encoding protein Photosynthesis Affected Mutant 68, which is essential for photosystem II assembly, is a candidate gene affecting PH and TGW. In addition, the QTL for HD on chromosome 3A was narrowed down to a 2.5 Mb interval, and a gene encoding an R3H domain-containing protein was speculated to be the causal gene influencing HD. The linked KASP markers developed in this study will be useful for marker-assisted selection in wheat breeding, and the candidate genes provide new insight into genetic study for those traits in wheat.

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.

New resources for genetic studies in maize (Zea mays L.): a genome-wide Maize6H-60K single nucleotide polymorphism array and its application

Despite the availability of numerous molecular markers in maize, effective evaluation of all types of germplasm resources, accurate identification of varieties and analysis of a large number of materials in a timely, low-cost manner is challenging. Here, we present Maize6H-60K, a genome-wide single nucleotide polymorphism (SNP) array to facilitate maize genotyping. We first identified 160 million variants by sequencing data of 388 representative inbreds and then tiled 200 000 high-quality variants on a screening array. These variants were further narrowed down to 61 282 using stringent filtering criteria. Among the 60 000 markers, 21 460 SNPs (35%) were within genic regions and 12 835 (21%) were located in coding regions. To assess their effectiveness, 329 inbreds, 221 hybrids, 34 parent-offspring sets and six breeding samples were genotyped. Overall, 48 972 SNPs (80%) were categorized into the highest quality class, that of 'poly high resolution'. A total of 54 658 (89.29%) and 53 091 (86.73%) SNPs had minor allele frequency values ≥ 0.20 in inbreds and hybrids respectively. A linkage disequilibrium (LD) analysis revealed that LD decline was in equilibrium when r² was between 0.10 and 0.15, which corresponds to a physical distance of 400-600 kb. UPGMA clustering analysis divided the 329 inbred lines into nine groups that were consistent with known pedigrees. A background analysis of breeding materials indicated that the 60 000 markers were suitable for evaluation of breeding populations constructed by materials between or within heterotic groups. The developed Maize6H-60K array should be an important tool in maize genetic studies, variety identification and molecular breeding.

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.

Genome-wide identification and characterization of long non-coding RNAs conferring resistance to Colletotrichum gloeosporioides in walnut (Juglans regia)

BACKGROUND

Walnut anthracnose caused by Colletotrichum gloeosporioides (Penz.) Penz. and Sacc. is an important walnut production problem in China. Although the long non-coding RNAs (lncRNAs) are important for plant disease resistance, the molecular mechanisms underlying resistance to C. gloeosporioides in walnut remain poorly understood.

RESULTS

The anthracnose-resistant F26 fruits from the B26 clone and the anthracnose-susceptible F423 fruits from the 4-23 clone of walnut were used as the test materials. Specifically, we performed a comparative transcriptome analysis of F26 and F423 fruit bracts to identify differentially expressed LncRNAs (DELs) at five time-points (tissues at 0 hpi, pathological tissues at 24 hpi, 48 hpi, 72 hpi, and distal uninoculated tissues at 120 hpi). Compared with F423, a total of 14,525 DELs were identified, including 10,645 upregulated lncRNAs and 3846 downregulated lncRNAs in F26. The number of upregulated lncRNAs in F26 compared to in F423 was significantly higher at the early stages of C. gloeosporioides infection. A total of 5 modules related to disease resistance were screened by WGCNA and the target genes of lncRNAs were obtained. Bioinformatic analysis showed that the target genes of upregulated lncRNAs were enriched in immune-related processes during the infection of C. gloeosporioides, such as activation of innate immune response, defense response to bacterium, incompatible interaction and immune system process, and enriched in plant hormone signal transduction, phenylpropanoid biosynthesis and other pathways. And 124 known target genes for 96 hub lncRNAs were predicted, including 10 known resistance genes. The expression of 5 lncRNAs and 5 target genes was confirmed by qPCR, which was consistent with the RNA-seq data.

CONCLUSIONS

The results of this study provide the basis for future functional characterizations of lncRNAs regarding the C. gloeosporioides resistance of walnut fruit bracts.

Pedigree analysis of 220 almond genotypes reveals two world mainstream breeding lines based on only three different cultivars

Loss of genetic variability is an increasing challenge in tree breeding programs due to the repeated use of a reduced number of founder genotypes. However, in almond, little is known about the genetic variability in current breeding stocks, although several cases of inbreeding depression have been reported. To gain insights into the genetic structure in modern breeding programs worldwide, marker-verified pedigree data of 220 almond cultivars and breeding selections were analyzed. Inbreeding coefficients, pairwise relatedness, and genetic contribution were calculated for these genotypes. The results reveal two mainstream breeding lines based on three cultivars: "Tuono", "Cristomorto", and "Nonpareil". Descendants from "Tuono" or "Cristomorto" number 76 (sharing 34 descendants), while "Nonpareil" has 71 descendants. The mean inbreeding coefficient of the analyzed genotypes was 0.041, with 14 genotypes presenting a high inbreeding coefficient, over 0.250. Breeding programs from France, the USA, and Spain showed inbreeding coefficients of 0.075, 0.070, and 0.037, respectively. According to their genetic contribution, modern cultivars from Israel, France, the USA, Spain, and Australia trace back to a maximum of six main founding genotypes. Among the group of 65 genotypes carrying the S_f allele for self-compatibility, the mean relatedness coefficient was 0.125, with "Tuono" as the main founding genotype (24.7% of total genetic contribution). The results broaden our understanding about the tendencies followed in almond breeding over the last 50 years and will have a large impact into breeding decision-making process worldwide. Increasing current genetic variability is required in almond breeding programs to assure genetic gain and continuing breeding progress.

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.

Development of a highly efficient 50K single nucleotide polymorphism genotyping array for the large and complex genome of Norway spruce (Picea abies L. Karst) by whole genome resequencing and its transferability to other spruce species

Norway spruce (Picea abies L. Karst) is one of the most important forest tree species with significant economic and ecological impact in Europe. For decades, genomic and genetic studies on Norway spruce have been challenging due to the large and repetitive genome (19.6 Gb with more than 70% being repetitive). To accelerate genomic studies, including population genetics, genome‐wide association studies (GWAS) and genomic selection (GS), in Norway spruce and related species, we here report on the design and performance of a 50K single nucleotide polymorphism (SNP) genotyping array for Norway spruce. The array is developed based on whole genome resequencing (WGS), making it the first WGS‐based SNP array in any conifer species so far. After identifying SNPs using genome resequencing data from 29 trees collected in northern Europe, we adopted a two‐step approach to design the array. First, we built a 450K screening array and used this to genotype a population of 480 trees sampled from both natural and breeding populations across the Norway spruce distribution range. These samples were then used to select high‐confidence probes that were put on the final 50K array. The SNPs selected are distributed over 45,552 scaffolds from the P. abies version 1.0 genome assembly and target 19,954 unique gene models with an even coverage of the 12 linkage groups in Norway spruce. We show that the array has a 99.5% probe specificity, >98% Mendelian allelic inheritance concordance, an average sample call rate of 96.30% and an SNP call rate of 98.90% in family trios and haploid tissues. We also observed that 23,797 probes (50%) could be identified with high confidence in three other spruce species (white spruce [Picea glauca], black spruce [P. mariana] and Sitka spruce [P. sitchensis]). The high‐quality genotyping array will be a valuable resource for genetic and genomic studies in Norway spruce as well as in other conifer species of the same genus.

Comparative Proteomic Analysis of Walnut (Juglans regia L.) Pellicle Tissues Reveals the Regulation of Nut Quality Attributes

Walnuts (Juglans regia L.) are a valuable dietary source of polyphenols and lipids, with increasing worldwide consumption. California is a major producer, with 'Chandler' and 'Tulare' among the cultivars more widely grown. 'Chandler' produces kernels with extra light color at a higher frequency than other cultivars, gaining preference by growers and consumers. Here we performed a deep comparative proteome analysis of kernel pellicle tissue from these two valued genotypes at three harvest maturities, detecting a total of 4937 J. regia proteins. Late and early maturity stages were compared for each cultivar, revealing many developmental responses common or specific for each cultivar. Top protein biomarkers for each developmental stage were also selected based on larger fold-change differences and lower variance among replicates, including proteins for biosynthesis of lipids and phenols, defense-related proteins and desiccation stress-related proteins. Comparison between the genotypes also revealed the common and specific protein repertoires, totaling 321 pellicle proteins with differential abundance at harvest stage. The proteomics data provides clues on antioxidant, secondary, and hormonal metabolism that could be involved in the loss of quality in the pellicles during processing for commercialization.

Major QTL with pleiotropic effects controlling time of leaf budburst and flowering-related traits in walnut (Juglans regia L.)

Breeding studies in walnut (Juglans regia L.) are usually time consuming due to the long juvenile period and therefore, this study aimed to determine markers associated with time of leaf budburst and flowering-related traits by performing a genome-wide association study (GWAS). We investigated genotypic variation and its association with time of leaf budburst and flowering-related traits in 188 walnut accessions. Phenotypic data was obtained from 13 different traits during 3 consecutive years. We used DArT-seq for genotyping with a total of 33,519 (14,761 SNP and 18,758 DArT) markers for genome-wide associations to identify marker underlying these traits. Significant correlations were determined among the 13 different traits. Linkage disequilibrium decayed very quickly in walnut in comparison with other plants. Sixteen quantitative trait loci (QTL) with major effects (R² between 0.08 and 0.23) were found to be associated with a minimum of two phenotypic traits each. Of these QTL, QTL05 had the maximum number of associated traits (seven). Our study is GWAS for time of leaf budburst and flowering-related traits in Juglans regia L. and has a strong potential to efficiently implement the identified QTL in walnut breeding programs.

STAGdb: a 30K SNP genotyping array and Science Gateway for Acropora corals and their dinoflagellate symbionts

Standardized identification of genotypes is necessary in animals that reproduce asexually and form large clonal populations such as coral. We developed a high-resolution hybridization-based genotype array coupled with an analysis workflow and database for the most speciose genus of coral, Acropora, and their symbionts. We designed the array to co-analyze host and symbionts based on bi-allelic single nucleotide polymorphisms (SNP) markers identified from genomic data of the two Caribbean Acropora species as well as their dominant dinoflagellate symbiont, Symbiodinium ‘fitti’. SNPs were selected to resolve multi-locus genotypes of host (called genets) and symbionts (called strains), distinguish host populations and determine ancestry of coral hybrids between Caribbean acroporids. Pacific acroporids can also be genotyped using a subset of the SNP loci and additional markers enable the detection of symbionts belonging to the genera Breviolum, Cladocopium, and Durusdinium. Analytic tools to produce multi-locus genotypes of hosts based on these SNP markers were combined in a workflow called the Standard Tools for Acroporid Genotyping (STAG). The STAG workflow and database are contained within a customized Galaxy environment (https://coralsnp.science.psu.edu/galaxy/), which allows for consistent identification of host genet and symbiont strains and serves as a template for the development of arrays for additional coral genera. STAG data can be used to track temporal and spatial changes of sampled genets necessary for restoration planning and can be applied to downstream genomic analyses. Using STAG, we uncover bi-directional hybridization between and population structure within Caribbean acroporids and detect a cryptic Acroporid species in the Pacific.

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.

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.

Comparative genomics of six Juglans species reveals disease-associated gene family contractions

Juglans (walnuts), the most speciose genus in the walnut family (Juglandaceae), represents most of the family's commercially valuable fruit and wood-producing trees. It includes several species used as rootstock for their resistance to various abiotic and biotic stressors. We present the full structural and functional genome annotations of six Juglans species and one outgroup within Juglandaceae (Juglans regia, J. cathayensis, J. hindsii, J. microcarpa, J. nigra, J. sigillata and Pterocarya stenoptera) produced using BRAKER2 semi-unsupervised gene prediction pipeline and additional tools. For each annotation, gene predictors were trained using 19 tissue-specific J. regia transcriptomes aligned to the genomes. Additional functional evidence and filters were applied to multi-exonic and mono-exonic putative genes to yield between 27 000 and 44 000 high-confidence gene models per species. Comparison of gene models to the BUSCO embryophyta dataset suggested that, on average, genome annotation completeness was 85.6%. We utilized these high-quality annotations to assess gene family evolution within Juglans, and among Juglans and selected Eurosid species. We found notable contractions in several gene families in J. hindsii, including disease resistance-related wall-associated kinase (WAK), Catharanthus roseus receptor-like kinase (CrRLK1L) and others involved in abiotic stress response. Finally, we confirmed an ancient whole-genome duplication that took place in a common ancestor of Juglandaceae using site substitution comparative analysis.

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.

An Axiom SNP genotyping array for Douglas-fir

BACKGROUND

In forest trees, genetic markers have been used to understand the genetic architecture of natural populations, identify quantitative trait loci, infer gene function, and enhance tree breeding. Recently, new, efficient technologies for genotyping thousands to millions of single nucleotide polymorphisms (SNPs) have finally made large-scale use of genetic markers widely available. These methods will be exceedingly valuable for improving tree breeding and understanding the ecological genetics of Douglas-fir, one of the most economically and ecologically important trees in the world.

RESULTS

We designed SNP assays for 55,766 potential SNPs that were discovered from previous transcriptome sequencing projects. We tested the array on ~ 2300 related and unrelated coastal Douglas-fir trees (Pseudotsuga menziesii var. menziesii) from Oregon and Washington, and 13 trees of interior Douglas-fir (P. menziesii var. glauca). As many as ~ 28 K SNPs were reliably genotyped and polymorphic, depending on the selected SNP call rate. To increase the number of SNPs and improve genome coverage, we developed protocols to 'rescue' SNPs that did not pass the default Affymetrix quality control criteria (e.g., 97% SNP call rate). Lowering the SNP call rate threshold from 97 to 60% increased the number of successful SNPs from 20,669 to 28,094. We used a subset of 395 unrelated trees to calculate SNP population genetic statistics for coastal Douglas-fir. Over a range of call rate thresholds (97 to 60%), the median call rate for SNPs in Hardy-Weinberg equilibrium ranged from 99.2 to 99.7%, and the median minor allele frequency ranged from 0.198 to 0.233. The successful SNPs also worked well on interior Douglas-fir.

CONCLUSIONS

Based on the original transcriptome assemblies and comparisons to version 1.0 of the Douglas-fir reference genome, we conclude that these SNPs can be used to genotype about 10 K to 15 K loci. The Axiom genotyping array will serve as an excellent foundation for studying the population genomics of Douglas-fir and for implementing genomic selection. We are currently using the array to construct a linkage map and test genomic selection in a three-generation breeding program for coastal Douglas-fir.

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 Axiom^TM 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.

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 (H_O = 0.34 and H_E = 0.38). The overall level of genetic differentiation among populations was moderate (F_ST = 0.07). However, the Semnan population showed high divergence from the other Iranian populations (on average F_ST = 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.

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 ² = 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.