The walnut (Juglans regia) genome sequence reveals diversity in genes coding for the biosynthesis of non-structural polyphenols

Genomic prediction in Persian walnut: Optimization levers according to genetic architecture of complex traits

Persian walnut (Juglans regia L.) is a widespread cultivated nut tree species in temperate regions. Advances in genomic tools, such as the high-density Axiom J. regia 700K single nucleotide polymorphism (SNP) genotyping array, enable the exploration of genomic prediction (GP) for this crop. This study is the first to evaluate GP accuracy and several influencing factors in walnut for traits related to phenology and nut quality. A core-collection of 170 accessions was phenotyped for 25 traits over 1 or 2 years. Highly heritable traits, such as budbreak date and female flowering date, were predicted with high accuracy (∼0.75) using ridge regression best linear unbiased prediction (rrBLUP). Three key factors influencing GP performance were examined: marker density, prediction model, and training set size. Selecting the top 1% of 364,275 SNPs based on their variance (∼3600 SNPs) was sufficient to achieve accurate predictions. Bayesian models slightly improved prediction accuracy for some traits when using this reduced SNP set, but rrBLUP provided robust results, balancing accuracy, simplicity, and computational efficiency. Training population size also influenced accuracy, with a subset comprising 50% of the population still yielding reliable predictions. Optimization of training set was assessed using coefficient of determination mean, prediction error variance mean, and mean relatedness (MeanRel) parameters, with MeanRel performing best for shell traits. However, incorporating SNPs identified in genome-wide association study into the prediction models did not enhance accuracy. In summary, this study demonstrates the feasibility and potential of GPs for walnut breeding programs using a core-collection, offering valuable insights for optimizing GP approaches in this crop.

Integrated multi-omics analyses provide new insights into genomic variation landscape and regulatory network candidate genes associated with walnut endocarp

Persian walnut (Juglans regia) is an economically important nut oil tree; the fruit has a hard endocarp/shell to protect seeds, thus playing a key role in its evolution, and the shell thickness is an important trait for walnut breeding. However, the genomic landscape and the gene regulatory networks associated with walnut shell development remain to be systematically elucidated. Here, we report a high-quality genome assembly of the walnut cultivar 'Xiangling' and construct a graphic structure pan-genome of eight Juglans species to reveal the genetic variations at the genome level. We re-sequence 285 accessions to characterize the genomic variation landscape. Through genome-wide association studies (GWAS), we identified 19 loci associated with more than 268 loci that underwent selection during walnut domestication and improvement. Multi-omics analyses, including transcriptomics, metabolomics, DNA methylation, and spatial transcriptomics across eleven developmental stages, revealed several candidate genes related to secondary cell biosynthesis and lignin accumulation. This integrated multi-omics approach revealed several candidate genes associated with secondary cell biosynthesis and lignin accumulation, such as UGP, MYB308, MYB83, NAC043, NAC073, CCoAOMT1, CCoAOMT7, CHS2, CESA7, LAC7, COBL4, and IRX12. Overexpression of JrUGP and JrMYB308 in Arabidopsis thaliana confirmed their roles in lignin biosynthesis and cell wall thickening. Consequently, our comprehensive multi-omics findings offer novel insights into walnut genetic variation and network regulation of endocarp development and shell thickness, which enable further genome-informed breeding strategies for walnut cultivar improvement.

Multi-omics analysis reveals the allelic variation in JrWDRC2A9 and JrGPIAP conferring resistance against anthracnose (Colletotrichum gloeosporioides) in walnut (Juglans regia)

Walnut anthracnose induced by Colletotrichum gloeosporioides is a devastating disease that seriously threatens walnut cultivation. Screening novel resistance genes and exploring the molecular mechanisms are essential for disease-resistant genetic improvement of walnut. We conducted a genome-wide association studies of disease resistance traits based on the relative resistance index and single nucleotide polymorphisms (SNPs) obtained from 182 resequenced walnut accessions and 10 loci and corresponding candidate genes associated with resistance against C. gloeosporioides were identified. Then, through combined transcriptome analysis during C. gloeosporioides infection and qRT-PCR, we identified JrWDRC2A9 in SNP Chr13_36265784 loci and JrGPIAP in SNP Chr07_10106470 loci as two walnut anthracnose resistance genes. The validation of the disease resistance function of transgenic strains indicated that both JrWDRC2A9 and JrGPIAP promote walnut resistance to anthracnose. SNP Chr13_36265784 (A>G) is located in the coding region of JrWDRC2A9 causing a glutamine (JrWDRC2A9HapI) to arginine (JrWDRC2A9HapII). Allelic variation in the WD domain attenuates JrWDRC2A9-mediated resistance against C. gloeosporioides and the binding affinity of JrWDRC2A9 for JrTLP1. On the contrary, the allelic variation caused by SNP Chr07_10106470 (T>G) increased the walnut accessions resistance to C. gloeosporioides by promoting the expression level of JrGPIAP. Functional genomics revealed that JrGPIAP binds to the promoter of JrPR1L and activates its transcription, which is strengthened by the interaction between JrGPIAP and JrEMP24. These findings reveal the allelic variation in JrWDRC2A9 and JrGPIAP conferring resistance against C. gloeosporioides, providing a genetic basis for walnut disease resistance breeding in the future.

Leveraging Walnut Somatic Embryos as a Biomanufacturing Platform for Recombinant Proteins and Metabolites

Biomanufacturing enables novel sources of compounds with constant demand, such as food coloring and preservatives, as well as new compounds with peak demand, such as diagnostics and vaccines. The COVID-19 pandemic has highlighted the need for alternative sources of research materials, thrusting research on diversification of biomanufacturing platforms. Here, we show initial results exploring the walnut somatic embryogenic system expressing the recombinant receptor binding domain (RBD) and ectodomain of the spike protein (Spike) from the SARS-CoV-2 virus. Stably transformed walnut embryo lines were selected and propagated in vitro. Both recombinant proteins were detected at 3-14 µg/g dry weight of tissue culture material. Although higher yields of recombinant protein have been obtained using more conventional biomanufacturing platforms, we also report on the production of the red pigment betanin in somatic embryos, reaching yields of 650 mg/g, even higher than red beet Beta vulgaris. This first iteration shows the potential of biomanufacturing using somatic walnut embryos that can now be further optimized for different applications sourcing specialized proteins and metabolites.

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.

JrPPO1/2 play distinct roles in regulating walnut fruit browning by different spatiotemporal expression and enzymatic characteristics

Polyphenol oxidase (PPO) activity drives walnut fruit browning, but the roles of its only two-family genes, JrPPO1 and JrPPO2, remain unclear. This study explores the spatiotemporal expression and enzymatic characteristics of JrPPO1 and JrPPO2 in walnut. Treatment with the PPO activator CuSO4 and H2O2 accelerated fruit browning and up-regulated JrPPO1/2 expression, whereas treatment with the PPO inhibitor ascorbic acid delayed browning, down-regulating JrPPO1 and up-regulating JrPPO2 expression. Compared to mJrPPO1, mJrPPO2 can exhibited better enzyme activity at higher temperatures (47 °C) and in more acidic environments (pH 4.25). mJrPPO2 exhibited a higher substrate specificity over mJrPPO1, and the preferred substrates are catechol, chlorogenic acid, and epicatechin. Additionally, mJrPPO2 adapted better to low concentration of oxygen (as low as 1.0% O2) and slightly elevated CO2 levels compared to mJrPPO1. Subcellular localization and spatiotemporal expression patterns showed that JrPPO1 is only expressed in green tissues and located in chloroplasts, while JrPPO2 is also located in chloroplasts, partly associated with membranes, and is expressed in both green and non-green tissues. Silencing JrPPO1/2 with virus-induced gene silencing (VIGS) reduced fruit browning, maintained higher total phenols, and decreased MDA production. Notably, silencing JrPPO1 had a greater impact on browning than JrPPO2, indicating JrPPO1's greater contribution to PPO activity and fruit browning in walnut fruits. Consequently, JrPPO1 can be effectively regulated both at the molecular level and by manipulating environmental conditions, to achieve the objective of controlling fruit browning.

Loss of Lateral suppressor gene is associated with evolution of root nodule symbiosis in Leguminosae

BACKGROUND

Root nodule symbiosis (RNS) is a fascinating evolutionary event. Given that limited genes conferring the evolution of RNS in Leguminosae have been functionally validated, the genetic basis of the evolution of RNS remains largely unknown. Identifying the genes involved in the evolution of RNS will help to reveal the mystery.

RESULTS

Here, we investigate the gene loss event during the evolution of RNS in Leguminosae through phylogenomic and synteny analyses in 48 species including 16 Leguminosae species. We reveal that loss of the Lateral suppressor gene, a member of the GRAS-domain protein family, is associated with the evolution of RNS in Leguminosae. Ectopic expression of the Lateral suppressor (Ls) gene from tomato and its homolog MONOCULM 1 (MOC1) and Os7 from rice in soybean and Medicago truncatula result in almost completely lost nodulation capability. Further investigation shows that Lateral suppressor protein, Ls, MOC1, and Os7 might function through an interaction with NODULATION SIGNALING PATHWAY 2 (NSP2) and CYCLOPS to repress the transcription of NODULE INCEPTION (NIN) to inhibit the nodulation in Leguminosae. Additionally, we find that the cathepsin H (CTSH), a conserved protein, could interact with Lateral suppressor protein, Ls, MOC1, and Os7 and affect the nodulation.

CONCLUSIONS

This study sheds light on uncovering the genetic basis of the evolution of RNS in Leguminosae and suggests that gene loss plays an essential role.

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.

Transcriptomic and coexpression network analyses revealed the regulatory mechanism of Cydia pomonella infestation on the synthesis of phytohormones in walnut husks

The codling moth (Cydia pomonella) has a major effect on the quality and yield of walnut fruit. Plant defences respond to insect infestation by activating hormonal signalling and the flavonoid biosynthetic pathway. However, little is known about the role of walnut husk hormones and flavonoid biosynthesis in response to C. pomonella infestation. The phytohormone content assay revealed that the contents of salicylic acid (SA), abscisic acid (ABA), jasmonic acid (JA), jasmonic acid-isoleucine conjugate (JA-ILE), jasmonic acid-valine (JA-Val) and methyl jasmonate (MeJA) increased after feeding at different time points (0, 12, 24, 36, 48, and 72 h) of walnut husk. RNA-seq analysis of walnut husks following C. pomonella feeding revealed a temporal pattern in differentially expressed genes (DEGs), with the number increasing from 3,988 at 12 h to 5,929 at 72 h postfeeding compared with the control at 0 h postfeeding. Walnut husks exhibited significant upregulation of genes involved in various defence pathways, including flavonoid biosynthesis (PAL, CYP73A, 4CL, CHS, CHI, F3H, ANS, and LAR), SA (PAL), ABA (ZEP and ABA2), and JA (AOS, AOC, OPR, JAZ, and MYC2) pathways. Three gene coexpression networks that had a significant positive association with these hormonal changes were constructed based on the basis of weighted gene coexpression network analysis (WGCNA). We identified several hub transcription factors, including the turquoise module (AIL6, MYB4, PRE6, WRKY71, WRKY31, ERF003, and WRKY75), the green module (bHLH79, PCL1, APRR5, ABI5, and ILR3), and the magenta module (ERF27, bHLH35, bHLH18, TIFY5A, WRKY31, and MYB44). Taken together, these findings provide useful genetic resources for exploring the defence response mediated by phytohormones in walnut husks.

Evolution of the Correlated Genomic Variation Landscape Across a Divergence Continuum in the Genus Castanopsis

The heterogeneous landscape of genomic variation has been well documented in population genomic studies. However, disentangling the intricate interplay of evolutionary forces influencing the genetic variation landscape over time remains challenging. In this study, we assembled a chromosome-level genome for Castanopsis eyrei and sequenced the whole genomes of 276 individuals from 12 Castanopsis species, spanning a broad divergence continuum. We found highly correlated genomic variation landscapes across these species. Furthermore, variations in genetic diversity and differentiation along the genome were strongly associated with recombination rates and gene density. These results suggest that long-term linked selection and conserved genomic features have contributed to the formation of a common genomic variation landscape. By examining how correlations between population summary statistics change throughout the species divergence continuum, we determined that background selection alone does not fully explain the observed patterns of genomic variation; the effects of recurrent selective sweeps must be considered. We further revealed that extensive gene flow has significantly influenced patterns of genomic variation in Castanopsis species. The estimated admixture proportion correlated positively with recombination rate and negatively with gene density, supporting a scenario of selection against gene flow. Additionally, putative introgression regions exhibited strong signals of positive selection, an enrichment of functional genes, and reduced genetic burdens, indicating that adaptive introgression has played a role in shaping the genomes of hybridizing species. This study provides insights into how different evolutionary forces have interacted in driving the evolution of the genomic variation landscape.

Genetic diversity and population structure of Botryosphaeria dothidea and Neofusicoccum parvum on English walnut (Juglans regia L.) in France

Botryosphaeriaceae species are the major causal agents of walnut dieback worldwide, along with Diaporthe species. Botryosphaeria dothidea and Neofusicoccum parvum are the only two Botryosphaeriaceae species associated with this recently emergent disease in France, and little is known about their diversity, structure, origin and dispersion in French walnut orchards. A total of 381 isolates of both species were genetically typed using a sequence-based microsatellite genotyping (SSR-seq) method. This analysis revealed a low genetic diversity and a high clonality of these populations, in agreement with their clonal mode of reproduction. The genetic similarity among populations, regardless of the tissue type and the presence of symptoms, supports the hypothesis that these pathogens can move between fruits and twigs and display latent pathogen lifestyles. Contrasting genetic patterns between N. parvum populations from Californian and Spanish walnut orchards and the French ones suggested no conclusive evidence for pathogen transmission from infected materials. The high genetic similarity with French vineyards populations suggested instead putative transmission between these hosts, which was also observed with B. dothidea populations. Overall, this study provides critical insight into the epidemiology of two important pathogens involved in the emerging dieback of French walnut orchards, including their distribution, potential to mate, putative origin and disease pathways.

JrPHL8-JrWRKY4-JrSTH2L module regulates resistance to Colletotrichum gloeosporioides in walnut

Walnut anthracnose (Colletotrichum gloeosporioides) reduces walnut yield and quality and seriously threatens the healthy development of the walnut industry. WRKY transcription factors (TFs) are crucial regulatory factors involved in plant-pathogen interactions. Our previous transcriptome analysis results indicate that JrWRKY4 responds to infection by C. gloeosporioides, but its specific regulatory network and disease resistance mechanism are still unclear. Herein, the characteristics of JrWRKY4 as a transcription activator located in the nucleus were first identified. Gain-of-function and loss-of-function analyses showed that JrWRKY4 could enhance walnut resistance against C. gloeosporioides. A series of molecular experiments showed that JrWRKY4 directly interacted with the promoter region of JrSTH2L and positively regulated its expression. In addition, JrWRKY4 interacted with JrVQ4 to form the protein complex, which inhibited JrWRKY4 for the activation of JrSTH2L. Notably, a MYB TF JrPHL8 interacting with the JrWRKY4 promoter has also been identified, which directly bound to the MBS element in the promoter of JrWRKY4 and induced its activity. Our study elucidated a novel mechanism of the JrPHL8-JrWRKY4-JrSTH2L in regulating walnut resistance to anthracnose. This mechanism improves our understanding of the molecular mechanism of WRKY TF mediated resistance to anthracnose in walnut, which provides new insights for molecular breeding of disease-resistant walnuts in the future.

Genomic basis of the distinct biosynthesis of β-glucogallin, a biochemical marker for hydrolyzable tannin production, in three oak species

Hydrolyzable tannins (HTs), predominant polyphenols in oaks, are widely used in grape wine aging, feed additives, and human healthcare. However, the limited availability of a high-quality reference genome of oaks greatly hampered the recognition of the mechanism of HT biosynthesis. Here, high-quality reference genomes of three Asian oak species (Quercus variabilis, Quercus aliena, and Quercus dentata) that have different HT contents were generated. Multi-omics studies were carried out to identify key genes regulating HT biosynthesis. In vitro enzyme activity assay was also conducted. Dual-luciferase and yeast one-hybrid assays were used to reveal the transcriptional regulation. Our results revealed that β-glucogallin was a biochemical marker for HT production in the cupules of the three Asian oaks. UGT84A13 was confirmed as the key enzyme for β-glucogallin biosynthesis. The differential expression of UGT84A13, rather than enzyme activity, was the main reason for different β-glucogallin and HT accumulation. Notably, sequence variations in UGT84A13 promoters led to different trans-activating activities of WRKY32/59, explaining the different expression patterns of UGT84A13 among the three species. Our findings provide three high-quality new reference genomes for oak trees and give new insights into different transcriptional regulation for understanding β-glucogallin and HT biosynthesis in closely related oak species.

Genome-Wide Analysis of Transcription Factor R2R3-MYB Gene Family and Gene Expression Profiles during Anthocyanin Synthesis in Common Walnut (Juglans regia L.)

The R2R3-MYB gene family, encoding plant transcriptional regulators, participates in many metabolic pathways of plant physiology and development, including flavonoid metabolism and anthocyanin synthesis. This study proceeded as follows: the JrR2R3-MYB gene family was analyzed genome-wide, and the family members were identified and characterized using the high-quality walnut reference genome "Chandler 2.0". All 204 JrR2R3-MYBs were established and categorized into 30 subgroups via phylogenetic analysis. JrR2R3-MYBs were unevenly distributed over 16 chromosomes. Most JrR2R3-MYBs had similar structures and conservative motifs. The cis-acting elements exhibit multiple functions of JrR2R3-MYBs such as light response, metabolite response, and stress response. We found that the expansion of JrR2R3-MYBs was mainly caused by WGD or segmental duplication events. Ka/Ks analysis indicated that these genes were in a state of negative purifying selection. Transcriptome results suggested that JrR2R3-MYBs were widely entangled in the process of walnut organ development and differentially expressed in different colored varieties of walnuts. Subsequently, we identified 17 differentially expressed JrR2R3-MYBs, 9 of which may regulate anthocyanin biosynthesis based on the results of a phylogenetic analysis. These genes were present in greater expression levels in 'Zijing' leaves than in 'Lvling' leaves, as revealed by the results of qRT-PCR experiments. These results contributed to the elucidation of the functions of JrR2R3-MYBs in walnut coloration. Collectively, this work provides a foundation for exploring the functional characteristics of the JrR2R3-MYBs in walnuts and improving the nutritional value and appearance quality of walnuts.

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.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

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.

A Multiplex PCR System of Novel Microsatellite Loci for Population Genetic Application in Walnuts

Multiplex polymerase chain reaction (PCR) of microsatellite loci allows for simultaneous amplification of two or more pairs of primers in a single PCR reaction; hence, it is cost and time effective. However, very few attempts have been reported in non-model species. In this study, by combining a genome-based de novo development and cross-species application approach, a multiplex PCR system comprising 5 PCR reactions of 33 microsatellites consisting of 26 novel genomic and 7 literature-sourced loci was tested for polymorphisms, cross-species transferability, and the ability to assess genetic diversity and population structure of three walnut species (Juglans spp.). We found that the genome-based approach is more efficient than other methods. An allelic ladder was developed for each locus to enhance consistent genotyping among laboratories. The population genetic analysis results showed that all 33 loci were successfully transferred across the three species, showing high polymorphism and a strong genetic structure. Hence, the multiplex PCR system is highly applicable in walnut species. Furthermore, we propose an efficient pipeline to characterize and genotype polymorphic microsatellite loci. The novel toolbox developed here will aid future ecology and evolution studies in walnut and could serve as a model for other plant species.

Identification and evolution of nsLTPs in the root nodule nitrogen fixation clade and molecular response of Frankia to AgLTP24

Non-specific lipid transfer proteins (nsLTPs) are antimicrobial peptides, involved in several plant biological processes including root nodule nitrogen fixation (RNF). Nodulating plants belonging to the RNF clade establish symbiosis with the nitrogen-fixing bacteria rhizobia (legumes symbiosis model) and Frankia (actinorhizal symbiosis model) leading to root nodule formation. nsLTPs are involved in processes active in early step of symbiosis and functional nodule in both models. In legumes, nsLTPs have been shown to regulate symbiont entry, promote root cortex infection, membrane biosynthesis, and improve symbiosis efficiency. More recently, a nsLTP, AgLTP24 has been described in the context of actinorhizal symbiosis between Alnus glutinosa and Frankia alni ACN14a. AgLTP24 is secreted at an early step of symbiosis on the deformed root hairs and targets the symbiont in the nitrogen-fixing vesicles in functional nodules. nsLTPs are involved in RNF, but their functions and evolutionary history are still largely unknown. Numerous putative nsLTPs were found up-regulated in functional nodules compared to non-infected roots in different lineages within the RNF clade. Here, results highlight that nodulating plants that are co-evolving with their nitrogen-fixing symbionts appear to have independently specialized nsLTPs for this interaction, suggesting a possible convergence of function, which opens perspectives to investigate nsLTPs functions in RNF.

DNA methylation role in subgenome expression dominance of Juglans regia and its wild relative J. mandshurica

Subgenome expression dominance plays a crucial role in the environmental adaptation of polyploids. However, the epigenetic molecular mechanism underlying this process has not been thoroughly investigated, particularly in perennial woody plants. Persian walnut (Juglans regia) and its wild relative, Manchurian walnut (Juglans mandshurica), are woody plants of great economic importance and are both paleopolyploids that have undergone whole-genome duplication events. In this study, we explored the characteristics of subgenome expression dominance in these 2 Juglans species and examined its epigenetic basis. We divided their genomes into dominant subgenome (DS) and submissive subgenome (SS) and found that the DS-specific genes might play critical roles in biotic stress response or pathogen defense. We comprehensively elucidated the characteristics of biased gene expression, asymmetric DNA methylation, transposable elements (TEs), and alternative splicing (AS) events of homoeologous gene pairs between subgenomes. The results showed that biased expression genes (BEGs) in 2 Juglans species were mainly related to external stimuli response, while non-BEGs were related to complexes that might be involved in signal transduction. DS genes had higher expression and more AS events while having less DNA methylation and TEs than homoeologous genes from the SS in the 2 Juglans species. Further studies showed that DNA methylation might contribute to the biased expression of gene pairs by modifying LTR/TIR/nonTIR TEs and improving the AS efficiency of corresponding precursor mRNAs in a particular context. Our study contributes to understanding the epigenetic basis of subgenome expression dominance and the environmental adaptation of perennial woody plants.

Molecular mechanism of thiamine in mitigating drought stress in Chinese wingnut (Pterocarya stenoptera): Insights from transcriptomics

Urban garden plants are frequently affected by drought, which can hinder their growth, development, and greening effect. Previous studies have indicated that Chinese wingnut (Pterocarya stenoptera) responds to drought stress by increasing the expression of thiamine synthesis genes. In this study, it was found that exogenous thiamine can effectively alleviate the negative effects of drought stress on plants. Forward transcriptome sequencing and physiological tests were further conducted to reveal the molecular mechanism of thiamine in alleviating drought stress. Results showed that exogenous thiamine activated the expression of eight chlorophyll synthesis genes in Chinese wingnut under drought stress. Moreover, physiological indicators proved that chlorophyll content increased in leaves of Chinese wingnut with thiamine treatment under drought stress. Photosynthesis genes were also activated in Chinese wingnut treated with exogenous thiamine under drought stress, as supported by photosynthetic indicators PIabs and PItotal. Additionally, exogenous thiamine stimulated the expression of genes in the auxin-activated signaling pathway, thus attenuating the effects of drought stress. This study demonstrates the molecular mechanism of thiamine in mitigating the effects of drought stress on non-model woody plants lacking transgenic systems. This study also provides an effective method to mitigate the negative impacts of drought stress on plants.

Genomic characterization of the NAC transcription factors, directed at understanding their functions involved in endocarp lignification of iron walnut (Juglans sigillata Dode)

The NAC (NAM, ATAF1/2, and CUC2) transcription factors (TF), one of the largest plant-specific gene families, play important roles in the regulation of plant growth and development, stress response and disease resistance. In particular, several NAC TFs have been identified as master regulators of secondary cell wall (SCW) biosynthesis. Iron walnut (Juglans sigillata Dode), an economically important nut and oilseed tree, has been widely planted in the southwest China. The thick and high lignified shell derived endocarp tissues, however, brings troubles in processing processes of products in industry. It is indispensable to dissect the molecular mechanism of thick endocarp formation for further genetic improvement of iron walnut. In the present study, based on genome reference of iron walnut, 117 NAC genes, in total, were identified and characterized in silico, which involves only computational analysis to provide insight into gene function and regulation. We found that the amino acids encoded by these NAC genes varied from 103 to 1,264 in length, and conserved motif numbers ranged from 2 to 10. The JsiNAC genes were unevenly distributed across the genome of 16 chromosomes, and 96 of these genes were identified as segmental duplication genes. Furthermore, 117 JsiNAC genes were divided into 14 subfamilies (A-N) according to the phylogenetic tree based on NAC family members of Arabidopsis thaliana and common walnut (Juglans regia). Furthermore, tissue-specific expression pattern analysis demonstrated that a majority of NAC genes were constitutively expressed in five different tissues (bud, root, fruit, endocarp, and stem xylem), while a total of 19 genes were specifically expressed in endocarp, and most of them also showed high and specific expression levels in the middle and late stages during iron walnut endocarp development. Our result provided a new insight into the gene structure and function of JsiNACs in iron walnut, and identified key candidate JsiNAC genes involved in endocarp development, probably providing mechanistic insight into shell thickness formation across nut species.

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.

Global Transcriptomic Analyses Provide New Insight into the Molecular Mechanisms of Endocarp Formation and Development in Iron Walnut (Juglans sigillata Dode)

Iron walnut (Juglans sigillata Dode) is a native species in southwestern China that exhibits variation in both fruit morphology and shell thickness. However, the underlying molecular processes controlling hardened endocarp development in walnut has not yet been reported. Here, we generated transcriptional profiles of iron walnut endocarp at three developmental stages using “Dapao”, the most common commercial variety. Using pairwise comparisons between these three stages, a total of 8555 non-redundant differentially expressed genes (DEGs) were identified, and more than one-half of the total DEGs exhibited significant differential expression in stage I as compared with stage II or stage III, suggesting that the first stage may ultimately determine the final characteristics of the mature walnut shell. Furthermore, in the clustering analysis of the above DEGs, 3682, 2349, and 2388 genes exhibited the highest expression in stages I, II, and III, respectively. GO enrichment analysis demonstrated that the major transcriptional variation among the three developmental stages was caused by differences in cell growth, plant hormones, metabolic process, and phenylpropanoid metabolism. Namely, using the tissue-specific expression analysis and a gene co-expression network, we identified MADS-box transcription factor JsiFBP2 and bHLH transcription factor JsibHLH94 as candidate regulators of endocarp formation in the early stage, and JsiNAC56 and JsiMYB78 might play key roles in regulating the lignification process of endocarp in the late stage. This study provides useful information for further research to dissect the molecular mechanisms governing the shell formation and development of iron walnut.

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.

Genome structure-based Juglandaceae phylogenies contradict alignment-based phylogenies and substitution rates vary with DNA repair genes

In lineages of allopolyploid origin, sets of homoeologous chromosomes may coexist that differ in gene content and syntenic structure. Presence or absence of genes and microsynteny along chromosomal blocks can serve to differentiate subgenomes and to infer phylogenies. We here apply genome-structural data to infer relationships in an ancient allopolyploid lineage, the walnut family (Juglandaceae), by using seven chromosome-level genomes, two of them newly assembled. Microsynteny and gene-content analyses yield identical topologies that place Platycarya with Engelhardia as did a 1980s morphological-cladistic study. DNA-alignment-based topologies here and in numerous earlier studies instead group Platycarya with Carya and Juglans, perhaps misled by past hybridization. All available data support a hybrid origin of Juglandaceae from extinct or unsampled progenitors nested within, or sister to, Myricaceae. Rhoiptelea chiliantha, sister to all other Juglandaceae, contains proportionally more DNA repair genes and appears to evolve at a rate 2.6- to 3.5-times slower than the remaining species.

Genome-Wide Characterization of the Mitogen-Activated Protein Kinase Gene Family and Their Expression Patterns in Response to Drought and Colletotrichum Gloeosporioides in Walnut (Juglans regia)

Mitogen-activated protein kinases (MAPKs) are a family of Ser/Thr (serine/threonine) protein kinases that play very important roles in plant responses to biotic and abiotic stressors. However, the MAPK gene family in the important crop walnut (Juglans regia L.) has been less well studied compared with other species. We discovered 25 JrMAPK members in the Juglans genome in this study. The JrMAPK gene family was separated into four subfamilies based on phylogenetic analysis, and members of the same subgroup had similar motifs and exons/introns. A variety of cis-acting elements, mainly related to the light response, growth and development, stress response, and hormone responses, were detected in the JrMAPK gene promoters. Collinearity analysis showed that purification selection was the main driving force in JrMAPK gene evolution, and segmental and tandem duplications played key roles in the expansion of the JrMAPK gene family. The RNA-Seq (RNA Sequencing) results indicated that many of the JrMAPK genes were expressed in response to different levels of Colletotrichum gloeosporioides infection. JrMAPK1, JrMAPK3, JrMAPK4, JrMAPK5, JrMAPK6, JrMAPK7, JrMAPK9, JrMAPK11, JrMAPK12, JrMAPK13, JrMAPK17, JrMAPK19, JrMAPK20, and JrMAPK21 were upregulated at the transcriptional level in response to the drought stress treatment. The results of this study will help in further investigations of the evolutionary history and biological functions of the MAPK gene family in walnut.

Advances in genomics for diversity studies and trait improvement in temperate fruit and nut crops under changing climatic scenarios

Genetic improvement of temperate fruit and nut crops through conventional breeding methods is not sufficient alone due to its extreme time-consuming, cost-intensive, and hard-to-handle approach. Again, few other constraints that are associated with these species, viz., their long juvenile period, high heterozygosity, sterility, presence of sexual incompatibility, polyploidy, etc., make their selection and improvement process more complicated. Therefore, to promote precise and accurate selection of plants based on their genotypes, supplement of advanced biotechnological tools, viz., molecular marker approaches along with traditional breeding methods, is highly required in these species. Different markers, especially the molecular ones, enable direct selection of genomic regions governing the trait of interest such as high quality, yield, and resistance to abiotic and biotic stresses instead of the trait itself, thus saving the overall time and space and helping screen fruit quality and other related desired traits at early stages. The availability of molecular markers like SNP (single-nucleotide polymorphism), DArT (Diversity Arrays Technology) markers, and dense molecular genetic maps in crop plants, including fruit and nut crops, led to a revelation of facts from genetic markers, thus assisting in precise line selection. This review highlighted several aspects of the molecular marker approach that opens up tremendous possibilities to reveal valuable information about genetic diversity and phylogeny to boost the efficacy of selection in temperate fruit crops through genome sequencing and thus cultivar improvement with respect to adaptability and biotic and abiotic stress resistance in temperate fruit and nut species.

The core triacylglycerol toolbox in woody oil plants reveals targets for oil production bioengineering

Woody oil plants are the most productive oil-bearing species that produce seeds with high levels of valuable triacylglycerols (TAGs). TAGs and their derivatives are the raw materials for many macromolecular bio-based products, such as nylon precursors, and biomass-based diesel. Here, we identified 280 genes encoding seven distinct classes of enzymes (i.e., G3PAT, LPAAT, PAP, DGAT, PDCT, PDAT, and CPT) involved in TAGs-biosynthesis. Several multigene families are expanded by large-scale duplication events, such as G3PATs, and PAPs. RNA-seq was used to survey the expression profiles of these TAG pathway-related genes in different tissues or development, indicating functional redundancy for some duplicated genes originated from the large-scale duplication events, and neo-functionalization or sub-functionalization for some of them. Sixty-two genes showed strong, preferential expression during the period of rapid seed lipid synthesis, suggesting that their might represented the core TAG-toolbox. We also revealed for the first time that there is no PDCT pathway in Vernicia fordii and Xanthoceras sorbifolium. The identification of key genes involved in lipid biosynthesis will be the foundation to plan strategies to develop woody oil plant varieties with enhanced processing properties and high oil content.

523, Beni Mellal, Morocco, Bouda S, Ennahli S, Hafida H, Laboratory of Biotechnologies and Valorization of Plant Gnetic Resources, University of Sultan Moulay Slimane, Faculty of Sciences and Techniques, P.B. 523, Beni Mellal, Morocco, Departement of Horticulture and Viticulture, National School of Agriculture, Meknes, Morocco, Laboratory of Development and Safety of Food Products, Faculty of Sciences and Techniques, University of Sultan Moulay Slimane Beni Mellal, Morocco. 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.

Domestication and selection footprints in Persian walnuts (Juglans regia)

Walnut (Juglans) species are economically important hardwood trees cultivated worldwide for both edible nuts and high-quality wood. Broad-scale assessments of species diversity, evolutionary history, and domestication are needed to improve walnut breeding. In this study, we sequenced 309 walnut accessions from around the world, including 55 Juglans relatives, 98 wild Persian walnuts (J. regia), 70 J. regia landraces, and 86 J. regia cultivars. The phylogenetic tree indicated that J. regia samples (section Dioscaryon) were monophyletic within Juglans. The core areas of genetic diversity of J. regia germplasm were southwestern China and southern Asia near the Qinghai-Tibet Plateau and the Himalayas, and the uplift of the Himalayas was speculated to be the main factor leading to the current population dynamics of Persian walnut. The pattern of genomic variation in terms of nucleotide diversity, linkage disequilibrium, single nucleotide polymorphisms, and insertions/deletions revealed the domestication and selection footprints in Persian walnut. Selective sweep analysis, GWAS, and expression analysis further identified two transcription factors, JrbHLH and JrMYB6, that influence the thickness of the nut diaphragm as loci under selection during domestication. Our results elucidate the domestication and selection footprints in Persian walnuts and provide a valuable resource for the genomics-assisted breeding of this important crop.

Genomic and Bioinformatic Resources for Perennial Fruit Species

In the post-genomic era, data management and development of bioinformatic tools are critical for the adequate exploitation of genomics data. In this review, we address the actual situation for the subset of crops represented by the perennial fruit species. The agronomical singularity of these species compared to plant and crop model species provides significant challenges on the implementation of good practices generally not addressed in other species. Studies are usually performed over several years in non-controlled environments, usage of rootstock is common, and breeders heavily rely on vegetative propagation. A reference genome is now available for all the major species as well as many members of the economically important genera for breeding purposes. Development of pangenome for these species is beginning to gain momentum which will require a substantial effort in term of bioinformatic tool development. The available tools for genome annotation and functional analysis will also be presented.

JrWRKY21 interacts with JrPTI5L to activate the expression of JrPR5L for resistance to Colletotrichum gloeosporioides in walnut

Walnut (Juglans regia L.) anthracnose, induced by Colletotrichum gloeosporioides, is a catastrophic disease impacting the walnut industry in China. Although WRKY transcription factors play a key role in plant immunity, the function of the WRKY gene family in walnut resistance to C. gloeosporioides is not clear. Here, through transcriptome sequencing and quantitative real-time polymerase chain reaction (qRT-PCR), we identified a differentially expressed gene, JrWRKY21, that was significantly upregulated upon C. gloeosporioides infection in walnut. JrWRKY21 positively regulated walnut resistance to C. gloeosporioides, as demonstrated by virus-induced gene silencing and transient gene overexpression. Additionally, JrWRKY21 directly interacted with the transcriptional activator of the pathogenesis-related (PR) gene JrPTI5L in vitro and in vivo, and could bind to the W-box in the JrPTI5L promoter for transcriptional activation. Moreover, JrPTI5L could induce the expression of the PR gene JrPR5L through binding to the GCCGAC motif in the promoter. Our data support that JrWRKY21 can indirectly activate the expression of the JrPR5L gene via the WRKY21-PTI5L protein complex to promote resistance against C. gloeosporioides in walnut. The results will enhance our understanding of the mechanism behind walnut disease resistance and facilitate the genetic improvement of walnut by molecular breeding for anthracnose-resistant varieties.

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.

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.

Genetic Diversity and Structure of Persian Walnut (Juglans regia L.) in Pakistan: Implications for Conservation

Persian (Common) walnut (Juglans regia L.) is a famous fruit tree species valued for its nutritious nuts and high-quality wood. Although walnut is widely distributed and plays an important role in the economy and culture of Pakistan, the genetic diversity and structure of its populations in the country remains poorly understood. Therefore, using 31 nuclear microsatellites, we assessed the genetic diversity and population structure of 12 walnut populations sampled across Pakistan. We also implemented the geostatistical IDW technique in ArcGIS to reveal "hotspots" of genetic diversity. Generally, the studied populations registered relatively low indices of genetic diversity (NA = 3.839, HO = 0.558, UHE = 0.580), and eight populations had positive inbreeding coefficient (FIS) values. Low among-population differentiation was indicated by AMOVA, pairwise FST and DC. STRUCTURE, PCoA and neighbor joining (NJ) analysis revealed a general lack of clear clustering in the populations except that one population in Upper Dir was clearly genetically distinct from the rest. Furthermore, the Mantel test showed no correlation between the geographic and genetic distance (r = 0.14, p = 0.22), while barrier analysis suggested three statistically significant genetic barriers. Finally, the spatial interpolation results indicated that populations in Ziarat, Kashmir, Dir, Swat, Chitral, and upper Dir had high intrapopulation genetic diversity, suggesting the need to conserve populations in those areas. The results from this study will be important for future breeding improvement and conservation of walnuts in Pakistan.

A Chromosome-Level Genome Assembly of the European Beech (Fagus sylvatica) Reveals Anomalies for Organelle DNA Integration, Repeat Content and Distribution of SNPs

The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here, we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome-wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate.

Identification and Characterization of Wall-Associated Kinase (WAK) and WAK-like (WAKL) Gene Family in Juglans regia and Its Wild Related Species Juglans mandshurica

Wall-associated kinase (WAK) and WAK-like kinase (WAKL) are receptor-like kinases (RLKs), which play important roles in signal transduction between the cell wall and the cytoplasm in plants. WAK/WAKLs have been studied in many plants, but were rarely studied in the important economic walnut tree. In this study, 27 and 14 WAK/WAKL genes were identified in Juglans regia and its wild related species Juglans mandshurica, respectively. We found tandem duplication might play a critical role in the expansion of WAK/WAKL gene family in J. regia, and most of the WAK/WAKL homologous pairs underwent purified selection during evolution. All WAK/WAKL proteins have the extracellular WAK domain and the cytoplasmic protein kinase domain, and the latter was more conserved than the former. Cis-acting elements analysis showed that WAK/WAKL might be involved in plant growth and development, plant response to abiotic stress and hormones. Gene expression pattern analysis further indicated that most WAK/WAKL genes in J. regia might play a role in the development of leaves and be involved in plant response to biotic stress. Our study provides a new perspective for the evolutionary analysis of gene families in tree species and also provides potential candidate genes for studying WAK/WAKL gene function in walnuts.

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.

Xanthomonas arboricola pv. juglandis and pv. corylina: Brothers or distant relatives? Genetic clues, epidemiology, and insights for disease management

BACKGROUND

The species Xanthomonas arboricola comprises up to nine pathovars, two of which affect nut crops: pv. juglandis, the causal agent of walnut bacterial blight, brown apical necrosis, and the vertical oozing canker of Persian (English) walnut; and pv. corylina, the causal agent of the bacterial blight of hazelnut. Both pathovars share a complex population structure, represented by different clusters and several clades. Here we describe our current understanding of symptomatology, population dynamics, epidemiology, and disease control.

TAXONOMIC STATUS

Bacteria; Phylum Proteobacteria; Class Gammaproteobacteria; Order Lysobacterales (earlier synonym of Xanthomonadales); Family Lysobacteraceae (earlier synonym of Xanthomonadaceae); Genus Xanthomonas; Species X. arboricola; Pathovars: pv. juglandis and pv. corylina.

HOST RANGE AND SYMPTOMS

The host range of each pathovar is not limited to a single species, but each infects mainly one plant species: Juglans regia (X. arboricola pv. juglandis) and Corylus avellana (X. arboricola. pv. corylina). Walnut bacterial blight is characterized by lesions on leaves and fruits, and cankers on twigs, branches, and trunks; brown apical necrosis symptoms consist of apical necrosis originating at the stigmatic end of the fruit. A peculiar symptom, the vertical oozing canker developing along the trunk, is elicited by a particular genetic lineage of the bacterium. Symptoms of hazelnut bacterial blight are visible on leaves and fruits as necrotic lesions, and on woody parts as cankers. A remarkable difference is that affected walnuts drop abundantly, whereas hazelnuts with symptoms do not.

DISTRIBUTION

Bacterial blight of walnut has a worldwide distribution, wherever Persian (English) walnut is cultivated; the bacterial blight of hazelnut has a more limited distribution, although disease outbreaks are currently more frequently reported. X. arboricola pv. juglandis is regulated almost nowhere, whereas X. arboricola pv. corylina is regulated in most European and Mediterranean Plant Protection Organization (EPPO) countries.

EPIDEMIOLOGY AND CONTROL

For both pathogens infected nursery material is the main pathway for their introduction and spread into newly cultivated areas; additionally, infected nursery material is the source of primary inoculum. X. arboricola pv. juglandis is also disseminated through pollen. Disease control is achieved through the phytosanitary certification of nursery material (hazelnut), although approved certification schemes are not currently available. Once the disease is present in walnut/hazelnut groves, copper compounds are widely used, mostly in association with dithiocarbamates; where allowed, antibiotics (preferably kasugamycin) are sprayed. The emergence of strains highly resistant to copper currently represents the major threat for effective management of the bacterial blight of walnut. USEFUL WEBSITES: https://gd.eppo.int/taxon/XANTJU, https://gd.eppo.int/taxon/XANTCY, https://www.euroxanth.eu, http://www.xanthomonas.org.

Transcriptome profiling based on Illumina- and SMRT-based RNA-seq reveals circadian regulation of key pathways in flower bud development in walnut

Flower bud development is a defining feature of walnut, which contributes to the kernel yield, yield stability, fruit quality and commodity value. However, little is known about the mechanism of the flower bud development in walnut. Here, the stages of walnut female flower bud development were divided into five period (P01-05) by using histological observation. They were further studied through PacBio Iso-Seq and RNA-seq analysis. Accordingly, we obtained 52,875 full-length transcripts, where 4,579 were new transcripts, 3,065 were novel genes, 1,437 were consensus lncRNAs and 20,813 were alternatively spliced isoforms. These transcripts greatly improved the current genome annotation and enhanced our understanding of the walnut transcriptome. Next, RNA sequencing of female flower buds at five periods revealed that circadian rhythm-plant was commonly enriched along with the flower bud developmental gradient. A total of 14 differentially expressed genes (DEGs) were identified, and six of them were confirmed by real-time quantitative analysis. Additionally, six and two differentially expressed clock genes were detected to be regulated by AS events and lncRNAs, respectively. All these detected plant circadian genes form a complex interconnected network to regulate the flower bud development. Thus, investigation of key genes associated with the circadian clock could clarify the process of flower bud development in walnut.

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.

Chromosome-scale genome assembly of Castanopsis tibetana provides a powerful comparative framework to study the evolution and adaptation of Fagaceae trees

Fagaceae species are increasingly used as models to elucidate the process and mechanism of adaptation and speciation by integrating ecology, evolution and genomics. The genus Castanopsis belongs to the family Fagaceae and is mainly distributed across subtropical and tropical Asia. In the present study, we reported the first chromosome-scale genome assembly of Castanopsis tibetana, a common species of evergreen broadleaved forests in subtropical China. The combination of Nanopore sequencing and Hi-C technologies enabled a high-quality genome assembly. The final assembled genome size of C. tibetana was 878.6 Mb (97.6% of the estimated genome size), consisting of 477 contigs with an N50 length of 3.3 Mb. The benchmarking universal single-copy orthologue (BUSCO) assessment indicated a completeness of 93.0%. Hi-C scaffolding generated 12 pseudochromosomes, representing 98.7% of the assembled genome. Subsequently, 40,937 protein-coding genes were predicted and 90.04% of them were functionally annotated. More than 476.9 Mb of repetitive sequences (54.3% of the genome) were identified, and the percentage of the genome covered by TE elements was 39.98%. Comparative genomics analysis revealed that C. tibetana was most closely related to Castanea mollissima and diverged at 18.48 Ma, and that C. tibetana has undergone considerable gene family expansion and contraction. Evidence of positive selection was detected in 53 genes, which showed different arrangement pattern compared to Quercus robur. The chromosome-scale genome assembly of C. tibetana will expand Fagaceae genome resources across the family and provide a powerful comparative framework to study the adaptation and evolution of Fagaceae trees.

Whole genome based insights into the phylogeny and evolution of the Juglandaceae

BACKGROUND

The walnut family (Juglandaceae) contains commercially important woody trees commonly called walnut, wingnut, pecan and hickory. Phylogenetic relationships and diversification within the Juglandaceae are classic and hot scientific topics that have been elucidated by recent fossil, morphological, molecular, and (paleo) environmental data. Further resolution of relationships among and within genera is still needed and can be achieved by analysis of the variation of chloroplast, mtDNA, and nuclear genomes.

RESULTS

We reconstructed the backbone phylogenetic relationships of Juglandaceae using organelle and nuclear genome data from 27 species. The divergence time of Juglandaceae was estimated to be 78.7 Mya. The major lineages diversified in warm and dry habitats during the mid-Paleocene and early Eocene. The plastid, mitochondrial, and nuclear phylogenetic analyses all revealed three subfamilies, i.e., Juglandoideae, Engelhardioideae, Rhoipteleoideae. Five genera of Juglandoideae were strongly supported. Juglandaceae were estimated to have originated during the late Cretaceous, while Juglandoideae were estimated to have originated during the Paleocene, with evidence for rapid diversification events during several glacial and geological periods. The phylogenetic analyses of organelle sequences and nuclear genome yielded highly supported incongruence positions for J. cinerea, J. hopeiensis, and Platycarya strobilacea. Winged fruit were the ancestral condition in the Juglandoideae, but adaptation to novel dispersal and regeneration regimes after the Cretaceous-Paleogene boundary led to the independent evolution of zoochory among several genera of the Juglandaceae.

CONCLUSIONS

A fully resolved, strongly supported, time-calibrated phylogenetic tree of Juglandaceae can provide an important framework for studying classification, diversification, biogeography, and comparative genomics of plant lineages. Our addition of new, annotated whole chloroplast genomic sequences and identification of their variability informs the study of their evolution in walnuts (Juglandaceae).

Whole genome based insights into the phylogeny and evolution of the Juglandaceae

BACKGROUND

The walnut family (Juglandaceae) contains commercially important woody trees commonly called walnut, wingnut, pecan and hickory. Phylogenetic relationships and diversification within the Juglandaceae are classic and hot scientific topics that have been elucidated by recent fossil, morphological, molecular, and (paleo) environmental data. Further resolution of relationships among and within genera is still needed and can be achieved by analysis of the variation of chloroplast, mtDNA, and nuclear genomes.

RESULTS

We reconstructed the backbone phylogenetic relationships of Juglandaceae using organelle and nuclear genome data from 27 species. The divergence time of Juglandaceae was estimated to be 78.7 Mya. The major lineages diversified in warm and dry habitats during the mid-Paleocene and early Eocene. The plastid, mitochondrial, and nuclear phylogenetic analyses all revealed three subfamilies, i.e., Juglandoideae, Engelhardioideae, Rhoipteleoideae. Five genera of Juglandoideae were strongly supported. Juglandaceae were estimated to have originated during the late Cretaceous, while Juglandoideae were estimated to have originated during the Paleocene, with evidence for rapid diversification events during several glacial and geological periods. The phylogenetic analyses of organelle sequences and nuclear genome yielded highly supported incongruence positions for J. cinerea, J. hopeiensis, and Platycarya strobilacea. Winged fruit were the ancestral condition in the Juglandoideae, but adaptation to novel dispersal and regeneration regimes after the Cretaceous-Paleogene boundary led to the independent evolution of zoochory among several genera of the Juglandaceae.

CONCLUSIONS

A fully resolved, strongly supported, time-calibrated phylogenetic tree of Juglandaceae can provide an important framework for studying classification, diversification, biogeography, and comparative genomics of plant lineages. Our addition of new, annotated whole chloroplast genomic sequences and identification of their variability informs the study of their evolution in walnuts (Juglandaceae).

Phenoloxidases in Plants-How Structural Diversity Enables Functional Specificity

The metabolism of polyphenolic polymers is essential to the development and response to environmental changes of organisms from all kingdoms of life, but shows particular diversity in plants. In contrast to other biopolymers, whose polymerisation is catalysed by homologous gene families, polyphenolic metabolism depends on phenoloxidases, a group of heterogeneous oxidases that share little beyond the eponymous common substrate. In this review, we provide an overview of the differences and similarities between phenoloxidases in their protein structure, reaction mechanism, substrate specificity, and functional roles. Using the example of laccases (LACs), we also performed a meta-analysis of enzyme kinetics, a comprehensive phylogenetic analysis and machine-learning based protein structure modelling to link functions, evolution, and structures in this group of phenoloxidases. With these approaches, we generated a framework to explain the reported functional differences between paralogs, while also hinting at the likely diversity of yet undescribed LAC functions. Altogether, this review provides a basis to better understand the functional overlaps and specificities between and within the three major families of phenoloxidases, their evolutionary trajectories, and their importance for plant primary and secondary metabolism.

Genome-Wide Identification and Transcriptional Expression Profiles of Transcription Factor WRKY in Common Walnut (Juglans regia L.)

The transcription factor WRKY is widely distributed in the plant kingdom, playing a significant role in plant growth, development and response to stresses. Walnut is an economically important temperate tree species valued for both its edible nuts and high-quality wood, and its response to various stresses is an important factor that determines the quality of its fruit. However, in walnut trees themselves, information about the WRKY gene family remains scarce. In this paper, we perform a comprehensive study of the WRKY gene family in walnut. In total, we identified 103 WRKY genes in the common walnut that are clustered into 4 groups and distributed on 14 chromosomes. The conserved domains all contained a WRKY domain, and motif 2 was observed in most WRKYs, suggesting a high degree of conservation and similar functions within each subfamily. However, gene structure was significantly differentiated between different subfamilies. Synteny analysis indicates that there were 56 gene pairs in J. regia and A. thaliana, 76 in J. regia and J. mandshurica, 75 in J. regia and J. microcarpa, 76 in J. regia and P. trichocarpa, and 33 in J. regia and Q. robur, indicating that the WRKY gene family may come from a common ancestor. GO and KEGG enrichment analysis showed that the WRKY gene family was involved in resistance traits and the plant-pathogen interaction pathway. In anthracnose-resistant F26 fruits (AR) and anthracnose-susceptible F423 fruits (AS), transcriptome and qPCR analysis results showed that JrWRKY83, JrWRKY73 and JrWRKY74 were expressed significantly more highly in resistant cultivars, indicating that these three genes may be important contributors to stress resistance in walnut trees. Furthermore, we investigate how these three genes potentially target miRNAs and interact with proteins. JrWRKY73 was target by the miR156 family, including 12 miRNAs; this miRNA family targets WRKY genes to enhance plant defense. JrWRKY73 also interacted with the resistance gene AtMPK6, showing that it may play a crucial role in walnut defense.