Haplotype-resolved Genome Research Articles

Gamete Binning to Achieve Haplotype-Resolved Genome Assembly.

Haplotype-resolved genome assemblies remain a challenge in practice. Here, we provide a step-by-step guide on gamete binning, a method to generate haplotype-resolved genome assemblies for diploid species. The protocol starts by phasing heterozygous variants to individual haplotypes of specific chromosomes using the genome information of individual haploid gametes of the focal individual. Using phased variants, the whole-genome sequencing reads from the diploid genome can be genotyped and assigned into groups, which represent the individual haplotypes of each of the chromosomes. Finally, haplotype-specific chromosomes can be assembled independently using standard assembly tools. First applications of gamete binning revealed a haplotyping accuracy over 99%, which outperformed sequence-only or Hi-C-based haplotype-resolved genome assemblies.Availability: github.com/schneebergerlab/GameteBinning_prac .

Population genomics of Puccinia graminis f.sp. tritici highlights the role of admixture in the origin of virulent wheat rust races

Puccinia graminis f.sp. tritici (Pgt) causes stem rust disease in wheat that can result in severe yield losses. The factors driving the evolution of its virulence and adaptation remain poorly characterized. We utilize long-read sequencing to develop a haplotype-resolved genome assembly of a U.S. isolate of Pgt. Using Pgt haplotypes as a reference, we characterize the structural variants (SVs) and single nucleotide polymorphisms in a diverse panel of isolates. SVs impact the repertoire of predicted effectors, secreted proteins involved in host-pathogen interaction, and show evidence of purifying selection. By analyzing global and local genomic ancestry we demonstrate that the origin of 8 out of 12 Pgt clades is linked with either somatic hybridization or sexual recombination between the diverged donor populations. Our study shows that SVs and admixture events appear to play an important role in broadening Pgt virulence and the origin of highly virulent races, creating a resource for studying the evolution of Pgt virulence and preventing future epidemic outbreaks.

Haplotype-resolved genome assembly of Bletilla striata (Thunb.) Reichb.f. to elucidate medicinal value.

Bletilla striata, commonly known as baiji, is a species used in traditional Chinese medicine; it is highly regarded for its medicinal applications and therefore has high economic value. Here, we report a high-quality haplotype-resolved genome of B. striata, haplotype A (2.37 Gb, with a scaffold N50 of 146.39 Mb and a contig N50 of 1.65 Mb) and haplotype B (2.43 Gb, with a scaffold N50 of 150.22 Mb and a contig N50 of 1.66 Mb), assembled from high-fidelity (HiFi) reads and chromosome conformation capture (Hi-C) reads. We find evidence that B. striata has undergone two whole-genome duplication (WGD) events: an ancient WGD event shared by most monocots and a recent WGD event unique to all orchids. We also reconstructed the ancestral orchid karyotype (AOK) of 18 ancient chromosomes and the evolutionary trajectories of 16 modern B. striata chromosomes. Comparative genomic analysis suggests that the expanded gene families of B. striata might play important roles in secondary metabolite biosynthesis and environmental adaptation. By combining genomic and transcriptomic data, we identified the 10 core members from nine gene families that were probably involved in B. striata polysaccharide (BSP) biosynthesis. Based on virus-induced gene silencing (VIGS) and yeast two-hybrid experiments, we present an MYB transcription factor (TF), BsMYB2, that can regulate BSP biosynthesis by directly interacting with eight key BSP-related genes: sacA1, HK1, scrK1, scrK2, GPI1, manA1, GMPP1 and UGP2_1. Our study will enhance the understanding of orchid evolution and accelerate the molecular-assisted breeding of B. striata for improving traits of medicinal value.

PtoNF-YC9-SRMT-PtoRD26 module regulates the high saline tolerance of a triploid poplar

BackgroundSensing and responding to stresses determine the tolerance of plants to adverse environments. The triploid Chinese white poplar is widely cultivated in North China because of its adaptation to a wide range of habitats including highly saline ones. However, its triploid genome complicates any detailed investigation of the molecular mechanisms underlying its adaptations.ResultsWe report a haplotype-resolved genome of this triploid poplar and characterize, using reverse genetics and biochemical approaches, a MYB gene, SALT RESPONSIVE MYB TRANSCRIPTION FACTOR (SRMT), which combines NUCLEAR FACTOR Y SUBUNIT C 9 (PtoNF-YC9) and RESPONSIVE TO DESICCATION 26 (PtoRD26), to regulate an ABA-dependent salt-stress response signaling. We reveal that the salt-inducible PtoRD26 is dependent on ABA signaling. We demonstrate that ABA or salt drives PtoNF-YC9 shuttling into the nucleus where it interacts with SRMT, resulting in the rapid expression of PtoRD26 which in turn directly regulates SRMT. This positive feedback loop of SRMT-PtoRD26 can rapidly amplify salt-stress signaling. Interference with either component of this regulatory module reduces the salt tolerance of this triploid poplar.ConclusionOur findings reveal a novel ABA-dependent salt-responsive mechanism, which is mediated by the PtoNF-YC9-SRMT-PtoRD26 module that confers salt tolerance to this triploid poplar. These genes may therefore also serve as potential and important modification targets in breeding programs.

Thousands of human mutation clusters are explained by short-range template switching.

Variation within human genomes is unevenly distributed, and variants show spatial clustering. DNA replication-related template switching is a poorly known mutational mechanism capable of causing major chromosomal rearrangements as well as creating short inverted sequence copies that appear as local mutation clusters in sequence comparisons. In this study, haplotype-resolved genome assemblies representing 25 human populations and multinucleotide variants aggregated from 140,000 human sequencing experiments were reanalyzed. Local template switching could explain thousands of complex mutation clusters across the human genome, the loci segregating within and between populations. During the study, computational tools were developed for identification of template switch events using both short-read sequencing data and genotype data, and for genotyping candidate loci using short-read data. The characteristics of template-switch mutations complicate their detection, and widely used analysis pipelines for short-read sequencing data, normally capable of identifying single nucleotide changes, were found to miss template-switch mutations of tens of base pairs, potentially invalidating medical genetic studies searching for a causative allele behind genetic diseases. Combined with the massive sequencing data now available for humans, the novel tools described here enable building catalogs of affected loci and studying the cellular mechanisms behind template switching in both healthy organisms and disease.

Genome architecture and tetrasomic inheritance of autotetraploid potato

Potato (Solanum tuberosum) is the most consumed non-cereal food crop. Most commercial potato cultivars are autotetraploids with highly heterozygous genomes, severely hampering genetic analyses and improvement. By leveraging the state-of-the-art sequencing technologies and polyploid graph binning, we achieved a chromosome-scale, haplotype-resolved genome assembly of a cultivated potato, Cooperation-88 (C88). Intra-haplotype comparative analyses revealed extensive sequence and expression differences in this tetraploid genome. We identified haplotype-specific pericentromeres on chromosomes, suggesting a distinct evolutionary trajectory of potato homologous centromeres. Furthermore, we detected double reduction events that are unevenly distributed on haplotypes in 1021 of 1034 selfing progeny, a feature of autopolyploid inheritance. By distinguishing maternal and paternal haplotype sets in C88, we simulated the origin of heterosis in cultivated tetraploid with a survey of 3110 tetra-allelic loci with deleterious mutations, which were masked in the heterozygous condition by two parents. This study provides insights into the genomic architecture of autopolyploids and will guide their breeding.

Chromosome-level and haplotype-resolved genome provides insight into the tetraploid hybrid origin of patchouli

Patchouli (Pogostemon cablin (Blanco) Benth.), a member of the Lamiaceae family, is an important aromatic plant that has been widely used in medicine and perfumery. Here, we report a 1.94 Gb chromosome-scale assembly of the patchouli genome (contig N50 = 7.97 Mb). The gene annotation reveals that tandem duplication of sesquiterpene biosynthetic genes may be a major contributor to the biosynthesis of patchouli bioactivity components. We further phase the genome into two distinct subgenomes (A and B), and identify a chromosome substitution event that have occurred between them. Further investigations show that a burst of universal LTR-RTs in the A subgenome lead to the divergence between two subgenomes. However, no significant subgenome dominance is detected. Finally, we track the evolutionary scenario of patchouli including whole genome tetraploidization, subgenome divergency, hybridization, and chromosome substitution, which are the key forces to determine the complexity of patchouli genome. Our work sheds light on the evolutionary history of patchouli and offers unprecedented genomic resources for fundamental patchouli research and elite germplasm development.

Chromosome-Scale, Haplotype-Resolved Genome Assembly of Non-Sex-Reversal Females of Swamp Eel Using High-Fidelity Long Reads and Hi-C Data.

The Asian swamp eel (Monopterus albus) is an excellent model species for studying sex change and chromosome evolution. M. albus is also widely reared in East Asia and South-East Asia because of its great nutritional value. The low fecundity of this species (about 300 eggs per fish) greatly hinders fries production and breeding programs. Interestingly, about 3–5% of the eels could remain as females for 3 years and lay more than 3,000 eggs per fish, which are referred to as non-sex-reversal (NSR) females. Here, we presented a new chromosome-level genome assembly of such NSR females using Illumina, HiFi, and Hi-C sequencing technologies. The new assembly (Mal.V2_NSR) is 838.39 Mb in length, and the N50 of the contigs is 49.8 Mb. Compared with the previous assembly obtained using the continuous long-read sequencing technology (Mal.V1_CLR), we found a remarkable increase of continuity in the new assembly Mal.V2_NSR with a 20-times longer contig N50. Chromosomes 2 and 12 were assembled into a single contig, respectively. Meanwhile, two highly contiguous haplotype assemblies were also obtained, with contig N50 being 14.54 and 12.13 Mb, respectively. BUSCO and Merqury analyses indicate completeness and high accuracy of these three assemblies. A comparative genomic analysis revealed substantial structural variations (SVs) between Mal.V2_NSR and Mal.V1_CLR and two phased haplotype assemblies, as well as whole chromosome fusion events when compared with the zig-zag eel. Additionally, our newly obtained assembly provides a genomic view of sex-related genes and a complete landscape of the MHC genes. Therefore, these high-quality genome assemblies would provide great help for future breeding works of the swamp eel, and it is a valuable new reference for genetic and genomic studies of this species.

TT-Mars: structural variants assessment based on haplotype-resolved assemblies

Variant benchmarking is often performed by comparing a test callset to a gold standard set of variants. In repetitive regions of the genome, it may be difficult to establish what is the truth for a call, for example, when different alignment scoring metrics provide equally supported but different variant calls on the same data. Here, we provide an alternative approach, TT-Mars, that takes advantage of the recent production of high-quality haplotype-resolved genome assemblies by providing false discovery rates for variant calls based on how well their call reflects the content of the assembly, rather than comparing calls themselves.Supplementary InformationThe online version contains supplementary material available at (10.1186/s13059-022-02666-2).

A High-Quality Haplotype-Resolved Genome of Common Bermudagrass (Cynodon dactylon L.) Provides Insights Into Polyploid Genome Stability and Prostrate Growth

Common bermudagrass (Cynodon dactylon L.) is an important perennial warm-season turfgrass species with great economic value. However, the reference genome is still deficient in C. dactylon, which severely impedes basic studies and breeding studies. In this study, a high-quality haplotype-resolved genome of C. dactylon cultivar Yangjiang was successfully assembled using a combination of multiple sequencing strategies. The assembled genome is approximately 1.01 Gb in size and is comprised of 36 pseudo chromosomes belonging to four haplotypes. In total, 76,879 protein-coding genes and 529,092 repeat sequences were annotated in the assembled genome. Evolution analysis indicated that C. dactylon underwent two rounds of whole-genome duplication events, whereas syntenic and transcriptome analysis revealed that global subgenome dominance was absent among the four haplotypes. Genome-wide gene family analyses further indicated that homologous recombination-regulating genes and tiller-angle-regulating genes all showed an adaptive evolution in C. dactylon, providing insights into genome-scale regulation of polyploid genome stability and prostrate growth. These results not only facilitate a better understanding of the complex genome composition and unique plant architectural characteristics of common bermudagrass, but also offer a valuable resource for comparative genome analyses of turfgrasses and other plant species.

Gene losses in the common vampire bat illuminate molecular adaptations to blood feeding.

Vampire bats are the only mammals that feed exclusively on blood. To uncover genomic changes associated with this dietary adaptation, we generated a haplotype-resolved genome of the common vampire bat and screened 27 bat species for genes that were specifically lost in the vampire bat lineage. We found previously unknown gene losses that relate to reduced insulin secretion (FFAR1 and SLC30A8), limited glycogen stores (PPP1R3E), and a unique gastric physiology (CTSE). Other gene losses likely reflect the biased nutrient composition (ERN2 and CTRL) and distinct pathogen diversity of blood (RNASE7) and predict the complete lack of cone-based vision in these strictly nocturnal bats (PDE6H and PDE6C). Notably, REP15 loss likely helped vampire bats adapt to high dietary iron levels by enhancing iron excretion, and the loss of CYP39A1 could have contributed to their exceptional cognitive abilities. These findings enhance our understanding of vampire bat biology and the genomic underpinnings of adaptations to blood feeding.

Haplotype-resolved assembly of diploid genomes without parental data.

Routine haplotype-resolved genome assembly from single samples remains an unresolved problem. Here we describe an algorithm that combines PacBio HiFi reads and Hi-C chromatin interaction data to produce a haplotype-resolved assembly without the sequencing of parents. Applied to human and other vertebrate samples, our algorithm consistently outperforms existing single-sample assembly pipelines and generates assemblies of similar quality to the best pedigree-based assemblies.

Gene prediction for leaf margin phenotype and fruit flesh color in pineapple (Ananas comosus) using haplotype-resolved genome sequencing.

Pineapple (Ananas comosus (L.) Merr.) is one of the most economically important tropical fruit species. The major aim of the breeding programs in several countries, including Japan, is quality improvement, mainly for the fresh market. ‘Yugafu’, a Japanese cultivar with distinctive pipe-type leaf margin phenotype and white flesh color, is popular for fresh consumption. Therefore, genome sequencing of ‘Yugafu’ is expected to assist pineapple breeding. Here, we developed a haplotype-resolved assembly for the heterozygous genome of ‘Yugafu’ using long-read sequencing technology and obtained a pair of 25 pseudomolecule sequences inherited from the parental accessions ‘Cream pineapple’ and ‘HI101’. The causative genes for leaf margin and fruit flesh color were identified. Fine mapping revealed a 162-kb region on CLG23 for the leaf margin phenotype. In this region, 20 kb of inverted repeat was specifically observed in the ‘Cream pineapple’ derived allele, and the WUSCHEL-related homeobox 3 (AcWOX3) gene was predicted as the key gene for leaf margin morphogenesis. Dominantly repressed AcWOX3 via RNAi was suggested to be the cause of the pipe-type leaf margin phenotype. Quantitative trait locus (QTL) analysis revealed that the terminal region of CLG08 contributed to white flesh and low carotenoid content. Carotenoid cleaved dioxygenase 4 (AcCCD4), a key gene for carotenoid degradation underlying this QTL, was predicted as the key gene for white flesh color through expression analysis. These findings could assist in modern pineapple breeding and facilitate marker-assisted selection for important traits.

Haplotype-resolved Chinese male genome assembly based on high-fidelity sequencing

The advantages of both the length and accuracy of high-fidelity (HiFi) reads enable chromosome-scale haplotype-resolved genome assembly. In this study, we sequenced a cell line named HJ, established from a Chinese Han male individual by using HiFi and Hi-C. We assembled two high-quality haplotypes of the HJ genome (haplotype 1 (H1): 3.1 Gb, haplotype 2 (H2): 2.9 Gb). The continuity (H1: contig N50 = 28.2 Mb, H2: contig N50 = 25.9 Mb) and completeness (BUSCO: H1 = 94.9%, H2 = 93.5%) are substantially better than those of other Chinese genomes, for example, HX1, NH1.0, and YH2.0. By comparing HJ genome with GRCh38, we reported the mutation landscape of HJ and found that 176 and 213 N-gaps were filled in H1 and H2, respectively. In addition, we detected 12.9 Mb and 13.4 Mb novel sequences containing 246 and 135 protein-coding genes in H1 and H2, respectively. Our results demonstrate the advantages of HiFi reads in haplotype-resolved genome assembly and provide two high-quality haplotypes of a potential Chinese genome as a reference for the Chinese Han population.

Chromosome-scale and haplotype-resolved genome assembly of a tetraploid potato cultivar

Potato is the most widely produced tuber crop worldwide. However, reconstructing the four haplotypes of its autotetraploid genome remained an unsolved challenge. Here, we report the 3.1 Gb haplotype-resolved (at 99.6% precision), chromosome-scale assembly of the potato cultivar ‘Otava’ based on high-quality long reads, single-cell sequencing of 717 pollen genomes and Hi-C data. Unexpectedly, ~50% of the genome was identical-by-descent due to recent inbreeding, which was contrasted by highly abundant structural rearrangements involving ~20% of the genome. Among 38,214 genes, only 54% were present in all four haplotypes with an average of 3.2 copies per gene. Taking the leaf transcriptome as an example, 11% of the genes were differently expressed in at least one haplotype, where 25% of them were likely regulated through allele-specific DNA methylation. Our work sheds light on the recent breeding history of potato, the functional organization of its tetraploid genome and has the potential to strengthen the future of genomics-assisted breeding.

The power of chromosome-scale, haplotype-resolved genomes

The assembly of chromosome-scale and haplotype-resolved reference genomes is now more easily attainable, largely due to various improvements in both assembly algorithms and long-read sequencing technologies (see recent review by Michael and VanBuren, 2020 Michael T.P. VanBuren R. Building near-complete plant genomes. Curr. Opin. Plant Biol. 2020; 54: 26-33 Google Scholar ). Due to these technological advancements, there has been a shift away from using highly inbred accessions, in some instances double haploids, to instead assembling the genomes of important reference genotypes. For example, a haplotype-resolved genome was recently assembled for an important cultivated garden strawberry (Fragaria x ananassa) cultivar Royal Royce, which is not only highly heterozygous but also octoploid (8×) ( Hardigan et al., 2021 Hardigan M.A. Feldmann M.J. Pincot D.D.A. Famula R.A. Vachev M.V. Madera M.A. Zerbe P. Mars K. Peluso P. Rank D. et al. Blueprint for phasing and assembling the genomes of heterozygous polyploids: application to the octoploid genome of strawberry. Preprint at bioRxiv. 2021; https://doi.org/10.1101/2021.11.03.467115 Google Scholar ). Similarly, haplotype-phased genomes have been assembled for other important crops, including cassava (Manihot esculenta) ( Hu et al., 2021 Hu W. Ji C. Shi H. Liang Z. Ding Z. Ye J. Ou W. Zhou G. Tie W. Yan Y. et al. Allele-defined genome reveals biallelic differentiation during cassava evolution. Mol. Plant. 2021; 14: 851-854 Google Scholar ) and potato (Solanum tuberosum) ( Hoopes et al., 2022 Hoopes G. Meng X. Hamilton J.P. Achakkagari S.R. de Alves Freitas Guesdes F. Bolger M.E. Coombs J.J. Esselink D. Kaiser N.R. Kodde L. et al. Phased, chromosome-scale genome assemblies of tetraploid potato reveals a complex genome, transcriptome, and predicted proteome landscape underpinning genetic diversity. Mol. Plant. 2022; (S1674-2052, 00003-X)https://doi.org/10.1016/j.molp.2022.01.003 Google Scholar ). This has resulted in the development of new genomic resources that can be directly used in guiding molecular breeding efforts, without worrying about the gene presence-absence variation that exists within a crop or species ( Golicz et al., 2016 Golicz A.A. Batley J. Edwards D. Towards plant pangenomics. Plant Biotechnol. J. 2016; 14: 1099-1105 Google Scholar ; Bayer et al., 2021 Bayer P.E. Scheben A. Golicz A.A. Yuan Y. Faure S. Lee H. Chawla H.S. Anderson R. Bancroft I. Raman H. et al. Modelling of gene loss propensity in the pangenomes of three Brassica species suggests different mechanisms between polyploids and diploids. Plant Biotechnol. J. 2021; 19: 2488-2500 Google Scholar ). An additional benefit of having a reference genome of an important cultivar and/or genotype, versus from a highly inbred line, is that valuable insights that are potentially gained from the analysis of genetic variants present among haplotypes are not lost, allowing for additional insights.

High quality haplotype-resolved genome assemblies of Populus tomentosa Carr., a stabilized interspecific hybrid species widespread in Asia.

Populus has a wide ecogeographical range spanning the Northern Hemisphere, and interspecific hybrids are common. Populus tomentosa Carr. is widely distributed and cultivated in the eastern region of Asia, where it plays multiple important roles in forestry, agriculture, conservation, and urban horticulture. Reference genomes are available for several Populus species, however, our goals were to produce a very high quality de novo chromosome-level genome assembly in P.tomentosa genome that could serve as a reference for evolutionary and ecological studies of hybrid speciation throughout the genus. Here, combining long-read sequencing and Hi-C scaffolding, we present a high-quality, haplotype-resolved genome assembly. The genome size was 740.2Mb, with a contig N50size of 5.47Mb and a scaffold N50size of 46.68Mb, consisting of 38 chromosomes, as expected with the known diploid chromosome number (2n=2x=38). A total of 59,124 protein-coding genes were identified. Phylogenomic analyses revealed that P.tomentosa is comprised of two distinct subgenomes, which we deomonstrate is likely to have resulted from hybridization between Populus adenopoda as the female parent and Populus alba var. pyramidalis as the male parent, with an origin of approximately 3.93Ma. Although highly colinear, significant structural variation was found between the two subgenomes. Our study provides a valuable resource for ecological genetics and forest biotechnology.

Genomic structural variants constrain and facilitate adaptation in natural populations of Theobroma cacao, the chocolate tree

Genomic structural variants (SVs) can play important roles in adaptation and speciation. Yet the overall fitness effects of SVs are poorly understood, partly because accurate population-level identification of SVs requires multiple high-quality genome assemblies. Here, we use 31 chromosome-scale, haplotype-resolved genome assemblies of Theobroma cacao-an outcrossing, long-lived tree species that is the source of chocolate-to investigate the fitness consequences of SVs in natural populations. Among the 31 accessions, we find over 160,000 SVs, which together cover eight times more of the genome than single-nucleotide polymorphisms and short indels (125 versus 15 Mb). Our results indicate that a vast majority of these SVs are deleterious: they segregate at low frequencies and are depleted from functional regions of the genome. We show that SVs influence gene expression, which likely impairs gene function and contributes to the detrimental effects of SVs. We also provide empirical support for a theoretical prediction that SVs, particularly inversions, increase genetic load through the accumulation of deleterious nucleotide variants as a result of suppressed recombination. Despite the overall detrimental effects, we identify individual SVs bearing signatures of local adaptation, several of which are associated with genes differentially expressed between populations. Genes involved in pathogen resistance are strongly enriched among these candidates, highlighting the contribution of SVs to this important local adaptation trait. Beyond revealing empirical evidence for the evolutionary importance of SVs, these 31 de novo assemblies provide a valuable resource for genetic and breeding studies in Tcacao.

Haplotype-resolved genome assembly and allele-specific gene expression in cultivated ginger

Ginger (Zingiber officinale) is one of the most valued spice plants worldwide; it is prized for its culinary and folk medicinal applications and is therefore of high economic and cultural importance. Here, we present a haplotype-resolved, chromosome-scale assembly for diploid ginger anchored to 11 pseudochromosome pairs with a total length of 3.1 Gb. Remarkable structural variation was identified between haplotypes, and two inversions larger than 15 Mb on chromosome 4 may be associated with ginger infertility. We performed a comprehensive, spatiotemporal, genome-wide analysis of allelic expression patterns, revealing that most alleles are coordinately expressed. The alleles that exhibited the largest differences in expression showed closer proximity to transposable elements, greater coding sequence divergence, more relaxed selection pressure, and more transcription factor binding site differences. We also predicted the transcription factors potentially regulating 6-gingerol biosynthesis. Our allele-aware assembly provides a powerful platform for future functional genomics, molecular breeding, and genome editing in ginger.

Haplotype-resolved genome of diploid ginger (Zingiber officinale) and its unique gingerol biosynthetic pathway

Ginger (Zingiber officinale), the type species of Zingiberaceae, is one of the most widespread medicinal plants and spices. Here, we report a high-quality, chromosome-scale reference genome of ginger ‘Zhugen’, a traditionally cultivated ginger in Southwest China used as a fresh vegetable, assembled from PacBio long reads, Illumina short reads, and high-throughput chromosome conformation capture (Hi-C) reads. The ginger genome was phased into two haplotypes, haplotype 1 (1.53 Gb with a contig N50 of 4.68 M) and haplotype 0 (1.51 Gb with a contig N50 of 5.28 M). Homologous ginger chromosomes maintained excellent gene pair collinearity. In 17,226 pairs of allelic genes, 11.9% exhibited differential expression between alleles. Based on the results of ginger genome sequencing, transcriptome analysis, and metabolomic analysis, we proposed a backbone biosynthetic pathway of gingerol analogs, which consists of 12 enzymatic gene families, PAL, C4H, 4CL, CST, C3’H, C3OMT, CCOMT, CSE, PKS, AOR, DHN, and DHT. These analyses also identified the likely transcription factor networks that regulate the synthesis of gingerol analogs. Overall, this study serves as an excellent resource for further research on ginger biology and breeding, lays a foundation for a better understanding of ginger evolution, and presents an intact biosynthetic pathway for species-specific gingerol biosynthesis.