A near complete genome assembly of the East Friesian sheep genome.

Advancements in sequencing have enabled the assembly of numerous sheep genomes, significantly advancing our understanding of the link between genetic variation and phenotypic traits. However, the genome of East Friesian sheep (Ostfriesisches Milchschaf), a key high-yield milk breed, remains to be fully assembled. Here, we constructed a near-complete and gap-free East Friesian genome assembly using PacBio HiFi, ultra-long ONT and Hi-C sequencing. The resulting genome assembly spans approximately 2.96 Gb, with a contig N50 length of 104.1 Mb and only 164 unplaced sequences. Remarkably, our assembly has captured 41 telomeres and 24 centromeres. The assembled sequence is of high quality on completeness (BUSCO score: 97.1%) and correctness (QV: 69.1). In addition, a total of 24,580 protein-coding genes were predicted, of which 97.2% (23,891) carried at least one conserved functional domain. Collectively, this assembly provides not only a near T2T gap-free genome, but also provides a valuable genetic resource for comparative genome studies of sheep and will serve as an important tool for the sheep research community.

Whole-genome Sequencing Reveals Autooctoploidy in Chinese Sturgeon and Its Evolutionary Trajectories.

The order Acipenseriformes, which includes sturgeons and paddlefishes, represents "living fossils" with complex genomes that are good models for understanding whole-genome duplication (WGD) and ploidy evolution in fishes. Here, we sequenced and assembled the first high-quality chromosome-level genome for the complex octoploid Acipenser sinensis (Chinese sturgeon), a critically endangered species that also represents a poorly understood ploidy group in Acipenseriformes. Our results show that A. sinensis is a complex autooctoploid species containing four kinds of octovalents (8n), a hexavalent (6n), two tetravalents (4n), and a divalent (2n). An analysis taking into account delayed rediploidization reveals that the octoploid genome composition of Chinese sturgeon results from two rounds of homologous WGDs, and further provides insights into the timing of its ploidy evolution. This study provides the first octoploid genome resource of Acipenseriformes for understanding ploidy compositions and evolutionary trajectories of polyploid fishes.

Two high-quality Prototheca zopfii genomes provide new insights into their evolution as obligate algal heterotrophs and their pathogenicity.

The majority of the nearly 10,000 described species of green algae are photoautotrophs; however, some species have lost their ability to photosynthesize and become obligate heterotrophs that rely on parasitism for survival. Two high-quality genomes of the heterotrophic algae Prototheca zopfii Pz20 and Pz23 were obtained using short- and long-read genomic as well as transcriptomic data. The genome sizes were 31.2 Mb and 31.3 Mb, respectively, and contig N50 values of 1.99 Mb and 1.26 Mb. Although P. zopfii maintained its plastid genome, the transition to heterotrophy led to a reduction in both plastid and nuclear genome size, including the loss of photosynthesis-related genes from both the nuclear and plastid genomes and the elimination of genes encoding for carotenoid oxygenase and pheophorbide an oxygenase. The loss of genes, including basic leucine-zipper (bZIP) transcription factors, flavin adenine dinucleotide-linked oxidase, and helicase, could have played a role in the transmission of autotrophy to heterotrophs and in the processes of abiotic stress resistance and pathogenicity. A total of 66 (1.37%) and 73 (1.49%) genes were identified as potential horizontal gene transfer events in the two P. zopfii genomes, respectively. Genes for malate synthase and isocitrate lyase, which are horizontally transferred from bacteria, may play a pivotal role in carbon and nitrogen metabolism as well as the pathogenicity of Prototheca and non-photosynthetic organisms. The two high-quality P. zopfii genomes provide new insights into their evolution as obligate heterotrophs and pathogenicity. IMPORTANCE: The genus Prototheca, characterized by its heterotrophic nature and pathogenicity, serves as an exemplary model for investigating pathobiology. The limited understanding of the protothecosis infectious disease is attributed to the lack of genomic resources. Using HiFi long-read sequencing, both nuclear and plastid genomes were generated for two strains of P. zopfii. The findings revealed a concurrent reduction in both plastid and nuclear genome size, accompanied by the loss of genes associated with photosynthesis, carotenoid oxygenase, basic leucine-zipper (bZIP) transcription factors, and others. The analysis of horizontal gene transfer revealed the presence of 1.37% and 1.49% bacterial genes, including malate synthase and isocitrate lyase, which play crucial roles in carbon and nitrogen metabolism, as well as pathogenicity and obligate heterotrophy. The two high-quality P. zopfii genomes represent valuable resources for investigating their adaptation and evolution as obligate heterotrophs, as well as for developing future prevention and treatment strategies against protothecosis.

A complete reference genome of broomcorn millet.

Broomcorn millet (Panicum miliaceum L.), known for its traits of drought resistance, adaptability to poor soil, short growth period, and high photosynthetic efficiency as a C4 plant, represents one of the earliest domesticated crops globally. This study reports the telomere-to-telomere (T2T) gap-free reference genome for broomcorn millet (AJ8) using PacBio high-fidelity (HiFi) long reads, Oxford Nanopore long-read technologies and high-throughput chromosome conformation capture (Hi-C) sequencing data. The size of AJ8 genome was approximately 834.7 Mb, anchored onto 18 pseudo-chromosomes. Notably, 18 centromeres and 36 telomeres were obtained. The assembled genome showed high quality in terms of completeness (BUSCO score: 99.6%, QV: 61.7, LAI value: 20.4). In addition, 63,678 protein-coding genes and 433.8 Mb (~52.0%) repetitive sequences were identified. The complete reference genome for broomcorn millet provides a valuable resource for genetic studies and breeding of this important cereal crop.

The telomere-to-telomere (T2T) genome of Peucedanum praeruptorum Dunn provides insights into the genome evolution and coumarin biosynthesis.

BACKGROUND: Traditional Chinese medicine has used Peucedanum praeruptorum Dunn (Apiaceae) for a long time. Various coumarins, including the significant constituents praeruptorin (A-E), are the active constituents in the dried roots of P. praeruptorum. Previous transcriptomic and metabolomic studies have attempted to elucidate the distribution and biosynthetic network of these medicinal-valuable compounds. However, the lack of a high-quality reference genome impedes an in-depth understanding of genetic traits and thus the development of better breeding strategies. RESULTS: A telomere-to-telomere (T2T) genome was assembled for P. praeruptorum by combining PacBio HiFi, ONT ultra-long, and Hi-C data. The final genome assembly was approximately 1.798 Gb, assigned to 11 chromosomes with genome completeness >98%. Comparative genomic analysis suggested that P. praeruptorum experienced 2 whole-genome duplication events. By the transcriptomic and metabolomic analysis of the coumarin metabolic pathway, we presented coumarins' spatial and temporal distribution and the expression patterns of critical genes for its biosynthesis. Notably, the COSY and cytochrome P450 genes showed tandem duplications on several chromosomes, which may be responsible for the high accumulation of coumarins. CONCLUSIONS: A T2T genome for P. praeruptorum was obtained, providing molecular insights into the chromosomal distribution of the coumarin biosynthetic genes. This high-quality genome is an essential resource for designing engineering strategies for improving the production of these valuable compounds.

The rise of baobab trees in Madagascar.

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

Multiscale X-ray phase-contrast CT unveils the evolution of bile infarct in obstructive biliary disease.

Bile infarct is a pivotal characteristic of obstructive biliary disease, but its evolution during the disease progression remains unclear. Our objective, therefore, is to explore morphological alterations of the bile infarct in the disease course by means of multiscale X-ray phase-contrast CT. Bile duct ligation is performed in mice to mimic the obstructive biliary disease. Intact liver lobes of the mice are scanned by phase-contrast CT at various resolution scales. Phase-contrast CT clearly presents three-dimensional (3D) images of the bile infarcts down to the submicron level with good correlation with histological images. The CT data illustrates that the infarct first appears on day 1 post-BDL, while a microchannel between the infarct and hepatic sinusoids is identified, the number of which increases with the disease progression. A 3D model of hepatic acinus is proposed, in which the infarct starts around the portal veins (zone I) and gradually progresses towards the central veins (zone III) during the disease process. Multiscale phase-contrast CT offers the comprehensive analysis of the evolutionary features of the bile infarct in obstructive biliary disease. During the course of the disease, the bile infarcts develop infarct-sinusoidal microchannels and gradually occupy the whole liver, promoting the disease progression.

Obliteration of portal venules contributes to portal hypertension in biliary cirrhosis.

The effects of the obliteration of portal venules (OPV) in cirrhotic portal hypertension are poorly understood. To investigate its contribution to portal hypertension in biliary cirrhosis and its underlying mechanism, we evaluated OPV using two-dimensional (2D) histopathology in liver explants from patients with biliary atresia (BA, n = 63), primary biliary cholangitis (PBC, n = 18), and hepatitis B-related cirrhosis (Hep-B-cirrhosis, n = 35). Then, three-dimensional (3D) OPV was measured by X-ray phase-contrast CT in two parallel models in rats following bile duct ligation (BDL) or carbon tetrachloride (CCl4) administration, representing biliary cirrhosis and post-necrotic cirrhosis, respectively. The portal pressure was also measured in the two models. Finally, the effects of proliferative bile ducts on OPV were investigated. We found that OPV was significantly more frequent in patients with biliary cirrhosis, including BA (78.57 ± 16.45%) and PBC (60.00 ± 17.15%), than that in Hep-B-cirrhotic patients (29.43 ± 14.94%, p < 0.001). OPV occurred earlier, evidenced by the paired liver biopsy at a Kasai procedure (KP), and was irreversible even after a successful KP in the patients with BA. OPV was also significantly more frequent in the BDL models than in the CCl4 models, as shown by 2D and 3D quantitative analysis. Portal pressure was significantly higher in the BDL model than that in the CCl4 model. With the proliferation of bile ducts, portal venules were compressed and irreversibly occluded, contributing to the earlier and higher portal pressure in biliary cirrhosis. OPV, as a pre-sinusoidal component, plays a key role in the pathogenesis of portal hypertension in biliary cirrhosis. The proliferated bile ducts and ductules gradually take up the 'territory' originally attributed to portal venules and compress the portal venules, which may lead to OPV in biliary cirrhosis. © 2024 The Pathological Society of Great Britain and Ireland.

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean.

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.

Gap-free genome assembly of Salangid icefish Neosalanx taihuensis.

Neosalanx taihuensis is widely distributed in freshwater and brackish water areas in China. Due to its high commercial value, it has been artificially introduced into many lakes and reservoirs, showing strong ecological adaptability. Here, a gap-free chromosome-level reference genome was constructed by combining short reads, PacBio HiFi long reads, Nanopore ultralong reads and Hi-C data. The reference genome of N. taihuensis was 397.29 Mb with a contig N50 of 15.61 Mb. The assembled sequences were anchored to 28 chromosomes. Furthermore, 20,024 protein-coding genes and 98.16% of the predicted genes were annotated in publicly available biological databases. This high-quality gap-free assembled genome will provide an essential reference for studying the evolution and ecological adaptability of N. taihuensis.

Gap-free genome assembly of anadromous Coilia nasus.

The Chinese tapertail anchovy, Coilia nasus, is a socioeconomically important anadromous fish that migrates from near ocean waters to freshwater to spawn every spring. The analysis of genomic architecture and information of C. nasus were hindered by the previously released versions of reference genomes with gaps. Here, we report the assembly of a chromosome-level gap-free genome of C. nasus by incorporating high-coverage and accurate long-read sequence data with multiple assembly strategies. All 24 chromosomes were assembled without gaps, representing the highest completeness and assembly quality. We assembled the genome with a size of 851.67 Mb and used BUSCO to estimate the completeness of the assembly as 92.5%. Using a combination of de novo prediction, protein homology and RNA-seq annotation, 21,900 genes were functionally annotated, representing 99.68% of the total predicted protein-coding genes. The availability of gap-free reference genomes for C. nasus will provide the opportunity for understanding genome structure and function, and will also lay a solid foundation for further management and conservation of this important species.

The first high-quality chromosome-level genome of the Sipuncula Sipunculus nudus using HiFi and Hi-C data.

Sipuncula is a class of exocoelomic unsegmented animals whose evolutionary relationships are unresolved. The peanut worm Sipunculus nudus is a globally distributed, economically important species belonging to the class Sipuncula. Herein, we present the first high-quality chromosome-level assembly of S. nudus based on HiFi reads and high-resolution chromosome conformation capture (Hi-C) data. The assembled genome was 1,427 Mb, with a contig N50 length of 29.46 Mb and scaffold N50 length of 80.87 Mb. Approximately 97.91% of the genome sequence was anchored to 17 chromosomes. A BUSCO assessment showed that 97.7% of the expectedly conserved genes were present in the genome assembly. The genome was composed of 47.91% repetitive sequences, and 28,749 protein-coding genes were predicted. A phylogenetic tree demonstrated that Sipuncula belongs to Annelida and diverged from the common ancestor of Polychaeta. The high-quality chromosome-level genome of S. nudus will serve as a valuable reference for studies of the genetic diversity and evolution of Lophotrochozoa.

Genome-wide association study analysis to resolve the key regulatory mechanism of biomineralization in Pinctada fucata martensii.

Biomineralization-controlled exo-/endoskeleton growth contributes to body growth and body size diversity. Molluscan shells undergo ectopic biomineralization to form the exoskeleton and biocalcified "pearl" involved in invading defence. Notably, exo-/endoskeletons have a common ancestral origin, but their regulation and body growth are largely unknown. This study employed the pearl oyster, Pinctada fucata marntensii, a widely used experimental model for biomineralization in invertebrates, to perform whole-genome resequencing of 878 individuals from wild and breeding populations. This study characterized the genetic architecture of biomineralization-controlled growth and ectopic biomineralization. The insulin-like growth factor (IGF) endocrine signal interacted with ancient single-copy transcription factors to form the regulatory network. Moreover, the "cross-phylum" regulation of key long noncoding RNA (lncRNA) in bivalves and mammals indicated the conserved genetic and epigenetic regulation in exo-/endoskeleton growth. Thyroid hormone signal and apoptosis regulation in pearl oysters affected ectopic biomineralization in pearl oyster. These findings provide insights into the mechanism underlying the evolution and regulation of biomineralization in exo-/endoskeleton animals and ectopic biomineralization.

Insight into microvascular adaptive alterations in the Glisson system of biliary atresia after Kasai portoenterostomy using X-ray phase-contrast CT.

OBJECTIVES: To investigate microvascular alterations in the Glisson system of biliary atresia (BA) patients after Kasai portoenterostomy (KP) using three-dimensional (3D) virtual histopathology based on X-ray phase-contrast CT (PCCT). METHODS: Liver explants from BA patients were imaged using PCCT, and 32 subjects were included and divided into two groups: KP (n = 16) and non-KP (n = 16). Combined with histological analysis and 3D visualization technology, 3D virtual histopathological assessment of the biliary, arterial, and portal venous systems was performed. According to loop volume ratio, 3D spatial density, relative surface area, tortuosity, and other parameters, pathological changes of microvasculature in the Glisson system were investigated. RESULTS: In the non-KP group, bile ducts mostly manifested as radial multifurcated hyperplasia and twisted into loops. In the KP group, the bile duct hyperplasia was less, and the loop volume ratio of bile ducts decreased by 13.89%. Simultaneously, the arterial and portal venous systems presented adaptive alterations in response to degrees of bile duct hyperplasia. Compared with the non-KP group, the 3D spatial density of arteries in the KP group decreased by 3.53%, and the relative surface area decreased from 0.088 ± 0.035 to 0.039 ± 0.015 (p < .01). Deformed portal branches gradually recovered after KP, with a 2.93% increase in 3D spatial density and a decrease in tortuosity from 1.17 ± 0.06 to 1.14 ± 0.04 (p < .01) compared to the non-KP group. CONCLUSION: 3D virtual histopathology via PCCT clearly reveals the microvascular structures in the Glisson system of BA patients and provides key insights into the morphological mechanism of microvascular adaptation induced by biliary tract dredging after KP in BA disease. KEY POINTS: • 3D virtual histopathology via X-ray phase-contrast computed tomography clearly presented the morphological structures and pathological changes of microvasculature in the Glisson system of biliary atresia patients. • The morphological alterations of microvasculature in the Glisson system followed the competitive occupancy mechanism in the process of biliary atresia.

A Chromosome-level assembly of the Japanese eel genome, insights into gene duplication and chromosomal reorganization.

Japanese eels (Anguilla japonica) are commercially important species, harvested extensively for food. Currently, this and related species (American and European eels) are challenging to breed on a commercial basis. As a result, the wild stock is used for aquaculture. Moreover, climate change, habitat loss, water pollution, and altered ocean currents affect eel populations negatively. Accordingly, the International Union for Conservation of Nature lists Japanese eels as endangered and on its red list. Here we presented a high-quality genome assembly for Japanese eels and demonstrated that large chromosome reorganizations occurred in the events of third-round whole-genome duplications (3R-WRDs). Several chromosomal fusions and fissions have reduced the ancestral protochromosomal number of 25 to 19 in the Anguilla lineage. A phylogenetic analysis of the expanded gene families showed that the olfactory receptors (group δ and ζ genes) and voltage-gated Ca2+ channels expanded significantly. Both gene families are crucial for olfaction and neurophysiology. Additional tandem and proximal duplications occurred following 3R-WGD to acquire immune-related genes for an adaptive advantage against various pathogens. The Japanese eel assembly presented here can be used to study other Anguilla species relating to evolution and conservation.

Gapless genome assembly of East Asian finless porpoise.

In recent years, conservation efforts have increased for rare and endangered aquatic wildlife, especially cetaceans. However, the East Asian finless porpoise (Neophocaena asiaeorientalis sunameri), which has a wide distribution in China, has received far less attention and protection. As an endangered small cetacean, the lack of a chromosomal-level reference for the East Asian finless porpoise limits our understanding of its population genetics and conservation biology. To address this issue, we combined PacBio HiFi long reads and Hi-C sequencing data to generate a gapless genome of the East Asian finless porpoise that is approximately 2.5 Gb in size over its 21 autosomes and two sex chromosomes (X and Y). A total of 22,814 protein-coding genes were predicted where ~97.31% were functionally annotated. This high-quality genome assembly of East Asian finless porpoise will not only provide new resources for the comparative genomics of cetaceans and conservation biology of threatened species, but also lay a foundation for more speciation, ecology, and evolutionary studies. Measurement(s) Neophocaena asiaeorientalis sunameri • Gapless genome assembly • sequence annotation Technology Type(s) MGISEQ. 2000 • PacBio HiFi Sequencing • Hi-C Sample Characteristic - Organism Neophocaena asiaeorientalis sunameri Sample Characteristic - Environment seawater Sample Characteristic - Location Yellow Sea near Lianyungang City, Jiangsu Province, China.

Genomic basis of the giga-chromosomes and giga-genome of tree peony Paeonia ostii.

Tree peony (Paeonia ostii) is an economically important ornamental plant native to China. It is also notable for its seed oil, which is abundant in unsaturated fatty acids such as α-linolenic acid (ALA). Here, we report chromosome-level genome assembly (12.28 Gb) of P. ostii. In contrast to monocots with giant genomes, tree peony does not appear to have undergone lineage-specific whole-genome duplication. Instead, explosive LTR expansion in the intergenic regions within a short period (~ two million years) may have contributed to the formation of its giga-genome. In addition, expansion of five types of histone encoding genes may have helped maintain the giga-chromosomes. Further, we conduct genome-wide association studies (GWAS) on 448 accessions and show expansion and high expression of several genes in the key nodes of fatty acid biosynthetic pathway, including SAD, FAD2 and FAD3, may function in high level of ALAs synthesis in tree peony seeds. Moreover, by comparing with cultivated tree peony (P. suffruticosa), we show that ectopic expression of class A gene AP1 and reduced expression of class C gene AG may contribute to the formation of petaloid stamens. Genomic resources reported in this study will be valuable for studying chromosome/genome evolution and tree peony breeding.

Synchrotron microtomography image restoration via regularization representation and deep CNN prior.

BACKGROUND AND OBJECTIVE: Synchrotron-based X-ray microtomography (S-µCT) is a promising imaging technique that plays an important role in modern medical science. S-µCT systems often cause various artifacts and noises in the reconstructed CT images, such as ring artifacts, quantum noise, and electronic noise. In most situations, such noise and artifacts occur simultaneously, which results in a deterioration in the image quality and affects subsequent research. Due to the complexity of the distribution of these mixed artifacts and noise, it is difficult to restore the corrupted images. To address this issue, we propose a novel algorithm to remove mixed artifacts and noise from S-µCT images simultaneously. METHODS: There are two important aspects of our method. Regarding ring artifacts, because of their specific structural characteristics, regularization-based methods are more suitable; thus, low-rank tensor decomposition and total variation are utilized to represent their directional and locally piecewise smoothness properties. Moreover, to determine the implicit prior of the random noise, a convolutional neural network (CNN) based method is used. The advantages of traditional regularization and the deep CNN are then combined and embedded in a plug-and-play framework. Hence, an efficient image restoration algorithm is proposed to address the problem of mixed artifacts and noise in S-µCT images. RESULTS: Our proposed method was assessed by utilizing simulations and real data experiments. The qualitative results showed that the proposed method could effectively remove ring artifacts as well as random noise. The quantitative results demonstrated that the proposed method achieved almost the best results in terms of PSNR, SSIM and MAE compared to other methods. CONCLUSIONS: The proposed method can serve as an effective tool for restoring corrupted S-µCT images, and it has the potential to promote the application of S-µCT.

Comparative genome anatomy reveals evolutionary insights into a unique amphitriploid fish.

Triploids are rare in nature because of difficulties in meiotic and gametogenic processes, especially in vertebrates. The Carassius complex of cyprinid teleosts contains sexual tetraploid crucian carp/goldfish (C. auratus) and unisexual hexaploid gibel carp/Prussian carp (C. gibelio) lineages, providing a valuable model for studying the evolution and maintenance mechanism of unisexual polyploids in vertebrates. Here we sequence the genomes of the two species and assemble their haplotypes, which contain two subgenomes (A and B), to the chromosome level. Sequencing coverage analysis reveals that C. gibelio is an amphitriploid (AAABBB) with two triploid sets of chromosomes; each set is derived from a different ancestor. Resequencing data from different strains of C. gibelio show that unisexual reproduction has been maintained for over 0.82 million years. Comparative genomics show intensive expansion and alterations of meiotic cell cycle-related genes and an oocyte-specific histone variant. Cytological assays indicate that C. gibelio produces unreduced oocytes by an alternative ameiotic pathway; however, sporadic homologous recombination and a high rate of gene conversion also exist in C. gibelio. These genomic changes might have facilitated purging deleterious mutations and maintaining genome stability in this unisexual amphitriploid fish. Overall, the current results provide novel insights into the evolutionary mechanisms of the reproductive success in unisexual polyploid vertebrates.

Chromosome-level genome of Entada phaseoloides provides insights into genome evolution and biosynthesis of triterpenoid saponins.

As a medicinal herbal plant, Entada phaseoloides has high levels of secondary metabolites, particularly triterpenoid saponins, which are important resources for scientific research and medical applications. However, the lack of a reference genome for this genus has limited research on its evolution and utilization of its medicinal potential. In this study, we report a chromosome-scale genome assembly for E. phaseoloides using Illumina, Nanopore long reads and high-throughput chromosome conformation capture technology. The assembled reference genome is 456.18 Mb (scaffold N50 = 30.9 Mb; contig N50 = 6.34 Mb) with 95.71% of the sequences anchored onto 14 pseudochromosomes. E. phaseoloides was estimated to have diverged from the Leguminosae lineage at ~72.0 million years ago. With the integration of transcriptomic and metabolomic data, gene expression patterns and metabolite profiling of E. phaseoloides were determined in different tissues. The pattern of gene expression and metabolic profile of the kernel were distinct from those of other tissues. Furthermore, the evolution of certain gene families involved in the biosynthesis of triterpenoid saponins and terpenes was analysed and offers new insights into the formation of these two metabolites. Four CYP genes, one UGT gene and related transcription factors were identified as candidate genes contributing to regulation of triterpenoid saponin biosynthesis. As the first high-quality assembled reference genome in the genus Entada, it will not only provide new information for the evolutionary study of this genus and conservation biology of E. phaseoloides but also lay a foundation for the formation and utilization of secondary metabolites in medicinal plants.