Thirteen Dipterocarpoideae genomes provide insights into their evolution and borneol biosynthesis.

Dipterocarpoideae, the largest subfamily of the Dipterocarpaceae, is a dominant component of Southeast Asian rainforests and is widely used as a source of wood, damar resin, medicine, and essential oil. However, many Dipterocarpoideae species are currently on the IUCN Red List owing to severe degradation of their habitats under global climate change and human disturbance. Genetic information regarding these taxa has only recently been reported with the sequencing of four Dipterocarp genomes, providing clues to the function and evolution of these species. Here, we report on 13 high-quality Dipterocarpoideae genome assemblies, ranging in size from 302.6 to 494.8 Mb and representing the five most species-rich genera in Dipterocarpoideae. Molecular dating analyses support the Western Gondwanaland origin of Dipterocarpaceae. Based on evolutionary analysis, we propose a three-step chromosome evolution scenario to describe the karyotypic evolution from an ancestor with six chromosomes to present-day species with 11 and 7 chromosomes. We discovered an expansion of genes encoding cellulose synthase (CesA), which is essential for cellulose biosynthesis and secondary cell-wall formation. We functionally identified five bornyl diphosphate synthase (BPPS) genes, which specifically catalyze the biosynthesis of borneol, a natural medicinal compound extracted from damar resin and oils, thus providing a basis for large-scale production of natural borneol in vitro.

Comparison of ONT and CCS sequencing technologies on the polyploid genome of a medicinal plant showed that high error rate of ONT reads are not suitable for self-correction.

BACKGROUND: Many medicinal plants are known for their complex genomes with high ploidy, heterozygosity, and repetitive content which pose severe challenges for genome sequencing of those species. Long reads from Oxford nanopore sequencing technology (ONT) or Pacific Biosciences Single Molecule, Real-Time (SMRT) sequencing offer great advantages in de novo genome assembly, especially for complex genomes with high heterozygosity and repetitive content. Currently, multiple allotetraploid species have sequenced their genomes by long-read sequencing. However, we found that a considerable proportion of these genomes (7.9% on average, maximum 23.7%) could not be covered by NGS (Next Generation Sequencing) reads (uncovered region by NGS reads, UCR) suggesting the questionable and low-quality of those area or genomic areas that can't be sequenced by NGS due to sequencing bias. The underlying causes of those UCR in the genome assembly and solutions to this problem have never been studied. METHODS: In the study, we sequenced the tetraploid genome of Veratrum dahuricum (Turcz.) O. Loes (VDL), a Chinese medicinal plant, with ONT platform and assembled the genome with three strategies in parallel. We compared the qualities, coverage, and heterozygosity of the three ONT assemblies with another released assembly of the same individual using reads from PacBio circular consensus sequencing (CCS) technology, to explore the cause of the UCR. RESULTS: By mapping the NGS reads against the three ONT assemblies and the CCS assembly, we found that the coverage of those ONT assemblies by NGS reads ranged from 49.15 to 76.31%, much smaller than that of the CCS assembly (99.53%). And alignment between ONT assemblies and CCS assembly showed that most UCR can be aligned with CCS assembly. So, we conclude that the UCRs in ONT assembly are low-quality sequences with a high error rate that can't be aligned with short reads, rather than genomic regions that can't be sequenced by NGS. Further comparison among the intermediate versions of ONT assemblies showed that the most probable origin of those errors is a combination of artificial errors introduced by "self-correction" and initial sequencing error in long reads. We also found that polishing the ONT assembly with CCS reads can correct those errors efficiently. CONCLUSIONS: Through analyzing genome features and reads alignment, we have found the causes for the high proportion of UCR in ONT assembly of VDL are sequencing errors and additional errors introduced by self-correction. The high error rates of ONT-raw reads make them not suitable for self-correction prior to allotetraploid genome assembly, as the self-correction will introduce artificial errors to > 5% of the UCR sequences. We suggest high-precision CCS reads be used to polish the assembly to correct those errors effectively for polyploid genomes.

Comparative and phylogenetic analyses of eleven complete chloroplast genomes of Dipterocarpoideae.

BACKGROUND: In South-east Asia, Dipterocarpoideae is predominant in most mature forest communities, comprising around 20% of all trees. As large quantity and high quality wood are produced in many species, Dipterocarpoideae plants are the most important and valuable source in the timber market. The d-borneol is one of the essential oil components from Dipterocarpoideae (for example, Dryobalanops aromatica or Dipterocarpus turbinatus) and it is also an important traditional Chinese medicine (TCM) formulation known as "Bingpian" in Chinese, with antibacterial, analgesic and anti-inflammatory effects and can enhance anticancer efficiency. METHODS: In this study, we analyzed 20 chloroplast (cp) genomes characteristics of Dipterocarpoideae, including eleven newly reported genomes and nine cp genomes previously published elsewhere, then we explored the chloroplast genomic features, inverted repeats contraction and expansion, codon usage, amino acid frequency, the repeat sequences and selective pressure analyses. At last, we constructed phylogenetic relationships of Dipterocarpoideae and found the potential barcoding loci. RESULTS: The cp genome of this subfamily has a typical quadripartite structure and maintains a high degree of consistency among species. There were slightly more tandem repeats in cp genomes of Dipterocarpus and Vatica, and the psbH gene was subjected to positive selection in the common ancestor of all the 20 species of Dipterocarpoideae compared with three outgroups. Phylogenetic tree showed that genus Shorea was not a monophyletic group, some Shorea species and genus Parashorea are placed in one clade. In addition, the rpoC2 gene can be used as a potential marker to achieve accurate and rapid species identification in subfamily Dipterocarpoideae. CONCLUSIONS: Dipterocarpoideae had similar cp genomic features and psbM, rbcL, psbH may function in the growth of Dipterocarpoideae. Phylogenetic analysis suggested new taxon treatment is needed for this subfamily indentification. In addition, rpoC2 is potential to be a barcoding gene to TCM distinguish.