Reference Genome Sequence Research Articles

Comparative Methods for Demystifying Spatial Transcriptomics.

Spatial Transcriptomics (ST), coined as the term for parallel RNA-Seq on cell populations ordered spatially on a histological tissue section, has recently become increasingly popular, especially in experiments where microfluidics-based single-cell sequencing fails, such as assays on neurons. ST platforms, like the 10x Visium technology investigated herein, therefore produce in a single experiment simultaneously thousands of RNA readouts, captured by an array of micrometer scale spots under the histological section. Therefore, a central challenge of analyzing ST experiments consists of analyzing the gene expression morphology of all spots to delineate clusters of similar cell mixtures, which are then compared to each other to identify up- or down-regulated marker genes. Moreover, another level of complexity in ST experiments, compared to traditional RNA-Seq, is imposed by staining the tissue section with protein markers of cells or cell components to identify spots providing relevant information afterward. The corresponding microscopy images need to be analyzed in addition to the RNA-Seq read mappings on the reference genome and transcriptome sequences. Focusing on the software suite provided by the Visium platform manufacturer, we break down the ST analysis pipeline into its four essential steps-the image analysis, the read alignment, the gene quantification, and the spot clustering-and compare results obtained when using reads from different subsets of spots and/or when employing alternative genome or transcriptome references. Our comparative analyses demonstrate the impact of spot selection and the choice of genome/transcriptome references on the analysis results when employing the manufacturer's pipeline.

Genomic predictions and candidate single nucleotide polymorphisms for growth, form, and wood properties of teak clones

Abstract Teak (Tectona grandis) is a premier hardwood timber, extensively cultivated throughout the tropics. The study aimed to undertake the first genomic selection of teak clones using single nucleotide polymorphism (SNP) markers and evaluate the prospects for prediction of stem volume (VOL), diameter (DBH), form Quality Score (QS), Specific Gravity (SG), and Heartwood Percentage (HW). Thirty-three elite teak clones from a genetic test in northern Costa Rica were genotyped at 10812 SNP loci and aligned to the teak reference genome sequence. The genomic relationship matrix revealed 7 cryptic groups of more highly related clones and within-group kinship values were consistent with full-sib families. Clone values for each of the 5 traits were calculated from 5-year field data for 26 of the geno-typed clones. Genomic prediction models were fitted to clone values using Bayesian Ridge Regression, BayesA, BayesB, BayesC, and Bayes Lasso. Predictive ability (PA) was greatest for DBH (up to 0.66), VOL (up to 0.63), SG (0.58), with 0.42 and 0.40 as the best PA for QS and HW, respectively. GWAS resulted in the discovery of 7 unique SNP markers for volume, 7 for DBH, 4 for QS, and 8 for SG. The small experimental size resulted in relatively large false discovery rates. Teak breeders are well positioned to benefit from the genomic tools and approaches to genotyping teak breeding and deployment populations on a large scale. Significant advantages are expected from widespread integration of genomic technologies into teak breeding systems to allow shortening of the breeding cycle.

Molecular tracking of Acidovorax citrulli: Unveiling pathogen dynamics and blotch disease outbreaks through specific markers

Acidovorax citrulli (Ac), a gram-negative bacterium, is the causal agent of bacterial fruit blotch (BFB), which poses a significant threat to cucurbit crop production worldwide. Understanding the genetic diversity of Ac is crucial for identifying sources of resistance and implementing effective disease management strategies. In this study, we conducted the first genetic characterization of Ac strains collected in Turkey using Inter-Simple Sequence Repeat (ISSR) markers. These markers were selected based on repetitive domains mapped on the complete reference genome sequence of Acidovorax citrulli strain NWB SC196. The identity of the Turkish strains was confirmed through molecular (PCR) and serological (Immunofluorescence test and ELISA) methods, while the selected ISSRs, which exhibited similarity to flanked regions in the pathogen's whole genome sequence, were employed to assess the genetic diversity among Ac strains. We compared the profiles of Turkish strains with those of a collection of Ac strains from various countries, including the US, to explore a possible common origin. Specifically, we considered the dissemination of these strains through rootstocks used for grafted seedling production (Cucurbita maxima ? Cucurbita moschata). The results demonstrated a shared genetic profile, suggesting a potential link between Ac strains collected in Turkey and foreign strains. The Mauve analysis, utilizing whole genome sequences of various Ac strains available in the NCBI database, displayed similar clustering patterns to those obtained using our selected molecular markers, confirming the discriminatory efficiency of our method. Based on the high discriminatory power of the selected markers, our proposed method offers a rapid and straightforward approach for genetic analysis of intraspecific variation and monitoring Ac gene flow across countries. The characterized strains and markers presented in this study serve as valuable resources and reference materials for further genetic investigations and tracking contamination sources associated with Ac.

The reference genome sequence of Artemisia argyi provides insights into secondary metabolism biosynthesis.

Artemisia argyi, a perennial herb of the genus Artemisia in the family Asteraceae, holds significant importance in Chinese traditional medicine, referred to as "Aicao". Here, we report a high-quality reference genome of Artemisia argyi L. cv. beiai, with a genome size up to 4.15 Gb and a contig N50 of 508.96 Kb, produced with third-generation Nanopore sequencing technology. We predicted 147,248 protein-coding genes, with approximately 68.86% of the assembled sequences comprising repetitive elements, primarily long terminal repeat retrotransposons(LTRs). Comparative genomics analysis shows that A. argyi has the highest number of specific gene families with 5121, and much more families with four or more members than the other 6 plant species, which is consistent with its more expanded gene families and fewer contracted gene families. Furthermore, through transcriptome sequencing of A. argyi in response to exogenous MeJA treatment, we have elucidated acquired regulatory insights into MeJA's impact on the phenylpropanoid, flavonoid, and terpenoid biosynthesis pathways of A. argyi. The whole-genome information obtained in this study serves as a valuable resource for delving deeper into the cultivation and molecular breeding of A. argyi. Moreover, it holds promise for enhancing genome assemblies across other members of the Asteraceae family. The identification of key genes establishes a solid groundwork for developing new varieties of Artemisia with elevated concentrations of active compounds.

A combination of long- and short-read genomics reveals frequent p-arm breakpoints within chromosome 21 complex genomic rearrangements

PurposeAlthough chromosome 21 is the smallest human chromosome, it is highly relevant in the pathogenicity of both cancer and congenital diseases, including Alzheimer disease and trisomy 21 (Down syndrome). In addition, cases with rare structural variants (SVs) of chromosome 21 have been reported. These events vary in size and include large chromosomal events, such as ring chromosomes and small partial aneuploidies. The p-arm of the acrocentric chromosome 21 was devoid of reference genomic sequence in GRCh37 and GRCh38, which hampered our ability to solve genomic rearrangements and find the mechanism of formation of disease-causing SVs. We hypothesize that conserved satellite structures and segmental duplications located on the p-arm play an important role in the formation of complex SVs involving chromosome 21. MethodsThree cases with complex chromosome 21 rearrangements were studied with a combination of short-read and long-read genome sequencing, as well as optical genome mapping. The data were aligned to the T2T-CHM13 assembly. ResultsWe were able to resolve all 3 complex chromosome 21 rearrangements in which 15, 8, and 26 breakpoints were identified, respectively. By comparing the identified SV breakpoints, we were able to pinpoint a region between 21p13 and 21p12 that appears to be frequently involved in chromosome 21 rearrangements. Importantly, we observed acrocentric satellite DNA at several breakpoint junctions suggesting an important role for those elements in the formation of complex SVs. ConclusionTaken together, our results provide further insights into the architecture and underlying mechanisms of complex rearrangements on acrocentric chromosomes.

Regulation of expression quantitative trait loci by SVA retrotransposons within the major histocompatibility complex.

Genomic and transcriptomic studies of expression quantitative trait loci (eQTL) revealed that SINE-VNTR-Alu (SVA) retrotransposon insertion polymorphisms (RIPs) within human genomes markedly affect the co-expression of many coding and noncoding genes by coordinated regulatory processes. This study examined the polymorphic SVA modulation of gene co-expression within the major histocompatibility complex (MHC) genomic region where more than 160 coding genes are involved in innate and adaptive immunity. We characterized the modulation of SVA RIPs utilizing the genomic and transcriptomic sequencing data obtained from whole blood of 1266 individuals in the Parkinson's Progression Markers Initiative (PPMI) cohort that included an analysis of human leukocyte antigen (HLA)-A regulation in a subpopulation of the cohort. The regulatory properties of eight SVAs located within the class I and class II MHC regions were associated with differential co-expression of 71 different genes within and 75 genes outside the MHC region. Some of the same genes were affected by two or more different SVA. Five SVA are annotated in the human genomic reference sequence GRCh38.p14/hg38, whereas the other three were novel insertions within individuals. We also examined and found distinct structural effects (long and short variants and the CT internal variants) for one of the SVA (R_SVA_24) insertions on the differential expression of the HLA-A gene within a subpopulation (550 individuals) of the PPMI cohort. This is the first time that many HLA and non-HLA genes (multilocus expression units) and splicing mechanisms have been shown to be regulated by eight structurally polymorphic SVA within the MHC genomic region by applying precise statistical analysis of RNA data derived from the blood samples of a human cohort population. This study shows that SVA within the MHC region are important regulators or rheostats of gene co-expression that might have potential roles in diversity, health, and disease.

Exome-wide benchmark of difficult-to-sequence regions using short-read next-generation DNA sequencing.

Next-generation DNA sequencing (NGS) in short-read mode has recently been used for genetic testing in various clinical settings. NGS data accuracy is crucial in clinical settings, and several reports regarding quality control of NGS data, primarily focusing on establishing NGS sequence read accuracy, have been published thus far. Variant calling is another critical source of NGS errors that remains unexplored at the single-nucleotide level despite its established significance. In this study, we used a machine-learning-based method to establish an exome-wide benchmark of difficult-to-sequence regions at the nucleotide-residue resolution using 10 genome sequence features based on real-world NGS data accumulated in The Genome Aggregation Database (gnomAD) of the human reference genome sequence (GRCh38/hg38). The newly acquired metric, designated the 'UNMET score,' along with additional lines of structural information from the human genome, allowed us to assess the sequencing challenges within the exonic region of interest using conventional short-read NGS. Thus, the UNMET score could provide a basis for addressing potential sequential errors in protein-coding exons of the human reference genome sequence GRCh38/hg38 in clinical sequencing.

Rapid mini-chromosome divergence among fungal isolates causing wheat blast outbreaks in Bangladesh and Zambia.

The fungal pathogen, Magnaporthe oryzae Triticum pathotype, causing wheat blast disease was first identified in South America and recently spread across continents to South Asia and Africa. Here, we studied the genetic relationship among isolates found on the three continents. Magnaporthe oryzae strains closely related to a South American field isolate B71 were found to have caused the wheat blast outbreaks in South Asia and Africa. Genomic variation among isolates from the three continents was examined using an improved B71 reference genome and whole-genome sequences. We found strong evidence to support that the outbreaks in Bangladesh and Zambia were caused by the introductions of genetically separated isolates, although they were all close to B71 and, therefore, collectively referred to as the B71 branch. In addition, B71 branch strains carried at least one supernumerary mini-chromosome. Genome assembly of a Zambian strain revealed that its mini-chromosome was similar to the B71 mini-chromosome but with a high level of structural variation. Our findings show that while core genomes of the multiple introductions are highly similar, the mini-chromosomes have undergone marked diversification. The maintenance of the mini-chromosome and rapid genomic changes suggest the mini-chromosomes may serve important virulence or niche adaptation roles under diverse environmental conditions.

Whole-genome sequences of Bipolaris bicolor, Curvularia hawaiiensis, Curvularia spicifera, and Exserohilum rostratum isolated from rice in Burkina Faso, France, Mali, and Pakistan.

Different fungal species of the Pleosporaceae family infect rice, causing similar symptoms. Reference genomic sequences are useful tools to study the evolution of these species and to develop accurate molecular diagnostic tools. Here, we report the complete genome sequences of Bipolaris bicolor, Curvularia hawaiiensis, Curvularia spicifera, and Exserohilum rostratum.

Chromosome-Level Genome Assembly for the Angiosperm Silene conica.

The angiosperm genus Silene has been the subject of extensive study in the field of ecology and evolution, but the availability of high-quality reference genome sequences has been limited for this group. Here, we report a chromosome-level assembly for the genome of Silene conica based on Pacific Bioscience HiFi, Hi-C, and Bionano technologies. The assembly produced 10 scaffolds (1 per chromosome) with a total length of 862 Mb and only ∼1% gap content. These results confirm previous observations that S. conica and its relatives have a reduced base chromosome number relative to the genus's ancestral state of 12. Silene conica has an exceptionally large mitochondrial genome (>11 Mb), predominantly consisting of sequence of unknown origins. Analysis of shared sequence content suggests that it is unlikely that transfer of nuclear DNA is the primary driver of this mitochondrial genome expansion. More generally, this assembly should provide a valuable resource for future genomic studies in Silene, including comparative analyses with related species that recently evolved sex chromosomes.

Reference-based genome compression using the longest matched substrings with parallelization consideration

BackgroundA large number of researchers have devoted to accelerating the speed of genome sequencing and reducing the cost of genome sequencing for decades, and they have made great strides in both areas, making it easier for researchers to study and analyze genome data. However, how to efficiently store and transmit the vast amount of genome data generated by high-throughput sequencing technologies has become a challenge for data compression researchers. Therefore, the research of genome data compression algorithms to facilitate the efficient representation of genome data has gradually attracted the attention of these researchers. Meanwhile, considering that the current computing devices have multiple cores, how to make full use of the advantages of the computing devices and improve the efficiency of parallel processing is also an important direction for designing genome compression algorithms.ResultsWe proposed an algorithm (LMSRGC) based on reference genome sequences, which uses the suffix array (SA) and the longest common prefix (LCP) array to find the longest matched substrings (LMS) for the compression of genome data in FASTA format. The proposed algorithm utilizes the characteristics of SA and the LCP array to select all appropriate LMSs between the genome sequence to be compressed and the reference genome sequence and then utilizes LMSs to compress the target genome sequence. To speed up the operation of the algorithm, we use GPUs to parallelize the construction of SA, while using multiple threads to parallelize the creation of the LCP array and the filtering of LMSs.ConclusionsExperiment results demonstrate that our algorithm is competitive with the current state-of-the-art algorithms in compression ratio and compression time.

Comparative and population genomics of buckwheat species reveal key determinants of flavor and fertility

Common buckwheat (Fagopyrum esculentum) is an ancient crop with a world-wide distribution. Due to its excellent nutritional quality and high economic and ecological value, common buckwheat is becoming increasingly important throughout the world. The availability of a high-quality reference genome sequence and population genomic data will accelerate the breeding of common buckwheat, but the high heterozygosity due to the outcrossing nature has greatly hindered the genome assembly. Here we report the assembly of a chromosome-scale high-quality reference genome of F. esculentum var. homotropicum, a homozygous self-pollinating variant of common buckwheat. Comparative genomics revealed that two cultivated buckwheat species, common buckwheat (F. esculentum) and Tartary buckwheat (F. tataricum), underwent metabolomic divergence and ecotype differentiation. The expansion of several gene families in common buckwheat, including FhFAR genes, is associated with its wider distribution than Tartary buckwheat. Copy number variation of genes involved in the metabolism of flavonoids is associated with the difference of rutin content between common and Tartary buckwheat. Furthermore, we present a comprehensive atlas of genomic variation based on whole-genome resequencing of 572 accessions of common buckwheat. Population and evolutionary genomics reveal genetic variation associated with environmental adaptability and floral development between Chinese and non-Chinese cultivated groups. Genome-wide association analyses of multi-year agronomic traits with the content of flavonoids revealed that Fh05G014970 is a potential major regulator of flowering period, a key agronomic trait controlling the yield of outcrossing crops, and that Fh06G015130 is a crucial gene underlying flavor-associated flavonoids. Intriguingly, we found that the gene translocation and sequence variation of FhS-ELF3 contribute to the homomorphic self-compatibility of common buckwheat. Collectively, our results elucidate the genetic basis of speciation, ecological adaptation, fertility, and unique flavor of common buckwheat, and provide new resources for future genomics-assisted breeding of this economically important crop.

A near-complete assembly of asparagus bean provides insights into anthocyanin accumulation in pods.

Asparagus bean (Vigna unguiculata ssp. sesquipedialis), a subspecies of V. unguiculata, is a vital legume crop widely cultivated in Asia for its tender pods consumed as vegetables. However, the existing asparagus bean assemblies still contain numerous gaps and unanchored sequences, which presents challenges to functional genomics research. Here, we present an improved reference genome sequence of an elite asparagus bean variety, Fengchan 6, achieved through the integration of nanopore ultra-long reads, PacBio high-fidelity reads, and Hi-C technology. The improved assembly is 521.3 Mb in length and demonstrates several enhancements, including a higher N50 length (46.4 Mb), an anchor ratio of 99.8%, and the presence of only one gap. Furthermore, we successfully assembled 14 telomeres and all 11 centromeres, including four telomere-to-telomere chromosomes. Remarkably, the centromeric regions cover a total length of 38.1 Mb, providing valuable insights into the complex architecture of centromeres. Among the 30 594 predicted protein-coding genes, we identified 2356 genes that are tandemly duplicated in segmental duplication regions. These findings have implications for defence responses and may contribute to evolutionary processes. By utilizing the reference genome, we were able to effectively identify the presence of the gene VuMYB114, which regulates the accumulation of anthocyanins, thereby controlling the purple coloration of the pods. This discovery holds significant implications for understanding the underlying mechanisms of color determination and the breeding process. Overall, the highly improved reference genome serves as crucial resource and lays a solid foundation for asparagus bean genomic studies and genetic improvement efforts.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs (HGs) both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

COPO - Managing sample metadata for biodiversity: considerations from the Darwin Tree of Life project.

Large-scale reference genome sequencing projects for all of biodiversity are underway and common standards have been in place for some years to enable the understanding and sharing of sequence data. However, the metadata that describes the collection, processing and management of samples, and link to the associated sequencing and genome data, are not yet adequately developed and standardised for these projects. At the time of writing, the Darwin Tree of Life (DToL) Project is over two years into its ten-year ambition to sequence all described eukaryotic species in Britain and Ireland. We have sought consensus from a wide range of scientists across taxonomic domains to determine the minimal set of metadata that we collectively deem as critically important to accompany each sequenced specimen. These metadata are made available throughout the subsequent laboratory processes, and once collected, need to be adequately managed to fulfil the requirements of good data management practice. Due to the size and scale of management required, software tools are needed. These tools need to implement rigorous development pathways and change management procedures to ensure that effective research data management of key project and sample metadata is maintained. Tracking of sample properties through the sequencing process is handled by Lab Information Management Systems (LIMS), so publication of the sequenced data is achieved via technical integration of LIMS and data management tools. Discussions with community members on how metadata standards need to be managed within large-scale programmes is a priority in the planning process. Here we report on the standards we developed with respect to a robust and reusable mechanism of metadata collection, in the hopes that other projects forthcoming or underway will adopt these practices for metadata.

Transcriptomic analysis of Andrias davidianus meat and experimental validation for exploring its bioactive components as functional foods

Andrias davidianus (Chinese giant salamander, CGS) is the largest and oldest extant amphibian species in the world and is a source of prospective functional food in China. However, the progress of functional peptides mining was slow due to lack of reference genome and protein sequence data. In this study, we illustrated full-length transcriptome sequencing to interpret the proteome of CGS meat and obtain 10703 coding DNA sequences. By functional annotation and amino acid composition analysis, we have discovered various genes related to signal transduction, and 16 genes related to longevity. We have also found vast variety of functional peptides through protein coding sequence (CDS) analysis by comparing the data obtained with the functional peptide database. Val-Pro-Ile predicted by the CDS analysis was released from the CGS meat through enzymatic hydrolysis, suggesting that our approach is reliable. This study suggested that transcriptomic analysis can be used as a reference to guide polypeptide mining in CGS meat, thereby providing a powerful mining strategy for the bioresources with unknown genomic and proteomic sequences.

The potentialities of omics resources for millet improvement.

Millets are nutrient-rich (nutri-rich) cereals with climate resilience attributes. However, its full productive potential is not realized due to the lack of a focused yield improvement approach, as evidenced by the available literature. Also, the lack of well-characterized genomic resources significantly limits millet improvement. But the recent availability of genomic data and advancement in omics tools has shown its enormous potential to enhance the efficiency and precision faced by conventional breeding in millet improvement. The development of high throughput genotyping platforms based on next-generation sequencing (NGS) has provided a low-cost method for genomic information, specifically for neglected nutri-rich cereals with the availability of a limited number of reference genome sequences. NGS has created new avenues for millet biotechnological interventions such as mutation-based study, GWAS, GS, and other omics technologies. The simultaneous discovery of high-throughput markers and multiplexed genotyping platform has aggressively aided marker-assisted breeding for millet improvement. Therefore, omics technology offers excellent opportunities to explore and combine useful variations for targeted traits that could impart high nutritional value to high-yielding cultivars under changing climatic conditions. In millet improvement, an in-depth account of NGS, integrating genomics data with different biotechnology tools, is reviewed in this context.

Analysis of coverage of Alu repeats by aligned genomic reads

Alu repeats occupy a notable part of the human genome and greatly affect processes related to genome integrity maintenance. One of the basic methods for studying variation in a genome, including Alu repeats is genome sequencing followed by mapping the sequenced reads to a reference genome sequence. The key feature of the read alignment is the depth of reference genome region coverage by mapped reads. In this paper, a new method is proposed for analyzing the coverage of Alu repeats and their flanking regions by whole-genome sequencing reads and the distribution of mean coverage in two aforementioned region types is explored.

Formation mechanism of glandular trichomes involved in the synthesis and storage of terpenoids in lavender

BackgroundLavender (genus Lavandula, family Lamiaceae) is an aromatic plant widely grown as an ornamental plant. The chemical composition of lavender is characterized by monoterpenoids, sesquiterpenoids, and other compounds, which are primarily synthesized and stored in epidermal secretory structures called glandular trichomes (GTs). Volatile organic compounds (VOCs) are responsible for the aroma characteristics of plant oil that drive consumer preference. Aroma is usually regarded as a characteristic trait for the classification of aromatic plants. Interestingly, VOCs are synthesized and stored in GTs. Lamiaceae species such as purple perilla, peppermint, basil, thyme, and oregano usually possess two types of GTs: peltate glandular trichomes (PGTs) and capitate glandular trichomes (CGTs). But the development process of PGTs in lavender has been reported in only a few studies to date.ResultsIn this study, we identified and quantified the VOCs in four lavender cultivars by headspace-solid phase micro extraction-gas chromatography mass spectrometry (HS–SPME–GC–MS). A total of 66 VOCs were identified in these four cultivars, the most prominent of which were linalyl acetate and linalool, and flowers were the main site of accumulation of these VOCs. Here, we examined the developmental process of PGTs, including the formation of their base, body, and apex. The apex cells contained secretory cavities, which produced VOCs. Based on the reference genome sequence of the lavender cultivar ‘Jingxun 2’, several R2R3-MYB subfamily genes related to GT formation were identified. These results will guide the engineering of GTs and molecular breeding of lavender for improving the VOC content.ConclusionsIn this study, we identified the VOCs in four lavender cultivars. We analyzed the formation of GTs, and compared the number and diameter size of PGTs among four lavender cultivars. Additionally, we identified four candidate genes belonging to the R2R3-MYB family.

High-resolution genome mapping and functional dissection of chlorogenic acid production in Lonicera maackii.

Amur honeysuckle (Lonicera maackii) is a widely used medicinal plant of the Caprifoliaceae family that produces chlorogenic acid. Research on this plant mainly focuses on its ornamental value and medicinal compounds, but a reference genome sequence and molecular resources for accelerated breeding are currently lacking. Herein, Nanopore sequencing and high-throughput chromosome conformation capture (Hi-C) allowed a chromosome-level genome assembly of Lonicera maackii (2n = 18). A global view of the gene regulatory network involved in the biosynthesis of chlorogenic acid and the dynamics of fruit coloration in L. maackii was established through metabolite profiling and transcriptome analyses. Moreover, we identified the genes encoding hydroxycinnamoyl-CoA quinate transferase (LmHQT) and hydroxycinnamoyl-CoA shikimic/quinate transferase (LmHCT), which localized to the cytosol and nucleus. Heterologous overexpression of these genes in Nicotiana benthamiana leaves resulted in elevated chlorogenic acid contents. Importantly, HPLC analyses revealed that LmHCT and LmHQTs recombinant proteins modulate the accumulation of chlorogenic acid (CGA) using quinic acid and caffeoyl CoA as substrates, highlighting the importance of LmHQT and LmHCT in chlorogenic acid biosynthesis. These results confirmed that LmHQTs and LmHCT catalyze the biosynthesis of CGA in vitro. The genomic data presented in this study will offer a valuable resource for the elucidation of CGA biosynthesis and facilitating selective molecular breeding.