Signatures of Selection in the Genomes of Chinese Chestnut (Castanea mollissima Blume): The Roots of Nut Tree Domestication

Inhibition of RNA degradation integrates the metabolic signals induced by osmotic stress into the Arabidopsis circadian system

The circadian clock system acts as an endogenous timing reference that coordinates many metabolic and physiological processes in plants. Previous studies have shown that the application of osmotic stress delays circadian rhythms via 3'-phospho-adenosine 5'-phosphate (PAP), a retrograde signalling metabolite that is produced in response to redox stress within organelles. PAP accumulation leads to the inhibition of exoribonucleases (XRNs), which are responsible for RNA degradation. Interestingly, we are now able to demonstrate that post-transcriptional processing is crucial for the circadian response to osmotic stress. Our data show that osmotic stress increases the stability of specific circadian RNAs, suggesting that RNA metabolism plays a vital role in circadian clock coordination during drought. Inactivation of XRN4 is sufficient to extend circadian rhythms as part of this response, with PRR7 and LWD1 identified as transcripts that are post-transcriptionally regulated to delay circadian progression.

Transcriptome analysis of differential sugar accumulation in the developing embryo of contrasting two Castanea mollissima cultivars

Chinese chestnut (Castanea mollissima) is an important nut tree species, and its embryo is rich in sugar. We combined metabolomic and transcriptomic data to analyze metabolites and genes related to sugar in two Chinese chestnut cultivars at 60, 70, 80, 90 and 100 days after flowering (DAF). The soluble sugar content of high-sugar cultivar at maturity is 1.5 times that of low-sugar cultivar. Thirty sugar metabolites were identified in embryo, with the most dominant being sucrose. Analysis of the gene expression patterns revealed that the high-sugar cultivar promoted the conversion of starch to sucrose by up-regulating genes related to starch degradation and sucrose synthesis at 90-100 DAF. It also strongly increased the enzyme activity of SUS-synthetic, which may promote sucrose synthesis. Gene co-expression network analysis showed that ABA and peroxide were related to starch decomposition during Chinese chestnut ripening. Our study analyzed the composition and molecular synthesis mechanism of sugar in Chinese chestnut embryos, and provided a new insight into the regulation pattern of high sugar accumulation in Chinese chestnut nuts.

Genetic Diversity and Population Structure of Chinese Chestnut (Castanea mollissima Blume) Cultivars Revealed by GBS Resequencing

Chinese chestnut (Castanea mollissima Bl.) is one of the earliest domesticated and cultivated fruit trees, and it is widely distributed in China. Because of the high quality of its nuts and its high resistance to abiotic and biotic stresses, Chinese chestnut could be used to improve edible chestnut varieties worldwide. However, the unclear domestication history and highly complex genetic background of Chinese chestnut have prevented the efficiency of breeding efforts. To explore the genetic diversity and structure of Chinese chestnut populations and generate new insights that could aid chestnut breeding, heterozygosity statistics, molecular variance analysis, ADMIXTURE analysis, principal component analysis, and phylogenetic analysis were conducted to analyze single nucleotide polymorphism data from 185 Chinese chestnut landraces from five geographical regions in China via genotyping by sequencing. Results showed that the genetic diversity level of the five populations from different regions was relatively high, with an observed heterozygosity of 0.2796-0.3427. The genetic diversity level of the population in the mid-western regions was the highest, while the population north of the Yellow River was the lowest. Molecular variance analysis showed that the variation among different populations was only 2.07%, while the intra-group variation reached 97.93%. The Chinese chestnut samples could be divided into two groups: a northern and southern population, separated by the Yellow River; however, some samples from the southern population were genetically closer to samples from the northern population. We speculate that this might be related to the migration of humans during the Han dynasty due to the frequent wars that took place during this period, which might have led to the introduction of chestnut to southern regions. Some samples from Shandong Province and Beijing City were outliers that did not cluster with their respective groups, and this might be caused by the special geographical, political, and economic significance of these two regions. The findings of our study showed the complex genetic relationships among Chinese chestnut landraces and the high genetic diversity of these resources.

Domestication and selection footprints in Persian walnuts (Juglans regia)

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

Inhibitory Effect and Mechanism of Dill Seed Essential Oil on Neofusicoccum parvum in Chinese Chestnut

The chestnut postharvest pathogen Neofusicoccum parvum (N. parvum) is an important postharvest pathogen that causes chestnut rot. Chestnut rot in postharvest reduces food quality and causes huge economic losses. This study aimed to evaluate the inhibitory effect of dill seed essential oil (DSEO) on N. parvum and its mechanism of action. The chemical characterization of DSEO by gas chromatography/mass spectrometry (GC/MS) showed that the main components of DSEO were apiole, carvone, dihydrocarvone, and limonene. DSEO inhibited the growth of mycelium in a dose-dependent manner. The antifungal effects are associated with destroying the fungal cell wall (cytoskeleton) and cell membrane. In addition, DSEO can induce oxidative damage and intracellular redox imbalance to damage cell function. Transcriptomics analysis showed DSEO treatment induced differently expressed genes most related to replication, transcription, translation, and lipid, DNA metabolic process. Furthermore, in vivo experiments showed that DSEO and DSEO emulsion can inhibit the growth of fungi and prolong the storage period of chestnuts. These results suggest that DSEO can be used as a potential antifungal preservative in food storage.

Pan-genome analysis of three main Chinese chestnut varieties

Chinese chestnut (Castanea mollissima Blume) is one of the earliest domesticated plants of high nutritional and ecological value, yet mechanisms of C. mollissima underlying its growth and development are poorly understood. Although individual chestnut species differ greatly, the molecular basis of the formation of their characteristic traits remains unknown. Though the draft genomes of chestnut have been previously released, the pan-genome of different variety needs to be studied. We report the genome sequence of three cultivated varieties of chestnut herein, namely Hei-Shan-Zhai-7 (H7, drought-resistant variety), Yan-Hong (YH, easy-pruning variety), and Yan-Shan-Zao-Sheng (ZS, early-maturing variety), to expedite convenience and efficiency in its genetics-based breeding. We obtained three chromosome-level chestnut genome assemblies through a combination of Oxford Nanopore technology, Illumina HiSeq X, and Hi-C mapping. The final genome assemblies are 671.99 Mb (YH), 790.99 Mb (ZS), and 678.90 Mb (H7), across 12 chromosomes, with scaffold N50 sizes of 50.50 Mb (YH), 65.05 Mb (ZS), and 52.16 Mb (H7). Through the identification of homologous genes and the cluster analysis of gene families, we found that H7, YH and ZS had 159, 131, and 91 unique gene families, respectively, and there were 13,248 single-copy direct homologous genes in the three chestnut varieties. For the convenience of research, the chestnut genome database was constructed. Based on the results of gene family identification, the presence/absence variations (PAVs) information of the three sample genes was calculated, and a total of 2,364, 2,232, and 1,475 unique genes were identified in H7, YH and ZS, respectively. Our results suggest that the GBSS II-b gene family underwent expansion in chestnut (relative to nearest source species). Overall, we developed high-quality and well-annotated genome sequences of three C. mollissima varieties, which will facilitate clarifying the molecular mechanisms underlying important traits, and shortening the breeding process.

Comparative Analysis of the Complete Chloroplast Genomes of Four Chestnut Species (Castanea)

Chloroplast (cp) DNA genomes are traditional workhorses for studying the evolution of species and reconstructing phylogenetic relationships in plants. Species of the genus Castanea (chestnuts and chinquapins) are valued as a source of nuts and timber wherever they grow, and chestnut species hybrids are common. We compared the cp genomes of C. mollissima, C. seguinii, C. henryi, and C. pumila. These cp genomes ranged from 160,805 bp to 161,010 bp in length, comprising a pair of inverted repeat (IR) regions (25,685 to 25,701 bp) separated by a large single-copy (LSC) region (90,440 to 90,560 bp) and a small single-copy (SSC) region (18,970 to 19,049 bp). Each cp genome encoded the same 113 genes; 82–83 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. There were 18 duplicated genes in the IRs. Comparative analysis of cp genomes revealed that rpl22 was absent in all analyzed species, and the gene ycf1 has been pseudo-genized in all Chinese chestnuts except C. pumlia. We analyzed the repeats and nucleotide substitutions in these plastomes and detected several highly variable regions. The phylogenetic analyses based on plastomes confirmed the monophyly of Castanea species.

Transcriptome analysis of genes involved in starch biosynthesis in developing Chinese chestnut (Castanea mollissima Blume) seed kernels

Chinese chestnut (Castanea mollissima Blume) seed kernels (CCSK) with high quality and quantity of starch has emerged as a potential raw material for food industry, but the molecular regulatory mechanism of starch accumulation in developing CCSK is still unclear. In this study, we firstly analyzed the fruit development, starch accumulation, and microscopic observation of dynamic accumulation of starch granules of developing CCSK from 10 days after flowering (DAF) to 100 DAF, of which six representative CCSK samples (50–100 DAF) were selected for transcriptome sequencing analysis. Approximately 40 million valid reads were obtained, with an average length of 124.95 bp, which were searched against a reference genome, returning 38,146 unigenes (mean size = 1164.19 bp). Using the DESeq method, 1968, 1573, 1187, 1274, and 1494 differentially expressed unigenes were identified at 60:50, 70:60, 80:70, 90:80 and 100:90 DAF, respectively. The relationship between the unigene transcriptional profiles and starch dynamic patterns in developing CCSK was comparatively analyzed, and the specific unigenes encoding for metabolic enzymes (SUSY2, PGM, PGI, GPT, NTT, AGP3, AGP2, GBSS1, SS1, SBE1, SBE2.1, SBE2.2, ISA1, ISA2, ISA3, and PHO) were characterized to be involved potentially in the biosynthesis of G-1-P, ADPG, and starch. Finally, the temporal transcript profiles of genes encoding key enzymes (susy2, pgi2, gpt1, agp2, agp3, gbss1, ss1, sbe1, sbe2.1, sbe2.2, isa1, isa2, isa3, and pho) were validated by quantitative real-time PCR (qRT-PCR). Our findings could help to reveal the molecular regulatory mechanism of starch accumulation in developing CCSK and may also provide potential candidate genes for increasing starch content in Chinese chestnut or other starchy crops.

Effects of GABA and Vigabatrin on the Germination of Chinese Chestnut Recalcitrant Seeds and Its Implications for Seed Dormancy and Storage

Recalcitrant chestnut seeds are rich in γ-aminobutyric acid (GABA), which negatively regulates adventitious root development by altering carbon/nitrogen metabolism. However, little is known regarding the role of this metabolite in chestnut seeds. In this study, we investigated the effects of GABA changes on the germination of chestnut seeds treated with exogenous GABA and vigabatrin (VGB, which inhibits GABA degradation). Both treatments significantly inhibited seed germination and primary root growth and resulted in the considerable accumulation of H₂O₂, but the endogenous GABA content decreased before germination at 48 h. Soluble sugar levels increased before germination, but subsequently decreased, whereas starch contents were relatively unchanged. Changes to organic acids were observed at 120 h after sowing, including a decrease and increase in citrate and malate levels, respectively. Similarly, soluble protein contents increased at 120 h, but the abundance of most free amino acids decreased at 48 h. Moreover, the total amino acid levels increased only in response to VGB at 0 h. Accordingly, GABA and VGB altered the balance of carbon and nitrogen metabolism, thereby inhibiting chestnut seed germination. These results suggested that changes to GABA levels in chestnut seeds might prevent seed germination. The study data may also help clarify the dormancy and storage of chestnut seeds, as well as other recalcitrant seeds.

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

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

Comparing Pool-seq, Rapture, and GBS genotyping for inferring weak population structure: The American lobster (Homarus americanus) as a case study

Unraveling genetic population structure is challenging in species potentially characterized by large population size and high dispersal rates, often resulting in weak genetic differentiation. Genotyping a large number of samples can improve the detection of subtle genetic structure, but this may substantially increase sequencing cost and downstream bioinformatics computational time. To overcome this challenge, alternative, cost-effective sequencing approaches, namely Pool-seq and Rapture, have been developed. We empirically measured the power of resolution and congruence of these two methods in documenting weak population structure in nonmodel species with high gene flow comparatively to a conventional genotyping-by-sequencing (GBS) approach. For this, we used the American lobster (Homarus americanus) as a case study. First, we found that GBS, Rapture, and Pool-seq approaches gave similar allele frequency estimates (i.e., correlation coefficient over 0.90) and all three revealed the same weak pattern of population structure. Yet, Pool-seq data showed F ST estimates three to five times higher than GBS and Rapture, while the latter two methods returned similar F ST estimates, indicating that individual-based approaches provided more congruent results than Pool-seq. We conclude that despite higher costs, GBS and Rapture are more convenient approaches to use in the case of species exhibiting very weak differentiation. While both GBS and Rapture approaches provided similar results with regard to estimates of population genetic parameters, GBS remains more cost-effective in project involving a relatively small numbers of genotyped individuals (e.g., <1,000). Overall, this study illustrates the complexity of estimating genetic differentiation and other summary statistics in complex biological systems characterized by large population size and migration rates.