High-quality assembly of the T2T genome for Isodon rubescens f. lushanensis reveals genomic structure variations between 2 typical forms of Isodon rubescens.

Background: Rabdosiae rubescentis herba (Isodon rubescens) is widely used as a folk medicine to treat esophageal cancer and sore throat in China. Its germplasm resources are abundant in China, with I. rubescens (Hemsl.) Hara and I. rubescens f. lushanensis as 2 typical forms. I. rubescens (Hemsl.) Hara is featured by biosynthesis of the diterpenoid oridonin with strong anticancer activity, while I. rubescens f. lushanensis produces another diterpenoid with anticancer activity, lushanrubescensin. However, the biosynthetic pathways of both still need to be fully understood. In particular, little is known about the genetic background of I. rubescens f. lushanensis.

Findings: We used Pacific Biosciences (PacBio) single-molecule real-time and Nanopore Ultra-long sequencing platforms, respectively, and obtained 139.07 Gb of high-quality data, with a sequencing depth of about 328×. We also obtained a high-quality reference genome for I. rubescens f. lushanensis, with a genome size of 349 Mb and a contig N50 of 28.8 Mb. The heterozygosity of the genome is 1.7% and the repeatability is 83.43%. In total, 34,865 protein-coding genes were predicted. Moreover, we found that most of the variant or unique genes in the diterpenoid synthesis pathways of I. rubescens f. lushanensis and I. rubescens (Hemsl.) Hara were enriched in diterpene synthases.

Conclusions: We provide the first genome sequence and gene annotation for the I. rubescens f. lushanensis, which provides molecular evidence for understanding the chemotypic differences of I. rubescens.