The Manchurian Walnut Genome: Insights into Juglone and Lipid Biosynthesis

Arbuscular mycorrhizal symbiosis with Rhizophagus irregularis DAOM197198 modifies the root transcriptome of walnut trees

Walnut trees are cultivated and exploited worldwide for commercial timber and nut production. They are heterografted plants, with the rootstock selected to grow in different soil types and conditions and to provide the best anchorage, vigor, and resistance or tolerance to soil borne pests and diseases. However, no individual rootstock is tolerant of all factors that impact walnut production. In Europe, Juglans regia is mainly used as a rootstock. Like most terrestrial plants, walnut trees form arbuscular mycorrhizal symbioses, improving water and nutrient uptake and providing additional ecosystem services. Effects of arbuscular mycorrhizal symbiosis on root gene regulation, however, has never been assessed. We analyzed the response of one rootstock of J. regia to colonization by the arbuscular mycorrhizal fungus Rhizophagus irregularis DAOM197198. Plant growth as well as the nitrogen and phosphorus concentrations in roots and shoots were significantly increased in mycorrhizal plants versus non-colonized plants. In addition, we have shown that 1,549 genes were differentially expressed, with 832 and 717 genes up- and down-regulated, respectively. The analysis also revealed that some rootstock genes involved in plant nutrition through the mycorrhizal pathway, are regulated similarly as in other mycorrhizal woody species: Vitis vinifera and Populus trichocarpa. In addition, an enrichment analysis performed on GO and KEGG pathways revealed some regulation specific to J. regia (i.e., the juglone pathway). This analysis reinforces the role of arbuscular mycorrhizal symbiosis on root gene regulation and on the need to finely study the effects of diverse arbuscular mycorrhizal fungi on root gene regulation, but also of the scion on the functioning of an arbuscular mycorrhizal fungus in heterografted plants such as walnut tree.

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.

Integrated metabolomic and transcriptomic dynamic profiles of endopleura coloration during fruit maturation in three walnut cultivars

BACKGROUND

The color of endopleura is a vital factor in determining the economic value and aesthetics appeal of nut. Walnuts (Juglans) are a key source of edible nuts, high in proteins, amino acids, lipids, carbohydrates. Walnut had a variety endopleura color as yellow, red, and purple. However, the regulation of walnut endopleura color remains little known.

RESULTS

To understand the process of coloration in endopleura, we performed the integrative analysis of transcriptomes and metabolomes at two developmental stages of walnut endopleura. We obtained total of 4,950 differentially expressed genes (DEGs) and 794 metabolites from walnut endopleura, which are involved in flavonoid and phenolic biosynthesis pathways. The enrichment analysis revealed that the cinnamic acid, coniferyl alcohol, naringenin, and naringenin-7-O-glucoside were important metabolites in the development process of walnut endopleura. Transcriptome and metabolome analyses revealed that the DEGs and differentially regulated metabolites (DRMs) were significantly enriched in flavonoid biosynthesis and phenolic metabolic pathways. Through co-expression analysis, CHS (chalcone synthase), CHI (chalcone isomerase), CCR (cinnamoyl CoA reductase), CAD (cinnamyl alcohol dehydrogenase), COMT (catechol-Omethyl transferase), and 4CL (4-coumaroyl: CoA-ligase) may be the key genes that potentially regulate walnut endopleura color in flavonoid biosynthesis and phenolic metabolic pathways.

CONCLUSIONS

This study illuminates the metabolic pathways and candidate genes that underlie the endopleura coloration in walnuts, lay the foundation for further study and provides insights into controlling nut's colour.

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.

Progress and prospect: Biosynthesis of plant natural products based on plant chassis

Plant-derived natural products are a specific class of active substances with numerous applications in the medical, energy, and industrial fields. Many of these substances are in high demand and have become the fundamental materials for various purposes. Recently, the use of synthetic biology to produce plant-derived natural products has become a significant trend. Plant chassis, in particular, offer unique advantages over microbial chassis in terms of cell structure, product affinity, safety, and storage. The development of the plant hairy root tissue culture system has accelerated the commercialization and industrialization of synthetic biology in the production of plant-derived natural products. This paper will present recent progress in the synthesis of various plant natural products using plant chassis, organized by the types of different structures. Additionally, we will summarize the four primary types of plant chassis used for synthesizing natural products from plant sources and review the enabling technologies that have contributed to the development of synthetic biology in recent years. Finally, we will present the role of isolated and combined use of different optimization strategies in breaking the upper limit of natural product production in plant chassis. This review aims to provide practical references for synthetic biologists and highlight the great commercial potential of plant chassis biosynthesis, such as hairy roots.

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.

Widely Targeted Metabolomic Profiling Combined with Transcriptome Analysis Provides New Insights into Lipid Biosynthesis in Seed Kernels of Pinus koraiensis

Lipid-rich Pinus koraiensis seed kernels are highly regarded for their nutritional and health benefits. To ascertain the molecular mechanism of lipid synthesis, we conducted widely targeted metabolomic profiling together with a transcriptome analysis of the kernels in P. koraiensis cones at various developmental stages. The findings reveal that 148 different types of lipid metabolites, or 29.6% of total metabolites, are present in kernels. Among those metabolites, the concentrations of linoleic acid, palmitic acid, and α-linolenic acid were higher, and they steadily rose as the kernels developed. An additional 10 hub genes implicated in kernel lipid synthesis were discovered using weighted gene co-expression network analysis (WGCNA), gene interaction network analysis, oil body biosynthesis, and transcriptome analysis. This study used lipid metabolome and transcriptome analyses to investigate the mechanisms of key regulatory genes and lipid synthesis molecules during kernel development, which served as a solid foundation for future research on lipid metabolism and the creation of P. koraiensis kernel food.

Genomic insights into biased allele loss and increased gene numbers after genome duplication in autotetraploid Cyclocarya paliurus

BACKGROUND

Autopolyploidy is a valuable model for studying whole-genome duplication (WGD) without hybridization, yet little is known about the genomic structural and functional changes that occur in autopolyploids after WGD. Cyclocarya paliurus (Juglandaceae) is a natural diploid-autotetraploid species. We generated an allele-aware autotetraploid genome, a chimeric chromosome-level diploid genome, and whole-genome resequencing data for 106 autotetraploid individuals at an average depth of 60 × per individual, along with 12 diploid individuals at an average depth of 90 × per individual.

RESULTS

Autotetraploid C. paliurus had 64 chromosomes clustered into 16 homologous groups, and the majority of homologous chromosomes demonstrated similar chromosome length, gene numbers, and expression. The regions of synteny, structural variation and nonalignment to the diploid genome accounted for 81.3%, 8.8% and 9.9% of the autotetraploid genome, respectively. Our analyses identified 20,626 genes (69.18%) with four alleles and 9191 genes (30.82%) with one, two, or three alleles, suggesting post-polyploid allelic loss. Genes with allelic loss were found to occur more often in proximity to or within structural variations and exhibited a marked overlap with transposable elements. Additionally, such genes showed a reduced tendency to interact with other genes. We also found 102 genes with more than four copies in the autotetraploid genome, and their expression levels were significantly higher than their diploid counterparts. These genes were enriched in enzymes involved in stress response and plant defense, potentially contributing to the evolutionary success of autotetraploids. Our population genomic analyses suggested a single origin of autotetraploids and recent divergence (~ 0.57 Mya) from diploids, with minimal interploidy admixture.

CONCLUSIONS

Our results indicate the potential for genomic and functional reorganization, which may contribute to evolutionary success in autotetraploid C. paliurus.

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 Identification of AP2/ERF Superfamily Genes in Juglans mandshurica and Expression Analysis under Cold Stress

Juglans mandshurica has strong freezing resistance, surviving temperatures as low as -40 °C, making it an important freeze tolerant germplasm resource of the genus Juglans. APETALA2/ethylene responsive factor (AP2/ERF) is a plant-specific superfamily of transcription factors that regulates plant development, growth, and the response to biotic and abiotic stress. In this study, phylogenetic analysis was used to identify 184 AP2/ERF genes in the J. mandshurica genome, which were classified into five subfamilies (JmAP2, JmRAV, JmSoloist, JmDREB, and JmERF). A significant amount of discordance was observed in the 184 AP2/ERF genes distribution of J. mandshurica throughout its 16 chromosomes. Duplication was found in 14 tandem and 122 segmental gene pairs, which indicated that duplications may be the main reason for JmAP2/ERF family expansion. Gene structural analysis revealed that 64 JmAP2/ERF genes contained introns. Gene evolution analysis among Juglandaceae revealed that J. mandshurica is separated by 14.23 and 15 Mya from Juglans regia and Carya cathayensis, respectively. Based on promoter analysis in J. mandshurica, many cis-acting elements were discovered that are related to light, hormones, tissues, and stress response processes. Proteins that may contribute to cold resistance were selected for further analysis and were used to construct a cold regulatory network based on GO annotation and JmAP2/ERF protein interaction network analysis. Expression profiling using qRT-PCR showed that 14 JmAP2/ERF genes were involved in cold resistance, and that seven and five genes were significantly upregulated under cold stress in female flower buds and phloem tissues, respectively. This study provides new light on the role of the JmAP2/ERF gene in cold stress response, paving the way for further functional validation of JmAP2/ERF TFs and their application in the genetic improvement of Juglans and other tree species.

Comprehensive Identification and Analyses of the GRF Gene Family in the Whole-Genome of Four Juglandaceae Species

The GRF gene family plays an important role in plant growth and development as regulators involved in plant hormone signaling and metabolism. However, the Juglandaceae GRF gene family remains to be studied. Here, we identified 15, 15, 19, and 20 GRF genes in J. regia, C. illinoinensis, J. sigillata, and J. mandshurica, respectively. The phylogeny shows that the Juglandaceae family GRF is divided into two subfamilies, the ε-group and the non-ε-group, and that selection pressure analysis did not detect amino acid loci subject to positive selection pressure. In addition, we found that the duplications of the Juglandaceae family GRF genes were all segmental duplication events, and a total of 79 orthologous gene pairs and one paralogous homologous gene pair were identified in four Juglandaceae families. The Ka/KS ratios between these homologous gene pairs were further analyzed, and the Ka/KS values were all less than 1, indicating that purifying selection plays an important role in the evolution of the Juglandaceae family GRF genes. The codon bias of genes in the GRF family of Juglandaceae species is weak, and is affected by both natural selection pressure and base mutation, and translation selection plays a dominant role in the mutation pressure in codon usage. Finally, expression analysis showed that GRF genes play important roles in pecan embryo development and walnut male and female flower bud development, but with different expression patterns. In conclusion, this study will serve as a rich genetic resource for exploring the molecular mechanisms of flower bud differentiation and embryo development in Juglandaceae. In addition, this is the first study to report the GRF gene family in the Juglandaceae family; therefore, our study will provide guidance for future comparative and functional genomic studies of the GRF gene family in the Juglandaceae specie.

Integrated metabolomic and transcriptomic analyses reveal different metabolite biosynthesis profiles of Juglans mandshurica in shade

Light is not only a very important source of energy for the normal growth and development of plants, but also a regulator of many development and metabolic processes. The mechanism of plant growth and development under low light conditions is an important scientific question. With the promulgation of the law to stop natural forest cutting, understory regeneration is an important method for artificial forest afforestation. Here, the growth and physiological indexes of Juglans mandshurica, an important hardwood species in Northeast China, were measured under different shade treatments. In addition, transcriptome and metabolome were compared to analyze the molecular mechanism of shade tolerance in J. mandshurica. The results showed that the seedling height of the shade treatment group was significantly higher than that of the control group, and the 50% light (L50) treatment was the highest. Compared with the control group, the contents of gibberellin, abscisic acid, brassinolide, chlorophyll a, and chlorophyll b in all shade treatments were significantly higher. However, the net photosynthetic rate and water use efficiency decreased with increasing shade. Furthermore, the transcriptome identified thousands of differentially expressed genes in three samples. Using enrichment analysis, we found that most of the differentially expressed genes were enriched in photosynthesis, plant hormone signal transduction and chlorophyll synthesis pathways, and the expression levels of many genes encoding transcription factors were also changed. In addition, analysis of differentially accumulated metabolites showed that a total of 470 differential metabolites were identified, and flavonoids were the major differential metabolites of J. mandshurica under light stress. These results improved our understanding of the molecular mechanism and metabolite accumulation under light stress in J. mandshurica.

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.