MARS: a motif-based autoregressive model for retrosynthesis prediction.

MOTIVATION: Retrosynthesis is a critical task in drug discovery, aimed at finding a viable pathway for synthesizing a given target molecule. Many existing approaches frame this task as a graph-generating problem. Specifically, these methods first identify the reaction center, and break a targeted molecule accordingly to generate the synthons. Reactants are generated by either adding atoms sequentially to synthon graphs or by directly adding appropriate leaving groups. However, both of these strategies have limitations. Adding atoms results in a long prediction sequence that increases the complexity of generation, while adding leaving groups only considers those in the training set, which leads to poor generalization. RESULTS: In this paper, we propose a novel end-to-end graph generation model for retrosynthesis prediction, which sequentially identifies the reaction center, generates the synthons, and adds motifs to the synthons to generate reactants. Given that chemically meaningful motifs fall between the size of atoms and leaving groups, our model achieves lower prediction complexity than adding atoms and demonstrates superior performance than adding leaving groups. We evaluate our proposed model on a benchmark dataset and show that it significantly outperforms previous state-of-the-art models. Furthermore, we conduct ablation studies to investigate the contribution of each component of our proposed model to the overall performance on benchmark datasets. Experiment results demonstrate the effectiveness of our model in predicting retrosynthesis pathways and suggest its potential as a valuable tool in drug discovery. AVAILABILITY AND IMPLEMENTATION: All code and data are available at https://github.com/szu-ljh2020/MARS.

Genomic Consequences of and Demographic Response to Pervasive Hybridization Over Time in Climate-Sensitive Pikas.

Rare and geographically restricted species may be vulnerable to genetic effects from inbreeding depression in small populations or from genetic swamping through hybridization with common species, but a third possibility is that selective gene flow can restore fitness (genetic rescue). Climate-sensitive pikas (Ochotona spp.) of the Qinghai-Tibetan Plateau (QHTP) and its vicinity have been reduced to residual populations through the movement of climatic zones during the Pleistocene and recent anthropogenic disturbance, whereas the plateau pika (O. curzoniae) remains common. Population-level whole-genome sequencing (n = 142) of six closely related species in the subgenus Ochotona revealed several phases of ancient introgression, lineage replacement, and bidirectional introgression. The strength of gene flow was the greatest from the dominant O. curzoniae to ecologically distinct species in areas peripheral to the QHTP. Genetic analyses were consistent with environmental reconstructions of past population movements. Recurrent periods of introgression throughout the Pleistocene revealed an increase in genetic variation at first but subsequent loss of genetic variation in later phases. Enhanced dispersion of introgressed genomic regions apparently contributed to demographic recovery in three peripheral species that underwent range shifts following climate oscillations on the QHTP, although it failed to drive recovery of northeastern O. dauurica and geographically isolated O. sikimaria. Our findings highlight differences in timescale and environmental background to determine the consequence of hybridization and the unique role of the QHTP in conserving key evolutionary processes of sky island species.

MEANGS: an efficient seed-free tool for de novo assembling animal mitochondrial genome using whole genome NGS data.

Advances in next-generation sequencing (NGS) technologies have led to an exponential increase in the number of whole genome sequences (WGS) in databases. This wealth of WGS data has greatly facilitated the recovery of full mitochondrial genomes (mitogenomes), which are vital for phylogenetic, evolutionary and ecological studies. Unfortunately, most existing software cannot easily assemble mitogenome reference sequences conveniently or efficiently. Therefore, we developed a seed-free de novo assembly tool, MEANGS, which applies the trie-search method to extend contigs from self-discovery seeds and assemble a mitogenome from animal WGS data. We then used data from 16 species with different qualities to compare the performance of MEANGS with three other available programs. MEANGS exhibited the best overall performance since it was the only one that completed all tests, and it assembled full or partial mitogenomes for all of the tested samples while the others failed. Furthermore, MEANGS selects superior assembly sequences and annotates protein-coding genes. Thus, MEANGS can be one of the most efficient software for generating high-quality mitogenomes so far, the further use of it will benefit the study on mitogenome based on whole genome NGS data. MEANGS is available at https://github.com/YanCCscu/meangs.

Genomic evidence reveals intraspecific divergence of the hot-spring snake (Thermophis baileyi), an endangered reptile endemic to the Qinghai-Tibet plateau.

Understanding how and why species evolve requires knowledge on intraspecific divergence. In this study, we examined intraspecific divergence in the endangered hot-spring snake (Thermophis baileyi), an endemic species on the Qinghai-Tibet Plateau (QTP). Whole-genome resequencing of 58 sampled individuals from 15 populations was performed to identify the drivers of intraspecific divergence and explore the potential roles of genes under selection. Our analyses resolved three groups, with major intergroup admixture occurring in regions of group contact. Divergence probably occurred during the Pleistocene as a result of glacial climatic oscillations, Yadong-Gulu rift, and geothermal fields differentiation, while complex gene flow between group pairs reflected a unique intraspecific divergence pattern on the QTP. Intergroup fixed loci involved selected genes functionally related to divergence and local adaptation, especially adaptation to hot spring microenvironments in different geothermal fields. Analysis of structural variants, genetic diversity, inbreeding, and genetic load indicated that the hot-spring snake population has declined to a low level with decreased diversity, which is important for the conservation management of this endangered species. Our study demonstrated that the integration of demographic history, gene flow, genomic divergence genes, and other information is necessary to distinguish the evolutionary processes involved in speciation.

Effects of climate change and anthropogenic activity on ranges of vertebrate species endemic to the Qinghai-Tibet Plateau over 40 years.

Over the past 40 years, the climate has been changing and human disturbance has increased in the vast Qinghai-Tibet Plateau (QTP). These 2 factors are expected to affect the distribution of a large number of endemic vertebrate species. However, quantitative relationships between range shifts and climate change and human disturbance of these species in the QTP have rarely been evaluated. We used occurrence records of 19 terrestrial vertebrate species (birds, mammals, amphibians, and reptiles) occurring in the QTP from 1980 to 2020 to quantify the effects of climate change and anthropogenic impacts on the distribution of these 4 taxonomic groups and estimated species range changes in each species. The trend in distribution changes differed among the taxonomic groups, although, generally, ranges shifted to central QTP. Climate change contributed more to range variation than human disturbance (the sum of the 4 climatic variables contributed more than the sum of the 4 human disturbance variables for all 4 taxonomic groups). Suitable geographic range increased for most mammals, amphibians, and reptiles (+27.6%, +18.4%, and +27.8% on average, respectively), whereas for birds range decreased on average by 0.9%. Quantitative evidence for climate change and human disturbance associations with range changes for endemic vertebrate species in the QTP can provide useful insights into biodiversity conservation under changing environments.

En los últimos 40 años, el clima ha cambiado y las perturbaciones humanas han aumentado en la vasta meseta Qinghai-Tíbet (MQT). Se espera que estos dos factores afecten la distribución de un gran número de especies de vertebrados endémicos. Sin embargo, las relaciones cuantitativas entre los cambios del área de distribución y el cambio climático y las perturbaciones humanas en estas especies de la MQT han sido poco evaluadas. Utilizamos los registros de presencia de 19 especies de vertebrados terrestres (aves, mamíferos, anfibios y reptiles) de la MQT tomados entre 1980 y 2020 para cuantificar los efectos del cambio climático y los impactos antropogénicos sobre la distribución de estos cuatro grupos taxonómicos y estimar los cambios en el área de distribución de cada especie. La tendencia en los cambios de distribución difirió entre los grupos taxonómicos, aunque, en general, las áreas de distribución se desplazaron hacia el centro de la MQT. El cambio climático contribuyó más a la variación del área de distribución que las perturbaciones humanas (la suma de las cuatro variables climáticas contribuyó más que la suma de las cuatro variables de perturbaciones humanas para los cuatro grupos taxonómicos). El área de distribución geográfica adecuada aumentó para la mayoría de los mamíferos, anfibios y reptiles (+27.6%, +18.4% y +27.8% en promedio, respectivamente), mientras que para las aves disminuyó en promedio un 0.9%. Las pruebas cuantitativas de la asociación del cambio climático y las perturbaciones humanas con los cambios en el área de distribución de las especies vertebradas endémicas de la MQT pueden aportar información útil para la conservación de la biodiversidad en entornos cambiantes.

Chiral Adaptive Induction of an Achiral Cucurbit[8]uril-Based Supramolecular Organic Framework by Dipeptides in Water.

Chiral induction by natural biomolecules can reveal the indispensable role of chiral structures in life and can be used to develop the chirality-sensing biomolecular recognition. Here, we present the synthesis and characterization of an achiral supramolecular organic framework (SOF-1) constructed from cucurbit[8]uril (CB[8]) and hexaphenylbenzene (HPB) derivative (1) in water. Due to the propeller-like rotational chiral conformation of HPB units and the specific recognition properties of CB[8], SOF-1 demonstrates chiral adaptive induction in water when interacting with the N-terminal Trp-/Phe-containing dipeptides including L-TrpX and L-PheX (X is an amino acid residue), respectively, exhibiting contrasting circular dichroism (CD) and circularly polarized luminescence (CPL) spectra. Consequently, SOF-1 has been developed as a supramolecular host and chiroptical sensor capable of recognizing and distinguishing the sequence-opposite Trp-/Phe-containing dipeptide pairs including L-TrpX/L-XTrp and L-PheX/L-XPhe based on the sequence-selective CD responses.

Genome of Laudakia sacra Provides New Insights into High-Altitude Adaptation of Ectotherms.

Anan's rock agama (Laudakia sacra) is a lizard species endemic to the harsh high-altitude environment of the Qinghai-Tibet Plateau, a region characterized by low oxygen tension and high ultraviolet (UV) radiation. To better understand the genetic mechanisms underlying highland adaptation of ectotherms, we assembled a 1.80-Gb L. sacra genome, which contained 284 contigs with an N50 of 20.19 Mb and a BUSCO score of 93.54%. Comparative genomic analysis indicated that mutations in certain genes, including HIF1A, TIE2, and NFAT family members and genes in the respiratory chain, may be common adaptations to hypoxia among high-altitude animals. Compared with lowland reptiles, MLIP showed a convergent mutation in L. sacra and the Tibetan hot-spring snake (Thermophis baileyi), which may affect their hypoxia adaptation. In L. sacra, several genes related to cardiovascular remodeling, erythropoiesis, oxidative phosphorylation, and DNA repair may also be tailored for adaptation to UV radiation and hypoxia. Of note, ERCC6 and MSH2, two genes associated with adaptation to UV radiation in T. baileyi, exhibited L. sacra-specific mutations that may affect peptide function. Thus, this study provides new insights into the potential mechanisms underpinning high-altitude adaptation in ectotherms and reveals certain genetic generalities for animals' survival on the plateau.

The Genome of Shaw's Sea Snake (Hydrophis curtus) Reveals Secondary Adaptation to Its Marine Environment.

The transition of terrestrial snakes to marine life ∼10 Ma is ideal for exploring adaptive evolution. Sea snakes possess phenotype specializations including laterally compressed bodies, paddle-shaped tails, valvular nostrils, cutaneous respiration, elongated lungs, and salt glands, yet, knowledge on the genetic underpinnings of the transition remains limited. Herein, we report the first genome of Shaw's sea snake (Hydrophis curtus) and use it to investigate sea snake secondary marine adaptation. A hybrid assembly strategy obtains a high-quality genome. Gene family analyses date a pulsed coding-gene expansion to ∼20 Ma, and these genes associate strongly with adaptations to marine environments. Analyses of selection pressure and convergent evolution discover the rapid evolution of protein-coding genes, and some convergent features. Additionally, 108 conserved noncoding elements appear to have evolved quickly, and these may underpin the phenotypic changes. Transposon elements may contribute to adaptive specializations by inserting into genomic regions around functionally related coding genes. The integration of genomic and transcriptomic analyses indicates independent origins and different components in sea snake and terrestrial snake venom; the venom gland of the sea snake harbors the highest PLA2 (17.23%) expression in selected elapids and these genes may organize tandemly in the genome. These analyses provide insights into the genetic mechanisms that underlay the secondary adaptation to marine and venom production of this sea snake.

Age-related gene expression and DNA methylation changes in rhesus macaque.

Aging is a very complicated biological process that can change gene expressions. The Chinese rhesus macaque (Macaca mulatta lasiota; CR) is closely related to humans. We explored gene expression with increasing age and DNA methylation changes in young and old CRs. Results showed blood transcriptome and DNA methylome significantly changed from young to old CRs. The age-associated differentially expressed genes (DEGs) and differentially methylated regions (DMRs) were associated with age-related biological features, such as immunity, blood coagulation, and biosynthetic process. The measurements of coagulation indicators confirmed old CRs had shorter coagulation time than young CRs, and the activities of coagulation factor II (FII) and factor VIII (FVIII) were enhanced in old CRs. Humans and CRs exhibited the same enhanced blood coagulation with age phenotype. Our study found aging is a critical factor affecting gene expression in CRs, and also provided new insights into the blood coagulation changes in non-human primates.

The first draft genome of Lophophorus: A step forward for Phasianidae genomic diversity and conservation.

The monal genus (Lophophorus) is a branch of Phasianidae and its species inhabit the high-altitude mountains of the Qinghai-Tibet Plateau. The Chinese monal, L. lhuysii, is a threatened endemic bird of China that possesses high-altitude adaptability, diversity of plumage color and potentially low reproductive life history. This is the first study to describe the monal genome using next generation sequencing technology. The Chinese monal genome size is 1.01 Gb, with 16,940 protein-coding genes. Gene annotation yielded 100.93 Mb (9.97%) repeat elements, 785 ncRNA, 5,465,549 bp (0.54%) SSR and 15,550 (92%) genes in public databases. Compared to other birds and mammals, the genome evolution analysis showed numerous expanded gene families and positive selected genes involved in high-altitude adaptation, especially related to the adaptation of low temperature and hypoxia. Consequently, this gene data can be used to investigate the molecular evolution of high-altitude adaptation in future bird research. Our first published genome of the genus Lophophorus will be integral for the study of monal population genetic diversity and conservation, genomic evolution and Galliformes species differentiation in the Qinghai-Tibetan Plateau.

Unraveling the structure-chirality sensing relationship between achiral anthracene-based tetracationic nanotubes and nucleosides in aqueous host-guest complexation.

In biological systems, nucleosides play crucial roles in various physiological processes. In this study, we designed and synthesized four achiral anthracene-based tetracationic nanotubes (1-4) as artificial hosts and chiroptical sensors for nucleosides in aqueous media. Notably, different nanotubes exhibit varied chirality sensing on circular dichroism (CD)/circularly polarized luminescence (CPL) spectra through the host-guest complexation, which prompted us to explore the factors influencing their chiroptical responses. Through systematic host-guest experiments, the structure-chirality sensing relationship between achiral anthracene-based tetracationic nanotubes and nucleosides in the host-guest complexation was unraveled. Firstly, the CD response originates from the anthracene rings situated at the side-wall position, resulting from the right-handed (P)- or left-handed (M)-twisted conformation of the macrocyclic structure. Secondly, the CPL signal is influenced by the presence of anthracene rings at the linking-wall position, which results from intermolecular chiral twisted stacking between these anthracene rings. Therefore, these nanotubes can serve as chiroptical sensor arrays to enhance the accuracy of nucleotide recognition through principal component analysis (PCA) analysis based on the diversified CD spectra. This study provides insights for the construction of adaptive chirality from achiral nanotubes with dynamic conformational nature and might facilitate further design of chiral functional materials for several applications.

"Gear-driven"-type chirality transfer of tetraphenylethene-based supramolecular organic frameworks for peptides in water.

Chirality transfer for natural chiral biomolecules can reveal the indispensable role of chiral structures in life and can be used to develop the chirality-sensing biomolecular recognition. Here, we report the synthesis and characterization of a series of achiral supramolecular organic frameworks (SOF-1, SOF-2, and SOF-3), constructed from cucurbit[8]uril (CB[8]) and tetraphenylethene (TPE) derivatives (1, 2, and 3), respectively, as chirality-sensing platforms to explore their chirality transfer mechanism for peptides in water. Given the right-handed (P) and left-handed (M) rotational conformation of TPE units and the selective binding of CB[8] to aromatic amino acids, these achiral SOFs can be selectively triggered in water by peptides containing N-terminal tryptophan (W) and phenylalanine (F) residues into their P- or M-rotational conformation, exhibiting significantly different circular dichroism (CD) spectra. Although various peptides have the same l-type chiral configuration, they can induce positive CD signals of SOF-1 and negative CD signals of SOF-2 and SOF-3, respectively. Based on the structural analysis of the linkage units between CB[8] and TPE units in these SOFs, a "gear-driven"-type chirality transfer mechanism has been proposed to visually illustrate the multiple-step chirality transfer process from the recognition site in the CB[8]'s cavity to TPE units. Furthermore, by utilizing the characteristic CD signals generated through the "gear-driven"-type chirality transfer, these SOFs can serve as chiroptical sensor arrays to effectively recognize and distinguish various peptides based on their distinctive CD spectra.

Molecule Sequence Generation with Rebalanced Variational Autoencoder Loss.

Molecule generation is the procedure to generate initial novel molecule proposals for molecule design. Molecules are first projected into continuous vectors in chemical latent space, and then, these embedding vectors are decoded into molecules under the variational autoencoder (VAE) framework. The continuous latent space of VAE can be utilized to generate novel molecules with desired chemical properties and further optimize the desired chemical properties of molecules. However, there is a posterior collapse problem with the conventional recurrent neural network-based VAEs for the molecule sequence generation, which deteriorates the generation performance. We investigate the posterior collapse problem and find that the underestimated reconstruction loss is the main factor in the posterior collapse problem in molecule sequence generation. To support our conclusion, we present both analytical and experimental evidence. What is more, we propose an efficient and effective solution to fix the problem and prevent posterior collapse. As a result, our method achieves competitive reconstruction accuracy and validity score on the benchmark data sets.

Comparison of spatiotemporal, kinematic, and kinetic gait characteristics in total and unicompartmental knee arthroplasty during level walking: A systematic review and meta-analysis.

PURPOSE: This meta-analysis was performed to compare the spatiotemporal, kinematic, and kinetic gait characteristics during level walking between total knee arthroplasty (TKA) and unicompartmental knee arthroplasty (UKA). METHODS: An electronic database literature search was performed to screen clinical trials which were included the studies evaluating not only spatiotemporal, kinematic, and kinetic gait parameters, but also knee range of motion and knee score (Knee Society Score and Oxford Knee Score, i.e., KSS and OKS). The data analysis was performed using statistical software Stata 14.0 and Review Manager 5.4. RESULTS: Thirteen studies (369 knees) that met the criteria were eventually included in this meta-analysis. The results revealed significant differences between UKA and TKA with regard to walking speed (P = 0.04), stride length (P = 0.02), maximum knee flexion at loading (P = 0.001), the 1st peak of vert-GRF (P = 0.006), the 1st valley of vert-GRF (P = 0.007), knee internal rotational moment (P = 0.04), knee extension (P < 0.00001), and KSS Function score (P = 0.05). In contrast, there were no statistical differences in the remaining spatiotemporal, kinematic, and kinetic gait parameters. CONCLUSION: Medial UKA design is superior to TKA design with regard to walking speed, stride length, maximum knee flexion at loading, the 1st peak and the 1st valley of vert-GRF, knee internal rotational moment, knee extension, and KSS Function score. And it could provide a stronger basis for physicians to make clinical decisions.

[Correlation between preoperative MRI measurement of cross-sectional area of hamstring tendon and graft in anterior cruciate ligament reconstruction].

OBJECTIVE: To study the corretation between the cross-sectional area of hamstring tenden measured by MRI and gragt in anterior cruciate ligament rexonstruction. METHODS: MRI data of 50 patients who planned to undergo anterior cruciate ligament reconstruction from November 2021 to March 2022 were collected, including 32 males and 18 females, aged from 19 to 48 years old with an average of(31.1±8.7) years. Before the operation, the semitendinosus and gracilis tendons were measured and recorded by MRI, and then the anterior cruciate ligament was reconstructed under arthroscope. During the operation, gracilis and semitendinosus tendons were taken to prepare the final tendon to be transplanted, and the diameter of the prepared final graft was measured during the operation. Finally, the data were analyzed by statistical software. RESULTS: The cross sectional areas of semitendinosus tendon, gracilis tendon, semitendinosus tendon and gracilis tendon measured by MRI were significantly and positively correlated with the diameter of grafts required in anterior cruciate ligament surgery, the r values were 0.858, 0.728, 0.842(P<0.001), respectively. The area under curre (AUC), sensitivity, and specificity of the sum of the cross sectional areas of semitendinosus tendon and gracilis tendon were 0.925, 90.48%, and 85.71%, respectively. CONCLUSION: In patients undergoing anterior cruciate ligament reconstruction, preoperative MRI measurement has a strong statistical correlation with the diameter of hamstring muscle transplantation during operation. The sum of the cross sectional areas of semitendinosus tendon and gracilis tendon has a high predictive value for the diameter of grafts during anterior cruciate ligament reconstruction, and can predict the size of grafts during operation.

Developmental regulation of conserved non-coding element evolution provides insights into limb loss in squamates.

Limb loss shows recurrent phenotypic evolution across squamate lineages. Here, based on three de novo-assembled genomes of limbless lizards from different lineages, we showed that divergence of conserved non-coding elements (CNEs) played an important role in limb development. These CNEs were associated with genes required for limb initiation and outgrowth, and with regulatory signals in the early stage of limb development. Importantly, we identified the extensive existence of insertions and deletions (InDels) in the CNEs, with the numbers ranging from 111 to 756. Most of these CNEs with InDels were lineage-specific in the limbless squamates. Nearby genes of these InDel CNEs were important to early limb formation, such as Tbx4, Fgf10, and Gli3. Based on functional experiments, we found that nucleotide mutations and InDels both affected the regulatory function of the CNEs. Our study provides molecular evidence underlying limb loss in squamate reptiles from a developmental perspective and sheds light on the importance of regulatory element InDels in phenotypic evolution.

Genetic mapping and molecular mechanism behind color variation in the Asian vine snake.

BACKGROUND: Reptiles exhibit a wide variety of skin colors, which serve essential roles in survival and reproduction. However, the molecular basis of these conspicuous colors remains unresolved. RESULTS: We investigate color morph-enriched Asian vine snakes (Ahaetulla prasina), to explore the mechanism underpinning color variations. Transmission electron microscopy imaging and metabolomics analysis indicates that chromatophore morphology (mainly iridophores) is the main basis for differences in skin color. Additionally, we assemble a 1.77-Gb high-quality chromosome-anchored genome of the snake. Genome-wide association study and RNA sequencing reveal a conservative amino acid substitution (p.P20S) in SMARCE1, which may be involved in the regulation of chromatophore development initiated from neural crest cells. SMARCE1 knockdown in zebrafish and immunofluorescence verify the interactions among SMARCE1, iridophores, and tfec, which may determine color variations in the Asian vine snake. CONCLUSIONS: This study reveals the genetic associations of color variation in Asian vine snakes, providing insights and important resources for a deeper understanding of the molecular and genetic mechanisms related to reptilian coloration.

Evolution of lasR mutants in polymorphic Pseudomonas aeruginosa populations facilitates chronic infection of the lung.

Chronic infection with the bacterial pathogen Pseudomonas aeruginosa often leads to coexistence of heterogeneous populations carrying diverse mutations. In particular, loss-of-function mutations affecting the quorum-sensing regulator LasR are often found in bacteria isolated from patients with lung chronic infection and cystic fibrosis. Here, we study the evolutionary dynamics of polymorphic P. aeruginosa populations using isolates longitudinally collected from patients with chronic obstructive pulmonary disease (COPD). We find that isolates deficient in production of different sharable extracellular products are sequentially selected in COPD airways, and lasR mutants appear to be selected first due to their quorum-sensing defects. Polymorphic populations including lasR mutants display survival advantages in animal models of infection and modulate immune responses. Our study sheds light on the multistage evolution of P. aeruginosa populations during their adaptation to host lungs.

RetroComposer: Composing Templates for Template-Based Retrosynthesis Prediction.

The main target of retrosynthesis is to recursively decompose desired molecules into available building blocks. Existing template-based retrosynthesis methods follow a template selection stereotype and suffer from limited training templates, which prevents them from discovering novel reactions. To overcome this limitation, we propose an innovative retrosynthesis prediction framework that can compose novel templates beyond training templates. As far as we know, this is the first method that uses machine learning to compose reaction templates for retrosynthesis prediction. Besides, we propose an effective reactant candidate scoring model that can capture atom-level transformations, which helps our method outperform previous methods on the USPTO-50K dataset. Experimental results show that our method can produce novel templates for 15 USPTO-50K test reactions that are not covered by training templates. We have released our source implementation.

Genome Sequencing of Amomum tsao-ko Provides Novel Insight Into Its Volatile Component Biosynthesis.

As an important economic and medicinal crop, Amomum tsao-ko is rich in volatile oils and widely used in food additives, essential oils, and traditional Chinese medicine. However, the lack of the genome remains a limiting factor for understanding its medicinal properties at the molecular level. Here, based on 288.72 Gb of PacBio long reads and 105.45 Gb of Illumina paired-end short reads, we assembled a draft genome for A. tsao-ko (2.70 Gb in size, contig N50 of 2.45 Mb). Approximately 90.07% of the predicted genes were annotated in public databases. Based on comparative genomic analysis, genes involved in secondary metabolite biosynthesis, flavonoid metabolism, and terpenoid biosynthesis showed significant expansion. Notably, the DXS, GGPPS, and CYP450 genes, which participate in rate-limiting steps for terpenoid backbone biosynthesis and modification, may form the genetic basis for essential oil formation in A. tsao-ko. The assembled A. tsao-ko draft genome provides a valuable genetic resource for understanding the unique features of this plant and for further evolutionary and agronomic studies of Zingiberaceae species.