The genomes of pecan and Chinese hickory provide insights into Carya evolution and nut nutrition

The Integrated Analysis of miRNome and Degradome Sequencing Reveals the Regulatory Mechanisms of Seed Development and Oil Biosynthesis in Pecan (Carya illinoinensis)

Pecan seed oil is a valuable source of essential fatty acids and various bioactive compounds; however, the functions of microRNAs and their targets in oil biosynthesis during seed development are still unknown. Here, we found that the oil content increased rapidly in the three early stages in three cultivars, and that oleic acid was the predominant fatty acid component in the mature pecan embryos. We identified, analyzed, and validated the expression levels of miRNAs related to seed development and oil biosynthesis, as well as their potential target genes, using small RNA sequencing data from three stages (120, 135, and 150 days after flowering). During the seed development process, 365 known and 321 novel miRNAs were discovered. In total, 91 known and 181 novel miRNAs were found to be differentially expressed, and 633 target genes were further investigated. The expression trend analysis revealed that the 91 known miRNAs were classified into eight groups, approximately two-thirds of which were up-regulated, whereas most novel miRNAs were down-regulated. The qRT-PCR and degradome sequencing data were used to identify five miRNA- target pairs. Overall, our study provides valuable insights into the molecular regulation of oil biosynthesis in pecan seeds.

Two Carya Species, Carya hunanensis and Carya illinoinensis, Used as Rootstocks Point to Improvements in the Heat Resistance of Carya cathayensis

Grafting as a crucial horticultural technique has been widely used in the cultivation of Carya cathayensis (Chinese hickory), which is a unique and important economic tree in the northeast of Zhejiang Province and the south of Anhui Province. However, the existing literature lacks research on the potential impact of various rootstocks on the thermal tolerance of Chinese hickory. The objectives of this study were to evaluate heat tolerance in four distinct groups of Chinese hickory, including C. cathayensis grafted onto Carya hunanensis and Carya illinoinensis, one self-grafted group (C. cathayensis grafted onto C. cathayensis), and one non-grafted group (C. cathayensis). We examined photosynthesis parameters, phytohormones, and differentially expressed genes in the four various hickory groups subjected to 25 °C, 35 °C, and 40 °C heat stress (HS). The results demonstrated that grafting onto C. hunanensis and C. illinoinensis exhibited a higher net photosynthetic rate and stomatal conductance, lower intercellular CO2 concentration, and smaller changes in plant hormone content compared to self-grafted and non-grafted group under HS. The transcriptome results revealed that the majority of differentially expressed genes (DEGs) associated with photosynthetic pathways exhibited downregulation under HS, while the degree of variation in grafted groups using C. hunanensis and C. illinoinensis as rootstocks was comparatively lower than that observed in self-grafted and non-grafted groups. The alteration in the expression patterns of DEGs involved in plant hormone synthesis and metabolism under HS corresponded to changes in plant hormone contents. Overall, Chinese hickory grafted onto C. hunanensis and C. illinoinensis exhibited enhanced resistance to high-temperature stress at the juvenile stage.

Integrated transcriptomic and lipidomic analysis provides key insights into lipid content changes during pecan (Carya illinoensis) fruit development

Pecans [Carya illinoinensis (Wangenh.) K. Koch] are highly valued for their abundance of quality healthy lipids, positively impacting human health and making themselves a preferred choice for nutritionally rich foods. However, a comprehensive understanding of the high-resolution characteristics of pecan fruit lipid composition and its dynamic changes, as well as the transfer between embryo and pericarp during development, remains incomplete. In this study, through integrated multi-omics analysis, we observed significant spatiotemporal heterogeneity in lipid changes between the pericarp and embryo. It showed smaller fluctuations and more stable lipid levels in the pericarp while exhibiting a dynamic pattern of initially increasing and then decreasing lipid content in the embryo. In this study, a total of 52 differentially expressed genes were identified, related to fatty acid synthesis and metabolism pathways in the two tissues, with changes in oleic acid and linoleic acid composition being the primary features of the embryo. This research lays the foundation for further understanding the differential regulation mechanisms of lipid metabolism between embryo and pericarp. Overall, this study filled the knowledge gap regarding dynamic changes in pericarp lipid metabolites, provided crucial insights into the lipid metabolism network during pecan fruit development, and established a scientific basis for the genetic improvement of pecan crops.

Polyphenols from hickory nut reduce the occurrence of atherosclerosis in mice by improving intestinal microbiota and inhibiting trimethylamine N-oxide production

BACKGROUND

Trimethylamine N-oxide (TMAO), a metabolite produced by intestinal microbiota through metabolizing phosphatidylcholine, choline, l-carnitine and betaine in the diet, has been implicated in the pathogenesis of atherosclerosis (AS). Concurrently, dietary polyphenols have garnered attention for their potential to ameliorate obesity, diabetes and atherosclerosis primarily by modulating the intestinal microbial structure. Hickory (Carya cathayensis) nut, a polyphenol-rich food product favored for its palatability, emerges as a candidate for exploration.

HYPOTHESIS/PURPOSE

The relationship between polyphenol of hickory nut and atherosclerosis prevention will be firstly clarified, providing theoretical basis for the discovery of natural products counteracting TMAO-induced AS process in hickory nut.

STUDY DESIGN AND METHODS

Employing Enzyme-linked Immunosorbent Assay (ELISA) and histological examination of aortic samples, the effects of total polyphenol extract on obesity index, inflammatory index and pathological changes of atherosclerosis in C57BL/6 J mice fed with high-fat and high choline diet were evaluated. Further, the composition, abundance, and function of mouse gut microbiota were analyzed through 16srDNA sequencing. Concurrently, the levels of TMAO and the expression of key enzymes (CutC and FMO3) involved in its synthesis are quantified using ELISA, Western Blot and Real-Time Quantitative PCR (RT-qPCR). Additionally, targeted metabolomic profiling of the hickory nut polyphenol extract was conducted, accompanied by molecular docking simulations to predict interactions between candidate polyphenols and the CutC/FMO3 using Autodock Vina. Finally, the docking prediction were verified by microscale thermophoresis (MST) .

RESULTS

Polyphenol extracts of hickory nut improved the index of obesity and inflammation, and alleviated the pathological changes of atherosclerosis in C57BL/6 J mice fed with high-fat and high-choline diet. Meanwhile, these polyphenol extracts also changed the composition and function of intestinal microbiota, and increased the abundance of microorganisms in mice. Notably, the abundance of intestinal microbiota endowed with CutC gene was significantly reduced, coherent with expression of CutC catalyzing TMA production. Moreover, polyphenol extracts also decreased the expression of FMO3 in the liver, contributing to the reduction of TMAO levels in serum. Furthermore, metabonomic profile analysis of these polyphenol extracts identified 647 kinds of polyphenols. Molecular docking predication further demonstrated that Casuariin and Cinnamtannin B2 had the most potential inhibition on the enzymatic activities of CutC or FMO3, respectively. Notably, MST analysis corroborated the potential for direct interaction between CutC enzyme and available polyphenols such as Corilagin, (-)-Gallocatechin gallate and Epigallocatechin gallate.

CONCLUSION

Hickory polyphenol extract can mitigate HFD-induced AS by regulating intestinal microflora in murine models. In addition, TMA-FMO3-TMAO pathway may play a key role in this process. This research unveils, for the inaugural time, the complex interaction between hickory nut-derived polyphenols and gut microbial, providing novel insights into the role of dietary polyphenols in AS prevention.

The Comparative Genomics of Botryosphaeriaceae Suggests Gene Families of Botryosphaeria dothidea Related to Pathogenicity on Chinese Hickory Tree

Trunk canker poses a major threat to the production of Chinese hickory tree (Carya cathayensis Sarg.), which is primarily determined by Botryosphaeriaceae. In our previous work, we identified Botryosphaeria dothidea as the predominant pathogen of this disease. However, it is still unclear about corresponding gene families and mechanisms associated with B. dothidea's pathogenicity on Chinese hickory tree. Here, we present a comparative analysis of high-quality genome assemblies of Botryosphaeria dothidea and other isolated pathogens, showing highly syntenic relationships between B. dothidea and its closely related species and the conservative evolution of the Botryosphaeriaceae family. Higher GC contents were found in the genomes of B. dothidea and three other isolated pathogens (Botryshaeria cortices, Botryshaeria fabicerciana, and Botryshaeria qingyuanensis) compared to Macrophomina phaseolina, Neofusicoccum parvum, Diplodia corticola, and Lasiodiplodia theobromae. An investigation of genes specific to or expanded in B. dothidea revealed that one secreted glucanase, one orsellinic acid biosynthesis enzyme, and two MFS transporters positively regulated B. dothidea's pathogenicity. We also observed an overrepresentation of viral integrase like gene and heterokaryon incompatibility proteins in the B. dothidea's genome. In addition, we observed one LRR-domain-containing protein and two Sec-domain-containing proteins (Sec_1 and Sec_7) that underwent positive selection. This study will help to understand B. dothidea's pathogenicity and potential influence on the infection of Chinese hickory, which will help in the development of disease control and ensure the security of Chinese hickory production.

Identification of ARF genes in Juglans Sigillata Dode and analysis of their expression patterns under drought stress

BACKGROUND

Auxin response factor (ARF), a transcription factors that controls the expression of genes responsive to auxin, plays a key role in the regulation of plant growth and development. Analyses aimed at identifying ARF family genes and characterizing their functions in Juglans sigillata Dode are lacking.

METHODS AND RESULTS

We used bioinformatic approaches to identify members of the J. sigillata ARF gene family and analyze their evolutionary relationships, collinearity, cis-acting elements, and tissue-specific expression patterns. The expression patterns of ARF gene family members under natural drought conditions were also analyzed. The J. sigillata ARF gene family contained 31 members, which were unevenly distributed across 16 chromosomes. We constructed a phylogenetic tree of JsARF genes and other plant ARF genes. Cis-acting elements in the promoters of JsARF were predicted. JsARF28 showed higher expressions in both the roots and leaves. A heat map of the transcriptome data of the cluster analysis under drought stress indicated that JsARF3/9/11/17/20/26 are responsive to drought. The expression of the 11 ARF genes varied under PEG treatment and JsARF18 and JsARF20 were significantly up-regulated.

CONCLUSIONS

The interactions between abiotic stresses and plant hormones are supported by our cumulative data, which also offers a theoretical groundwork for comprehending the ARF mechanism and drought resistance in J. sigillata.

Pecan secondary metabolites influenced the population of Zeuzera coffeae by affecting the structure and function of the larval gut microbiota

Background

The plant secondary metabolites (PSMs), as important plant resistance indicators, are important targets for screening plant insect resistance breeding. In this study, we aimed to investigate whether the population of Zeuzera coffeae (ZC) is affected by different varieties of Carya illinoinensis PSMs content. At the same time, the structure and function of the gut microbiome of ZC were also analyzed in relation to different pecan varieties.

Methods

We counted the populations of ZC larvae in four pecan varieties and determined the content of four types of PSMs. The structure and function of the larval gut microbiota were studied in connection to the number of larvae and the content of PSMs. The relationships were investigated between larval number, larval gut microbiota, and PSM content.

Results

We found that the tannins, total phenolics, and total saponins of 4 various pecans PSMs stifled the development of the ZC larval population. The PSMs can significantly affect the diversity and abundance of the larval gut microbiota. Enrichment of ASV46 (Pararhizobium sp.), ASV994 (Olivibacter sp.), ASV743 (Rhizobium sp.), ASV709 (Rhizobium sp.), ASV671 (Luteolibacter sp.), ASV599 (Agrobacterium sp.), ASV575 (Microbacterium sp.), and ASV27 (Rhizobium sp.) in the gut of larvae fed on high-resistance cultivars was positively associated with their tannin, total saponin, and total phenolic content. The results of the gut microbiome functional prediction for larvae fed highly resistant pecan varieties showed that the enriched pathways in the gut were related to the breakdown of hazardous chemicals.

Conclusion

Our findings provide further evidence that pecan PSMs influence the structure and function of the gut microbiota, which in turn affects the population stability of ZC. The study's findings can serve as a theoretical foundation for further work on selecting ZC-resistant cultivars and developing green management technology for ZC.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BACKGROUND

Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.

RESULTS

We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.

CONCLUSION

Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

Characteristic analysis of BZR genes family and their responses to hormone treatments and abiotic stresses in Carya illinoinensis

As the core of Brassinosteroids (BR) signaling pathway, BR-resistant (BZR) transcription factor regulates thousands of targeted genes mediating photomophogenesis, pollen sterility, cell expansion and stress response. Pecan (Carya illinoinensis) is a famous trees species of Carya, and its nut has high nutritional and economic values. However, there has no report on BZR genes family in pecan yet. Herein, totals of seven CiBZR members were identified in pecan genome, which were predicted to be hydrophilic unstable proteins and located in the nucleus. CiBZR genes had close evolutionary relationships with CcBZRs and JrBZRs in both Carya cathayensis and Juglans regia. These seven CiBZR genes were located independently on 7 chromosomes without doubling or tandem duplication. Based on the analysis of conserved motifs and gene structures, CiBZR genes were divided into three categories. More than 40 cis-acting elements were found in the 2 kb promoter regions of CiBZRs, which were mainly involved in hormone, light, and stress response, and plant growth and development. Notably, some of these CiBZR proteins were mainly located in the nucleus, had the self-activation ability and interaction relationship with BIN2 kinase, and negatively regulated the expression of CiCPD and CiDWF4. Gene expressions analysis further showed that CiBZR genes could express in many tissues and shared similar expression trends during embryo development. Moreover, most CiBZR genes responded to BR, Gibberellin (GA), Strigolactone (SL), salt, acid and osmotic stress. This study provides theoretical basis for the subsequent study on the role of CiBZR family genes in plant growth, development and stress responses.

Genome-wide association study of growth curve parameters reveals novel genomic regions and candidate genes associated with metatarsal bone traits in chickens

The growth and development of chicken bones have an enormous impact on the health and production performance of chickens. However, the development pattern and genetic regulation of the chicken skeleton are poorly understood. This study aimed to evaluate metatarsal bone growth and development patterns in chickens via non-linear models, and to identify the genetic determinants of metatarsal bone traits using a genome-wide association study (GWAS) based on growth curve parameters. Data on metatarsal length (MeL) and metatarsal circumference (MeC) were obtained from 471 F2 chickens (generated by crossing broiler sires, derived from a line selected for high abdominal fat, with Baier layer dams) at 4, 6, 8, 10, and 12 weeks of age. Four non-linear models (Gompertz, Logistic, von Bertalanffy, and Brody) were used to fit the MeL and MeC growth curves. Subsequently, the estimated growth curve parameters of the mature MeL or MeC (A), time-scale parameter (b), and maturity rate (K) from the non-linear models were utilized as substitutes for the original bone data in GWAS. The Logistic and Brody models displayed the best goodness-of-fit for MeL and MeC, respectively. Single-trait and multi-trait GWASs based on the growth curve parameters of the Logistic and Brody models revealed 4 618 significant single nucleotide polymorphisms (SNPs), annotated to 332 genes, associated with metatarsal bone traits. The majority of these significant SNPs were located on Gallus gallus chromosome (GGA) 1 (167.433-176.318 Mb), GGA2 (96.791-103.543 Mb), GGA4 (65.003-83.104 Mb) and GGA6 (64.685-95.285 Mb). Notably, we identified 12 novel GWAS loci associated with chicken metatarsal bone traits, encompassing 35 candidate genes. In summary, the combination of single-trait and multi-trait GWASs based on growth curve parameters uncovered numerous genomic regions and candidate genes associated with chicken bone traits. The findings benefit an in-depth understanding of the genetic architecture underlying metatarsal growth and development in chickens.

Distinct ancient structural polymorphisms control heterodichogamy in walnuts and hickories

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

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

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

Analysis of Delta(9) fatty acid desaturase gene family and their role in oleic acid accumulation in Carya cathayensis kernel

Carya cathayensis, commonly referred to as Chinese hickory, produces nuts that contain high-quality edible oils, particularly oleic acid (18:1). It is known that stearoyl-ACP desaturase (SAD) is the first key step converting stearic acid (C18:0, SA) to oleic acid (C18:1, OA) in the aminolevulinic acid (ALA) biosynthetic pathway and play an important role in OA accumulation. Thus far, there is little information about SAD gene family in C. cathayensis and the role of individual members in OA accumulation. This study searched the Chinese Hickory Genome Database and identified five members of SAD genes, designated as CcSADs, at the whole genome level through the comparison with the homologous genes from Arabidopsis. RNA-Seq analysis showed that CcSSI2-1, CcSSI2-2, and CcSAD6 were highly expressed in kernels. The expression pattern of CcSADs was significantly correlated with fatty acid accumulation during the kernel development. In addition, five full-length cDNAs encoding SADs were isolated from the developing kernel of C. cathayensis. CcSADs-green fluorescent protein (GFP) fusion construct was infiltrated into tobacco epidermal cells, and results indicated their chloroplast localization. The catalytic function of these CcSADs was further analyzed by heterologous expression in Saccharomyces cerevisiae, Nicotiana benthamiana, and walnut. Functional analysis demonstrated that all CcSADs had fatty acid desaturase activity to catalyze oleic acid biosynthesis. Some members of CcSADs also have strong substrate specificity for 16:0-ACP to synthesize palmitoleic acid (C16:1, PA). Our study documented SAD gene family in C. cathayensis and the role of CcSSI2-1, CcSSI2-2, and CcSAD6 in OA accumulation, which could be important for future improvement of OA content in this species via genetic manipulation.

Lipidomic and comparative transcriptomic analysis of fatty acid synthesis pathway in Carya illinoinensis embryo

Pecan (Carya illinoinensis (Wagenh.) K. Koch) is an important oilseed nut and is rich in fatty acids (FAs) and flavonols. Pecan FA has significantly edible, industrial and clinical value. To investigate the dynamic patterns and compositions of FA, and the molecular mechanism that controls FA accumulation in pecan, lipidomic and transcriptomic analyses were performed to determine lipid profiles and gene expression in pecan's FA biosynthesis pathway. In the present study, compared with cultivars 'Caddo' and 'Y-01', 'Mahan' formed larger and heavier embryos and accumulated higher oil content. Lipidomic analysis showed that FA and (O-acyl)-1-hydroxy FA contents were higher in 'Mahan' at the mature stage. Based on full-length and comparative RNA-Seq, differential expression gene enrichment analysis revealed that many functional genes participated in the pathways of 'fatty acid biosynthesis', 'fatty acid metabolism' and 'linoleic acid metabolism'. High FA accumulation model from 'Mahan' demonstrated that key enzyme-encoding genes played an important role in regulating FA biosynthesis. Co-expression module analysis indicated that several transcription factors (TFs; MYB, TCP, bHLH, Dof, ERF, NAC) were involved in FA accumulation by regulating the expression of functional genes, and real-time quantitative PCR verification proved that these TFs had a high correlation with the pecan FA accumulation pattern. These findings provided an insight into the molecular mechanism of FA accumulation in C. illinoinensis embryo, which contributes to pecan oil yielding and pecan molecular breeding.

Whole-genome identification and multiple abiotic stresses expression pattern profiling analysis of PLATZ transcription factor family members in Pecan (Carya illinoensis)

Plant AT-rich sequence and zinc-binding (PLATZ), as a plant-specific transcription factor, have been identified and studied in a variety of plants. However, there are no reports about PLATZ proteins in Carya illinoensis (pecan). Here, 24 C. illinoensis CiPLATZs have been identified and divided into 4 groups. Gene structure, motif composition, conserved domain and cis-acting elements analysis indicated that the PLATZ gene family was highly conserved. Transcriptome data combined with qRT-PCR analysis revealed that CiPLATZ6, CiPLATZ12, CiPLATZ13, CiPLATZ14 and CiPLATZ23 were highly expressed in multiple tissues of C. illinoensis and strongly responded to drought, salt and heat stress. Among them, CiPLATZ6, CiPLATZ12 and CiPLATZ23 were all located in the nucleus and had no transcriptional autoactivation ability in yeast cells, and acted as transcriptional suppressors in plants. In addition, the CiPLATZ23-overexpressing transgenic Arabidopsis thaliana showed enhanced tolerance to drought. Measurements of physiological indicators and analysis of stress-related genes expression levels in transgenic A. thaliana were used to support this conclusion. The results of this study are helpful to understand the structural feature and function of CiPLATZs, and provide candidate genes for molecular breeding of drought tolerance of C. illinoensis.

Exploring the Core Bacteria and Functional Traits in Pecan (Carya illinoinensis) Rhizosphere

Pecan (Carya illinoinensis) and Chinese hickory (Carya cathayensis) are important commercially cultivated nut trees. They are phylogenetically closely related plants; however, they exhibit significantly different phenotypes in response to abiotic stress and development. The rhizosphere selects core microorganisms from bulk soil, playing a pivotal role in the plant's resistance to abiotic stress and growth. In this study, we used metagenomic sequencing to compare the selection capabilities of seedling pecan and seedling hickory at taxonomic and functional levels in bulk soil and the rhizosphere. We observed that pecan has a stronger capacity to enrich rhizosphere plant-beneficial microbe bacteria (e.g., Rhizobium, Novosphingobium, Variovorax, Sphingobium, and Sphingomonas) and their associated functional traits than hickory. We also noted that the ABC transporters (e.g., monosaccharide transporter) and bacterial secretion systems (e.g., type IV secretion system) are the core functional traits of pecan rhizosphere bacteria. Rhizobium and Novosphingobium are the main contributors to the core functional traits. These results suggest that monosaccharides may help Rhizobium to efficiently enrich this niche. Novosphingobium may use a type IV secretion system to interact with other bacteria and thereby influence the assembly of pecan rhizosphere microbiomes. Our data provide valuable information to guide core microbial isolation and expand our knowledge of the assembly mechanisms of plant rhizosphere microbes. IMPORTANCE The rhizosphere microbiome has been identified as a fundamental factor in maintaining plant health, helping plants to fight the deleterious effects of diseases and abiotic stresses. However, to date, studies on the nut tree microbiome have been scarce. Here, we observed a significant "rhizosphere effect" on the seedling pecan. We furthermore demonstrated the core rhizosphere microbiome and function in the seedling pecan. Moreover, we deduced possible factors that help the core bacteria, such as Rhizobium, to efficiently enrich the pecan rhizosphere and the importance of the type IV system for the assembly of pecan rhizosphere bacterial communities. Our findings provide information for understanding the mechanism of the rhizosphere microbial community enrichment process.

Genome-Wide Identification, Characterization, and Expression Analysis of Long-Chain Acyl-CoA Synthetases in Carya illinoinensis under Different Treatments

As crucial enzymes in the lipid metabolic network, long-chain acyl-CoA synthases (LACSs) are members of the acyl-activated enzyme superfamily and play a crucial role in epidermal wax synthesis, plant lipid anabolic metabolism, and stress tolerance. In this study, 11 pecan LACS genes were identified and categorized into five groups and located on nine chromosomes. The significant degree of conservation in the AtLACS and CiLACS protein sequences was demonstrated by multiple sequence alignment and conserved motif analysis. Cis-acting element analysis identified numerous stress-responsive and hormone-inducible elements in the promoter regions of CiLACS genes. The expression levels of CiLACS9 and CiLACS9-1 were considerably up-regulated under salt and drought stress, according to the qRT-RCR study. Treatment with ABA also led to increased expression levels of CiLACS1, CiLACS1-1, CiLACS2, and CiLACS9-1. Notably, CiLACS4, CiLACS4-1, CiLACS9, and CiLACS9-1 exhibited peak expression levels at 135 days after anthesis and are likely to have been crucial in the accumulation of seed kernel oil. Moreover, the CiLACS9 gene was shown to be located in the cytoplasm. These findings offer a theoretical framework for clarifying the roles of LACS genes in the processes of pecan kernel oil synthesis and response to abiotic stressors.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae)

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes and its only congeneric species, Platycarya strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole-genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics and the driving forces underneath the incipient speciation of the two Platycarya lineages.

Chromosome-scale genome assembly and insights into the metabolome and gene regulation of leaf color transition in an important oak species, Quercus dentata

Quercus dentata Thunb., a dominant forest tree species in northern China, has significant ecological and ornamental value due to its adaptability and beautiful autumn coloration, with color changes from green to yellow into red resulting from the autumnal shifts in leaf pigmentation. However, the key genes and molecular regulatory mechanisms for leaf color transition remain to be investigated. First, we presented a high-quality chromosome-scale assembly for Q. dentata. This 893.54 Mb sized genome (contig N50 = 4.21 Mb, scaffold N50 = 75.55 Mb; 2n = 24) harbors 31 584 protein-coding genes. Second, our metabolome analyses uncovered pelargonidin-3-O-glucoside, cyanidin-3-O-arabinoside, and cyanidin-3-O-glucoside as the main pigments involved in leaf color transition. Third, gene co-expression further identified the MYB-bHLH-WD40 (MBW) transcription activation complex as central to anthocyanin biosynthesis regulation. Notably, transcription factor (TF) QdNAC (QD08G038820) was highly co-expressed with this MBW complex and may regulate anthocyanin accumulation and chlorophyll degradation during leaf senescence through direct interaction with another TF, QdMYB (QD01G020890), as revealed by our further protein-protein and DNA-protein interaction assays. Our high-quality genome assembly, metabolome, and transcriptome resources further enrich Quercus genomics and will facilitate upcoming exploration of ornamental values and environmental adaptability in this important genus.

WRINKLED1 Positively Regulates Oil Biosynthesis in Carya cathayensis

Hickory (Carya cathayensis Sarg.) is a kind of important woody oil tree species, and its nut has high nutritional value. Previous gene coexpression analysis showed that WRINKLED1 (WRI1) may be a core regulator during embryo oil accumulation in hickory. However, its specific regulatory mechanism on hickory oil biosynthesis has not been investigated. Herein, two hickory orthologs of WRI1 (CcWRI1A and CcWRI1B) containing two AP2 domains with AW-box binding sites and three intrinsically disordered regions (IDRs) but lacking the PEST motif in the C-terminus were characterized. They are nucleus-located and have self-activated ability. The expression of these two genes was tissue-specific and relatively high in the developing embryo. Notably, CcWRI1A and CcWRI1B can restore the low oil content, shrinkage phenotype, composition of fatty acid, and expression of oil biosynthesis pathway genes of Arabidopsis wri1-1 mutant seeds. Additionally, CcWRI1A/B were shown to modulate the expression of some fatty acid biosynthesis genes in the transient expression system of nonseed tissues. Transcriptional activation analysis further indicated that CcWRI1s directly activated the expression of SUCROSE SYNTHASE2 (SUS2), PYRUVATE KINASE β SUBUNIT 1 (PKP-β1), and BIOTIN CARBOXYL CARRIER PROTEIN2 (BCCP2) involved in oil biosynthesis. These results suggest that CcWRI1s can promote oil synthesis by upregulating some late glycolysis- and fatty acid biosynthesis-related genes. This work reveals the positive function of CcWRI1s in oil accumulation and provides a potential target for improving plant oil by bioengineering technology.

Transcriptome Analysis of Resistant and Susceptible Pecan (Carya illinoinensis) Reveals the Mechanism of Resistance to Black Spot Disease (Colletotrichum fioriniae)

Pecan, Carya illinoinensis (Wangenh.) K. Koch, is an important dried fruit and woody oil tree species grown worldwide. With continuous expansion of pecan cultivation, the frequency and scope of diseases, especially black spot disease, are increasing, damaging trees and reducing yields. In this study, the key factors in resistance to black spot disease (Colletotrichum fioriniae) were investigated between the high-resistance pecan variety "Kanza" and the low-resistance variety "Mahan". Leaf anatomy and antioxidase activities confirmed much stronger resistance to black spot disease in "Kanza" than in "Mahan". Transcriptome analysis indicated that the increased expression of genes associated with defense response, oxidation-reduction, and catalytic activity was involved in disease resistance. A connection network identified a highly expressed hub gene CiFSD2 (CIL1242S0042), which might participate in redox reactions to affect disease resistance. Overexpression of CiFSD2 in tobacco inhibited enlargement of necrotic spots and increased disease resistance. Overall, the expression of differentially expressed genes differed in pecan varieties with different levels of resistance to C. fioriniae infection. In addition, the hub genes associated with black spot resistance were identified and the functions clarified. The in-depth understanding of resistance to black spot disease provides new insights for early screening of resistant varieties and molecular-assisted breeding in pecan.

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

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

Transcriptional Regulation of Female and Male Flower Bud Initiation and Development in Pecan (Carya illinoensis)

Pecan (Carya illinoensis) nuts are delicious and rich in unsaturated fatty acids, which are beneficial for human health. Their yield is closely related to several factors, such as the ratio of female and male flowers. We sampled and paraffin-sectioned female and male flower buds for one year and determined the stages of initial flower bud differentiation, floral primordium formation, and pistil and stamen primordium formation. We then performed transcriptome sequencing on these stages. Our data analysis suggested that FLOWERING LOCUS T (FT) and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 play a role in flower bud differentiation. J3 was highly expressed in the early stage of female flower buds and may play a role in regulating flower bud differentiation and flowering time. Genes such as NF-YA1 and STM were expressed during male flower bud development. NF-YA1 belongs to the NF-Y transcription factor family and may initiate downstream events leading to floral transformation. STM promoted the transformation of leaf buds to flower buds. AP2 may have been involved in the establishment of floral meristem characteristics and the determination of floral organ characteristics. Our results lay a foundation for the control and subsequent regulation of female and male flower bud differentiation and yield improvement.

Integrative multi-omics analysis of three early diverged rosid species reveals an ancient hierarchical cold-responsive regulatory network

Elucidating regulators, including transcription factors (TFs) and RNA-binding proteins (RBPs), underlying gene transcriptional and post-transcriptional co-regulatory network is key to understand plant cold responses. Previous studies were mainly conducted on single species, and whether the regulators are conserved across different species remains elusive. Here, we selected three species that diverged at the early evolution of rosids (~99-113 million years ago), performed cold-responsive phylotranscriptome experiments, and integrated chromatin immunoprecipitation- and DNA affinity purification-sequencing (ChIP/DAP-seq) analysis to explore cold-responsive regulators and their regulatory networks. First, we detected over 10,000 cold-induced differentially expressed genes (DEGs) and alternative splicing genes (DASGs) in each species. Among the DEGs, a set of TFs and RBPs were conserved in rosid cold response. Compared to TFs, RBPs displayed a delayed cold-responsive pattern, implying a hierarchical regulation of DEGs and DASGs. By integrating DEGs and DASGs, we identified 259 overlapping DE-DASG orthogroups (closely-related homologs) that were cold-regulated at both transcriptional and post-transcriptional levels in all three studied species. Notably, pathway analysis on each of the DEGs, DASGs, and DE-DASGs in the three species showed a common enrichment connected to the circadian rhythm. Evidently, 226 cold-responsive genes were directly targeted by at least two circadian rhythm components (CCA1, LHY, RV4, RVE7, and RVE8). Finally, we revealed an ancient hierarchy of cold-responsive regulatory networks at transcriptional and post-transcriptional levels launched by circadian components in rosids. Altogether, this study sheds light on conserved regulators underlying cold-responsive regulatory networks across rosid species, despite a long evolutionary history after their divergence.

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.

Proteomic Analysis of Pecan (Carya illinoinensis) Nut Development

Pecan (Carya illinoinensis) nuts are an economically valuable crop native to the United States and Mexico. A proteomic summary from two pecan cultivars at multiple time points was used to compare protein accumulation during pecan kernel development. Patterns of soluble protein accumulation were elucidated using qualitative gel-free and label-free mass-spectrometric proteomic analyses and quantitative (label-free) 2-D gel electrophoresis. Two-dimensional (2-D) gel electrophoresis distinguished a total of 1267 protein spots and shotgun proteomics identified 556 proteins. Rapid overall protein accumulation occurred in mid-September during the transition to the dough stage as the cotyledons enlarge within the kernel. Pecan allergens Car i 1 and Car i 2 were first observed to accumulate during the dough stage in late September. While overall protein accumulation increased, the presence of histones diminished during development. Twelve protein spots accumulated differentially based on 2-D gel analysis in the weeklong interval between the dough stage and the transition into a mature kernel, while eleven protein spots were differentially accumulated between the two cultivars. These results provide a foundation for more focused proteomic analyses of pecans that may be used in the future to identify proteins that are important for desirable traits, such as reduced allergen content, improved polyphenol or lipid content, increased tolerance to salinity, biotic stress, seed hardiness, and seed viability.

Genome-wide identification and expression analysis of AUX/LAX family genes in Chinese hickory (Carya cathayensis Sarg.) Under various abiotic stresses and grafting

Auxin is essential for regulating plant growth and development as well as the response of plants to abiotic stresses. AUX/LAX proteins are auxin influx transporters belonging to the amino acid permease family of proton-driven transporters, and are involved in the transport of indole-3-acetic acid (IAA). However, how AUX/LAX genes respond to abiotic stresses in Chinese hickory is less studied. For the first time identification, structural characteristics as well as gene expression analysis of the AUX/LAX gene family in Chinese hickory were conducted by using techniques of gene cloning and real-time fluorescent quantitative PCR. Eight CcAUX/LAXs were identified in Chinese hickory, all of which had the conserved structural characteristics of AUX/LAXs. CcAUX/LAXs were most closely related to their homologous proteins in Populus trichocarpa , which was in consistence with their common taxonomic character of woody trees. CcAUX/LAXs exhibited different expression profiles in different tissues, indicating their varying roles during growth and development. A number of light-, hormone-, and abiotic stress responsive cis-acting regulatory elements were detected on the promoters of CcAUX/LAX genes. CcAUX/LAX genes responded differently to drought and salt stress treatments to varying degrees. Furthermore, CcAUX/LAX genes exhibited complex expression changes during Chinese hickory grafting. These findings not only provide a valuable resource for further functional validation of CcAUX/LAXs, but also contribute to a better understanding of their potential regulatory functions during grafting and abiotic stress treatments in Chinese hickory.

Comprehensive identification and analysis of circRNAs during hickory (Carya cathayensis Sarg.) flower bud differentiation

Flower bud differentiation represents a crucial transition from vegetative growth to reproductive development. Carya cathayensis (hickory) is an important economic species in China, with a long juvenile period that hinders its commercial development. In recent years, circular RNAs (circRNAs) have been widely studied and identified as sponges for miRNA regulation of mRNA expression. However, little is known regarding the role of circRNAs in flower buds. In this study, we sequenced circRNAs at three developmental stages (undifferentiated, differentiating, and fully differentiated) in both female and male buds. A total of 6,931 circRNAs were identified in the three developmental stages and 4,449 and 2,209 circRNAs were differentially expressed in female and male buds, respectively. Gene ontology demonstrated that many circRNA host genes participated in various processes, for example, cellular and intracellular pH regulation. Function annotation identified 46 differentially expressed circRNAs involved in flowering regulation, with 28 circRNAs found only in female buds, 4 found only in male buds, and 11 found in both female and male buds. A circRNA-miRNA-mRNA network was predicted based on 13 flowering-related circRNAs and their seven putative interacting miRNAs to describe the regulatory mechanism. Our preliminary results demonstrated a potential involvement of circRNA in bud differentiation. They provided a preliminary theoretical basis for how circRNA might participate in flower development in hickory, perhaps in woody plants.

Metabolomics and Transcriptomics Analyses Reveal Regulatory Networks Associated with Fatty Acid Accumulation in Pecan Kernels

Pecans are a globally important tree nut crop. Pecan nuts are rich in fatty acids (FAs), proteins, and flavonoids in addition to thiamine and numerous micronutrients. Although several of these nutriments have been studied in this plant, the comprehensive metabolite variations and molecular mechanisms associated with them have not been fully elucidated. In this study, untargeted metabolomics and transcriptomics were integrated to reveal the metabolite accumulation patterns and their associated molecular mechanisms during pecan kernel development. In total, 4260 (under positive mode) and 2726 (under negative mode) high quality features were retained. Overall, 163 differentially accumulated metabolites were identified. Most components were classified into the categories "organic acids and derivatives" and "lipids and lipid-like molecules." The accumulation patterns of amino acids, FAs, carbohydrates, organic acids, vitamins, flavonoids, and phenylpropanoids alongside embryo development were determined. Furthermore, transcriptomes from four pecan kernel developmental stages were used to assess transcript expression levels. Coexpression analyses were performed between FAs and their related genes. This study provides a comprehensive overview of the metabolic changes and regulations during pecan kernel development. We believe that the identification of nutriment accumulation trends and hub genes associated with the biosynthesis of the components will be valuable for genetically improving this plant.

Identification and Expression Analysis of MPK and MKK Gene Families in Pecan (Carya illinoinensis)

Mitogen-activated protein kinases consist of three kinase modules composed of MPKs, MKKs, and MPKKKs. As members of the protein kinase (PK) superfamily, they are involved in various processes, such as developmental programs, cell division, hormonal progression, and signaling responses to biotic and abiotic stresses. In this study, a total of 18 MPKs and 10 MKKs were annotated on the pecan genome, all of which could be classified into four subgroups, respectively. The gene structures and conserved sequences of family members in the same branch were relatively similar. All MPK proteins had a conserved motif TxY, and D(L/I/V)K and VGTxxYMSPER existed in all MKK proteins. Duplication events contributed largely to the expansion of the pecan MPK and MKK gene families. Phylogenetic analysis of protein sequences from six plants indicated that species evolution occurred in pecan. Organ-specific expression profiles of MPK and MKK showed functional diversity. Ka/Ks values indicated that all genes with duplicated events underwent strong negative selection. Seven CiPawMPK and four CiPawMKK genes with high expression levels were screened by transcriptomic data from different organs, and these candidates were validated by qRT-PCR analysis of hormone-treated and stressed samples.

Population Genetic Characteristics and Mating Type Frequency of Venturia effusa from Pecan in South America

Scab, caused by the plant-pathogenic fungus Venturia effusa, is a major disease of pecan in South America, resulting in loss of quantity and quality of nut yield. Characteristics of the populations of V. effusa in South America are unknown. We used microsatellites to describe the genetic diversity and population structure of V. effusa in South America, and determined the mating type status of the pathogen. The four hierarchically sampled orchard populations from Argentina (AR), Brazil (BRC and BRS), and Uruguay (UR) had moderate to high genotypic and gene diversity. There was evidence of population differentiation (F_st = 0.196) but the correlation between geographic distance and genetic distance was not statistically significant. Genetic differentiation was minimal between the UR, BRC, and BRS populations, and these populations were more clearly differentiated from the AR population. The MAT1-1 and MAT1-2 mating types occurred in all four orchards and their frequencies did not deviate from the 1:1 ratio expected under random mating; however, multilocus linkage equilibrium was rejected in three of the four populations. The population genetics of South American populations of V. effusa has many similarities to the population genetics of V. effusa previously described in the United States. Characterizing the populations genetics and reproductive systems of V. effusa are important to establish the evolutionary potential of the pathogen and, thus, its adaptability-and can provide a basis for informed approaches to utilizing available host resistance and determining phytosanitary needs.

Expression analysis of PIN family genes in Chinese hickory reveals their potential roles during grafting and salt stress

Grafting is an effective way to improve Chinese hickory while salt stress has caused great damage to the Chinese hickory industry. Grafting and salt stress have been regarded as the main abiotic stress types for Chinese hickory. However, how Chinese hickory responds to grafting and salt stress is less studied. Auxin has been proved to play an essential role in the stress response through its re-distribution regulation mediated by polar auxin transporters, including PIN-formed (PIN) proteins. In this study, the PIN gene family in Chinese hickory (CcPINs) was identified and structurally characterized for the first time. The expression profiles of the genes in response to grafting and salt stress were determined. A total of 11 CcPINs with the open reading frames (ORFs) of 1,026-1,983 bp were identified. Transient transformation in tobacco leaves demonstrated that CcPIN1a, CcPIN3, and CcPIN4 were localized in the plasma membrane. There were varying phylogenetic relationships between CcPINs and homologous genes in different species, but the closest relationships were with those in Carya illinoinensis and Juglans regia. Conserved N- and C-terminal transmembrane regions as well as sites controlling the functions of CcPINs were detected in CcPINs. Five types of cis-acting elements, including hormone- and stress-responsive elements, were detected on the promoters of CcPINs. CcPINs exhibited different expression profiles in different tissues, indicating their varied roles during growth and development. The 11 CcPINs responded differently to grafting and salt stress treatment. CcPIN1a might be involved in the regulation of the grafting process, while CcPIN1a and CcPIN8a were related to the regulation of salt stress in Chinese hickory. Our results will lay the foundation for understanding the potential regulatory functions of CcPIN genes during grafting and under salt stress treatment in Chinese hickory.

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.

Genome-Wide Analysis of the GDSL Genes in Pecan (Carya illinoensis K. Koch): Phylogeny, Structure, Promoter Cis-Elements, Co-Expression Networks, and Response to Salt Stresses

The Gly-Asp-Ser-Leu (GDSL)-lipase family is a large subfamily of lipolytic enzymes that plays an important role in plant growth and defense against environmental stress. However, little is known about their function in pecans (Carya illinoensis K. Koch). In this study, 87 CilGDSLs were identified and divided into 2 groups and 12 subgroups using phylogenetic analysis; members of the same sub-branch had conserved gene structure and motif composition. The majority of the genes had four introns and were composed of an α-helix and a β-strand. Subcellular localization analysis revealed that these genes were localized in the extracellular matrix, chloroplasts, cytoplasm, nucleus, vacuole, and endoplasmic reticulum, and were validated by transient expression in tobacco mesophyll cells. Furthermore, the analysis of the promoter cis-elements for the CilGDSLs revealed the presence of plant anaerobic induction regulatory, abscisic acid response, light response elements, jasmonic acid (JA) response elements, etc. The qRT-PCR analysis results in “Pawnee” with salt treatment showed that the CilGDSL42.93 (leaf) and CilGDSL39.88 (root) were highly expressed in different tissues. After salt stress treatment, isobaric tags for relative and absolute quantitation (iTRAQ) analysis revealed the presence of a total of ten GDSL proteins. Moreover, the weighted gene co-expression network analysis (WGCNA) showed that one set of co-expressed genes (module), primarily CilGDSL41.11, CilGDSL39.49, CilGDSL34.85, and CilGDSL41.01, was significantly associated with salt stress in leaf. In short, some of them were shown to be involved in plant defense against salt stress in this study.

Whole-Transcriptome Analysis Reveals Long Noncoding RNAs Involved in Female Floral Development of Hickory (Carya cathayensis Sarg.)

Hickory, an endemic woody oil and fruit tree species in China, is of great economic value. However, hickory has a long juvenile period and an inconsistent flowering of males and females, thus influencing the bearing rates and further limiting fruits yield. Currently, it is reported that long noncoding RNAs (lncRNAs) play critical regulatory roles in biological processes. However, the role of lncRNAs in the development of hickory female flowers remains unclear. In this study, a total of 6,862 putative lncRNAs were identified from the female flower transcriptomes in three different growth stages of hickory. We proposed that lncRNAs might play an important role in phytohormone signaling processes for flower formation, especially in the abscisic acid and jasmonic acid pathways, according to the results of our Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. Moreover, we predicted the interactions among four microRNAs (miRNAs), three lncRNAs, and four genes. We proposed that facing the changing environment, LNC_002115 competes with PHOSPHATE2 (PHO2) for the binding sites on cca-miR399f, and protects PHO2 from suppression. In addition, cis-acting LNC_002115 regulates the expression of the SHORT VEGETATIVE PHASE (SVP) by influencing ABRE-binding factor (ABF). In brief, LNC_002115 regulates hickory female floral development by impacting both PHO2 and SVP. This study was the first to identify lncRNAs involved in hickory female floral development, and provided new insight to elucidate how lncRNAs and their targets play a role in female floral development in hickory, thus unfolding the opportunities for functional characterization of blossom-related lncRNAs in further studies.

Insight into the CBL and CIPK gene families in pecan (Carya illinoinensis): identification, evolution and expression patterns in drought response

BACKGROUND

Calcium (Ca²⁺) serves as a ubiquitous second messenger and plays a pivotal role in signal transduction. Calcineurin B-like proteins (CBLs) are plant-specific Ca²⁺ sensors that interact with CBL-interacting protein kinases (CIPKs) to transmit Ca²⁺ signals. CBL-CIPK complexes have been reported to play pivotal roles in plant development and response to drought stress; however, limited information is available about the CBL and CIPK genes in pecan, an important nut crop.

RESULTS

In the present study, a total of 9 CBL and 30 CIPK genes were identified from the pecan genome and divided into four and five clades based on phylogeny, respectively. Gene structure and distribution of conserved sequence motif analysis suggested that family members in the same clade commonly exhibited similar exon-intron structures and motif compositions. The segmental duplication events contributed largely to the expansion of pecan CBL and CIPK gene families, and Ka/Ks values revealed that all of them experienced strong negative selection. Phylogenetic analysis of CIPK proteins from 14 plant species revealed that CIPKs in the intron-poor clade originated in seed plants. Tissue-specific expression profiles of CiCBLs and CiCIPKs were analysed, presenting functional diversity. Expression profiles derived from RNA-Seq revealed distinct expression patterns of CiCBLs and CiCIPKs under drought treatment in pecan. Moreover, coexpression network analysis helped to elucidate the relationships between these genes and identify potential candidates for the regulation of drought response, which were verified by qRT-PCR analysis.

CONCLUSIONS

The characterization and analysis of CBL and CIPK genes in pecan genome could provide a basis for further functional analysis of CiCBLs and CiCIPKs in the drought stress response of pecan.

Integrated transcriptome and proteome analysis of developing embryo reveals the mechanisms underlying the high levels of oil accumulation in Carya cathayensis Sarg

Hickory (Carya cathayensis Sarg.) is an extraordinary nut-bearing deciduous arbor with high content of oil in its embryo. However, the molecular mechanism underlying high oil accumulation is mostly unknown. Here, we reported that the lipid droplets and oil accumulation gradually increased with the embryo development and the oil content was up to ~76% at maturity. Furthermore, transcriptome and proteome analysis of developing hickory embryo identified 32,907 genes and 9857 proteins. Time-series analysis of gene expressions showed that these genes were divided into 12 clusters and lipid metabolism-related genes were enriched in Cluster 3, with the highest expression levels at 95 days after pollination (S2). Differentially expressed genes and proteins indicated high correlation, and both were enriched in the lipid metabolism. Notably, the genes involved in biosynthesis, transport of fatty acid/lipid and lipid droplets formation had high expression levels at S2, while the expression levels of other genes required for suberin/wax/cutin biosynthesis and lipid degradation were very low at all the sampling time points, ultimately promoting the accumulation of oil. Quantitative reverse-transcription PCR analysis also verified the results of RNA-seq. The co-regulatory networks of lipid metabolism were further constructed and WRINKLED1 (WRI1) was a core transcriptional factor located in the nucleus. Of note, CcWRI1A/B could directly activate the expression of some genes (CcBCCP2A, CcBCCP2B, CcFATA and CcFAD3) required for fatty acid synthesis. These results provided in-depth evidence for revealing the molecular mechanism of high oil accumulation in hickory embryo.

Application of titanium regulates the functional components of photosynthetic apparatus in grafted seedlings of Carya cathayensis Sarg. under shade

Light acts as a key environmental factor for normal growth and development of plants. Carya cathayensis Sarg. (hickory) faces low light conditions, especially those caused by cloudy or rainy days during the rapid growth period, which has caused adverse effects on its growth. In the current investigation, to alleviate the adverse effects of insufficient light on the cultivation of hickory, anti-hydrolyze stabilized ionic titanium (ASIT) was sprayed on the leaves of the three kinds of grafted seedlings and the non-grafted seedlings of hickory grown under different shade conditions. Results showed that the leaf mass per area and chlorophyll content of grafted hickory seedlings were increased after ASIT application. Rubisco content and photosynthetic rate (P_n) of seedlings grown under shading conditions were positively affected by ASIT treatment, especially on the 45th day of treatment, while the interaction effects of the two parameters between ASIT application and different shade treatments were significant. Titanium accumulation was the highest in roots, followed by leaves, and then in stems, while ASIT had the most significant effects on roots and leaves under 50 ± 5% shade. Severe shading inhibited growth and lead to serious destruction of chloroplast ultrastructure. In addition, the role of ASIT was rootstock-dependent, since ASIT had the weakest mitigation effect on the C/H grafted seedlings. To sum up, the application of ASIT to the grafted seedlings of hickory could improve its ability to resist shade stress.

Comparative Transcriptome Analysis Reveals Differential Regulation of Flavonoids Biosynthesis Between Kernels of Two Pecan Cultivars

Flavonoids influence the flavor and nutritional value of pecan nuts. However, limited information is available regarding the molecular mechanisms underlying pecan flavonoid biosynthesis. Here, we used a high ("YLC28") and a low ("Oconee") flavonoid content cultivar as the research objects. The changes in flavonoid content and the gene transcription patterns during kernel development were identified. Different accumulation patterns of total flavonoids (TF) and condensed tannins (CT) were observed between the two cultivars. The contents of TF and CT in "YLC28" were 1.76- and 2.67-fold higher levels than that of "Oconee" on 150 days after full bloom of female flowers, respectively. In total, 30 RNA-Seq libraries were constructed and sequenced. The upregulated genes in "YLC28" were highly enriched in flavonoid-related pathways. Thirty-three structural genes were identified, and the expression of two phenylalanine ammonia lyases, one chalcone synthase, one flavonoid 3',5'-hydroxylase, and one flavonol synthase exhibited high correlation (r ≥ 0.7, p < 0.01) with the condensed tannin content in "YLC28." A putative MYB transcription factor, CIL1093S0100, might act as a flavonoid biosynthesis repressor during kernel development. Altogether, these results will be useful for uncovering the molecular mechanisms of flavonoid biosynthesis and subsequently accelerating quality pecan breeding.

The Complete Chloroplast Genome of Carya cathayensis and Phylogenetic Analysis

Carya cathayensis, an important economic nut tree, is narrowly endemic to eastern China in the wild. The complete cp genome of C. cathayensis was sequenced with NGS using an Illumina HiSeq2500, analyzed, and compared to its closely related species. The cp genome is 160,825 bp in length with an overall GC content of 36.13%, presenting a quadripartite structure comprising a large single copy (LSC; 90,115 bp), a small single copy (SSC; 18,760 bp), and a pair of inverted repeats (IRs; 25,975 bp). The genome contains 129 genes, including 84 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. A total of 252 simple sequence repeats (SSRs) and 55 long repeats were identified. Gene selective pressure analysis showed that seven genes (rps15, rpoA, rpoB, petD, ccsA, atpI, and ycf1-2) were possibly under positive selection compared with the other Juglandaceae species. Phylogenetic relationships of 46 species inferred that Juglandaceae is monophyletic, and that C. cathayensis is sister to Carya kweichowensis and Carya illinoinensis. The genome comparison revealed that there is a wide variability of the junction sites, and there is higher divergence in the noncoding regions than in coding regions. These results suggest a great potential in phylogenetic research. The newly characterized cp genome of C. cathayensis provides valuable information for further studies of this economically important species.

Comparative plastomes of Carya species provide new insights into the plastomes evolution and maternal phylogeny of the genus

Carya, in the Juglandiodeae subfamily, is to a typical temperate-subtropical forest-tree genus for studying the phylogenetic evolution and intercontinental disjunction between eastern Asia (EA) and North America (NA). Species of the genus have high economic values worldwide for their high-quality wood and the rich healthy factors of their nuts. Although previous efforts based on multiple molecular markers or genome-wide SNPs supported the monophyly of Carya and its two EA and NA major subclades, the maternal phylogeny of Carya still need to be comprehensively evaluated. The variation of Carya plastome has never been thoroughly characterized. Here, we novelly present 19 newly generated plastomes of congeneric Carya species, including the recently rediscovered critically endangered C. poilanei. The overall assessment of plastomes revealed highly conservative in the general structures. Our results indicated that remarkable differences in several plastome features are highly consistent with the EA-NA disjunction and showed the relatively diverse matrilineal sources among EA Carya compared to NA Carya. The maternal phylogenies were conducted with different plastome regions and full-length plastome datasets from 30 plastomes, representing 26 species in six genera of Juglandoideae and Myrica rubra (as root). Six out of seven phylogenetic topologies strongly supported the previously reported relationships among genera of Juglandoideae and the two subclades of EA and NA Carya, but displayed significant incongruencies between species within the EA and NA subclades. The phylogenetic tree generated from full-length plastomes demonstrated the optimal topology and revealed significant geographical maternal relationships among Carya species, especially for EA Carya within overlapping distribution areas. The full-length plastome-based phylogenetic topology also strongly supported the taxonomic status of five controversial species as separate species of Carya. Historical and recent introgressive hybridization and plastid captures might contribute to plastome geographic patterns and inconsistencies between topologies built from different datasets, while incomplete lineage sorting could account for the discordance between maternal topology and the previous nuclear genome data-based phylogeny. Our findings highlight full-length plastomes as an ideal tool for exploring maternal relationships among the subclades of Carya, and potentially in other outcrossing perennial woody plants, for resolving plastome phylogenetic relationships.

Global Transcriptome Analysis Revealed the Molecular Regulation Mechanism of Pigment and Reactive Oxygen Species Metabolism During the Stigma Development of Carya cathayensis

Hickory (Carya cathayensis Sarg.) is a monoecious plant of the genus Carya of the Juglandaceae family. Its nuts contain a number of nutritional compounds and are deeply loved by consumers. Interestingly, it was observed that the color of hickory stigma changed obviously from blooming to mature. However, the molecular mechanism underlying color formation during stigma development and the biological significance of this phenomenon was mostly unknown. In this work, pigment content, reactive oxygen species (ROS) removal capacity, and transcriptome analysis of developing stigma of hickory at 4 differential sampling time points (S1, S2, S3, and S4) were performed to reveal the dynamic changes of related pigment, antioxidant capacity, and its internal molecular regulatory mechanism. It was found that total chlorophyll content was decreased slightly from S1 to S4, while total carotenoids content was increased from S1 to S3 but decreased gradually from S3 to S4. Total anthocyanin content continued to increase during the four periods of stigma development, reaching the highest level at the S4. Similarly, the antioxidant capacity of stigma was also gradually improved from S1 to S4. Furthermore, transcriptome analysis of developing hickory stigma identified 31,027 genes. Time-series analysis of gene expressions showed that these genes were divided into 12 clusters. Cluster 5 was enriched with some genes responsible for porphyrin and chlorophyll metabolism, carotenoid metabolism, and photosynthesis. Meanwhile, cluster 10 was enriched with genes related to flavonoid metabolism, including anthocyanin involved in ROS scavenging, and its related genes were mainly distributed in cluster 12. Based on the selected threshold values, a total of 10432 differentially expressed genes were screened out and enriched in the chlorophyll, carotenoid, anthocyanin, and ROS metabolism. The expression trends of these genes provided plausible explanations for the dynamic change of color and ROS level of hickory stigma with development. qRT-PCR analyses were basically consistent with the results of RNA-seq. The gene co-regulatory networks of pigment and ROS metabolism were further constructed and MYB113 (CCA0887S0030) and WRKY75 (CCA0573S0068) were predicted to be two core transcriptional regulators. These results provided in-depth evidence for revealing the molecular mechanism of color formation in hickory stigma and its biological significance.

Reference-free discovery of nuclear SNPs permits accurate, sensitive identification of Carya (hickory) species and hybrids

PREMISE

DNA-based species identification is critical when morphological identification is restricted, but DNA-based identification pipelines typically rely on the ability to compare homologous sequence data across species. Because many clades lack robust genomic resources, we present here a bioinformatics pipeline capable of generating genome-wide single-nucleotide polymorphism (SNP) data while circumventing the need for any reference genome or annotation data.

METHODS

Using the SISRS bioinformatics pipeline, we generated de novo ortholog data for the genus Carya, isolating sites where genetic variation was restricted to a single Carya species (i.e., species-informative SNPs). We leveraged these SNPs to identify both full-species and hybrid Carya specimens, even at very low sequencing depths.

RESULTS

We identified between 46,000 and 476,000 species-identifying SNPs for each of eight diploid Carya species, and all species identifications were concordant with the species of record. For all putative F1 hybrid specimens, both parental species were correctly identified in all cases, and more punctate patterns of introgression were detectable in more cryptic crosses.

DISCUSSION

Bioinformatics pipelines that use only short-read sequencing data provide vital new tools enabling rapid expansion of DNA identification assays for model and non-model clades alike.

Characterization and expression analysis of the SPL gene family during floral development and abiotic stress in pecan (Carya illinoinensis)

SQUAMOSA promoter binding protein-like (SPL) genes are a type of plant-specific transcription factors that play crucial roles in the regulation of phase transition, floral transformation, fruit development, and various stresses. Although SPLs have been characterized in several model species, no systematic analysis has been studied in pecans, an important woody oil tree species. In this study, a total of 32 SPL genes (CiSPLs) were identified in the pecan genome. After conducting phylogenetic analysis of the conserved SBP proteins from Arabidopsis, rice, and poplar, the CiSPLs were separated into eight subgroups. The CiSPL genes within the same subgroup contained very similar exon-intron structures and conserved motifs. Nine segmentally duplicated gene pairs in the pecan genome and 16 collinear gene pairs between the CiSPL and AtSPL genes were identified. Cis-element analysis showed that CiSPL genes may regulate plant meristem differentiation and seed development, participate in various biological processes, and respond to plant hormones and environmental stresses. Therefore, we focused our study on the expression profiles of CiSPL genes during flower and fruit development. Most of the CiSPL genes were predominantly expressed in buds and/or female flowers. Additionally, quantitative real time PCR (qRT-PCR) analyses confirmed that CiSPL genes showed distinct spatiotemporal expression patterns in response to drought and salt treatments. The study provides foundation for the further exploration of the function and evolution of SPL genes in pecan.

Reassessment of Annamocarya sinesis (Carya sinensis) Taxonomy through Concatenation and Coalescence Phylogenetic Analysis

Due to its peculiar morphological characteristics, there is dispute as to whether the genus of Annamocarya sinensis, a species of Juglandaceae, is Annamocarya or Carya. Most morphologists believe it should be distinguished from the Carya genus while genomicists suggest that A. sinensis belongs to the Carya genus. To explore the taxonomic status of A. sinensis using chloroplast genes, we collected chloroplast genomes of 16 plant species and assembled chloroplast genomes of 10 unpublished Carya species. We analyzed all 26 species' chloroplast genomes through two analytical approaches (concatenation and coalescence), using the entire and unique chloroplast coding sequence (CDS) and entire and protein sequences. Our results indicate that the analysis of the CDS and protein sequences or unique CDS and unique protein sequence of chloroplast genomes shows that A. sinensis indeed belongs to the Carya genus. In addition, our analysis shows that, compared to single chloroplast genes, the phylogeny trees constructed using numerous genes showed higher consistency. Moreover, the phylogenetic analysis calculated with the coalescence method and unique gene sequences was more robust than that done with the concatenation method, particularly for analyzing phylogenetically controversial species. Through the analysis, our results concluded that A. sinensis should be called C. sinensis.

The monoploid chromosome complement of reconstructed ancestral genomes in a phylogeny

Using RACCROCHE, a method for reconstructing gene content and order of ancestral chromosomes from a phylogeny of extant genomes represented by the gene orders on their chromosomes, we study the evolution of three orders of woody plants. The method retrieves the monoploid complement of each Ancestor in a phylogeny, consisting a complete set of distinct chromosomes, despite some of the extant genomes being recently or historically polyploidized. The three orders are the Sapindales, the Fagales and the Malvales. All of these are independently estimated to have ancestral monoploid number [Formula: see text].

Chromosome-scale assembly reveals asymmetric paleo-subgenome evolution and targets for the acceleration of fungal resistance breeding in the nut crop, pecan

Pecan (Carya illinoinensis) is a tree nut crop of worldwide economic importance that is rich in health-promoting factors. However, pecan production and nut quality are greatly challenged by environmental stresses such as the outbreak of severe fungal diseases. Here, we report a high-quality, chromosome-scale genome assembly of the controlled-cross pecan cultivar 'Pawnee' constructed by integrating Nanopore sequencing and Hi-C technologies. Phylogenetic and evolutionary analyses reveal two whole-genome duplication (WGD) events and two paleo-subgenomes in pecan and walnut. Time estimates suggest that the recent WGD event and considerable genome rearrangements in pecan and walnut account for expansions in genome size and chromosome number after the divergence from bayberry. The two paleo-subgenomes differ in size and protein-coding gene sets. They exhibit uneven ancient gene loss, asymmetrical distribution of transposable elements (especially LTR/Copia and LTR/Gypsy), and expansions in transcription factor families (such as the extreme pecan-specific expansion in the far-red impaired response 1 family), which are likely to reflect the long evolutionary history of species in the Juglandaceae. A whole-genome scan of resequencing data from 86 pecan scab-associated core accessions identified 47 chromosome regions containing 185 putative candidate genes. Significant changes were detected in the expression of candidate genes associated with the chitin response pathway under chitin treatment in the scab-resistant and scab-susceptible cultivars 'Excell' and 'Pawnee'. These findings enable us to identify key genes that may be important susceptibility factors for fungal diseases in pecan. The high-quality sequences are valuable resources for pecan breeders and will provide a foundation for the production and quality improvement of tree nut crops.

Genome-wide identification and expression profiles of ABCB gene family in Chinese hickory (Carya cathayensis Sarg.) during grafting

Chinese hickory (Carya cathayensis Sarg.) is an important nut tree species native to China. Excessive plant height and long juvenile phase has restricted development of its industry. Recently, grafting has been used increasingly in production practice of this species to solve the problems above. Previous studies have proved the importance of auxin during Chinese hickory grafting. However, the function of ATP-binding cassette subfamily B (ABCB) genes during Chinese hickory grafting is less studied. In this study, 23 ABCB genes were identified and characterized in Chinese hickory (CcABCBs). The expression profiles of auxin-related ABCBs among tissues, under auxin-related phytohormone treatments and during grafting were determined. CcABCB proteins were divided into half-size and full-size transporters. Many phytohormone-related cis-acting regulatory elements were detected on the promoters of CcABCB genes. Four CcABCB genes homologous to auxin-related AtABCB1, 6, 19 and 20 in Arabidopsis were selected for expression analysis. The four genes displayed varying expression patterns in different tissues of Chinese hickory. Expressions of the four CcABCB genes were regulated by auxin-related phytohormones to varying degrees. Expression levels of the four genes were significantly changed at different stages of grafting, especially 7 days after grafting, indicating their involvement of auxin homeostasis regulation during grafting. In addition, the expressions of CcABCB1 were regulated by IAA and NPA treatments during grafting in comparison with CK treatment, while expressions of the other 3 CcABCB genes were slightly affected. This study will lay the foundation for understanding the potential regulatory roles of CcABCB genes during Chinese hickory grafting.

Comparing whole-genome shotgun sequencing and DNA metabarcoding approaches for species identification and quantification of pollen species mixtures

Molecular identification of mixed-species pollen samples has a range of applications in various fields of research. To date, such molecular identification has primarily been carried out via amplicon sequencing, but whole-genome shotgun (WGS) sequencing of pollen DNA has potential advantages, including (1) more genetic information per sample and (2) the potential for better quantitative matching. In this study, we tested the performance of WGS sequencing methodology and publicly available reference sequences in identifying species and quantifying their relative abundance in pollen mock communities. Using mock communities previously analyzed with DNA metabarcoding, we sequenced approximately 200Mbp for each sample using Illumina HiSeq and MiSeq. Taxonomic identifications were based on the Kraken k-mer identification method with reference libraries constructed from full-genome and short read archive data from the NCBI database. We found WGS to be a reliable method for taxonomic identification of pollen with near 100% identification of species in mixtures but generating higher rates of false positives (reads not identified to the correct taxon at the required taxonomic level) relative to rbcL and ITS2 amplicon sequencing. For quantification of relative species abundance, WGS data provided a stronger correlation between pollen grain proportion and sequence read proportion, but diverged more from a 1:1 relationship, likely due to the higher rate of false positives. Currently, a limitation of WGS-based pollen identification is the lack of representation of plant diversity in publicly available genome databases. As databases improve and costs drop, we expect that eventually genomics methods will become the methods of choice for species identification and quantification of mixed-species pollen samples.

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

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