Atractylodes lancea (Thunb.) DC. is a traditional Chinese medicine rich in sesquiterpenes that has been widely used in China and Japan for the treatment of viral infections. Despite its important pharmacological value, genomic information regarding A. lancea is currently unavailable. In this study, the whole genome sequence of A. lancea was obtained using an Illumina sequencing platform. The results revealed an estimated genome size for A. lancea of 4,159.24 Mb, with 2.28% heterozygosity, and a repeat rate of 89.2%, all of which indicate a highly heterozygous genome. Based on the genomic data of A. lancea, 27,582 simple sequence repeat (SSR) markers were identified. The differences in representation among nucleotide repeat types were large, e.g., the mononucleotide repeat type was the most abundant (54.74%) while the pentanucleotide repeats were the least abundant (0.10%), and sequence motifs GA/TC (31.17%) and TTC/GAA (7.23%) were the most abundant among the dinucleotide and trinucleotide repeat motifs, respectively. A total of 93,434 genes matched known genes in common databases including 48,493 genes in the GO database and 34,929 genes in the KEGG database. This is the first report to sequence and characterize the whole genome of A. lancea and will provide a theoretical basis and reference for further genome-wide deep sequencing and SSR molecular marker development of A. lancea.

Introduction: The uridine diphosphate (UDP)-glycosyltransferase (UGT) family is the largest glycosyltransferase family, which is involved in the biosynthesis of natural plant products and response to abiotic stress. UGT has been studied in many medicinal plants, but there are few reports on Platycodon grandiflorus. This study is devoted to genome-wide analysis of UGT family and identification of UGT genes involved in drought stress of Platycodon grandiflorus (PgUGTs).Methods: The genome data of Platycodon grandiflorus was used for genome-wide identification of PgUGTs, online website and bioinformatics analysis software was used to conduct bioinformatics analysis of PgUGT genes and the genes highly responsive to drought stress were screened out by qRT-PCR, these genes were cloned and conducted bioinformatics analysis.Results: A total of 75 PgUGT genes were identified in P.grandiflorus genome and clustered into 14 subgroups. The PgUGTs were distributed on nine chromosomes, containing multiple cis-acting elements and 22 pairs of duplicate genes were identified. Protein-protein interaction analysis was performed to predict the interaction between PgUGT proteins. Additionally, six genes were upregulated after 3d under drought stress and three genes (PGrchr09G0563, PGrchr06G0523, PGrchr06G1266) responded significantly to drought stress, as confirmed by qRT-PCR. This was especially true for PGrchr06G1266, the expression of which increased 16.21-fold after 3d of treatment. We cloned and conducted bioinformatics analysis of three candidate genes, both of which contained conserved motifs and several cis-acting elements related to stress response, PGrchr06G1266 contained the most elements.Discussion: PgGT1 was confirmed to catalyze the C-3 position of platycodin D and only eight amino acids showed differences between gene PGr008G1527 and PgGT1, which means PGr008G1527 may be able to catalyze the C-3 position of platycodin D in the same manner as PgGT1. Seven genes were highly expressed in the roots, stems, and leaves, these genes may play important roles in the development of the roots, stems, and leaves of P. grandiflorus. Three genes were highly responsive to drought stress, among which the expression of PGrchr06G1266 was increased 16.21-fold after 3d of drought stress treatment, indicating that PGrchr06G1266 plays an important role in drought stress tolerance. To summarize, this study laied the foundation to better understand the molecular bases of responses to drought stress and the biosynthesis of platycodin.