The Corylus mandshurica genome provides insights into the evolution of Betulaceae genomes and hazelnut breeding

Abstract

Hazelnut is popular for its flavor, and it has also been suggested that hazelnut is beneficial to cardiovascular health because it is rich in oleic acid. Here, we report the first high-quality chromosome-scale genome for the hazelnut species Corylus mandshurica (2n = 22), which has a high concentration of oleic acid in its nuts. The assembled genome is 367.67 Mb in length, and the contig N50 is 14.85 Mb. All contigs were assembled into 11 chromosomes, and 28,409 protein-coding genes were annotated. We reconstructed the evolutionary trajectories of the genomes of Betulaceae species and revealed that the 11 chromosomes of the hazelnut genus were derived from the most ancestral karyotype in Betula pendula, which has 14 protochromosomes, by inferring homology among five Betulaceae genomes. We identified 96 candidate genes involved in oleic acid biosynthesis, and 10 showed rapid evolution or positive selection. These findings will help us to understand the mechanisms of lipid synthesis and storage in hazelnuts. Several gene families related to salicylic acid metabolism and stress responses experienced rapid expansion in this hazelnut species, which may have increased its stress tolerance. The reference genome presented here constitutes a valuable resource for molecular breeding and genetic improvement of the important agronomic properties of hazelnut.

References

[1] Wu, Z. Vegetation in China (1995).

[2] Wang, G. Progress in cultivation and utilization of Corylus L. Resources in China (I)-Corylus germplasm resources. For. Res. 31, 105-112 (2018).

[3] Ji, J. M., Ge, Z. F., Feng, Y. S. & Wang, X. D. Lipid characterization of Chinese Wild Hazelnuts (Corylus mandshurica Maxim.) . J. Oleo Sci. 68, 13-20 (2019).

[4] Alasalvar, C., Amaral, J. S. & Shahidi, F . Functional lipid characteristics of Turkish Tombul hazelnut (Corylus avellana L.) . J. Agr. Food Chem. 54, 10177-10183 (2006).

[5] Tufekci, F. & Karatas, S. Determination of geographical origin Turkish hazelnuts according to fatty acid composition. Food Sci. Nutr. 6, 557-562 (2018).

[6] Teres, S. et al. Oleic acid content is responsible for the reduction in blood pressure induced by olive oil. Proc. Natl Acad. Sci. USA 105, 13811-13816 (2008).

[7] Cerain, L. D. & Adela. in Mutagenic Act. Meat Samples Deep-Fry Olive Oil Vol. 106, 989-996 (2010).

[8] Perdomo, L. et al. Protective role of oleic acid against cardiovascular insulin resistance and in the early and late cellular atherosclerotic process. Cardiovasc. Diabetol. 14, 75 (2015).

[9] Boccacci, P. & Botta, R. Investigating the origin of hazelnut (Corylus avellana L.) cultivars using chloroplast microsatellites . Genet. Resour. Crop Evol. 56, 851-859 (2009).

[10] Johnson, K. B. et al. Eastern filbert blight of European hazelnut: It’s becoming a manageable disease. Plant Dis. 80, 1308-1316 (1996).

[11] Bhattarai, G., Mehlenbacher, S. & Smith, D. C. Novel sources of resistance to Eastern Filbert blight in Hazelnut. Hortscience 50, S401- S401 (2015).

[12] Kask, K. Nut quality of wild European hazelnut in Estonia and attempts at hazelnut breeding. Acta Hortic. 556, 37-40 (2001).

[13] Li, T. D. et al. Domestication of wild tomato is accelerated by genome editing. Nat. Biotechnol. 36, 1160-1165 (2018).

[14] Chen, Y. Z. et al. High oleic acid content, nontransgenic allotetraploid cotton (Gossypium hirsutum L.) generated by knockout of GhFAD2 genes with CRISPR/Cas9 system . Plant Biotechnol. J. (2020).

[15] Dangl, J. L. & Jones, J. D. G. Plant pathogens and integrated defence responses to infection. Nature 411, 826-833 (2001).

[16] Meyers, B. C. et al. Genome-wide analysis of NBS-LRR-encoding genes in Arabidopsis. Plant Cell 15, 809-834 (2003).

[17] Clarice, J. C., Shawn, A. M. & David, C. S. Sources of resistance to Eastern Filbert Blight in Hazelnut. J. Am. Soc. Hortic. Sci. 123, 253-257 (1998).

[18] Erdogan, V. & Mehlenbacher, S. A. Interspecific hybridization in hazelnut (Corylus) . J. Am. Soc. Hortic. Sci. 125, 489-497 (2000).

[19] Robert, H. W. Cytological studies on the Betulaceae. II. Coryolus Alnus. Chic. J. 88, 383-399 (1929).

[20] Robert, H. W. Cytological studies on the Betulaceae. IV. Betula, Carpinus, Ostrya, Ostryopsis. Chic. J. 90, 108-115 (1930).

[21] Chen, Z. Phylogeny ang phytogeography of the Betulaceae. Acta Phyytotax. Sin. 32, 1-31 (1994).

[22] Rowley, E. R. et al. A draft genome and high-density genetic map of European Hazelnut (Corylus avellana L.) . Preprint at https://www.biorxiv.org/content/10.1101/469015v1 (2018).

[23] Hu, J., Fan, J. P., Sun, Z. Y. & Liu, S. L. NextPolish: a fast and efficient genome polishing tool for long-read assembly. Bioinformatics 36, 2253-2255 (2020).

[24] Simao, F. A. et al. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210-3212 (2015).

[25] Salojarvi, J. et al. Genome sequencing and population genomic analyses provide insights into the adaptive landscape of silver birch. Nat. Genet. 49, 904-915 (2017).

[26] Yang, X. Y. et al. A chromosome-level reference genome of the hornbeam, Carpinus fangiana . Sci. Data 7, 1-24 (2020).

[27] Yang, Y. Z. et al. Genomic effects of population collapse in a critically endangered ironwood tree Ostrya rehderiana . Nat. Commun. 9, 5449-5457 (2018).

[28] Wang, Z. F. et al. Hybrid speciation via inheritance of alternate alleles of parental isolating genes. Mol. Plant 14, 208-222 (2020).

[29] Chen, F. et al. The sequenced angiosperm genomes and genome databases. Front. Plant Sci. 9, 418 (2018).

[30] Bairoch, A. & Apweiler, R. The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res. 28, 45-48 (2000).

[31] Kanehisa, M. & Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res 28, 27-30 (2000).

[32] Xu, Z. S., Chen, M., Li, L. C. & Ma, Y. Z. Functions and application of the AP2/ERF transcription factor family in crop improvement. J. Integr. Plant Biol. 53, 570-585 (2011).

[33] Doherty, C. J., Van Buskirk, H. A., Myers, S. J. & Thomashow, M. F. Roles for Arabidopsis CAMTA transcription factors in cold-regulated gene expression and freezing tolerance. Plant Cell 21, 972-984 (2009).

[34] Roy, S. Function of MYB domain transcription factors in abiotic stress and epigenetic control of stress response in plant genome. Plant Signal Behav. 11, 1-7 (2016).

[35] Toledo-Ortiz, G., Huq, E. & Quail, P. H. The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell 15, 1749-1770 (2003).

[36] Ye, G. F. et al. De novo genome assembly of the stress tolerant forest species Casuarina equisetifolia provides insight into secondary growth. Plant J. 97, 779-794 (2019).

[37] Martinez-Garcia, P. J. et al. The walnut (Juglans regia) genome sequence reveals diversity in genes coding for the biosynthesis of non-structural polyphenols . Plant J. 87, 507-532 (2016).

[38] Plomion, C. et al. Oak genome reveals facets of long lifespan. Nat. Plants 4, 440-452 (2018).

[39] Zapata, L. et al. Chromosome-level assembly of Arabidopsis thaliana Ler reveals the extent of translocation and inversion polymorphisms. Proc. Natl Acad. Sci. USA 113, E4052-E4060 (2016).

[40] Buti, M. et al. The genome sequence and transcriptome of Potentilla micrantha and their comparison to Fragaria vesca (the woodland strawberry) . Gigascience 7, 1-41 (2018).

[41] Jaillon, O. et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature 449, 463-465 (2007).

[42] Chen, Z. D., Manchester, S. R. & Sun, H. Y. Phylogeny and evolution of the Betulaceae as inferred from DNA sequences, morphology, and paleobotany. Am. J. Bot. 86, 1168-1181 (1999).

[43] Chen, Z. D. & Lu, A. M. Phylogeny and evolution od Betulaceae. China Acad. J. 3, 189-191 (2001).

[44] Kosztarab, M ., Roane, M. K. & Drake, C. R. Reduction of Eastern Filbert Blight on Corylus Avellana . Phytopathology 70, 690- 690 (1980).

[45] Sathuvalli, V., Mehlenbacher, S. A. & Smith, D. C. High-resolution genetic and physical mapping of the eastern filbert blight resistance region in ‘Jefferson’ Hazelnut (Corylus avellana L.) . Plant Genome 10, 1-12 (2017).

[46] Bates, P. D., Stymne, S. & Ohlrogge, J. Biochemical pathways in seed oil synthesis. Curr. Opin. Plant Biol. 16, 358-364 (2013).

[47] Thelen, J. J. & Ohlrogge, J. B. Metabolic engineering of fatty acid biosynthesis in plants. Metab. Eng. 4, 12-21 (2002).

[48] Ohlrogge, J. B. Design of new plant-products - engineering of fatty-acid metabolism. Plant Physiol. 104, 821-826 (1994).

[49] Nikolau, B. J., Ohlrogge, J. B. & Wurtele, E. S. Plant biotin-containing carboxylases. Arch. Biochem. Biophys. 414, 211-222 (2003).

[50] Yang, Z. H. PAML: a program package for phylogenetic analysis by maximum likelihood. Comput. Appl. Biosci. 13, 555-556 (1997).

[51] Alasalvar, C., Shahidi, F. & Cadwallader, K. R. Comparison of natural and roasted turkish tombul hazelnut (corylus avellana l.) volatiles and flavor by dha/gc/ms and descriptive sensory analysis . J. Agric. Food Chem. 51, 5067-5072 (2003).

[52] Dormann, P., Voelker, T. A. & Ohlrogge, J. B. Cloning and expression in Escherichia coli of a novel Thioesterase from Arabidopsis thaliana specific for long-chain acyl-acyl carrier proteins. Arch. Biochem. Biophys. 316, 612-618 (1995).

[53] Roesler, K. et al. Targeting of the Arabidopsis homomeric acetyl-coenzyme A carboxylase to plastids of rapeseeds. Plant Physiol. 113, 75-81 (1997).

[54] Kim, Y. S. et al. CAMTA-mediated regulation of salicylic acid immunity pathway genes in Arabidopsis exposed to low temperature and pathogen infection. Plant Cell 29, 2465-2477 (2017).

[55] Wang, Z. Y. & Wang, X. Y. Evolutionary genomics model chromosome number reduction B chromosome production. Sci. Sin. Vitae 50, 524-537 (2020).

[56] Chen, S. F., Zhou, Y. Q., Chen, Y. R. & Gu, J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, 884-890 (2018).

[57] Li, R. Q. et al. The sequence and de novo assembly of the giant panda genome. Nature 463, 311-317 (2010).

[58] Marcais, G. & Kingsford, C. A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics 27, 764-770 (2011).

[59] Vurture, G. W. et al. GenomeScope: fast reference-free genome profiling from short reads. Bioinformatics 33, 2202-2204 (2017).

[60] Roach, M. J., Schmidt, S. A. & Borneman, A. R. Purge Haplotigs: allelic contig reassignment for third-gen diploid genome assemblies. BMC Bioinformatics 19, 460 (2018).

[61] Connell, L. W., Sexton, F. W. & Prinja, A. K. Further development of the heavy ion cross section for single event UPset: model (HICUP). IEEE Trans. Nucl. Sci. 42, 2026-2034 (1995).

[62] Zhang, X. T. et al. Assembly of allele-aware, chromosomal-scale autopolyploid genomes based on Hi-C data. Nat. Plants 5, 833-845 (2019).

[63] Chen, N. in Using RepeatMasker to Identify Repetitive Elements in Genomic Sequences Ch. 4 (ed. Andreas, D. B.) (2004).

[64] Bao, W. D., Kojima, K. K. & Kohany, O. Repbase Update, a database of repetitive elements in eukaryotic genomes. Mob. DNA-Uk 6, 1-6 (2015).

[65] Xu, Z. & Wang, H. LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 35, W265-W268 (2007).

[66] tanke, M. & Morgenstern, B. AUGUSTUS: a web server for gene prediction in eukaryotes that allows user-defined constraints. Nucleic Acids Res. 33, W465-W467 (2005).

[67] ajoros, W. H., Pertea, M. & Salzberg, S. L. TigrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders. Bioinformatics 20, 2878-2879 (2004).

[68] Haas, B. J. et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments. Genome Biol. 9, R7 (2008).

[69] Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389-3402 (1997).

[70] Zdobnov, E. M. & Apweiler, R. InterProScan-an integration platform for the signature-recognition methods in InterPro. Bioinformatics 17, 847-848 (2001).

[71] Reagan, R. L. & Bernstein, R. L. Data mining of signaling proteins using the HMMER method: high selectivity for protein sequence homology searches. Proc. Int. Conf. Math. Eng. Tech. Med. Biol. Sci. sI, II, 185-191 (2000).

[72] Zheng, Y. et al. iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. Mol. Plant 9, 1667-1670 (2016).

[73] Tian, F. et al. PlantRegMap: charting functional regulatory maps in plants. Nucleic Acids Res. 48, D1104-D1113 (2020).

[74] Li, L., Stoeckert, C. J. & Roos, D. S. OrthoMCL: Identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178-2189 (2003).

[75] Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312-1313 (2014).

[76] Puttick, M. N. MCMCtreeR: functions to prepare MCMCtree analyses and visualize posterior ages on trees. Bioinformatics 35, 5321-5322 (2019).

[77] De Bie, T., Cristianini, N., Demuth, J. P. & Hahn, M. W. CAFE: a computational tool for the study of gene family evolution. Bioinformatics 22, 1269-1271 (2006).

[78] Salvatore, C. & Wataru, I. Sonicparanoid: fast, accurate, and easy orthology inference. Bioinformatics 1, 1-3 (2018).

[79] Wang, Y. P. et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40, 1-14 (2012).