Chromosome-scale Genome Assembly Research Articles

The brittle star genome illuminates the genetic basis of animal appendage regeneration

Species within nearly all extant animal lineages are capable of regenerating body parts. However, it remains unclear whether the gene expression programme controlling regeneration is evolutionarily conserved. Brittle stars are a species-rich class of echinoderms with outstanding regenerative abilities, but investigations into the genetic bases of regeneration in this group have been hindered by the limited genomic resources. Here we report a chromosome-scale genome assembly for the brittle star Amphiura filiformis. We show that the brittle star genome is the most rearranged among echinoderms sequenced so far, featuring a reorganized Hox cluster reminiscent of the rearrangements observed in sea urchins. In addition, we performed an extensive profiling of gene expression during brittle star adult arm regeneration and identified sequential waves of gene expression governing wound healing, proliferation and differentiation. We conducted comparative transcriptomic analyses with other invertebrate and vertebrate models for appendage regeneration and uncovered hundreds of genes with conserved expression dynamics, particularly during the proliferative phase of regeneration. Our findings emphasize the crucial importance of echinoderms to detect long-range expression conservation between vertebrates and classical invertebrate regeneration model systems.

A wound-induced differentiation trajectory for neurons

Animals capable of whole-body regeneration can replace any missing cell type and regenerate fully functional new organs, including new brains, de novo. The regeneration of a new brain requires the formation of diverse neural cell types and their assembly into an organized structure with correctly wired circuits. Recent work in various regenerative animals has revealed transcriptional programs required for the differentiation of distinct neural subpopulations, however, how these transcriptional programs are initiated in response to injury remains unknown. Here, we focused on the highly regenerative acoel worm, Hofstenia miamia, to study wound-induced transcriptional regulatory events that lead to the production of neurons and subsequently a functional brain. Footprinting analysis using chromatin accessibility data on a chromosome-scale genome assembly revealed that binding sites for the Nuclear Factor Y (NFY) transcription factor complex were significantly bound during regeneration, showing a dynamic increase in binding within one hour upon amputation specifically in tail fragments, which will regenerate a new brain. Strikingly, NFY targets were highly enriched for genes with neuronal function. Single-cell transcriptome analysis combined with functional studies identified soxC+ stem cells as a putative progenitor population for multiple neural subtypes. Further, we found that wound-induced soxC expression is likely under direct transcriptional control by NFY, uncovering a mechanism for the initiation of a neural differentiation pathway by early wound-induced binding of a transcriptional regulator.

The complex hexaploid oil-Camellia genome traces back its phylogenomic history and multi-omics analysis of Camellia oil biosynthesis.

Oil-Camellia (Camellia oleifera), belonging to the Theaceae family Camellia, is an important woody edible oil tree species. The Camellia oil in its mature seed kernels, mainly consists of more than 90% unsaturated fatty acids, tea polyphenols, flavonoids, squalene and other active substances, which is one of the best quality edible vegetable oils in the world. However, genetic research and molecular breeding on oil-Camellia are challenging due to its complex genetic background. Here, we successfully report a chromosome-scale genome assembly for a hexaploid oil-Camellia cultivar Changlin40. This assembly contains 8.80 Gb genomic sequences with scaffold N50 of 180.0 Mb and 45 pseudochromosomes comprising 15 homologous groups with three members each, which contain 135 868 genes with an average length of 3936 bp. Referring to the diploid genome, intragenomic and intergenomic comparisons of synteny indicate homologous chromosomal similarity and changes. Moreover, comparative and evolutionary analyses reveal three rounds of whole-genome duplication (WGD) events, as well as the possible diversification of hexaploid Changlin40 with diploid occurred approximately 9.06 million years ago (MYA). Furthermore, through the combination of genomics, transcriptomics and metabolomics approaches, a complex regulatory network was constructed and allows to identify potential key structural genes (SAD, FAD2 and FAD3) and transcription factors (AP2 and C2H2) that regulate the metabolism of Camellia oil, especially for unsaturated fatty acids biosynthesis. Overall, the genomic resource generated from this study has great potential to accelerate the research for the molecular biology and genetic improvement of hexaploid oil-Camellia, as well as to understand polyploid genome evolution.

Chromosome-scale genome assemblies of sexually dimorphic male and female Acrossocheilus fasciatus

Acrossocheilus fasciatus is a stream-dwelling fish species of the Barbinae subfamily. It is valued for its colorfully striped appearance and delicious meat. This species is also characterized by apparent sexual dimorphism and toxic ovum. Biology and aquaculture researches of A. fasciatus are hindered by the lack of a high-quality reference genome. Here, we report chromosome-level genome assemblies of the male and female A. fasciatus. The HiFi-only genome assemblies for both female and male individuals were 899.13 Mb (N50 length of 32.58 Mb) and 885.68 Mb (N50 length of 33.06 Mb), respectively. Notably, a substantial proportion of the assembled sequences, accounting for 96.15% and 98.35% for female and male genomes, respectively, were successfully anchored onto 25 chromosomes utilizing Hi-C data. We annotated the female assembly as a reference genome and identified a total of 400.62 Mb (44.56%) repetitive sequences, 27,392 protein-coding genes, and 35,869 ncRNAs. The high-quality male and female reference genomes will provide genomic resources for developing sex-specific molecular markers, inform single-sex breeding, and elucidate genetic mechanisms of sexual dimorphism.

Chromosome-level genome assembly of the cereal cyst nematode Heterodera flipjevi

As an economically important plant parasitic nematode (PPN), Heterodera filipjevi causes great damage on wheat, and now it was widely recorded in many countries. While multiple genomes of PPNs have been published, high-quality genome assembly and annotation on H. filipjevi have yet to be performed. This study presents a chromosome-scale genome assembly and annotation for H. filipjevi, utilizing a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing technologies. The genome consists of 9 pseudo-chromosomes that contain 134.19 Mb of sequence, with a scaffold N50 length of 11.88 Mb. In total, 10,036 genes were annotated, representing 75.20% of the total predicted protein-coding genes. Our study provides the first chromosome-scale genome for H. filipjevi, which is also the inaugural high-quality genome of cereal cyst nematodes (CCNs). It provides a valuable genomic resource for further biological research and pest management of cereal cyst nematodes disease.

An improved genome assembly of Chrysanthemum nankingense reveals expansion and functional diversification of terpene synthase gene family

BackgroundTerpenes are important components of plant aromas, and terpene synthases (TPSs) are the key enzymes driving terpene diversification. In this study, we characterized the volatile terpenes in five different Chrysanthemum nankingense tissues. In addition, genome-wide identification and expression analysis of TPS genes was conducted utilizing an improved chromosome-scale genome assembly and tissue-specific transcriptomes. The biochemical functions of three representative TPSs were also investigated.ResultsWe identified tissue-specific volatile organic compound (VOC) and volatile terpene profiles. The improved Chrysanthemum nankingense genome assembly was high-quality, including a larger assembled size (3.26 Gb) and a better contig N50 length (3.18 Mb) compared to the old version. A total of 140 CnTPS genes were identified, with the majority representing the TPS-a and TPS-b subfamilies. The chromosomal distribution of these TPS genes was uneven, and 26 genes were included in biosynthetic gene clusters. Closely-related Chrysanthemum taxa were also found to contain diverse TPS genes, and the expression profiles of most CnTPSs were tissue-specific. The three investigated CnTPS enzymes exhibited versatile activities, suggesting multifunctionality.ConclusionsWe systematically characterized the structure and diversity of TPS genes across the Chrysanthemum nankingense genome, as well as the potential biochemical functions of representative genes. Our results provide a basis for future studies of terpene biosynthesis in chrysanthemums, as well as for the breeding of improved chrysanthemum varieties.

Haplotype-resolved chromosome-level genome assembly of Ehretia macrophylla

Ehretia macrophylla Wall, known as wild loquat, is an ecologically, economically, and medicinally significant tree species widely grown in China, Japan, Vietnam, and Nepal. In this study, we have successfully generated a haplotype-resolved chromosome-scale genome assembly of E. macrophylla by integrating PacBio HiFi long-reads, Illumina short-reads, and Hi-C data. The genome assembly consists of two haplotypes, with sizes of 1.82 Gb and 1.58 Gb respectively, and contig N50 lengths of 28.11 Mb and 21.57 Mb correspondingly. Additionally, 99.41% of the assembly was successfully anchored into 40 pseudo-chromosomes. We predicted 58,886 protein-coding genes, of which 99.60% were functionally annotated from databases. We furthermore detected 2.65 Gb repeat sequences, 659,290 rRNAs, 4,931 tRNAs and 4,688 other ncRNAs. The high-quality assembly of the genome offers a solid basis for furthering the fields of molecular breeding and functional genomics of E. macrophylla.

A chromosome-scale genome assembly of the nipa palm hispid beetle Octodonta nipae

Nipa palm hispid beetle (Octodonta nipae) is an insect species that is native to Malaysia but has spread to southern China and beyond, seriously threatening palm production. A lack of high-quality genome resources has hindered understanding of the insect’s invasive characteristics and ecological adaptations. Here, we combined Illumina short read, PacBio long-read, and high-throughput chromosome conformation capture (Hi-C) sequencing technologies to generate a high-quality, chromosome-scale genome assembly of nipa palm hispid beetle. The genome assembly was 1.31 Gb in size, consisting of nine chromosomes. The contig and scaffold N50 values were 1.022 Mb and 148.6 Mb, respectively. The genome assembly completeness was estimated at 99.1%. Annotation revealed 16,305 protein-coding genes and 62.16% repeat sequences. This high-quality genome assembly is a valuable resource that will contribute to understanding of the genetic factors underlying the invasive characteristics of nipa palm hispid beetle, ultimately promoting development of efficient control policies.

Chromosome-Scale Genome Assembly and Characterization of Top-Quality Japanese Green Tea Cultivar 'Seimei'.

Japanese green tea, an essential beverage in Japanese culture, is characterized by the initial steaming of freshly harvested leaves during production. This process efficiently inactivates endogenous enzymes such as polyphenol oxidases, resulting in the production of sencha, gyokuro and matcha that preserves the vibrant green color of young leaves. Although genome sequences of several tea cultivars and germplasms have been published, no reference genome sequences are available for Japanese green tea cultivars. Here, we constructed a reference genome sequence of the cultivar 'Seimei', which is used to produce high-quality Japanese green tea. Using the PacBio HiFi and Hi-C technologies for chromosome-scale genome assembly, we obtained 15 chromosome sequences with a total genome size of 3.1 Gb and an N50 of 214.9 Mb. By analyzing the genomic diversity of 23 Japanese tea cultivars and lines, including the leading green tea cultivars 'Yabukita' and 'Saemidori', it was revealed that several candidate genes could be related to the characteristics of Japanese green tea. The reference genome of 'Seimei' and information on genomic diversity of Japanese green tea cultivars should provide crucial information for effective breeding of such cultivars in the future.

A fully haplotype-resolved and nearly gap-free genome assembly of wheat stripe rust fungus

Stripe rust fungus Puccinia striiformis f. sp. tritici (Pst) is a destructive pathogen of wheat worldwide. Pst has a macrocyclic-heteroecious lifecycle, in which one-celled urediniospores are dikaryotic, each nucleus containing one haploid genome. We successfully generated the first fully haplotype-resolved and nearly gap-free chromosome-scale genome assembly of Pst by combining PacBio HiFi sequencing and trio-binning strategy. The genome size of the two haploid assemblies was 75.59 Mb and 75.91 Mb with contig N50 of 4.17 Mb and 4.60 Mb, and both had 18 pseudochromosomes. The high consensus quality values of 55.57 and 59.02 for both haplotypes confirmed the correctness of the assembly. Of the total 18 chromosomes, 15 and 16 were gapless while there were only five and two gaps for the remaining chromosomes of the two haplotypes, respectively. In total, 15,046 and 15,050 protein-coding genes were predicted for the two haplotypes, and the complete BUSCO scores achieved 97.7% and 97.9%, respectively. The genome will lay the foundation for further research on genetic variations and the evolution of rust fungi.

An integrated multi-omics approach reveals polymethoxylated flavonoid biosynthesis in Citrus reticulata cv. Chachiensis

Citrus reticulata cv. Chachiensis (CRC) is an important medicinal plant, its dried mature peels named “Guangchenpi”, has been used as a traditional Chinese medicine to treat cough, indigestion, and lung diseases for several hundred years. However, the biosynthesis of the crucial natural products polymethoxylated flavonoids (PMFs) in CRC remains unclear. Here, we report a chromosome-scale genome assembly of CRC with the size of 314.96 Mb and a contig N50 of 16.22 Mb. Using multi-omics resources, we discover a putative caffeic acid O-methyltransferase (CcOMT1) that can transfer a methyl group to the 3-hydroxyl of natsudaidain to form 3,5,6,7,8,3’,4’-heptamethoxyflavone (HPMF). Based on transient overexpression and virus-induced gene silencing experiments, we propose that CcOMT1 is a candidate enzyme in HPMF biosynthesis. In addition, a potential gene regulatory network associated with PMF biosynthesis is identified. This study provides insights into PMF biosynthesis and may assist future research on mining genes for the biosynthesis of plant-based medicines.

A chromosome-scale assembly of Brassica carinata (BBCC) accession HC20 containing resistance to multiple pathogens and an early generation assessment of introgressions into B.juncea (AABB).

Brassica carinata (BBCC) commonly referred to as Ethiopian mustard is a natural allotetraploid containing the genomes of Brassica nigra (BB) and Brassica oleracea (CC). It is an oilseed crop endemic to the northeastern regions of Africa. Although it is under limited cultivation, B. carinata is valuable as it is resistant/highly tolerant to most of the pathogens affecting widely cultivated Brassica species of the U's triangle. We report a chromosome-scale genome assembly of B. carinata accession HC20 using long-read Oxford Nanopore sequencing and Bionano optical maps. The assembly has a scaffold N50 of ~39.8 Mb and covers ~1.11 Gb of the genome. We compared the long-read genome assemblies of the U's triangle species and found extensive gene collinearity between the diploids and allopolyploids with no evidence of major gene losses. Therefore, B. juncea (AABB), B. napus (AACC), and B. carinata can be regarded as strict allopolyploids. We cataloged the nucleotide-binding and leucine-rich repeat immune receptor (NLR) repertoire of B.carinata and, identified 465 NLRs, and compared these with the NLRs in the other Brassica species. We investigated the extent and nature of early-generation genomic interactions between the constituent genomes of B. carinata and B. juncea in interspecific crosses between the two species. Besides the expected recombination between the constituent B genomes, extensive homoeologous exchanges were observed between the A and C genomes. Interspecific crosses, therefore, can be used for transferring disease resistance from B. carinata to B. juncea and broadening the genetic base of the two allotetraploid species.

A high-quality chromosome-scale genome assembly of blood orange, an important pigmented sweet orange variety

Blood orange (BO) is a rare red-fleshed sweet orange (SWO) with a high anthocyanin content and is associated with numerous health-related benefits. Here, we reported a high-quality chromosome-scale genome assembly for Neixiu (NX) BO, reaching 336.63 Mb in length with contig and scaffold N50 values of 30.6 Mb. Furthermore, 96% of the assembled sequences were successfully anchored to 9 pseudo-chromosomes. The genome assembly also revealed the presence of 37.87% transposon elements and 7.64% tandem repeats, and the annotation of 30,395 protein-coding genes. A high level of genome synteny was observed between BO and SWO, further supporting their genetic similarity. The speciation event that gave rise to the Citrus species predated the duplication event found within them. The genome-wide variation between NX and SWO was also compared. This first high-quality BO genome will serve as a fundamental basis for future studies on functional genomics and genome evolution.

Chromosome-scale genome assembly clarifies the mechanism of flooding tolerance and evolutionary history of Myricaria laxiflora

Myricaria laxiflora, belonging to the Tamaricaceae family, is a valuable plant resource naturally occurring along the riverbanks of the Yangtze River Basin. In this study, we constructed a chromosome-level M. laxiflora genome which was 1.34 Gb in size, with 12 chromosomes and a contig N50 of 11.69 Mb. In total, 23,719 protein-coding genes were discovered across the M. laxiflora genome. Phylogenetic investigation indicated that Tamaricaceae and Polygonaceae diverged between 65 and 81 million years ago (Mya). M. laxiflora underwent one whole genome duplication (WGD) between 19.74 and 24.60 Mya, and experienced a long terminal repeat (LTR) insertion around 0.15 Mya. These genomic changes may be related to environmental upheavals (such as the uplift of the Qinghai-Tibet Plateau) experienced by M. laxiflora during its evolution, and also may have provided a foundation for its evolution and adaptation. To analyze the flooding adaption mechanism of M. laxiflora, we compared morphological, physiological, and transcriptomic changes under control and submerged conditions. Overall, M. laxiflora was found to adopt a "quiescence" strategy when completely submerged. Specifically, in response to submersion, photosynthesis weakens, growth stagnates, glycolysis is activated, and the antioxidant enzyme system is enhanced. The “ethylene response pathway” was ascertained as a likely regulator of this process, and a MlERF-VII gene (Mla11G026380) may be the key gene regulating this process. Our findings shed light on the underlying mechanisms behind the flooding resistance displayed by M. laxiflora, while also establishing a basis for the conservation of its genetic resources.

Identification of a viral gene essential for the genome replication of a domesticated endogenous virus in ichneumonid parasitoid wasps.

Thousands of endoparasitoid wasp species in the families Braconidae and Ichneumonidae harbor "domesticated endogenous viruses" (DEVs) in their genomes. This study focuses on ichneumonid DEVs, named ichnoviruses (IVs). Large quantities of DNA-containing IV virions are produced in ovary calyx cells during the pupal and adult stages of female wasps. Females parasitize host insects by injecting eggs and virions into the body cavity. After injection, virions rapidly infect host cells which is followed by expression of IV genes that promote the successful development of wasp offspring. IV genomes consist of two components: proviral segment loci that serve as templates for circular dsDNAs that are packaged into capsids, and genes from an ancestral virus that produce virions. In this study, we generated a chromosome-scale genome assembly for Hyposoter didymator that harbors H. didymator ichnovirus (HdIV). We identified a total of 67 HdIV loci that are amplified in calyx cells during the wasp pupal stage. We then focused on an HdIV gene, U16, which is transcribed in calyx cells during the initial stages of replication. Sequence analysis indicated that U16 contains a conserved domain in primases from select other viruses. Knockdown of U16 by RNA interference inhibited virion morphogenesis in calyx cells. Genome-wide analysis indicated U16 knockdown also inhibited amplification of HdIV loci in calyx cells. Altogether, our results identified several previously unknown HdIV loci, demonstrated that all HdIV loci are amplified in calyx cells during the pupal stage, and showed that U16 is required for amplification and virion morphogenesis.

Chromosome-level genome assembly of the predatory stink bug Arma custos

The stink bug Arma custos (Hemiptera: Pentatomidae) is a predatory enemy successfully used for biocontrol of lepidopteran and coleopteran pests in notorious invasive species. In this study, a high-quality chromosome-scale genome assembly of A. custos was achieved through a combination of Illumina sequencing, PacBio HiFi sequencing, and Hi-C scaffolding techniques. The final assembled genome was 969.02 Mb in size, with 935.94 Mb anchored to seven chromosomes, and a scaffold N50 length of 135.75 Mb. This genome comprised 52.78% repetitive elements. The detected complete BUSCO score was 99.34%, indicating its completeness. A total of 13,708 protein-coding genes were predicted in the genome, and 13219 of them were annotated. This genome provides an invaluable resource for further research on various aspects of predatory bugs, such as biology, genetics, and functional genomics.

The genomic and enzymatic basis for iridoid biosynthesis in cat thyme (Teucrium marum).

Iridoids are non-canonical monoterpenoids produced by both insects and plants. An example is the cat-attracting and insect-repelling volatile iridoid nepetalactone, produced by Nepeta sp. (catmint) and aphids. Recently, both nepetalactone biosynthetic pathways were elucidated, showing a remarkable convergent evolution. The iridoid, dolichodial, produced by Teucrium marum (cat thyme) and multiple insect species, has highly similar properties to nepetalactone but its biosynthetic origin remains unknown. We set out to determine the genomic, enzymatic, and evolutionary basis of iridoid biosynthesis in T. marum. First, we generated a de novo chromosome-scale genome assembly for T. marum using Oxford Nanopore Technologies long reads and proximity-by-ligation Hi-C reads. The 610.3 Mb assembly spans 15 pseudomolecules with a 32.9 Mb N50 scaffold size. This enabled identification of iridoid biosynthetic genes, whose roles were verified via activity assays. Phylogenomic analysis revealed that the evolutionary history of T. marum iridoid synthase, the iridoid scaffold-forming enzyme, is not orthologous to typical iridoid synthases but is derived from its conserved paralog. We discovered an enzymatic route from nepetalactol to diverse iridoids through the coupled activity of an iridoid oxidase cytochrome P450 and acetyltransferases, via an inferred acylated intermediate. This work provides a genomic resource for specialized metabolite research in mints and demonstration of the role of acetylation in T. marum iridoid diversity. This work will enable future biocatalytic or biosynthetic production of potent insect repellents, as well as comparative studies into iridoid biosynthesis in insects.

Genome assembly and multi-omic analyses reveal the mechanisms underlying flower color formation in Torenia fournieri.

Torenia fournieri Lind. is an ornamental plant that is popular for its numerous flowers and variety of colors. However, its genomic evolutionary history and the genetic and metabolic bases of flower color formation remain poorly understood. Here, we report the first T. fournieri reference genome, which was resolved to the chromosome scale and was 164.4 Mb in size. Phylogenetic analyses clarified relationships with other plant species, and a comparative genomic analysis indicated that the shared ancestor of T. fournieri and Antirrhinum majus underwent a whole genome duplication event. Joint transcriptomic and metabolomic analyses identified many metabolites related to pelargonidin, peonidin, and naringenin production in rose (TfR)-colored flowers. Samples with blue (TfB) and deep blue (TfD) colors contained numerous derivatives of petunidin, cyanidin, quercetin, and malvidin; differences in the abundances of these metabolites and expression levels of the associated genes were hypothesized to be responsible for variety-specific differences in flower color. Furthermore, the genes encoding flavonoid 3-hydroxylase, anthocyanin synthase, and anthocyanin reductase were differentially expressed between flowers of different colors. Overall, we successfully identified key genes and metabolites involved in T. fournieri flower color formation. The data provided by the chromosome-scale genome assembly establish a basis for understanding the differentiation of this species and will facilitate future genetic studies and genomic-assisted breeding.

Chromosome-scale genome assembly of Poa trivialis and population genomics reveal widespread gene flow in a cool-season grass seed production system.

Poa trivialis (L.) is a cool-season grass species found in various environments worldwide. In addition to being a desired turfgrass species, it is a common weed of agricultural systems and natural areas. As a weed, it is an important contaminant of commercial cool-season grass seed lots, resulting in widespread gene flow facilitated by human activities and causing significant economic losses to farmers. To better understand and manage infestations, we assembled and annotated a haploid genome of P. trivialis and studied troublesome field populations from Oregon, the largest cool-season grass seed producing region in the United States. The genome assembly resulted in 1.35 Gb of DNA sequence distributed among seven chromosome-scale scaffolds, revealing a high content of transposable elements, conserved synteny with Poa annua, and a close relationship with other C3 grasses. A reduced-representation sequencing analysis of field populations revealed limited genetic diversity and suggested potential gene flow and human-assisted dispersal in the region. The genetic resources and insights into P. trivialis provided by this study will improve weed management strategies and enable the development of molecular detection tests for contaminated seed lots to limit seed-mediated gene flow. These resources should also be beneficial for turfgrass breeders seeking to improve desirable traits of commercial P. trivialis varieties and help to guide breeding efforts in other crops to enhance the resiliency of agricultural ecosystems under climate change. Significance Statement: The chromosome-scale assembly of Poa trivialis and population genomic analyses provide crucial insights into the gene flow of weedy populations in agricultural systems and contribute a valuable genomic resource for the plant science community.

Genome assembly of autotetraploid Actinidia arguta highlights adaptive evolution and enables dissection of important economic traits

Actinidia arguta, the most widely distributed Actinidia species and the second cultivated species in the genus, can be distinguished from the currently cultivated Actinidia chinensis on the basis of its small and smooth fruit, rapid softening, and excellent cold tolerance. Adaptive evolution of tetraploid Actinidia species and the genetic basis of their important agronomic traits are still unclear. Here, we generated a chromosome-scale genome assembly of an autotetraploid male A. arguta accession. The genome assembly was 2.77 Gb in length with a contig N50 of 9.97 Mb and was anchored onto 116 pseudo-chromosomes. Resequencing and clustering of 101 geographically representative accessions showed that they could be divided into two geographic groups, Southern and Northern, which first diverged 12.9 million years ago. A. arguta underwent two prominent expansions and one demographic bottleneck from the mid-Pleistocene climate transition to the late Pleistocene. Population genomics studies using paleoclimate data enabled us to discern the evolution of the species’ adaptation to different historical environments. Three genes (AaCEL1, AaPME1, and AaDOF1) related to flesh softening were identified by multi-omics analysis, and their ability to accelerate flesh softening was verified through transient expression assays. A set of genes that characteristically regulate sexual dimorphism located on the sex chromosome (Chr3) or autosomal chromosomes showed biased expression during stamen or carpel development. This chromosome-level assembly of the autotetraploid A. arguta genome and the genes related to important agronomic traits will facilitate future functional genomics research and improvement of A. arguta.