Resequencing of durian genomes reveals large genetic variations among different cultivars

Durian (Durio zibethinus), which yields the fruit known as the “King of Fruits,” is an important economic crop in Southeast Asia. Several durian cultivars have been developed in this region. In this study, we resequenced the genomes of three popular durian cultivars in Thailand, including Kradumthong (KD), Monthong (MT), and Puangmanee (PM) to investigate genetic diversities of cultivated durians. KD, MT, and PM genome assemblies were 832.7, 762.6, and 821.6 Mb, and their annotations covered 95.7, 92.4, and 92.7% of the embryophyta core proteins, respectively. We constructed the draft durian pangenome and analyzed comparative genomes with related species in Malvales. Long terminal repeat (LTR) sequences and protein families in durian genomes had slower evolution rates than that in cotton genomes. However, protein families with transcriptional regulation function and protein phosphorylation function involved in abiotic and biotic stress responses appeared to evolve faster in durians. The analyses of phylogenetic relationships, copy number variations (CNVs), and presence/absence variations (PAVs) suggested that the genome evolution of Thai durians was different from that of the Malaysian durian, Musang King (MK). Among the three newly sequenced genomes, the PAV and CNV profiles of disease resistance genes and the expressions of methylesterase inhibitor domain containing genes involved in flowering and fruit maturation in MT were different from those in KD and PM. These genome assemblies and their analyses provide valuable resources to gain a better understanding of the genetic diversity of cultivated durians, which may be useful for the future development of new durian cultivars.

[1]  M. Zaynab,et al.  Genome wide study of cysteine rich receptor like proteins in Gossypium sp. , 2022, Scientific Reports.

[2]  Badruzsaufari,et al.  Molecular diversity of exotic durian (Durio spp.) germplasm: a case study of Kalimantan, Indonesia , 2022, Journal of Genetic Engineering and Biotechnology.

[3]  Cris Q. Cortaga,et al.  Identification and characterization of genome-wide resistance gene analogs (RGAs) of durian (Durio zibethinus L.) , 2022, Journal of Genetic Engineering and Biotechnology.

[4]  P. Bayer,et al.  The pangenome of banana highlights differences between genera and genomes , 2021, The plant genome.

[5]  Robert C. Edgar,et al.  MUSCLE v5 enables improved estimates of phylogenetic tree confidence by ensemble bootstrapping , 2021 .

[6]  D. Choi,et al.  Comparative analysis of de novo genomes reveals dynamic intra-species divergence of NLRs in pepper , 2021, BMC plant biology.

[7]  Rafael Della Coletta,et al.  How the pan-genome is changing crop genomics and improvement , 2021, Genome Biology.

[8]  Steven L Salzberg,et al.  Liftoff: accurate mapping of gene annotations , 2020, Bioinform..

[9]  Ben Fulton,et al.  CAFE 5 models variation in evolutionary rates among gene families , 2020, Bioinform..

[10]  Evgeny M. Zdobnov,et al.  OrthoDB in 2020: evolutionary and functional annotations of orthologs , 2020, Nucleic Acids Res..

[11]  Maojun Wang,et al.  Comparative Genome Analyses Highlight Transposon-Mediated Genome Expansion and the Evolutionary Architecture of 3D Genomic Folding in Cotton , 2020, Molecular biology and evolution.

[12]  K. Mayer,et al.  European maize genomes highlight intraspecies variation in repeat and gene content , 2020, Nature Genetics.

[13]  J. Batley,et al.  Plant pan-genomes are the new reference , 2020, Nature Plants.

[14]  Cédric Feschotte,et al.  RepeatModeler2 for automated genomic discovery of transposable element families , 2020, Proceedings of the National Academy of Sciences.

[15]  P. Bayer,et al.  Plant pangenomics: approaches, applications and advancements. , 2020, Current opinion in plant biology.

[16]  Chaochun Wei,et al.  Computational Strategies for Eukaryotic Pangenome Analyses , 2020, The Pangenome.

[17]  R. Paull,et al.  The Durian , 2020 .

[18]  Shujun Ou,et al.  LTR_FINDER_parallel: parallelization of LTR_FINDER enabling rapid identification of long terminal repeat retrotransposons , 2019, Mobile DNA.

[19]  Marie E. Bolger,et al.  MapMan4: A Refined Protein Classification and Annotation Framework Applicable to Multi-Omics Data Analysis. , 2019, Molecular plant.

[20]  Z. Fei,et al.  The tomato pan-genome uncovers new genes and a rare allele regulating fruit flavor , 2019, Nature Genetics.

[21]  Alexey M. Kozlov,et al.  RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference , 2019, Bioinform..

[22]  Alexey M. Kozlov,et al.  ModelTest-NG: A New and Scalable Tool for the Selection of DNA and Protein Evolutionary Models , 2019, bioRxiv.

[23]  D. Neale,et al.  Comparative Genomics of Six Juglans Species Reveals Disease-associated Gene Family Contractions , 2019, bioRxiv.

[24]  S. Kelly,et al.  OrthoFinder: phylogenetic orthology inference for comparative genomics , 2019, Genome Biology.

[25]  Brook T. Moyers,et al.  Sunflower pan-genome analysis shows that hybridization altered gene content and disease resistance , 2018, Nature Plants.

[26]  Xiyin Wang,et al.  Recursive Paleohexaploidization Shaped the Durian Genome , 2018, Plant Physiology.

[27]  Kenneth L. McNally,et al.  Genomic variation in 3,010 diverse accessions of Asian cultivated rice , 2018, Nature.

[28]  Heng Li,et al.  Minimap2: pairwise alignment for nucleotide sequences , 2017, Bioinform..

[29]  Shujun Ou,et al.  LTR_retriever: A Highly Accurate and Sensitive Program for Identification of Long Terminal Repeat Retrotransposons1[OPEN] , 2017, Plant Physiology.

[30]  Mark A. Arick,et al.  Comparative Genomics of an Unusual Biogeographic Disjunction in the Cotton Tribe (Gossypieae) Yields Insights into Genome Downsizing , 2017, Genome biology and evolution.

[31]  N. Nagarajan,et al.  The draft genome of tropical fruit durian (Durio zibethinus) , 2017, Nature Genetics.

[32]  Yue Zhao,et al.  EUPAN enables pan‐genome studies of a large number of eukaryotic genomes , 2017, Bioinform..

[33]  C. K. Chan,et al.  The pangenome of an agronomically important crop plant Brassica oleracea , 2016, Nature Communications.

[34]  Jeffrey T Leek,et al.  Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown , 2016, Nature Protocols.

[35]  J. Roustan,et al.  Ethylene Control of Fruit Ripening: Revisiting the Complex Network of Transcriptional Regulation1 , 2015, Plant Physiology.

[36]  Evgeny M. Zdobnov,et al.  BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs , 2015, Bioinform..

[37]  S. Somsri.  Current status of durian breeding program in Thailand , 2014 .

[38]  Kemal Kazan,et al.  MYC2: the master in action. , 2013, Molecular plant.

[39]  Diqiu Yu,et al.  Arabidopsis WRKY46 coordinates with WRKY70 and WRKY53 in basal resistance against pathogen Pseudomonas syringae. , 2012, Plant science : an international journal of experimental plant biology.

[40]  S. Gupta,et al.  Isolation and characterization of ripening related pectin methylesterase inhibitor gene from banana fruit , 2012, Physiology and Molecular Biology of Plants.

[41]  Mark Yandell,et al.  MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects , 2011, BMC Bioinformatics.

[42]  Toni Gabaldón,et al.  trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses , 2009, Bioinform..

[43]  Jeffery P. Demuth,et al.  The life and death of gene families , 2009, BioEssays : news and reviews in molecular, cellular and developmental biology.

[44]  Stefan Kurtz,et al.  LTRharvest, an efficient and flexible software for de novo detection of LTR retrotransposons , 2008, BMC Bioinformatics.

[45]  Ulrike Zentgraf,et al.  The Antagonist Function of Arabidopsis WRKY53 and ESR/ESP in Leaf Senescence Is Modulated by the Jasmonic and Salicylic Acid Equilibrium , 2007, The Plant Cell Online.

[46]  T. Mengiste,et al.  Arabidopsis WRKY33 transcription factor is required for resistance to necrotrophic fungal pathogens. , 2006, The Plant journal : for cell and molecular biology.

[47]  Burkhard Morgenstern,et al.  Gene prediction in eukaryotes with a generalized hidden Markov model that uses hints from external sources , 2006, BMC Bioinformatics.

[48]  Juan Miguel García-Gómez,et al.  BIOINFORMATICS APPLICATIONS NOTE Sequence analysis Manipulation of FASTQ data with Galaxy , 2005 .

[49]  J. D. Faveri,et al.  A detached leaf bioassay to screen Durian cultivars for susceptibility to Phytophthora palmivora , 2005, Australasian Plant Pathology.

[50]  Robert C. Edgar,et al.  MUSCLE: a multiple sequence alignment method with reduced time and space complexity , 2004, BMC Bioinformatics.

[51]  Ian Korf,et al.  Gene finding in novel genomes , 2004, BMC Bioinformatics.

[52]  D. Guest,et al.  Screening for resistance to Phytophthora. , 2004 .

[53]  J. Beynon,et al.  RPP13 is a simple locus in Arabidopsis thaliana for alleles that specify downy mildew resistance to different avirulence determinants in Peronospora parasitica. , 2000, The Plant journal : for cell and molecular biology.