Genomic Variation Among and Within Six Juglans Species

Ancient structural variants control sex-specific flowering time morphs in walnuts and hickories

Balanced mating type polymorphisms offer a distinct window into the forces shaping sexual reproduction strategies. Multiple hermaphroditic genera in Juglandaceae, including walnuts (Juglans) and hickories (Carya), show a 1:1 genetic dimorphism for male versus female flowering order (heterodichogamy). We map two distinct Mendelian inheritance mechanisms to ancient (>37 million years old) genus-wide structural DNA polymorphisms. The dominant haplotype for female-first flowering in Juglans contains tandem repeats of the 3' untranslated region of a gene putatively involved in trehalose-6-phosphate metabolism and is associated with increased cis gene expression in developing male flowers, possibly mediated by small RNAs. The Carya locus contains ~20 syntenic genes and shows molecular signatures of sex chromosome-like evolution. Inheritance mechanisms for heterodichogamy are deeply conserved, yet may occasionally turn over, as in sex determination.

Comprehensive insights on association mapping in perennial fruit crops breeding - Its implications, current status and future perspectives

In order to provide food and nutritional security for the world's rapidly expanding population, fruit crop researchers have identified two critical priorities: increasing production and preserving fruit quality during the pre- and post-harvest periods. The genetic basis of these complex, commercially important fruit traits which are uniquely regulated by polygenes or multi-allelic genes that interact with one another and the environment can be analyzed with the aid of trait mapping tools. The most interesting trait mapping approach that offers the genetic level investigation for marker-trait associations (MTAs) for these complex fruit traits, without the development of mapping population, is association mapping. This approach was used during the genetic improvement program, emphasizing the obstacles (breeding strategies adopted, generation interval, and their genomic status) pertaining to perennial fruit crops. This method of studying population diversity and linkage disequilibrium in perennial fruit crops has been made possible by recent developments in genotyping, phenotyping, and statistical analysis. Thus, the purpose of this review is to provide an overview of different trait mapping techniques, with a focus on association mapping (method, essential components, viability, constraints, and future perspective) and its advantages, disadvantages, and possibilities for breeding perennial fruit crops.

Comprehensive Evaluation of Cryptic Juglans Genotypes: Insight from Molecular Markers and Phylogenetic Analysis

Background/Objectives: The current research work aimed to evaluate the cryptic walnut genotypes of the Hazara region in Pakistan by using DNA barcoding and phylogenetic analysis. Methods: Based on morphological traits such as nut size, nut shape, and the number of leaflets, five genotypes were chosen and samples were collected for the current study. For molecular analysis, gDNA was isolated from the fresh leaves, and the five most effective angiosperm-specific markers, ITS2, rbcLa, rbcLc, rpoC1, and UBE3, were utilized. Based on amplification, sequencing, and identification success rates, ITS2 and UBE3 were recorded as the most efficient markers followed by rbcLa, rbcLc, and rpoC1. Results: During phylogenetic analysis, the query genotype-1 based on ITS2 and genotype-2 based on UBE3 clustered with (KF454101.1-Juglans regia) and (KC870919.1-J. regia) with bootstraps of 56 and 100, respectively. Genotype-3 based on rbcla clustered in a major clade with J. regia L., cultivars (MN397935.1 J. regia 'Vina') and (MN397934.1-J. regia 'Serr'), (MN397933.1 J. regia 'Pedro'), (MN397932.1 J. regia 'Lara'), (MN397931.1 J. regia 'Howard'), and (MN397930.1 J. regia 'Hartley') with bootstrap of 100. Meanwhile, genotype-4 and genotype-5 based on rbclc and rpoC1 clustered with (MN397935.1 J. regia 'Vina') and (MN397934.1 J. regia 'Serr'), across the database sequences. To clarify the taxonomic status of cryptic walnut genotypes, it is necessary to combine diverse DNA barcodes. The results of ITS2 and UBE3, followed by rbcL barcoding markers, are promising taxonomic tools for cryptic walnut genotypes in Pakistan. Conclusions: It has been determined that the genotypes of walnuts in the study area are both J. regia L. and its cultivars and that the accuracy of discrimination regarding the genus Juglans L. is greater than 90%. The reported DNA barcodes are recommended for the correct identification and genetic evaluation of Juglans taxa and its population.

Genome-wide identification analysis of the 4-Coumarate: CoA ligase (4CL) gene family expression profiles in Juglans regia and its wild relatives J. Mandshurica resistance and salt stress

Persian walnut (Juglans regia) and Manchurian walnut (Juglans mandshurica) belong to Juglandaceae, which are vulnerable, temperate deciduous perennial trees with high economical, ecological, and industrial values. 4-Coumarate: CoA ligase (4CL) plays an essential function in plant development, growth, and stress. Walnut production is challenged by diverse stresses, such as salinity, drought, and diseases. However, the characteristics and expression levels of 4CL gene family in Juglans species resistance and under salt stress are unknown. Here, we identified 36 Jr4CL genes and 31 Jm4CL genes, respectively. Based on phylogenetic relationship analysis, all 4CL genes were divided into three branches. WGD was the major duplication mode for 4CLs in two Juglans species. The phylogenic and collinearity analyses showed that the 4CLs were relatively conserved during evolution, but the gene structures varied widely. 4CLs promoter region contained multiply cis-acting elements related to phytohormones and stress responses. We found that Jr4CLs may be participated in the regulation of resistance to anthracnose. The expression level and some physiological of 4CLs were changed significantly after salt treatment. According to qRT-PCR results, positive regulation was found to be the main mode of regulation of 4CL genes after salt stress. Overall, J. mandshurica outperformed J. regia. Therefore, J. mandshurica can be used as a walnut rootstock to improve salt tolerance. Our results provide new understanding the potential functions of 4CL genes in stress tolerance, offer the theoretical genetic basis of walnut varieties adapted to salt stress, and provide an important reference for breeding cultivated walnuts for stress tolerance.

A haplotype-resolved genome provides insight into allele-specific expression in wild walnut (Juglans regia L.)

Wild germplasm resources are crucial for gene mining and molecular breeding because of their special trait performance. Haplotype-resolved genome is an ideal solution for fully understanding the biology of subgenomes in highly heterozygous species. Here, we surveyed the genome of a wild walnut tree from Gongliu County, Xinjiang, China, and generated a haplotype-resolved reference genome of 562.99 Mb (contig N50 = 34.10 Mb) for one haplotype (hap1) and 561.07 Mb (contig N50 = 33.91 Mb) for another haplotype (hap2) using PacBio high-fidelity (HiFi) reads and Hi-C technology. Approximately 527.20 Mb (93.64%) of hap1 and 526.40 Mb (93.82%) of hap2 were assigned to 16 pseudochromosomes. A total of 41039 and 39744 protein-coding gene models were predicted for hap1 and hap2, respectively. Moreover, 123 structural variations (SVs) were identified between the two haplotype genomes. Allele-specific expression genes (ASEGs) that respond to cold stress were ultimately identified. These datasets can be used to study subgenome evolution, for functional elite gene mining and to discover the transcriptional basis of specific traits related to environmental adaptation in wild walnut.

Distinct ancient structural polymorphisms control heterodichogamy in walnuts and hickories

The maintenance of stable mating type polymorphisms is a classic example of balancing selection, underlying the nearly ubiquitous 50/50 sex ratio in species with separate sexes. One lesser known but intriguing example of a balanced mating polymorphism in angiosperms is heterodichogamy - polymorphism for opposing directions of dichogamy (temporal separation of male and female function in hermaphrodites) within a flowering season. This mating system is common throughout Juglandaceae, the family that includes globally important and iconic nut and timber crops - walnuts (Juglans), as well as pecan and other hickories (Carya). In both genera, heterodichogamy is controlled by a single dominant allele. We fine-map the locus in each genus, and find two ancient (>50 Mya) structural variants involving different genes that both segregate as genus-wide trans-species polymorphisms. The Juglans locus maps to a ca. 20 kb structural variant adjacent to a probable trehalose phosphate phosphatase (TPPD-1), homologs of which regulate floral development in model systems. TPPD-1 is differentially expressed between morphs in developing male flowers, with increased allele-specific expression of the dominant haplotype copy. Across species, the dominant haplotype contains a tandem array of duplicated sequence motifs, part of which is an inverted copy of the TPPD-1 3' UTR. These repeats generate various distinct small RNAs matching sequences within the 3' UTR and further downstream. In contrast to the single-gene Juglans locus, the Carya heterodichogamy locus maps to a ca. 200-450 kb cluster of tightly linked polymorphisms across 20 genes, some of which have known roles in flowering and are differentially expressed between morphs in developing flowers. The dominant haplotype in pecan, which is nearly always heterozygous and appears to rarely recombine, shows markedly reduced genetic diversity and is over twice as long as its recessive counterpart due to accumulation of various types of transposable elements. We did not detect either genetic system in other heterodichogamous genera within Juglandaceae, suggesting that additional genetic systems for heterodichogamy may yet remain undiscovered.

Genome-Wide Identification of the CBF Gene Family and ICE Transcription Factors in Walnuts and Expression Profiles under Cold Conditions

Cold stress impacts woody tree growth and perennial production, especially when the temperature rapidly changes in late spring. To address this issue, we conducted the genome-wide identification of two important transcription factors (TFs), CBF (C-repeat binding factors) and ICE (inducers of CBF expression), in three walnut (Juglans) genomes. Although the CBF and ICE gene families have been identified in many crops, very little systematic analysis of these genes has been carried out in J. regia and J. sigillata. In this study, we identified a total of 16 CBF and 12 ICE genes in three Juglans genomes using bioinformatics analysis. Both CBF and ICE had conserved domains, motifs, and gene structures, which suggests that these two TFs were evolutionarily conserved. Most ICE genes are located at both ends of the chromosomes. The promoter cis-regulatory elements of CBF and ICE genes are largely involved in light and phytohormone responses. Based on 36 RNA sequencing of leaves from four walnut cultivars ('Zijing', 'Lvling', 'Hongren', and 'Liao1') under three temperature conditions (8 °C, 22 °C, and 5 °C) conditions in late spring, we found that the ICE genes were expressed more highly than CBFs. Both CBF and ICE proteins interacted with cold-related proteins, and many putative miRNAs had interactions with these two TFs. These results determined that CBF1 and ICE1 play important roles in the tolerance of walnut leaves to rapid temperature changes. Our results provide a useful resource on the function of the CBF and ICE genes related to cold tolerance in walnuts.

Reference genome of California walnut, Juglans californica, and resemblance with other genomes in the order Fagales

Juglans californica, California walnut, is a vulnerable small tree that is locally abundant but restricted to woodland and chaparral habitats of Southern California threatened by urbanization and land use change. This species is the dominant species in a unique woodland ecosystem in California. It is one of 2 endemic California walnut species (family Juglandaceae). The other species, Northern California black walnut (J. hindsii), has been suggested controversially to be a variety of J. californica. Here, we report a new, chromosome-level assembly of J. californica as part of the California Conservation Genomics Project (CCGP). Consistent with the CCGP common methodology across ~150 genomes, we used Pacific Biosciences HiFi long reads and Omni-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises 137 scaffolds spanning 551,065,703 bp, has a contig N50 of 30 Mb, a scaffold N50 of 37 Mb, and BUSCO complete score of 98.9%. Additionally, the mitochondrial genome has 701,569 bp. In addition, we compare this genome with other existing high-quality Juglans and Quercus genomes, which are in the same order (Fagales) and show relatively high synteny within the Juglans genomes. Future work will utilize the J. californica genome to determine its relationship with the Northern California walnut and assess the extent to which these 2 endemic trees might be at risk from fragmentation and/or climate warming.

Genome assembly of two diploid and one auto-tetraploid Cyclocarya paliurus genomes

Cyclocarya paliurus, an endemic species in the genus Juglandaceae with the character of heterodichogamy, is one of triterpene-rich medicinal plants in China. To uncover the genetic mechanisms behind the special characteristics, we sequenced the genomes of two diploid (protandry, PA-dip and protogyny, PG-dip) and one auto-tetraploid (PA-tetra) C. paliurus genomes. Based on 134.9 (~225x), 75.5 (~125x) and 271.8 Gb (~226x) subreads of PacBio platform sequencing data, we assembled 586.62 Mb (contig N50 = 1.9 Mb), 583.45 Mb (contig N50 = 1.4 Mb), and 2.38 Gb (contig N50 = 430.9 kb) for PA-dip, PG-dip and PA-tetra genome, respectively. Furthermore, 543.53, 553.87, and 2168.65 Mb in PA-dip, PG-dip, and PA-tetra, were respectively anchored to 16, 16, and 64 pseudo-chromosomes using over 65.4 Gb (~109x), 68 Gb (~113x), and 264 (~220x) Hi-C sequencing data. Annotation of PA-dip, PG-dip, and PA-tetra genome assembly identified 34,699, 35,221, and 34,633 protein-coding genes (90,752 gene models) or allele-defined genes, respectively. In addition, 45 accessions from nine locations were re-sequenced, and more than 10 × coverage reads were generated.

Haplotyping interspecific hybrids by dual alignment to both parental genomes

Sequencing-based genotyping of heterozygous diploids requires sufficient depth to accurately call heterozygous genotypes. In interspecific hybrids, alignment of reads to both parental genomes simultaneously can generate haploid data, potentially eliminating the problem of heterozygosity. Two populations of interspecific hybrid rootstocks of walnut (Juglans) and pistachio (Pistacia) were genotyped using alignment to the maternal genome, paternal genome, and dual alignment to both genomes simultaneously. Downsampling was used to examine concordance of imputed genotype calls as a function of sequencing depth. Dual alignment resulted in datasets essentially free of heterozygous genotypes, simplifying the identification and removal of cross-contaminated samples. Concordance between full and downsampled genotype calls was always highest after dual alignment. Nearly all single nucleotide polymorphisms (SNPs) in dual alignment datasets were shared with the corresponding single-parent datasets, but 60%-90% of single-parent SNPs were private to that dataset. Private SNPs in single-parent datasets had higher rates of heterozygosity, lower levels of concordance, and were enriched in fixed differences between parental genomes ("homeo-SNPs") compared to shared SNPs in the same dataset. In multi-parental walnut hybrids, the paternal-aligned dataset was ineffective at resolving population structure in the maternal parent. Overall, the dual alignment strategy effectively produced phased, haploid data, increasing data quality and reducing cost.

Pan-genome and transcriptome analyses provide insights into genomic variation and differential gene expression profiles related to disease resistance and fatty acid biosynthesis in eastern black walnut (Juglans nigra)

Walnut (Juglans) species are used as nut crops worldwide. Eastern black walnut (EBW, Juglans nigra), a diploid, horticultural important woody species is native to much of eastern North America. Although it is highly valued for its wood and nut, there are few resources for understanding EBW genetics. Here, we present a high-quality genome assembly of J. nigra based on Illumina, Pacbio, and Hi-C technologies. The genome size was 540.8 Mb, with a scaffold N50 size of 35.1 Mb, and 99.0% of the assembly was anchored to 16 chromosomes. Using this genome as a reference, the resequencing of 74 accessions revealed the effective population size of J. nigra declined during the glacial maximum. A single whole-genome duplication event was identified in the J. nigra genome. Large syntenic blocks among J. nigra, Juglans regia, and Juglans microcarpa predominated, but inversions of more than 600 kb were identified. By comparing the EBW genome with those of J. regia and J. microcarpa, we detected InDel sizes of 34.9 Mb in J. regia and 18.3 Mb in J. microcarpa, respectively. Transcriptomic analysis of differentially expressed genes identified five presumed NBS-LRR (NUCLEOTIDE BINDING SITE-LEUCINE-RICH REPEAT) genes were upregulated during the development of walnut husks and shells compared to developing embryos. We also identified candidate genes with essential roles in seed oil synthesis, including FAD (FATTY ACID DESATURASE) and OLE (OLEOSIN). Our work advances the understanding of fatty acid bioaccumulation and disease resistance in nut crops, and also provides an essential resource for conducting genomics-enabled breeding in walnut.

K-Mer-Based Genome Size Estimation in Theory and Practice

Recent advances in sequencing technologies have made genome sequencing of non-model organisms with very large and complex genomes possible. The data can be used to estimate diverse genome characteristics, including genome size, repeat content, and levels of heterozygosity. K-mer analysis is a powerful biocomputational approach with a wide range of applications, including estimation of genome sizes. However, interpretation of the results is not always straightforward. Here, I review k-mer-based genome size estimation, focusing specifically on k-mer theory and peak calling in k-mer frequency histograms. I highlight common pitfalls in data analysis and result interpretation, and provide a comprehensive overview on current methods and programs developed to conduct these analyses.

The Manchurian Walnut Genome: Insights into Juglone and Lipid Biosynthesis

BACKGROUND

Manchurian walnut (Juglans mandshurica Maxim.) is a tree with multiple industrial uses and medicinal properties in the Juglandaceae family (walnuts and hickories). J. mandshurica produces juglone, which is a toxic allelopathic agent and has potential utilization value. Furthermore, the seed of J. mandshurica is rich in various unsaturated fatty acids and has high nutritive value.

FINDINGS

Here, we present a high-quality chromosome-scale reference genome assembly and annotation for J. mandshurica (n = 16) with a contig N50 of 21.4 Mb by combining PacBio high-fidelity reads with high-throughput chromosome conformation capture data. The assembled genome has an estimated sequence size of 548.7 Mb and consists of 657 contigs, 623 scaffolds, and 40,453 protein-coding genes. In total, 60.99% of the assembled genome consists of repetitive sequences. Sixteen super-scaffolds corresponding to the 16 chromosomes were assembled, with a scaffold N50 length of 33.7 Mb and a BUSCO complete gene percentage of 98.3%. J. mandshurica displays a close sequence relationship with Juglans cathayensis, with a divergence time of 13.8 million years ago. Combining the high-quality genome, transcriptome, and metabolomics data, we constructed a gene-to-metabolite network and identified 566 core and conserved differentially expressed genes, which may be involved in juglone biosynthesis. Five CYP450 genes were found that may contribute to juglone accumulation. NAC, bZip, NF-YA, and NF-YC are positively correlated with the juglone content. Some candidate regulators (e.g., FUS3, ABI3, LEC2, and WRI1 transcription factors) involved in the regulation of lipid biosynthesis were also identified.

CONCLUSIONS

Our genomic data provide new insights into the evolution of the walnut genome and create a new platform for accelerating molecular breeding and improving the comprehensive utilization of these economically important tree species.

Advances in genomics and genome editing for breeding next generation of fruit and nut crops

Fruit tree crops are an essential part of the food production systems and are key to achieve food and nutrition security. Genetic improvement of fruit trees by conventional breeding has been slow due to the long juvenile phase. Advancements in genomics and molecular biology have paved the way for devising novel genetic improvement tools like genome editing, which can accelerate the breeding of these perennial crops to a great extent. In this article, advancements in genomics of fruit trees covering genome sequencing, transcriptome sequencing, genome editing technologies (GET), CRISPR-Cas system based genome editing, potential applications of CRISPR-Cas9 in fruit tree crops improvement, the factors influencing the CRISPR-Cas editing efficiency and the challenges for CRISPR-Cas9 applications in fruit tree crops improvement are reviewed. Besides, base editing, a recently emerging more precise editing system, and the future perspectives of genome editing in the improvement of fruit and nut crops are covered.

Transcriptomic and Metabolic Analysis of Fruit Development and Identification of Genes Involved in Raffinose and Hydrolysable Tannin Biosynthesis in Walnuts

Walnut (Juglans regia L.) is an important fruit tree with high nutrition in its nuts. Here, the development of walnut fruits was monitored, and nine biological samples at five developmental stages were collected and analyzed by transcriptomic and metabolic assays. Many phenolic metabolites accumulated in the peel of mature fruits, while lipids, carbohydrates, and amino acids and their derivatives mainly accumulated in the kernel. Fatty acid biosynthesis occurred at 13 weeks after pollination, and photosynthesis might occur in the exocarp of walnuts. By coexpression analysis of the transcriptome and metabolome, genes responsible for some metabolic pathways were predicted. Three genes encoding shikimate dehydrogenases (SDHs) that convert 3-dehydroshikimic acid to gallic acid (GA) and four genes encoding UDP-glycosyltransferase (UGT) that convert GA to β-glucogallin in the biosynthesis of hydrolysable tannins (HTs) were selected for functional confirmation. These three SDH genes were then expressed in Escherichia coli, and their recombinant proteins showed GA formation activity. Moreover, heterologous expression of the three SDH and four UGT genes in poplar hairy roots also showed a significant increase in GA and β-glucogallin accumulation, respectively. Taken together, we have provided an overview of walnut fruit development and uncovered genes involved in HT biosynthesis.

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.

Plant NLR diversity: the known unknowns of pan-NLRomes

Plants and pathogens constantly adapt to each other. As a consequence, many members of the plant immune system, and especially the intracellular nucleotide-binding site leucine-rich repeat receptors, also known as NOD-like receptors (NLRs), are highly diversified, both among family members in the same genome, and between individuals in the same species. While this diversity has long been appreciated, its true extent has remained unknown. With pan-genome and pan-NLRome studies becoming more and more comprehensive, our knowledge of NLR sequence diversity is growing rapidly, and pan-NLRomes provide powerful platforms for assigning function to NLRs. These efforts are an important step toward the goal of comprehensively predicting from sequence alone whether an NLR provides disease resistance, and if so, to which pathogens.

Co-located quantitative trait loci mediate resistance to Agrobacterium tumefaciens, Phytophthora cinnamomi, and P. pini in Juglans microcarpa × J. regia hybrids

Soil-borne plant pathogens represent a serious threat that undermines commercial walnut (Juglans regia) production worldwide. Crown gall, caused by Agrobacterium tumefaciens, and Phytophthora root and crown rots, caused by various Phytophthora spp., are among the most devastating walnut soil-borne diseases. A recognized strategy to combat soil-borne diseases is adoption of resistant rootstocks. Here, resistance to A. tumefaciens, P. cinnamomi, and P. pini is mapped in the genome of Juglans microcarpa, a North American wild relative of cultivated walnut. Half-sib J. microcarpa mother trees DJUG 31.01 and DJUG 31.09 were crossed with J. regia cv. Serr, producing 353 and 400 hybrids, respectively. Clonally propagated hybrids were genotyped by sequencing to construct genetic maps for the two populations and challenged with the three pathogens. Resistance to each of the three pathogens was mapped as a major QTL on the long arm of J. microcarpa chromosome 4D and was associated with the same haplotype, designated as haplotype b, raising the possibility that the two mother trees were heterozygous for a single Mendelian gene conferring resistance to all three pathogens. The deployment of this haplotype in rootstock breeding will facilitate breeding of a walnut rootstock resistant to both crown gall and Phytophthora root and crown rots.

Improved de novo chromosome-level genome assembly of the vulnerable walnut tree Juglans mandshurica reveals gene family evolution and possible genome basis of resistance to lesion nematode

Manchurian walnut (Juglans mandshurica Maxim.) is a synonym of J. cathayensis, a diploid, vulnerable, temperate deciduous tree valued for its wood and nut. It is also valued as a rootstock for Juglans regia because of its reported tolerance of lesion nematode. Reference genomes are available for several Juglans species, our goal was to produce a de novo, chromosome-level assembly of the J. mandshurica genome. Here, we reported an improved assembly of J. mandshurica with a contig N50 size of 6.49 Mb and a scaffold N50 size of 36.1 Mb. The total genome size was 548 Mb encoding 29,032 protein coding genes which were annotated. The collinearity analysis showed that J. mandshurica and J. regia originated from a common ancestor, with both species undergoing two WGD events. A genomic comparison showed that J. mandshurica was missing 1657 genes found in J. regia, and J. mandshurica includes 2827 genes not found in of the J. regia genome. The J. mandshurica contained 1440 unique paralogues that were highly enriched for flavonoid biosynthesis, phenylpropanoid biosynthesis, and plant-pathogen interaction. Four gene families related to disease resistance notable contraction (rapidly evolving; LEA, WAK, PPR, and PR) in J. mandshurica compared to eight species. JmaPR10 and JmaPR8 contained three orthologous gene pairs with J. regia that were highly expressed in root bark. JmaPR10 is a strong candidate gene for lesion nematodes resistance in J. mandshurica. The J. mandshurica genome should be a useful resource for study of the evolution, breeding, and genetic variation in walnuts (Juglans).

Comparative analysis of drought-responsive and -adaptive genes in Chinese wingnut (Pterocarya stenoptera C. DC)

Background

Drought is the main stress factor for the cultivation of Pterocarya stenoptera in urban areas, and this factor will cause its dehydration and affect its growth. Identifying drought-related genes will be useful for understanding the drought adaptation mechanism of P. stenoptera.

Results

We used physiological indicator detection, comparative transcriptome sequencing, and reanalysis on the results of previous landscape genomics studies to investigate the drought adaptation mechanism in P. stenoptera. The changes in malondialdehyde content showed that P. stenoptera was remarkably affected by drought stress, and the increase in soluble sugar content suggested its important role in response to drought stress. Results of comparative transcriptome sequencing showed that P. stenoptera initiated a series of programs, such as increasing the gene expression of unsaturated fatty acids, tyrosine, and plant pathogen resistance, to deal with the transient drought stress. According to the annotated results in a previous study, P. stenoptera adapts to the long-term differential drought stress by regulating the thickness of cell walls and expressing upper or lower limits of the downstream genes in the hormone signaling pathway. Through the comparative analysis of drought-responsive and -adaptive genes in P. stenoptera, this study supports the hypothesis that the environment-responsive genes (ERGs) introduced by the transient environmental stresses will be substantially more than the environment-adaptive genes (EAGs) in response to long-term differential environmental stresses, and the EAGs are not necessarily ERGs.

Conclusions

Our study identified drought-responsive and -adaptive genes in P. stenoptera and revealed that P. stenoptera increased the gene expression of unsaturated fatty acids, tyrosine, and plant pathogen resistance in response to transient drought stress. This study reveals the different adaptation mechanism of P. stenoptera under the transient and long-term differential drought stresses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07470-z.

Convergent evolution of plant specialized 1,4-naphthoquinones: metabolism, trafficking, and resistance to their allelopathic effects

Plant 1,4-naphthoquinones encompass a class of specialized metabolites known to mediate numerous plant-biotic interactions. This class of compounds also presents a remarkable case of convergent evolution. The 1,4-naphthoquinones are synthesized by species belonging to nearly 20 disparate orders spread throughout vascular plants, and their production occurs via one of four known biochemically distinct pathways. Recent developments from large-scale biology and genetic studies corroborate the existence of multiple pathways to synthesize plant 1,4-naphthoquinones and indicate that extraordinary events of metabolic innovation and links to respiratory and photosynthetic quinone metabolism probably contributed to their independent evolution. Moreover, because many 1,4-naphthoquinones are excreted into the rhizosphere and they are highly reactive in biological systems, plants that synthesize these compounds also needed to independently evolve strategies to deploy them and to resist their effects. In this review, we highlight new progress made in understanding specialized 1,4-naphthoquinone biosynthesis and trafficking with a focus on how these discoveries have shed light on the convergent evolution and diversification of this class of compounds in plants. We also discuss how emerging themes in metabolism-based herbicide resistance may provide clues to mechanisms plants employ to tolerate allelopathic 1,4-naphthoquinones.

Computational Identification and Comparative Analysis of Conserved miRNAs and Their Putative Target Genes in the Juglans regia and J. microcarpa Genomes

MicroRNAs (miRNAs) are important factors for the post-transcriptional regulation of protein-coding genes in plants and animals. They are discovered either by sequencing small RNAs or computationally. We employed a sequence-homology-based computational approach to identify conserved miRNAs and their target genes in Persian (English) walnut, Juglans regia, and its North American wild relative, J. microcarpa. A total of 119 miRNA precursors (pre-miRNAs) were detected in the J. regia genome and 121 in the J. microcarpa genome and miRNA target genes were predicted and their functional annotations were performed in both genomes. In the J. regia genome, 325 different genes were targets; 87.08% were regulated by transcript cleavage and 12.92% by translation repression. In the J. microcarpa genome, 316 different genes were targets; 88.92% were regulated by transcript cleavage and 11.08% were regulated by translation repression. Totals of 1.3% and 2.0% of all resistance gene analogues (RGA) and 2.7% and 2.6% of all transcription factors (TFs) were regulated by miRNAs in the J. regia and J. microcarpa genomes, respectively. Juglans genomes evolved by a whole genome duplication (WGD) and consist of eight pairs of fractionated homoeologous chromosomes. Within each pair, the chromosome that has more genes with greater average transcription also harbors more pre-miRNAs and more target genes than its homoeologue. While only minor differences were detected in pre-miRNAs between the J. regia and J. microcarpa genomes, about one-third of the pre-miRNA loci were not conserved between homoeologous chromosome within each genome. Pre-miRNA and their corresponding target genes showed a tendency to be collocated within a subgenome.

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.

Identification of the amino acid position controlling the different enzymatic activities in walnut tyrosinase isoenzymes (jrPPO1 and jrPPO2)

Polyphenol oxidases (PPOs) are ubiquitously distributed among plants, bacteria, fungi and animals. They catalyze the hydroxylation of monophenols (monophenolase activity) and the oxidation of o-diphenols (diphenolase activity) to o-quinones. PPOs are commonly present as an isoenzyme family. In walnut (Juglans regia), two different genes (jrPPO1 and jrPPO2) encoding PPOs have been identified. In this study, jrPPO2 was, for the first time, heterologously expressed in E. coli and characterized as a tyrosinase (TYR) by substrate scope assays and kinetic investigations, as it accepted tyramine and L-tyrosine as substrates. Moreover, the substrate acceptance and kinetic parameters (k_cat and K_m values) towards 16 substrates naturally present in walnut were assessed for jrPPO2 (TYR) and its isoenzyme jrPPO1 (TYR). The two isoenzymes prefer different substrates, as jrPPO1 shows a higher activity towards monophenols, whereas jrPPO2 is more active towards o-diphenols. Molecular docking studies performed herein revealed that the amino acid residue in the position of the 1st activity controller (His_B1 + 1; in jrPPO1 Asn240 and jrPPO2 Gly240) is responsible for the different enzymatic activities. Additionally, interchanging the 1st activity controller residue of the two enzymes in two mutants (jrPPO1-Asn240Gly and jrPPO2-Gly240Asn) proved that the amino acid residue located in this position allows plants to selectively target or dismiss substrates naturally present in walnut.

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.

Combining independent de novo assemblies to optimize leaf transcriptome of Persian walnut

Transcriptome resources can facilitate to increase yield and quality of walnuts. Finding the best transcriptome assembly has not been the subject of walnuts research as yet. This research generated 240,179,782 reads from 11 walnut leaves according to cDNA libraries. The reads provided a complete de novo transcriptome assembly. Fifteen different transcriptome assemblies were constructed from five different well-known assemblers used in scientific literature with different k-mer lengths (Bridger, BinPacker, SOAPdenovo-Trans, Trinity and SPAdes) as well as two merging approaches (EvidentialGene and Transfuse). Based on the four quality metrics of assembly, the results indicated an efficiency in the process of merging the assemblies after being generated by de novo assemblers. Finally, EvidentialGene was recognized as the best assembler for the de novo assembly of the leaf transcriptome in walnut. Among a total number of 183,191 transcripts which were generated by EvidentialGene, there were 109,413 transcripts capable of protein potential (59.72%) and 104,926 were recognized as ORFs (57.27%). In addition, 79,185 transcripts were predicted to exist with at least one hit to the Pfam database. A number of 3,931 transcription factors were identified by BLAST searching against PlnTFDB. Furthermore, 6,591 of the predicted peptide sequences contained signaling peptides, while 92,704 contained transmembrane domains. Comparison of the assembled transcripts with transcripts of the walnut and published genome assembly for the 'Chandler' cultivar using the BLAST algorithm led to identify a total number of 27,304 and 19,178 homologue transcripts, respectively. De novo transcriptomes in walnut leaves can be developed for the future studies in functional genomics and genetic studies of walnuts.

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.

Portal of Juglandaceae: A comprehensive platform for Juglandaceae study

Juglandaceae species are plants of great economic value and have been cultivated, domesticated, and utilized by human society for a long time. Their edible, nutrient-rich nuts and tough, durable wood have attracted the attention of botanists and breeders. With the advent of the genomics era, genome sequencing of the Juglandaceae family has been greatly accelerated, and a large amount of data has been generated. In this paper, we introduce the Portal of Juglandaceae (PJU), a tool to bring all these data together. The PJU contains genomes, gene-coding sequences, protein sequences, various types of annotation information, expression data, and miRNA data, which are configured with BLAST, JBrowse, and our self-developed synteny analysis tool. The PJU has a user-friendly and straightforward interface that performs a variety of query tasks with a few simple operations. In the future, we hope that the PJU will serve as a hub for the study of the Juglandaceae family.

Morpho-Physiological and Genomic Evaluation of Juglans Species Reveals Regional Maladaptation to Cold Stress

Climate change may have unpredictable effects on the cold hardiness of woody species planted outside of their range of origin. Extreme undulations in temperatures may exacerbate susceptibility to cold stress, thereby interfering with productivity and ecosystem functioning. Juglans L. and their naturally occurring interspecific F1 hybrids, are distributed natively across many temperate regions, and J. regia has been extensively introduced. Cold hardiness, an environmental and genetic factor yet to be evaluated in many native and introduced Juglans species, may be a limiting factor under future climate change and following species introductions. We evaluated cold hardiness of native North American and Eastern Asian Juglans along with J. regia genotypes using field data from the Midwestern United States (Indiana), controlled freezing tests, and genome sequencing with close assessment of Juglans cold hardy genes. Many Juglans species previously screened for cold-hardiness were genotypes derived from the Midwest, California, and Europe. In 2014, despite general climate adaptation, Midwestern winter temperatures of -30°C killed J. regia originating from California; however, naturalized Midwestern J. regia survived and displayed low damage. Hybridization of J. regia with black walnut (J. nigra) and butternut (J. cinerea) produced F1s displaying greater cold tolerance than pure J. regia. Cold hardiness and growth are variable in Midwestern J. regia compared to native Juglans, East Asian Juglans, and F1 hybrids. Phylogeny analyses revealed that J. cinerea sorted with East Asian species using the nuclear genome but with North American species using the organellar genome. Investigation of selected cold hardy genes revealed that J. regia was distinct from other species and exhibited less genetic diversity than native Juglans species Average whole genome heterozygosity and Tajima's D for cold hardy genes was low within J. regia samples and significantly higher for hybrid as well as J. nigra. We confirmed that molecular and morpho-physiological data were highly correlated and thus can be used effectively to characterize cold hardiness in Juglans species. We conclude that the genetic diversity within local J. regia populations is low and additional germplasm is needed for development of more regionally adapted J. regia varieties.

Conversion of walnut tyrosinase into a catechol oxidase by site directed mutagenesis

Polyphenol oxidases (PPOs) comprise tyrosinases (TYRs) and catechol oxidases (COs), which catalyse the initial reactions in the biosynthesis of melanin. TYRs hydroxylate monophenolic (monophenolase activity) and oxidize diphenolic (diphenolase activity) substrates, whereas COs react only with diphenols. In order to elucidate the biochemical basis for the different reactions in PPOs, cDNA from walnut leaves was synthesized, the target gene encoding the latent walnut tyrosinase (jrPPO1) was cloned, and the enzyme was heterologously expressed in Escherichia coli. Mutations targeting the two activity controller residues (Asn240 and Leu244) as well as the gatekeeper residue (Phe260) were designed to impair monophenolase activity of jrPPO1. For the first time, monophenolase activity of jrPPO1 towards L-tyrosine was blocked in two double mutants (Asn240Lys/Leu244Arg and Asn240Thr/Leu244Arg) while its diphenolase activity was partially preserved, thereby converting jrPPO1 into a CO. Kinetic data show that recombinant jrPPO1 resembles the natural enzyme, and spectrophotometric investigations proved that the copper content remains unaffected by the mutations. The results presented herein provide experimental evidence that a precisely tuned interplay between the amino acids located around the active center controls the substrate specificity and therewith the mono- versus diphenolase activity in the type-III copper enzyme jrPPO1.

Genomics: cracking the mysteries of walnuts

The Juglans plants are economically important as they provide nuts, wood and garden trees. They also play an important ecological role by supplying food for wild insects and animals. The decoding of genome sequences has fundamental values for understanding the evolution of Juglans plants and molecules, and is also a prerequisite for molecular breeding. During the last three years, the rapid development of sequencing technology has made walnut research into the genome era. Here, we reviewed the progress of genome sequencing of six Juglans species, the resequencing of four Juglans populations as well as the genome sequencing of the closely related species Pterocarya stenoptera. The analysis of the Juglans regia genome uncovers a whole genome duplication (WGD) event. Based on the molecular dating of the divergence time of six Juglans species, we proposed this WGD event was associated with Cretaceous-Palaeogene (K-Pg) boundary occurred ∼65 million years ago. Genomic sequences also provide clear details for understanding the evolution and development of GGT and PPO genes involved in fruit development. The decoding of these genomes has made it easier for us to understand and enhance the use of walnuts. We expect that the functional genomics research of walnut will also develop rapidly in the near future.

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.

Novel assembly strategy cracks open the mysteries of walnut genome evolution

High quality chromosome-scale assemblies from an interspecific hybrid between walnut and a wild relative reveal the persistence of asymmetric fractionation between the sub-genomes and suggest a late-Miocene origin for the genus Juglans.

Sequencing a Juglans regia × J. microcarpa hybrid yields high-quality genome assemblies of parental species

Members of the genus Juglans are monecious wind-pollinated trees in the family Juglandaceae with highly heterozygous genomes, which greatly complicates genome sequence assembly. The genomes of interspecific hybrids are usually comprised of haploid genomes of parental species. We exploited this attribute of interspecific hybrids to avoid heterozygosity and sequenced an interspecific hybrid Juglans microcarpa × J. regia using a novel combination of single-molecule sequencing and optical genome mapping technologies. The resulting assemblies of both genomes were remarkably complete including chromosome termini and centromere regions. Chromosome termini consisted of arrays of telomeric repeats about 8 kb long and heterochromatic subtelomeric regions about 10 kb long. The centromeres consisted of arrays of a centromere-specific Gypsy retrotransposon and most contained genes, many of them transcribed. Juglans genomes evolved by a whole-genome-duplication dating back to the Cretaceous-Paleogene boundary and consist of two subgenomes, which were fractionated by numerous short gene deletions evenly distributed along the length of the chromosomes. Fractionation was shown to be asymmetric with one subgenome exhibiting greater gene loss than the other. The asymmetry of the process is ongoing and mirrors an asymmetry in gene expression between the subgenomes. Given the importance of J. microcarpa × J. regia hybrids as potential walnut rootstocks, we catalogued disease resistance genes in the parental genomes and studied their chromosomal distribution. We also estimated the molecular clock rates for woody perennials and deployed them in estimating divergence times of Juglans genomes and those of other woody perennials.

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.