Eustoma grandiflorum (Raf.) Shinn. is an annual herbaceous plant native to the southern United States, Mexico, and the Greater Antilles. It has a large flower with a variety of colors and is an important flower crop. In this study, we established a chromosome-scale de novo assembly of E. grandiflorum genome sequences by integrating four genomic and genetic approaches: (1) Pacific Biosciences (PacBio) Sequel deep sequencing, (2) error correction of the assembly by Illumina short reads, (3) scaffolding by chromatin conformation capture sequencing (Hi-C), and (4) genetic linkage maps derived from an F2 mapping population. Thirty-six pseudomolecules and 64 unplaced scaffolds were created, with a total length of 1,324.8 Mb. A total of 36,619 genes were predicted on the genome as high-confidence genes. A comparison of genome structure between E. grandiflorum and C. canephora or O. pumila suggested whole genome duplication after divergence between the families Gentianaceae and Rubiaceae. Phylogenetic analysis with single-copy genes suggested that the divergence time between Gentianaceae and Rubiaceae was 74.94 MYA. Genetic diversity analysis was performed for nine commercial E. grandiflorum varieties bred in Japan, from which 254,205 variants were identified. This first report of the construction of a reference genome sequence in genus Eustoma is expected to contribute to genetic and genomic studies in this genus and in the family Gentianaceae.

Eustoma grandiflorum (Raf.) Shinn., is an annual herbaceous plant native to the southern United States, Mexico, and the Greater Antilles. It has a large flower with a variety of colors and an important flower crop. In this study, we established a chromosome-scale de novo assembly of E. grandiflorum by integrating four genomic and genetic approaches: (1) Pacific Biosciences (PacBio) Sequel deep sequencing, (2) error correction of the assembly by Illumina short reads, (3) scaffolding by chromatin conformation capture sequencing (Hi-C), and (4) genetic linkage maps derived from an F2 mapping population. The 36 pseudomolecules and unplaced 64 scaffolds were created with total length of 1,324.8 Mb. Full-length transcript sequencing was obtained by PacBio Iso-Seq sequencing for gene prediction on the assembled genome, Egra_v1. A total of 36,619 genes were predicted on the genome as high confidence HC) genes. Of the 36,619, 25,936 were annotated functions by ZenAnnotation. Genetic diversity analysis was also performed for nine commercial E. grandiflorum varieties bred in Japan, and 254,205 variants were identified. This is the first report of the construction of reference genome sequences in E. grandiflorum as well as in the family Gentianaceae.