De novo assembly and analysis of Sonneratia ovata genome and population analysis.

Mangroves are an important part of coastal and estuarine ecosystems where they serve as nurseries for marine species and prevent coastal erosion. Here we report the genome of Sonneratia ovata, which is a true mangrove that grows in estuarine environments and can tolerate moderate salt exposure. We sequenced the S. ovata genome and assembled it into chromosome-level scaffolds through the use of Hi-C. The genome is 212.3 Mb and contains 12 chromosomes that range in size from 12.2 to 23.2 Mb. Annotation identified 29,829 genes with a BUSCO completeness of 95.9%. We identified salt genes and found copy number expansion of salt genes such as ADP-ribosylation factor 1, and elongation factor 1-alpha. Population analysis identified a low level of genetic variation and a lack of population structure within S. ovata.

Rice height QTLs in KDML105 chromosome segment substitution lines.

Rice is an important crop that is consumed by approximately half of the world's population on a regular basis. Plant height is an important characteristic with shorter rice often having higher lodging resistance and better soil nutrient utilization allowing for lower fertilizer use. We used a Chromosome Segment Substitution Line (CSSL) population generated by introgressing segments of CT9993 and IR62266 into KDML 105. We identified height QTLs on chromosomes 1 and 4. We performed whole genome sequencing of the parental lines and found that IR62266 has the deletion in Gibberellin 20-oxidase 2 corresponding to the semi-dwarf 1 locus. However, short height on chromosome 1 came from CT9993 with no mutation in Gibberellin 20-oxidase 2, or any known height genes. The height QTL on chromosome 4 contains mutations in Peroxisome biogenesis protein 6, which has been linked to a reduced growth phenotype in A. thaliana, making this a good candidate height gene.

Assembly of a hybrid mangrove, Bruguiera hainesii, and its two ancestral contributors, Bruguiera cylindrica and Bruguiera gymnorhiza.

Mangroves are plants that live in tropical and subtropical coastal regions of the world, they are adapted to high salt environments and cyclic tidal flooding. Mangroves play important ecological roles, including acting as breeding grounds for many fish species and to prevent coastal erosion. The genomes of three mangrove species, Bruguiera gymnorhiza, Bruguiera cylindrica, and a hybrid of the two, Bruguiera hainesii were sequenced, assembled and annotated. The two progenitor species, B. gymnorhiza and B. cylindrica, were found to be highly similar to each other and sufficiently similar to B. parviflora to allow it to be used for reference based scaffolding to generate chromosome level scaffolds. The two subgenomes of B. hainesii were independently assembled and scaffolded. Analysis of B. hainesii confirms that it is a hybrid and the hybridisation event was estimated at 2.4 to 3.5 million years ago using a Bayesian Relaxed Molecular Clock approach.

De novo chromosome-level assembly of the Centella asiatica genome.

Centella asiatica is a herbaceous, perennial species indigenous to India and Southeast Asia. C. asiatica possesses several medicinal properties: anti-aging, anti-inflammatory, wound healing and memory enhancing. The lack of available genomics resources significantly impedes the improvement of C. asiatica varieties through molecular breeding. Here, we combined the 10× Genomics linked-read technology and the long-range HiC technique to obtain the genome assembly. The final assembly contained nine pseudomolecules, corresponding to the haploid chromosome number in C. asiatica. These nine chromosomes covered 402,536,584 bases or 93.6% of the 430-Mb assembly. Comparative genomics analyses based on single-copy orthologous genes showed that C. asiatica and the common ancestor of Coriandrum sativum (coriander) and Daucus carota (carrot) diverged about 48 million years ago. This assembly provides a valuable reference genome for future molecular studies, varietal development through marker-assisted breeding and comparative genomics studies in C. asiatica.

Taxonomic profiling of Symbiodiniaceae and bacterial communities associated with Indo-Pacific corals in the Gulf of Thailand using PacBio sequencing of full-length ITS and 16S rRNA genes.

Corals live with complex assemblages of microbes including bacteria, the dinoflagellate Symbiodiniaceae, fungi and viruses in a coral holobiont. These coral-associated microorganisms play an important role in their host fitness and survival. Here, we investigated the structure and diversity of algal and bacterial communities associated with five Indo-Pacific coral species, using full-length 16S rRNA and internal transcribed spacer sequences. While the dinoflagellate communities associated with Poriteslutea were dominated with Symbiodiniaceae genus Cladocopium, the other four coral hosts were associated mainly with members of the Durusdinium genus, suggesting that host species was one of the underlying factors influencing the structure and composition of dinoflagellate communities associated with corals in the Gulf of Thailand. Alphaproteobacteria dominated the microbiomes of Pocillopora spp. while Pavonafrondifera and P. lutea were associated primarily with Gammaproteobacteria. Finally, we demonstrated a superior performance of full-length 16S rRNA sequences in achieving species-resolution taxonomic classification of coral-associated microbiota.

Quantitative analysis of methoxyflavones discriminates between the two types of Kaempferia parviflora.

INTRODUCTION: Kaempferia parviflora or black ginger is abundantly cultivated because its rhizomes contain methoxyflavones that have many pharmacological properties. K. parviflora can be divided into two types, based on morphological characteristics, but differences in their chemical compositions have never been explored. OBJECTIVES: This research aims to find chemical markers that can be used to differentiate between the two types of K. parviflora, the red-leaf and green-leaf types, by quantifying the amounts of methoxyflavones. MATERIAL AND METHODS: K. parviflora samples were collected from 39 locations in Thailand. Their genetic diversity was assessed by a genotyping-by-sequencing (GBS) technique to construct the population structure. Their chemical compositions were analyzed by high performance liquid chromatography-photodiode array detection to determine the methoxyflavone contents. RESULTS: The population structure based on >3,000 single nucleotide polymorphism (SNP) markers showed that the samples can be divided into two groups, which were consistent with the classification by leaf margin color (red-leaf and green-leaf types). HPLC analysis revealed 3,5,7,3',4'-pentamethoxyflavone (PMF), 5,7-dimethoxyflavone (DMF), 5,7,4'-trimethoxyflavone (TMF), 3,5,7-trimethoxyflavone and 3,5,7,4'-tetramethoxyflavone as major methoxyflavones that can be used as chemical markers. The red-leaf type showed higher amounts of PMF, TMF and 3,5,7,4'-tetramethoxyflavone than the green-leaf type, while the green-leaf type showed higher amounts of DMF and 3,5,7-trimethoxyflavone than the red-leaf type. CONCLUSION: These results provide another approach to discriminate the two types of K. parviflora using chemical profiles alongside genetic and morphological analyses. Therefore, a specific type of K. parviflora can be selected over the other based on preferences for a certain methoxyflavone.

ACC oxidase and miRNA 159a, and their involvement in fresh fruit bunch yield (FFB) via sex ratio determination in oil palm.

Oil palm (Elaeis guineesis Jacq.) is the most productive oil-bearing crop, yielding more oil per area than any other oil-bearing crops. However, there are still efforts to improve oil palm yield, in order to serve consumer and manufacturer demand. Oil palm produces female and male inflorescences in an alternating cycle. So, high sex ratio (SR), the ratio of female inflorescences to the total inflorescences, is a favorable trait in term of increasing yields in oil palm. This study aims to understand the genetic control for SR related traits, such as fresh fruit bunch yield (FFB), by characterizing genes at FFB quantitative trait loci (QTLs) on linkage 10 (chromosome 6) and linkage 15 (chromosome 10). Published oil palm sequences at the FFB QTLs were used to develop gene-based and simple sequence repeat (SSR) markers. We used the multiple QTL analysis model (MQM) to characterize the relationship of new markers with the SR traits in the oil palm population. The RNA expression of the most linked QTL genes was also evaluated in various tissues of oil palm. We identified EgACCO1 (encoding aminocyclopropane carboxylate (ACC) oxidase) at chromosome 10 and EgmiR159a (microRNA 159a) at chromosome 6 to be the most linked QTL genes or determinants for FFB yield and/or female inflorescence number with a phenotype variance explained (PVE) from 10.4 to 15 % and suggest that these play the important roles in sex determination and differentiation in oil palm. The strongest expression of EgACCO1 and the predicted precursor of EgmiR159a was found in ovaries and, to a lesser extent, fruit development. In addition, highly normalized expression of EgmiR159a was found in female flowers. In summary, the QTL analysis and the RNA expression reveal that EgACCO1 and EgmiR159a are the potential genetic factors involved in female flower determination and hence would affect yield in oil palm. However, to clarify how these genetic factors regulate female flower determination, more investigation of their down regulation or target may be essential. Additionally, if more sex determination genes controlled by plant hormones are identified, it may possible to reveal a crosstalk of sex determination genes with hormones and environment factors.

Genome-wide SNP discovery and identification of QTL associated with agronomic traits in oil palm using genotyping-by-sequencing (GBS).

Oil palm has become one of the most important oil crops in the world. Marker-assisted selections have played a pivotal role in oil palm breeding programs. Here, we report the use of genotyping-by-sequencing (GBS) approach for a large-scale SNP discovery and genotyping of a mapping population. Reduced representation libraries of 108 F2 progeny were sequenced and a total of 524 million reads were obtained. We detected 21,471 single nucleotide substitutions, most of which (62.6%) represented transition events. Of 3417 fully informative SNP markers, we were able to place 1085 on a linkage map, which spanned 1429.6 cM and had an average of one marker every 1.26 cM. Three QTL affecting trunk height were detected on LG 10, 14 and 15, whereas a single QTL associated with fruit bunch weight was identified on LG 3. The use of GBS approach proved to be rapid, cost-effective and highly reproducible in this species.

The AKR gene family and modifying sex ratios in palms through abiotic stress responsiveness.

Sex ratio (SR), the ratio of female inflorescences to total inflorescences, is one of the main yield components of oil palm (Elaeis guineensis Jacq). The SR quantitative trait locus (QTL) was recently identified on linkage (LG) 8 with a phenotype variance explained (PVE) of 11.3 %. The use of both genetic and physical mapping is one strategy for uncovering the genetic basis of the traits. Here, we report the construction of bacterial artificial chromosome (BAC) and fosmid libraries, and their use for physical mapping in oil palm. Combined, the libraries consist of more than 200,000 clones, representing 6.35 genome equivalents. Physical mapping at the SR locus was implemented by incorporating the published oil palm genome sequence and positive BAC/fosmid clones as identified by colony PCR screening. Based on the previously published sequences, the interval (about 184 kb) was comprised of 19 contigs of the known sequences (~117 kb, 64 %). After, combining the 454 pyrosequences of 15 positive clones and the previously published sequences, the known sequences were revealed to cover about 82 % of the interval (~150 kb), and were used for identifying the new markers by designing 35 gene-based and 23 simple sequence repeat (SSR)-amplified primers. As a result, a putative aldo-keto reductase gene (named EgAKR1) was revealed to be a promising candidate for sex ratio determination, via controlling female inflorescence number (11 % of PVE). This was predicted from the two newly identified polymorphic marker loci (mEgSSRsr8-21LB and mEgAKR1-9) designing from EgAKR1. The functions of AKR gene families in other plant species and our promoter analysis suggested that EgAKR1 may contribute to the sex ratio through abiotic stress responsiveness.

The complete mitochondrial genome of the Hipposideros pendleburyi (Pendlebury's leaf-nosed bat) an endemic species in Thailand.

This study presents the first complete mitochondrial genome of the Hipposideros pendleburyi (Pendlebury's leaf-nosed bat), an endemic species in Thailand. The mitochondrial genome was 16,820 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. The overall base composition was 31.5% A, 26.2% T, 28.3% C, and 14.0% G. A maximum-likelihood tree revealed that H. pendleburyi was grouped with Hipposideros armiger within the Hipposideridae clade, which has Rhinolophidae as a sister clade.

De Novo Reference Assembly of the Upriver Orange Mangrove (Bruguiera sexangula) Genome.

Upriver orange mangrove (Bruguiera sexangula) is a member of the most mangrove-rich taxon (Rhizophoraceae family) and is commonly distributed in the intertidal zones in tropical and subtropical latitudes. In this study, we employed the 10× Genomics linked-read technology to obtain a preliminary de novo assembly of the B. sexangula genome, which was further scaffolded to a pseudomolecule level using the Bruguiera parviflora genome as a reference. The final assembly of the B. sexangula genome contained 260 Mb with an N50 scaffold length of 11,020,310 bases. The assembly comprised 18 pseudomolecules (corresponding to the haploid chromosome number in B. sexangula), covering 204,645,832 bases or 78.6% of the 260-Mb assembly. We predicted a total of 23,978 protein-coding sequences, 17,598 of which were associated with gene ontology terms. Our gene prediction recovered 96.6% of the highly conserved orthologs based on the Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. The chromosome-level assembly presented in this work provides a valuable genetic resource to help strengthen our understanding of mangroves' physiological and morphological adaptations to the intertidal zones.

The complete chloroplast genome sequence of Intsia bijuga (Colebr.) Kuntze (Fabaceae: Detaroideae: Afzelieae).

Intsia bijuga (Colebr.) Kuntze. (1891) is a threatened mangrove species, belonging to the Fabaceae family and is native to the western Pacific coast and Southeast Asia. Here, we applied short-read Illumina technology to sequence and assemble its chloroplast genome. The complete chloroplast genome is 158,363 bp in length, composed of one large single-copy (LSC) region of 87,489 bp, one small single-copy (SSC) region of 19,438 bp, and a pair of inverted repeats (IRs) of 25,719 bp. A total of 129 unique genes were annotated, comprising 84 protein-coding genes, eight rRNA genes, and 37 tRNA genes. Our phylogenetic analysis showed the placement of I. bijuga (OL699920.1) with Afzelia species within Fabaceae family.

A Chromosome-Scale Genome Assembly of Mitragyna speciosa (Kratom) and the Assessment of Its Genetic Diversity in Thailand.

Mitragyna speciosa (Kratom) is a tropical narcotic plant native to Southeast Asia with unique pharmacological properties. Here, we report the first chromosome-scale assembly of the M. speciosa genome. We employed PacBio sequencing to obtain a preliminary assembly, which was subsequently scaffolded using the chromatin contact mapping technique (Hi-C) into 22 pseudomolecules. The final assembly was 692 Mb with a scaffold N50 of 26 Mb. We annotated a total of 39,708 protein-coding genes, and our gene predictions recovered 98.4% of the highly conserved orthologs based on the BUSCO analysis. The phylogenetic analysis revealed that M. speciosa diverged from the last common ancestors of Coffea arabica and Coffea canephora approximately 47.6 million years ago. Our analysis of the sequence divergence at fourfold-degenerate sites from orthologous gene pairs provided evidence supporting a genome-wide duplication in M. speciosa, agreeing with the report that members of the genus Mitragyna are tetraploid. The STRUCTURE and principal component analyses demonstrated that the 85 M. speciosa accessions included in this study were an admixture of two subpopulations. The availability of our high-quality chromosome-level genome assembly and the transcriptomic resources will be useful for future studies on the alkaloid biosynthesis pathway, as well as comparative phylogenetic studies in Mitragyna and related species.

A SNP variation in an expansin (EgExp4) gene affects height in oil palm.

Oil palm (Elaeis guineensis Jacq.), an Aracaceae family plant, is utilized for both consumable and non-consumable products, including cooking oil, cosmetics and biodiesel production. Oil palm is a perennial tree with 25 years of optimal harvesting time and a height of up to 18 m. However, harvesting of oil palm fruit bunches with heights of more than 2-3 meters is challenging for oil palm farmers. Thus, understanding the genetic control of height would be beneficial for using gene-based markers to speed up oil palm breeding programs to select semi-dwarf oil palm varieties. This study aims to identify Insertion/Deletions (InDels) and single nucleotide polymorphisms (SNPs) of five height-related genes, including EgDELLA1, EgGRF1, EgGA20ox1, EgAPG1 and EgExp4, in short and tall oil palm groups by PacBio SMRT sequencing technology. Then, the SNP variation's association with height was validated in the Golden Tenera (GT) population. All targeted genes were successfully amplified by two rounds of PCR amplification with expected sizes that ranged from 2,516 to 3,015 base pair (bp), covering 5' UTR, gene sequences and 3' UTR from 20 short and 20 tall oil palm trees. As a result, 1,166, 909, 1,494, 387 and 5,384 full-length genomic DNA sequences were revealed by PacBio SMRT sequencing technology, from EgDELLA1, EgGRF1, EgGA20ox1, EgAPG1 and EgExp4 genes, respectively. Twelve variations, including eight InDels and four SNPs, were identified from EgDELLA1, EgGRF1, EgGA20ox1 and EgExp4. No variation was found for EgAPG1. After SNP through-put genotyping of 4 targeted SNP markers was done by PACE™ SNP genotyping, the association with height was determined in the GT population. Only the mEgExp4_SNP118 marker, designed from EgExp4 gene, was found to associate with height in 2 of 4 height-recordings, with p values of 0.0383 for height (HT)-1 and 0.0263 for HT-4. In conclusion, this marker is a potential gene-based marker that may be used in oil palm breeding programs for selecting semi-dwarf oil palm varieties in the near future.

A draft chromosome-scale genome assembly of a commercial sugarcane.

Sugarcane accounts for a large portion of the worlds sugar production. Modern commercial cultivars are complex hybrids of S. officinarum, S. spontaneum, and several other Saccharum species, resulting in an auto-allopolyploid with 8-12 copies of each chromosome. The current genome assembly gold standard is to generate a long read assembly followed by chromatin conformation capture sequencing to scaffold. We used the PacBio RSII and chromatin conformation capture sequencing to sequence and assemble the genome of a South East Asian commercial sugarcane cultivar, known as Khon Kaen 3. The Khon Kaen 3 genome assembled into 104,477 contigs totalling 7 Gb, which scaffolded into 56 pseudochromosomes containing 5.2 Gb of sequence. Genome annotation produced 242,406 genes from 30,927 orthogroups. Aligning the Khon Kaen 3 genome sequence to S. officinarum and S. spontaneum revealed a high level of apparent recombination, indicating a chimeric assembly. This assembly error is explained by high nucleotide identity between S. officinarum and S. spontaneum, where 91.8% of S. spontaneum aligns to S. officinarum at 94% identity. Thus, the subgenomes of commercial sugarcane are so similar that using short reads to correct long PacBio reads produced chimeric long reads. Future attempts to sequence sugarcane must take this information into account.

Single nucleotide polymorphism marker development in the rubber tree, Hevea brasiliensis (Euphorbiaceae).

PREMISE OF THE STUDY: We demonstrated the application of high-throughput 454 sequencing technology in the identification of single nucleotide polymorphism (SNP) markers in Hevea brasiliensis. METHODS AND RESULTS: A total of 5883 putative SNP positions were discovered in silico, and 10 biallelic SNP markers were validated from 454-derived EST sequences. The polymorphism information content (PIC) and the observed heterozygosity (H(o)) ranged from 0.0963-0.5135 and 0.1071-0.4643, respectively. CONCLUSIONS: These markers can be useful for the construction of genetic maps, the identification of quantitative trait loci linked to commercially desirable traits, and the study of genetic structure in H. brasiliensis.

Complete chloroplast genome sequences of five Bruguiera species (Rhizophoraceae): comparative analysis and phylogenetic relationships.

Bruguiera is a genus of true mangroves that are mostly distributed in the Indo-West Pacific region. However, the number of published whole chloroplast genome sequences of Bruguiera species are limited. Here, the complete chloroplast sequences of five Bruguiera species were sequenced and assembled using Illumina data. The chloroplast genomes of B. gymnorhiza, B. hainesii, B. cylindrica, B. parviflora and B. sexangula were assembled into 161,195, 164,295, 164,297, 163,228 and 164,170 bp, respectively. All chloroplast genomes contain 37 tRNA and eight rRNA genes, with either 84 or 85 protein-coding genes. A comparative analysis of these genomes revealed high similarity in gene structure, gene order and boundary position of the LSC, SSC and two IR regions. Interestingly, B. gymnorhiza lost a rpl32 gene in the SSC region. In addition, a ndhF gene in B. parviflora straddles both the SSC and IRB boundary regions. These genes reveal differences in chloroplast evolution among Bruguiera species. Repeats and SSRs in the chloroplast genome sequences were found to be highly conserved between B. cylindrica and B. hainesii as well as B. gymnorhiza and B. sexangula indicating close genetic relationships based on maternal inheritance. Notably, B. hainesii, which is considered a hybrid between B. gymnorhiza and B. cylindrica, appears to have inherited the chloroplast from B. cylindrica. Investigating the effects of selection events on shared protein-coding genes showed a positive selection in rps7 and rpl36 genes in all species compared to land-plant species. A phylogenetic analysis, based on 59 conserved chloroplast protein-coding genes, showed strong support that all Bruguiera species are in the clade Rhizophoraceae. This study provides valuable genetic information for the study of evolutionary relationships and population genetics in Bruguiera and other mangrove species.

De novo assemblies of Luffa acutangula and Luffa cylindrica genomes reveal an expansion associated with substantial accumulation of transposable elements.

Luffa spp. (sponge gourd or ridge gourd) is an economically important vegetable crop widely cultivated in China, India and Southeast Asia. Here, we employed PacBio long-read single-molecule real-time (SMRT) sequencing to perform de novo genome assemblies of two commonly cultivated Luffa species, L. acutangula and L. cylindrica. We obtained preliminary draft genomes of 734.6 Mb and 689.8 Mb with scaffold N50 of 786,130 and 578,616 bases for L. acutangula and L. cylindrica, respectively. We also applied long-range Chicago and HiC techniques to obtain the first chromosome-scale whole-genome assembly of L. acutangula. The final assembly contained 13 pseudomolecules, corresponding to the haploid chromosome number in Luffa spp. (1n = 13, 2n = 26). The sizes of the assembled Luffa genomes are approximately twice as large as the genome assemblies of related Cucurbitaceae. A large proportion of L. acutangula (62.17%; 456.69 Mb) and L. cylindrica (56.78%; 391.65 Mb) genome assemblies contained repetitive elements. Phylogenetic analyses revealed that the substantial accumulation of transposable elements likely contributed to the expansion of the Luffa genomes. We also investigated alternative splicing events in Luffa using full-length transcript sequences obtained from PacBio Isoform Sequencing (Iso-seq). While the predominant form of alternative splicing in most plant species examined was intron retention, alternative 3' acceptor site selection appeared to be a major event observed in Luffa. High-quality genome assemblies for L. acutangula and L. cylindrica reported here provide valuable resources for Luffa breeding and future genetics and comparative genomics studies in Cucurbitaceae.

A chromosome-scale assembly of the black gram (Vigna mungo) genome.

Black gram (Vigna mungo) is an important short duration grain legume crop. Black gram seeds provide an inexpensive source of dietary protein. Here, we applied the 10X Genomics linked-read technology to obtain a de novo whole genome assembly of V. mungo cultivated variety Chai Nat 80 (CN80). The preliminary assembly contained 12,228 contigs and had an N50 length of 5.2 Mb. Subsequent scaffolding using the long-range Chicago and HiC techniques yielded the first high-quality, chromosome-level assembly of 499 Mb comprising 11 pseudomolecules. Comparative genomics analyses based on sequence information from single-copy orthologous genes revealed that black gram and mungbean (Vigna radiata) diverged about 2.7 million years ago . The transversion rate (4DTv) analysis in V. mungo revealed no evidence supporting a recent genome-wide duplication event observed in the tetraploid créole bean (Vigna reflexo-pilosa). The proportion of repetitive elements in the black gram genome is slightly lower than the numbers reported for related Vigna species. The majority of long terminal repeat retrotransposons appeared to integrate into the genome within the last five million years. We also examined alternative splicing events in V. mungo using full-length transcript sequences. While intron retention was the most prevalent mode of alternative splicing in several plant species, alternative 3' acceptor site selection represented the majority of events in black gram. Our high-quality genome assembly along with the genomic variation information from the germplasm provides valuable resources for accelerating the development of elite varieties through marker-assisted breeding and for future comparative genomics and phylogenetic studies in legume species.

The complete mitochondrial genome of Luffa acutangula.

Based on PacBio de novo assembly, we report the first complete mitochondrial genome of Luffa acutangula (460,333 bp) containing nine large chloroplast-derived sequences (1.9-17.3 kb) across the mitogenome. The base composition of the mitogenome in descending order is A: 28.02%, C: 22.04%, G: 21.83% and T: 28.10%, and the G + C content is 43.87%. There are 63 mitochondrial genes including 40 protein-coding genes, 3 rRNA genes and 20 tRNA genes. Additionally, a total of 288 repeats ranging from 31 to 5,301 bp were identified, accounting for 5.7% of the mitogenome. Two large direct repeats (5,301 and 405 bp) within the mitogenome were found for the formation of four subgenomic molecules. A phylogenetic analysis showed that L. acutangula was closely related to other species in Cucurbiaceae. This mitogenome provides useful genetic information for evolutionary studies.