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

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.

Genetic Mechanisms and Adaptive Benefits of Anthocyanin Red Stigmas in a Wind-Pollinated Tree

Anthocyanin accumulation in leaves or flowers mitigates photooxidation damage from reactive oxygen species (ROS) and functions in plant/animal interactions. Among the most conspicuously anthocyanin-accumulating tissues are stigmas, especially in wind-pollinated trees. In the walnut genus (Juglans), yellow stigmas are ancestral, but a few species have dark red stigmas. We have used a natural F1 hybrid resulting from crosses between yellow stigma and red stigma species to investigate the genetic basis of the red stigmas. We found that a Copia transposable element (TE) insertion in the ubiquitin-protein ligase gene MIEL1 suppresses its expression in stigmas through RNA-directed DNA methylation and has gone to fixation in red stigma species. A younger Gypsy TE insertion fully inhibits MIEL1 expression, but is not fixed, explaining the color segregation in hybrid populations. Based on reference genomes and whole-genome sequencing data representing 20 of the 22 species of Juglans, we traced the evolution of MIEL1, finding the insertions in all consistently red stigma species. Red stigmas had lower levels of ROS than yellow stigmas, and population genetic data reveal strong positive selection on the TE-bearing MIEL1 allele. In combination, these results suggest that anthocyanin-accumulating stigma tissues support pollen germination and growth by protecting cells from ROS.

Genome-wide association study dissects the genetic architecture of progesterone content in Persian walnut leaves (Juglans regia L.)

Progesterone (P4) is an endogenous sex steroid hormone involved in the ovulatory cycle and pregnancy of animal species. In sheep and goats, P4 analogues are used to induce synchronized ovulations and oestrus behavior of the females. In humans, P4 from chemical synthesis is used to treat peri-menopausal disorders. However, such molecules are released into aquatic environment and can be a source of pollution, are prohibited in organic farms and go against the trend of "naturality" in animal production as well as in human health. A natural alternative may consist in the extraction and use of P4 in plants. Mammalian hormones were discovered in an increasing number of plant species, including walnut leaves that contain high levels of P4. We compared the content of P4 in leaves of 170 accessions of Juglans regia from the walnut germplasm collection of INRAE Prunus-Juglans Biological Resources Center previously genotyped using the Axiom™ J. regia 700 K SNP array. We conducted a genome-wide association study (GWAS) using multi-locus models. When collected in October, P4 content goes from 34,1 to 287,5 mg/kg dry weight of leaves. The two laciniate accessions have the largest P4 content. We identified seven significant marker-trait associations on chromosomes 1, 2, 3, 6, 7, 15 and 16, and a candidate gene involved in the metabolism of sterols, precursors of plant steroid hormones. Our results raise the huge variability of P4 content within J. regia and propose a candidate gene which may have a role in the control of this variability, opening the way to a potential use of walnut P4 by the pharmaceutical industry towards more natural source of chemical compounds.

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.

Integrated analysis of the transcriptome and hormone metabolome elucidates the regulatory mechanisms governing walnut bud germination

The walnut (Juglans regia) is an important oilseed tree species characterized by its extensive distribution, high oil yield, and nutrient-dense kernels, which provide substantial economic benefits. However, the rising incidence of late-spring frosts, exacerbated by global climate change, has adversely affected walnut yields. A comprehensive understanding of the regulatory mechanisms involved in bud dormancy, germination, and development is essential for developing strategies to mitigate the effects of late-spring frosts and for breeding frost-resistant cultivars. This study focused on W13, a protogynous walnut variety with early germination of dormant buds in spring, employing a combination of transcriptomic and hormone metabolomic analyses. Our results emphasized four key biological processes-cellular response to ethylene stimulus, phenylpropanoid metabolic process, ethylene-activated signaling pathway, and monooxygenase activity-along with several relevant pathways, including plant hormone signal transduction, flavone and flavonol biosynthesis, biosynthesis of secondary metabolites, and MAPK signaling pathway, all crucial for walnut bud germination. Additionally, bud germination is closely associated with alterations in various hormone signaling pathways, including abscisic acid, auxin, cytokinin, ethylene, gibberellins, jasmonic acid, and salicylic acid. By assessing hormone levels and gene expression at different developmental stages, we pinpointed potential regulatory genes and critical hormones associated with bud germination. Furthermore, through weighted correlation network analysis, we constructed a co-expression network, identifying gene modules specifically expressed during dormancy, germination, budding, and leafing phases. The hub genes within these modules are likely pivotal in regulating walnut bud germination. Our analysis also revealed that genes from various transcription factor families are central within the co-expression network, indicating their significant roles in the bud germination process. Correlation network analysis of hormone and gene further illuminated the mechanisms through which genes and hormones jointly influence walnut bud germination. These findings establish a crucial molecular basis for a more comprehensive understanding of the mechanisms governing germination and development in dormant walnut buds.

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.

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.

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.

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.

Unveiling the evolutionary history of lingonberry (Vaccinium vitis-idaea L.) through genome sequencing and assembly of European and North American subspecies

Lingonberry (Vaccinium vitis-idaea L.) produces tiny red berries that are tart and nutty in flavor. It grows widely in the circumpolar region, including Scandinavia, northern parts of Eurasia, Alaska, and Canada. Although cultivation is currently limited, the plant has a long history of cultural use among indigenous communities. Given its potential as a food source, genomic resources for lingonberry are significantly lacking. To advance genomic knowledge, the genomes for 2 subspecies of lingonberry (V. vitis-idaea ssp. minus and ssp. vitis-idaea var. 'Red Candy') were sequenced and de novo assembled into contig-level assemblies. The assemblies were scaffolded using the bilberry genome (Vaccinium myrtillus) to generate a chromosome-anchored reference genome consisting of 12 chromosomes each with a total length of 548.07 Mb [contig N50 = 1.17 Mb, BUSCO (C%) = 96.5%] for ssp. vitis-idaea and 518.70 Mb [contig N50 = 1.40 Mb, BUSCO (C%) = 96.9%] for ssp. minus. RNA-seq-based gene annotation identified 27,243 and 25,718 genes on the respective assembly, and transposable element detection methods found that 45.82 and 44.58% of the genome were repeats. Phylogenetic analysis confirmed that lingonberry was most closely related to bilberry and was more closely related to blueberries than cranberries. Estimates of past effective population size suggested a continuous decline over the past 1-3 MYA, possibly due to the impacts of repeated glacial cycles during the Pleistocene leading to frequent population fragmentation. The genomic resource created in this study can be used to identify industry-relevant genes (e.g. anthocyanin production), infer phylogeny, and call sequence-level variants (e.g. SNPs) in future research.

Chromosome-level genome assembly of Platycarya strobilacea

Platycarya strobilacea belongs to the walnut family (Juglandaceae), is commonly known as species endemic to East Asia, and is an ecologically important, wind pollinated, woody deciduous tree. To facilitate this ancient tree for the ecological value and conservation of this ancient tree, we report a new high-quality genome assembly of P. strobilacea. The genome size was 677.30 Mb, with a scaffold N50 size of 45,791,698 bp, and 98.43% of the assembly was anchored to 15 chromosomes. We annotated 32,246 protein-coding genes in the genome, of which 96.30% were functionally annotated in six databases. This new high-quality assembly of P. strobilacea provide valuable resource for the phylogenetic and evolutionary analysis of the walnut family and angiosperm.

T2T reference genome assembly and genome-wide association study reveal the genetic basis of Chinese bayberry fruit quality

Chinese bayberry (Myrica rubra or Morella rubra; 2n = 16) produces fruit with a distinctive flavor, high nutritional, and economic value. However, previous versions of the bayberry genome lack sequence continuity. Moreover, to date, no large-scale germplasm resource association analysis has examined the allelic and genetic variations determining fruit quality traits. Therefore, in this study, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar 'Zaojia' using PacBio HiFi long reads. The resulting 292.60 Mb T2T genome, revealed 8 centromeric regions, 15 telomeres, and 28 345 genes. This represents a substantial improvement in the genome continuity and integrity of Chinese bayberry. Subsequently, we re-sequenced 173 accessions, identifying 6 649 674 single nucleotide polymorphisms (SNPs). Further, the phenotypic analyses of 29 fruit quality-related traits enabled a genome-wide association study (GWAS), which identified 1937 SNPs and 1039 genes significantly associated with 28 traits. An SNP cluster pertinent to fruit color was identified on Chr6: 3407532 to 5 153 151 bp region, harboring two MYB genes (MrChr6G07650 and MrChr6G07660), exhibiting differential expression in extreme phenotype transcriptomes, linked to anthocyanin synthesis. An adjacent, closely linked gene, MrChr6G07670 (MLP-like protein), harbored an exonic missense variant and was shown to increase anthocyanin production in tobacco leaves tenfold. This SNP cluster, potentially a quantitative trait locus (QTL), collectively regulates bayberry fruit color. In conclusion, our study presented a complete reference genome, uncovered a suite of allelic variations related to fruit-quality traits, and identified functional genes that could be harnessed to enhance fruit quality and breeding efficiency of bayberries.

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.

Conserving a threatened North American walnut: a chromosome-scale reference genome for butternut (Juglans cinerea)

With the advent of affordable and more accurate third-generation sequencing technologies, and the associated bioinformatic tools, it is now possible to sequence, assemble, and annotate more species of conservation concern than ever before. Juglans cinerea, commonly known as butternut or white walnut, is a member of the walnut family, native to the Eastern United States and Southeastern Canada. The species is currently listed as Endangered on the IUCN Red List due to decline from an invasive fungus known as Ophiognomonia clavigignenti-juglandacearum (Oc-j) that causes butternut canker. Oc-j creates visible sores on the trunks of the tree which essentially starves and slowly kills the tree. Natural resistance to this pathogen is rare. Conserving butternut is of utmost priority due to its critical ecosystem role and cultural significance. As part of an integrated undergraduate and graduate student training program in biodiversity and conservation genomics, the first reference genome for Juglans cinerea is described here. This chromosome-scale 539 Mb assembly was generated from over 100 × coverage of Oxford Nanopore long reads and scaffolded with the Juglans mandshurica genome. Scaffolding with a closely related species oriented and ordered the sequences in a manner more representative of the structure of the genome without altering the sequence. Comparisons with sequenced Juglandaceae revealed high levels of synteny and further supported J. cinerea's recent phylogenetic placement. Comparative assessment of gene family evolution revealed a significant number of contracting families, including several associated with biotic stress response.

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.

Transcriptomic Evidence of a Link between Cell Wall Biogenesis, Pathogenesis, and Vigor in Walnut Root and Trunk Diseases

Crown gall disease (Agrobacterium tumefaciens), crown/root rot disease (Phytophthora spp.), root lesion disease (Pratylenchus vulnus) and tree vigor are key traits affecting the productivity and quality of walnuts in California. Unchallenged hybrid rootstocks were analyzed by RNA-seq to examine pre-formed factors affecting these traits. Enrichment analysis of the differentially expressed genes revealed that the increased expression of cell wall biogenesis-related genes plays a key role in susceptibility to A. tumefaciens, susceptibility to Phytophthora spp. and increased vigor. Analysis of the predicted subcellular loci of the encoded proteins revealed that many gene products associated with vigor and susceptibility were targeted to the plasma membrane and extracellular space, connecting these traits to sustaining barrier function. We observed that RNA processing and splicing, along with predicted nuclear targeting, were associated with resistance to A. tumefaciens, resistance to Phytophthora spp. and low vigor. Four genes within the J. microcarpa QTL region for resistance to A. tumefaciens and Phytophthora spp. were represented among our transcripts, with two of the genes being differentially expressed in association with resistance to A. tumefaciens and decreased vigor. No differential expression related to Phytophthora spp. or P. vulnus resistance was observed in this region. Additionally, the J. microcarpa haplotype expressed more transcripts associated with resistance to A. tumefaciens, Phytophthora spp. and low vigor, but not P. vulnus, than the J. regia haplotype. We also report unique and shared hormone and defense responses associated with each trait. This research suggests a link between cell wall biogenesis, vigor and critical root diseases of walnut.

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 reference transcriptome for walnut anthracnose pathogen, Ophiognomonia leptostyla, guides the discovery of candidate virulence genes

Despite the economic losses due to the walnut anthracnose, Ophiognomonia leptostyla is an orphan fungus with respect to genomic resources. In the present study, the transcriptome of O. leptostyla was assembled for the first time. RNA sequencing was conducted for the fungal mycelia grown in a liquid media, and the inoculated leaf samples of walnut with the fungal conidia sampled at 48, 96 and 144 h post inoculation (hpi). The completeness, correctness, and contiguity of the de novo transcriptome assemblies generated with Trinity, Oases, SOAPdenovo-Trans and Bridger were compared to identify a single superior reference assembly. In most of the assessment criteria including N50, Transrate score, number of ORFs with known description in gene bank, the percentage of reads mapped back to the transcript (RMBT), BUSCO score, Swiss-Prot coverage bin and RESM-EVAL score, the Bridger assembly was the superior and thus used as a reference for profiling the O. leptostyla transcriptome in liquid media vs. during walnut infection. The k-means clustering of transcripts resulted in four distinct transcription patterns across the three sampling time points. Most of the detected CAZy transcripts had elevated transcription at 96 hpi that is hypothetically concurrent with the start of intracellular growth. The in-silico analysis revealed 103 candidate effectors of which six were members of Necrosis and Ethylene Inducing Like Protein (NLP) gene family belonging to three distinct k-means clusters. This study provided a complex and temporal pattern of the CAZys and candidate effectors transcription during six days post O. leptostyla inoculation on walnut leaves, introducing a list of candidate virulence genes for validation in future studies.

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.

LncRNA109897-JrCCR4-JrTLP1b forms a positive feedback loop to regulate walnut resistance against anthracnose caused by Colletotrichum gloeosporioides

Walnut anthracnose induced by Colletotrichum gloeosporioides is a disastrous disease that severely restricts the development of the walnut industry in China. Long non-coding RNAs (lncRNAs) are involved in adaptive responses to disease, but their roles in the regulation of walnut anthracnose resistance response are not well defined. In this study, transcriptome analysis demonstrated that a C. gloeosporioides-induced lncRNA, lncRNA109897, located upstream from the target gene JrCCR4, upregulated the expression of JrCCR4. JrCCR4 interacted with JrTLP1b and promoted its transcriptional activity. In turn, JrTLP1b induced the transcription of lncRNA109897 to promote its expression. Meanwhile, transient expression in walnut leaves and stable transformation of Arabidopsis thaliana further proved that lncRNA, JrCCR4, and JrTLP1b improve the resistance of C. gloeosporioides. Collectively, these findings provide insights into the mechanism by which the lncRNA109897-JrCCR4-JrTLP1b transcriptional cascade regulates the resistance of walnut to anthracnose.

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.

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.

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.

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.

Chromosome-level genome assembly provides new insights into Japanese chestnut (Castanea crenata) genomes

Japanese chestnut (Castanea crenata Sieb. et Zucc) is an economically and ecologically important chestnut species in East Asia. Here, we presented a high-quality chromosome-level reference genome of the Japanese chestnut cultivar 'Tsukuba' by combining Nanopore long reads and Hi-C sequencing. The final assembly has a size of 718.30 Mb and consists of 12 pseudochromosomes ranging from 41.03 to 92.03 Mb, with a BUSCO complete gene percentage of 97.6%. A total of 421.37 Mb repetitive sequences and 46,744 gene models encoding 46,463 proteins were predicted in the genome. Genome evolution analysis showed that Japanese chestnut is closely related to Chinese chestnut and these species shared a common ancestor ~6.5 million years ago. This high-quality Japanese chestnut genome represents an important resource for the chestnut genomics community and will improve our understanding of chestnut biology and evolution.

Isolating Linum usitatissimum L. Nuclear DNA Enabled Assembling High-Quality Genome

High-quality genome sequences help to elucidate the genetic basis of numerous biological processes and track species evolution. For flax (Linum usitatissimum L.)-a multifunctional crop, high-quality assemblies from Oxford Nanopore Technologies (ONT) data were unavailable, largely due to the difficulty of isolating pure high-molecular-weight DNA. This article proposes a scheme for gaining a contiguous L. usitatissimum assembly using Nanopore data. We developed a protocol for flax nuclei isolation with subsequent DNA extraction, which allows obtaining about 5 μg of pure high-molecular-weight DNA from 0.5 g of leaves. Such an amount of material can be collected even from a single plant and yields more than 30 Gb of ONT data in two MinION runs. We performed a comparative analysis of different genome assemblers and polishers on the gained data and obtained the final 447.1-Mb assembly of L. usitatissimum line 3896 genome using the Canu-Racon (two iterations)-Medaka combination. The genome comprised 1695 contigs and had an N50 of 6.2 Mb and a completeness of 93.8% of BUSCOs from eudicots_odb10. Our study highlights the impact of the chosen genome construction strategy on the resulting assembly parameters and its eligibility for future genomic studies.

Integrated metabolomic and transcriptomic analysis of the anthocyanin and proanthocyanidin regulatory networks in red walnut natural hybrid progeny leaves

Background

Walnuts are among the most important dry fruit crops worldwide, typically exhibiting green leaves and yellow-brown or gray-yellow seed coats. A specific walnut accession with red leaves and seed coats, 'RW-1', was selected for study because of its high anthocyanin and proanthocyanidin (PA) contents. Anthocyanins and PAs are important secondary metabolites and play key roles in plant responses to biotic and abiotic stresses. However, few studies have focused on the molecular mechanism of anthocyanin biosynthesis in walnuts.

Methods

In this study, we determined the anthocyanin and PA components and their contents in different color leaves of 'RW-1' natural hybrid progenies at various developmental stages. Integrated transcriptome and metabolome analyses were used to identify the differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs). We also performed conjoint analyses on DEGs and DAMs to ascertain the degree pathways, and explore the regulation of anthocyanin and PA biosynthesis.

Results

The results of widely targeted metabolome profiling and anthocyanin detection revealed 395 substances, including four PAs and 26 anthocyanins, in red (SR) and green leaves (SG) of 'RW-1' natural hybrid progenies. From the research, the contents of all anthocyanin components in SR were higher than that in SG. Among them, the contents of delphinidin 3-O-galactoside, cyanidin 3-O-galactoside, delphinidin 3-O-arabinoside and cyanidin 3-O-glucoside were significantly higher than others, and they were considered as the main types of anthocyanins. However, nine anthocyanins were detected only in SR. For PAs, the content of procyanidin C1 was higher in SR compared with SG, while procyanidin B1 and procyanidin B3 were higher in SR-1 and SR-3 but downregulated in SR-2 compared with the controls. Furthermore, transcriptome analysis revealed that the expressions of structural genes (C4H, F3H, F3'5'H, UFGT, LAR and ANR), three MYBs predicted as the activators of anthocyanin and PA biosynthesis, two MYBs predicted as the repressors of anthocyanin biosynthesis, and five WD40s in the anthocyanin and PA biosynthetic pathways were significantly higher in the SR walnuts. Gene-metabolite correlation analyses revealed a core set of 31 genes that were strongly correlated with four anthocyanins and one PA metabolites. The alteration of gene coding sequence altered the binding or regulation of regulatory factors to structural genes in different color leaves, resulting in the effective increase of anthocyanins and PAs accumulation in red walnut.

Conclusions

This study provides valuable information on anthocyanin and PA metabolites and candidate genes for anthocyanin and PA biosynthesis, yielding new insights into anthocyanin and PA biosynthesis in walnuts.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Transcriptome-wide identification of walnut PP2C family genes in response to external stimulus

Walnut is an important economic tree species while confronting with global environmental stress, resulting in decline in quality and yield. Therefore, it is urgent to elucidate the molecular mechanism for the regulation of walnut response to adversity. The protein phosphatase 2C (PP2C) gene family participates in cellular processes in eukaryotes through reversible phosphorylation of proteins and signal transduction regulation. However, the stress response function of PP2C genes was far to be clarified. Therefore, to understand the stress response mechanism of walnut tree, in this study, a total of 41 PP2C genes with complete ORFs were identified from Juglans regia, whose basic bio-information and expression patterns in response to multiple stresses and ABA were confirmed. The results showed that the ORFs of JrPP2Cs were 495 ~ 3231 bp in length, the predicted JrPP2C proteins contained 164 to 1076 amino acids and the molecular weights were 18,581.96 ~ 118,853.34 Da, the pI was 4.55 ~ 9.58. These JrPP2C genes were unevenly distributed on 14 chromosomes, among which Chr11 and Chr13 contained the most genes. Phylogenetic analysis found that these JrPP2C proteins were classed into 9 subfamilies, among which group F covered most JrPP2Cs. The JrPP2Cs in the same subfamily exhibited similarities in the composition of conserved domains, amino acid sequences of motifs and exon/intron organization in DNA sequences. Each JrPP2C includes 4 ~ 10 motifs and each motif contained 15 ~ 37 amino acids. Among the motifs, motif1, motif2, motif3 and motif8 were most abundant. Most of the JrPP2C genes diversely response to osmotic, cadmium, and Colletotrichum gloeosporioide stress as well as ABA treatments, among which JrPP2C28, JrPP2C17, JrPP2C09, JrPP2C36 were more obvious and deserves further attention. All these results indicated that JrPP2C genes play potential vital roles in plant response to multiple stimulus, and are possibly involved in ABA-dependent signaling pathway.

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.

Structural and functional leaf diversity lead to variability in photosynthetic capacity across a range of Juglans regia genotypes

Similar to other cropping systems, few walnut cultivars are used as scion in commercial production. Germplasm collections can be used to diversify cultivar options and hold potential for improving crop productivity, disease resistance and stress tolerance. In this study, we explored the anatomical and biochemical bases of photosynthetic capacity and response to water stress in 11 Juglans regia accessions in the U.S. department of agriculture, agricultural research service (USDA-ARS) National Clonal Germplasm. Net assimilation rate (An ) differed significantly among accessions and was greater in lower latitudes coincident with higher stomatal and mesophyll conductances, leaf thickness, mesophyll porosity, gas-phase diffusion, leaf nitrogen and lower leaf mass and stomatal density. High CO2 -saturated assimilation rates led to increases in An under diffusional and biochemical limitations. Greater An was found in lower-latitude accessions native to climates with more frost-free days, greater precipitation seasonality and lower temperature seasonality. As expected, water stress consistently impaired photosynthesis with the highest % reductions in lower-latitude accessions (A3, A5 and A9), which had the highest An under well-watered conditions. However, An for A3 and A5 remained among the highest under dehydration. J. regia accessions, which have leaf structural traits and biochemistry that enhance photosynthesis, could be used as commercial scions or breeding parents to enhance productivity.

A chromosome-scale genome assembly of the okapi (Okapia johnstoni)

The okapi (Okapia johnstoni), or forest giraffe, is the only species in its genus and the only extant sister group of the giraffe within the family Giraffidae. The species is one of the remaining large vertebrates surrounded by mystery because of its elusive behavior as well as the armed conflicts in the region where it occurs, making it difficult to study. Deforestation puts the okapi under constant anthropogenic pressure, and it is currently listed as “Endangered” on the IUCN Red List. Here, we present the first annotated de novo okapi genome assembly based on PacBio continuous long reads, polished with short reads, and anchored into chromosome-scale scaffolds using Hi-C proximity ligation sequencing. The final assembly (TBG_Okapi_asm_v1) has a length of 2.39 Gbp, of which 98% are represented by 28 scaffolds > 3.9 Mbp. The contig N50 of 61 Mbp and scaffold N50 of 102 Mbp, together with a BUSCO score of 94.7%, and 23 412 annotated genes, underline the high quality of the assembly. This chromosome-scale genome assembly is a valuable resource for future conservation of the species and comparative genomic studies among the giraffids and other ruminants.

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.

De Novo Arginine Synthesis Is Required for Full Virulence of Xanthomonas arboricola pv. juglandis During Walnut Bacterial Blight Disease

Walnut blight (WB) disease caused by Xanthomonas arboricola pv. juglandis (Xaj) threatens orchards worldwide. Nitrogen metabolism in this bacterial pathogen is dependent on arginine, a nitrogen-enriched amino acid that can either be synthesized or provided by the plant host. The arginine biosynthetic pathway uses argininosuccinate synthase (argG), associated with increased bacterial virulence. We examined the effects of bacterial arginine and nitrogen metabolism on the plant response during WB by proteomic analysis of the mutant strain Xaj argG^-. Phenotypically, the mutant strain produced 42% fewer symptoms and survived in the plant tissue with 2.5-fold reduced growth compared with wild type, while showing itself to be auxotrophic for arginine in vitro. Proteomic analysis of infected tissue enabled the profiling of 676 Xaj proteins and 3,296 walnut proteins using isobaric labeling in a data-dependent acquisition approach. Comparative analysis of differentially expressed proteins revealed distinct plant responses. Xaj wild type (WT) triggered processes of catabolism and oxidative stress in the host under observed disease symptoms, while most of the host biosynthetic processes triggered by Xaj WT were inhibited during Xaj argG^- infection. Overall, the Xaj proteins revealed a drastic shift in carbon and energy management induced by disruption of nitrogen metabolism while the top differentially expressed proteins included a Fis transcriptional regulator and a peptidyl-prolyl isomerase. Our results show the critical role of de novo arginine biosynthesis to sustain virulence and minimal growth during WB. This study is timely and critical as copper-based control methods are losing their effectiveness, and new sustainable methods are urgently needed in orchard environments.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

A chromosome-level reference genome of Ensete glaucum gives insight into diversity and chromosomal and repetitive sequence evolution in the Musaceae

BACKGROUND

Ensete glaucum (2n = 2x = 18) is a giant herbaceous monocotyledonous plant in the small Musaceae family along with banana (Musa). A high-quality reference genome sequence assembly of E. glaucum is a resource for functional and evolutionary studies of Ensete, Musaceae, and the Zingiberales.

FINDINGS

Using Oxford Nanopore Technologies, chromosome conformation capture (Hi-C), Illumina and RNA survey sequence, supported by molecular cytogenetics, we report a high-quality 481.5 Mb genome assembly with 9 pseudo-chromosomes and 36,836 genes. A total of 55% of the genome is composed of repetitive sequences with predominantly LTR-retroelements (37%) and DNA transposons (7%). The single 5S ribosomal DNA locus had an exceptionally long monomer length of 1,056 bp, more than twice that of the monomers at multiple loci in Musa. A tandemly repeated satellite (1.1% of the genome, with no similar sequence in Musa) was present around all centromeres, together with a few copies of a long interspersed nuclear element (LINE) retroelement. The assembly enabled us to characterize in detail the chromosomal rearrangements occurring between E. glaucum and the x = 11 species of Musa. One E. glaucum chromosome has the same gene content as Musa acuminata, while others show multiple, complex, but clearly defined evolutionary rearrangements in the change between x= 9 and 11.

CONCLUSIONS

The advance towards a Musaceae pangenome including E. glaucum, tolerant of extreme environments, makes a complete set of gene alleles, copy number variation, and a reference for structural variation available for crop breeding and understanding environmental responses. The chromosome-scale genome assembly shows the nature of chromosomal fusion and translocation events during speciation, and features of rapid repetitive DNA change in terms of copy number, sequence, and genomic location, critical to understanding its role in diversity and evolution.

Genomic Prediction of Complex Traits in Perennial Plants: A Case for Forest Trees

This chapter provides an overview of the genomic selection progress in long-lived forest tree species. Factors affecting the prediction accuracy in genomic prediction are assessed with examples from empirical studies. Infrastructure and resources required for the implementation of genomic selection are evaluated. Some general guidelines are provided for the successful application of genomic selection in forest tree breeding programs.

Walnut N-Acetylserotonin Methyltransferase Gene Family Genome-Wide Identification and Diverse Functions Characterization During Flower Bud Development

Melatonin widely mediates multiple developmental dynamics in plants as a vital growth stimulator, stress protector, and developmental regulator. N-acetylserotonin methyltransferase (ASMT) is the key enzyme that catalyzes the final step of melatonin biosynthesis in plants and plays an essential role in the plant melatonin regulatory network. Studies of ASMT have contributed to understanding the mechanism of melatonin biosynthesis in plants. However, AMST gene is currently uncharacterized in most plants. In this study, we characterized the JrASMT gene family using bioinformatics in a melatonin-rich plant, walnut. Phylogenetic, gene structure, conserved motifs, promoter elements, interacting proteins and miRNA analyses were also performed. The expansion and differentiation of the ASMT family occurred before the onset of the plant terrestrialization. ASMT genes were more differentiated in dicotyledonous plants. Forty-six ASMT genes were distributed in clusters on 10 chromosomes of walnut. Four JrASMT genes had homologous relationships both within walnut and between species. Cis-regulatory elements showed that JrASMT was mainly induced by light and hormones, and targeted cleavage of miRNA172 and miR399 may be an important pathway to suppress JrASMT expression. Transcriptome data showed that 13 JrASMT were differentially expressed at different periods of walnut bud development. WGCNA showed that JrASMT1/10/13/23 were coexpressed with genes regulating cell fate and epigenetic modifications during early physiological differentiation of walnut female flower buds. JrASMT12/28/37/40 were highly expressed during morphological differentiation of flower buds, associated with altered stress capacity of walnut flower buds, and predicted to be involved in the regulatory network of abscisic acid, salicylic acid, and cytokinin in walnut. The qRT-PCR validated the results of differential expression analysis and further provided three JrASMT genes with different expression profiles in walnut flower bud development. Our study explored the evolutionary relationships of the plant ASMT gene family and the functional characteristics of walnut JrASMT. It provides a valuable perspective for further understanding the complex melatonin mechanisms in plant developmental regulation.

Genomic evidence for homoploid hybrid speciation between ancestors of two different genera

Homoploid hybrid speciation (HHS) has been increasingly recognized as occurring widely during species diversification of both plants and animals. However, previous studies on HHS have mostly focused on closely-related species while it has been rarely reported or tested between ancestors of different genera. Here, we explore the likely HHS origin of Carpinus sect. Distegocarpus between sect. Carpinus and Ostrya in the family Betulaceae. We generate a chromosome-level reference genome for C. viminea of sect. Carpinus and re-sequence genomes of 44 individuals from the genera Carpinus and Ostrya. Our integrated analyses of all genomic data suggest that sect. Distegocarpus, which has three species, likely originates through HHS during the early divergence between Carpinus and Ostrya. Our study highlights the likelihood of an HHS event between ancestors of the extant genera during their initial divergences, which may have led to reticulate phylogenies at higher taxonomic levels.

Genome-Wide Identification of TLP Gene Family and Their Roles in Carya cathayensis Sarg in Response to Botryosphaeria dothidea

Hickory (Carya cathayensis) is a critical tree species of the genus Carya from the Juglandaceae family that contains nutrient-rich nuts. Due to large-scale soil degradation, the pests and diseases of hickory are becoming more and more serious. Thaumatin-like proteins (TLPs) are vital proteins involved in the complex defense process of plant pathogens. In this study, 40 CcTLP genes were identified genome-widely and phylogenetically grouped into three subfamilies. The sequence of CcTLPs had a conservative pattern, such as eight stable disulfide bonds, REDDD, and G-X-[GF]-X-C-X-T-[GA]-D-C-X(1,2)-G-X-(2,3)-C structure. In total, 57 cis-elements related to stress-responsive, light-responsive, phytohormone-responsive, and plant-responsive were discovered. Under salicylate (SA), methyl jasmonate (MeJA), and ethephon (ETH) treatments, the expressions of CcTLP28, CcTLP29, CcTLP30, CcTLP31, CcTLP32, CcTLP33, CcTLP37, CcTLP38, and CcTLP39 had different patterns. This is an indication that most of the TLP genes were upregulated by SA and downregulated by MeJA. Notably, seven TLP genes were significantly upregulated under the Botryosphaeria dothidea inoculation, especially CcTLP31, with an over 20-fold change. Nine genes were shown by subcellular localization analysis to be located at the plasma membrane and cytoplasm. The knowledge of the disease-resistant function of the CcTLP family in hickory is promoted by these results. A foundation reference for the molecular breeding of this plant in the future is provided by our findings.

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 Comparative Analysis of Genes and MicroRNAs Involved in the Floral Transition of the Xinjiang Early-Flowering Walnut (Juglans regia L.)

For tree crops, shortening the juvenile phase is a vital strategy to advance fruit bearing and enhance the breeding efficiency. Walnut (Juglans regia L.) seedlings usually take at least eight to 10 years to flower, but early-flowering (EF) types can flower one or two years after planting. In this study, RNA sequencing (RNA-Seq) and microRNA sequencing (miRNA-Seq) were used for a transcriptome-wide analysis of gene and miRNA expression in hybrids of the Xinjiang EF walnut variety ‘Xinwen 81’ and later-flowering (LF) walnut. Based on a high-quality chromosome-scale reference genome, a total of 3009 differentially expressed genes (DEGs) were identified, of which 933 were upregulated (accounting for 31%) and 2076 were downregulated (accounting for 69%). DEGs were functionally annotated, and the flowering-related genes FLOWERING LOCUS T (FT), SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), and LEAFY (LFY) showed remarkable upregulation in EF compared with in the LF walnut. In addition, miRNAs associated with floral transition were screened as candidates for flowering time regulation in the walnut. This work provides new insights into walnut floral transition, which may ultimately contribute to genetic improvement of the walnut.

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.

Dead-End Hybridization in Walnut Trees Revealed by Large-Scale Genomic Sequence Data

Although hybridization plays a large role in speciation, some unknown fraction of hybrid individuals never reproduces, instead remaining as genetic dead-ends. We investigated a morphologically distinct and culturally important Chinese walnut, Juglans hopeiensis, suspected to have arisen from hybridization of Persian walnut (J. regia) with Asian butternuts (J. cathayensis, J. mandshurica, and hybrids between J. cathayensis and J. mandshurica). Based on 151 whole-genome sequences of the relevant taxa, we discovered that all J. hopeiensis individuals are first-generation hybrids, with the time for the onset of gene flow estimated as 370,000 years, implying both strong postzygotic barriers and the presence of J. regia in China by that time. Six inversion regions enriched for genes associated with pollen germination and pollen tube growth may be involved in the postzygotic barriers that prevent sexual reproduction in the hybrids. Despite its long-recurrent origination and distinct traits, J. hopeiensis does not appear on the way to speciation.