Taxonomy, comparative genomics of Mullein (Verbascum, Scrophulariaceae), with implications for the evolution of Verbascum and Lamiales

The chloroplast genomes of two medicinal species (Veronica anagallis-aquatica L. and Veronica undulata Wall.) and its comparative analysis with related Veronica species

Veronica anagallis-aquatica L. and Veronica undulata Wall. are widely used ethnomedicinal plants in China. The two species have different clinical efficacies, while their extremely similar morphology and unclear interspecific relationship make it difficult to accurately identify them, leading to increased instances of mixed usage. This article reports on the complete chloroplast genomes sequence of these two species and their related Veronica species to conduct a comparative genomics analysis and phylogenetic construction. The results showed that the chloroplast (cp) genomes of Veronica exhibited typical circular quadripartite structures, with total lengths of 149,386 to 152,319 base pairs (bp), and GC content of 37.9 to 38.1%, and the number of genes was between 129-134. The total number of simple sequence repeats (SSRs) in V. anagallis-aquatica and V. undulata is 37 and 36, while V. arvensis had the highest total number of 56, predominantly characterized by A/T single bases. The vast majority of long repeat sequence types are forward repeats and palindromic repeats. Selective Ka/Ks values showed that three genes were under positive selection. Sequence differences often occur in the non-coding regions of the large single-copy region (LSC) and small single-copy region (SSC), with the lowest sequence variation in the inverted repeat regions (IR). Seven highly variable regions (trnT-GGU-psbD, rps8-rpl16, trnQ-UUG, trnN-GUU-ndhF, petL, ycf3, and ycf1) were detected, which may be potential molecular markers for identifying V. anagallis-aquatica and V. undulata. The phylogenetic tree indicates that there is a close genetic relationship between the genera Veronica and Neopicrorhiza, and V. anagallis-aquatica and V. undulata are sister groups. The molecular clock analysis results indicate that the divergence time of Veronica may occur at ∼ 9.09 Ma, and the divergence time of these two species occurs at ∼ 0.48 Ma. It is speculated that climate change may be the cause of Veronica species diversity and promote the radiation of the genus. The chloroplast genome data of nine Veronica specie provides important insights into the characteristics and evolution of the chloroplast genome of this genus, as well as the phylogenetic relationships of the genus Veronica.

Genome assembly of M. spongiola and comparative genomics of the genus Morchella provide initial insights into taxonomy and adaptive evolution

Morchella spongiola is a highly prized mushroom for its delicious flavor and medical value and is one of the most flourishing, representative, and dominant macrofungi in the Qilian Mountains of the Qinghai-Tibet Plateau subkingdoms (QTPs). However, the understanding of M. spongiola remains largely unknown, and its taxonomy is ambiguous. In this study, we redescribed a unique species of M. spongiola, i.e., micromorphology, molecular data, genomics, and comparative genomics, and the historical biogeography of M. spongiola were estimated for 182 single-copy homologous genes. A high-quality chromosome-level reference genome of M. spongiola M12-10 was obtained by combining PacBio HiFi data and Illumina sequencing technologies; it was approximately 57.1 Mb (contig N50 of 18.14 Mb) and contained 9775 protein-coding genes. Comparative genome analysis revealed considerable conservation and unique characteristics between M. spongiola M12-10 and 32 other Morchella species. Molecular phylogenetic analysis indicated that M. spongiola M12-10 is similar to the M. prava/Mes-7 present in sandy soil near rivers, differentiating from black morels ~ 43.06 Mya (million years ago), and diverged from M. parva/Mes-7 at approximately 12.85 Mya (in the Miocene epoch), which is closely related to the geological activities in the QTPs (in the Neogene). Therefore, M. spongiola is a unique species rather than a synonym of M. vulgaris/Mes-5, which has a distinctive grey-brown sponge-like ascomata. This genome of M. spongiola M12-10 is the first published genome sequence of the species in the genus Morchella from the QTPs, which could aid future studies on functional gene identification, germplasm resource management, and molecular breeding efforts, as well as evolutionary studies on the Morchella taxon in the QTPs.

Verbascum species as a new source of saffron apocarotenoids and molecular tools for the biotechnological production of crocins and picrocrocin

Crocins are glucosylated apocarotenoids present in flowers and fruits of a few plant species, including saffron, gardenia, and Buddleja. The biosynthesis of crocins in these plants has been unraveled, and the enzymes engineered for the production of crocins in heterologous systems. Mullein (Verbascum sp.) has been identified as a new source of crocins and picrocrocin. In this work, we have identified eight enzymes involved in the cleavage of carotenoids in two Verbascum species, V. giganteum and V. sinuatum. Four of them were homologous to the previously identified BdCCD4.1 and BdCCD4.3 from Buddleja, involved in the biosynthesis of crocins. These enzymes were analyzed for apocarotenogenic activity in bacteria and Nicotiana benthamiana plants using a virus-driven system. Metabolic analyses of bacterial extracts and N. benthamiana leaves showed the efficient activity of these enzymes to produce crocins using β-carotene and zeaxanthin as substrates. Accumulations of 0.17% of crocins in N. benthamiana dry leaves were reached in only 2 weeks using a recombinant virus expressing VgCCD4.1, similar to the amounts previously produced using the canonical saffron CsCCD2L. The identification of these enzymes, which display a particularly broad substrate spectrum, opens new avenues for apocarotenoid biotechnological production.

Phylogenomics sheds new light on the drivers behind a long-lasting systematic riddle: the figwort family Scrophulariaceae

The figwort family, Scrophulariaceae, comprises c. 2000 species whose evolutionary relationships at the tribal level have proven difficult to resolve, hindering our ability to understand their origin and diversification. We designed a specific probe kit for Scrophulariaceae, targeting 849 nuclear loci and obtaining plastid regions as by-products. We sampled c. 87% of the genera described in the family and use the nuclear dataset to estimate evolutionary relationships, timing of diversification, and biogeographic patterns. Ten tribes, including two new tribes, Androyeae and Camptolomeae, are supported, and the phylogenetic positions of Androya, Camptoloma, and Phygelius are unveiled. Our study reveals a major diversification at c. 60 million yr ago in some Gondwanan landmasses, where two different lineages diversified, one of which gave rise to nearly 81% of extant species. A Southern African origin is estimated for most modern-day tribes, with two exceptions, the American Leucophylleae, and the mainly Australian Myoporeae. The rapid mid-Eocene diversification is aligned with geographic expansion within southern Africa in most tribes, followed by range expansion to tropical Africa and multiple dispersals out of Africa. Our robust phylogeny provides a framework for future studies aimed at understanding the role of macroevolutionary patterns and processes that generated Scrophulariaceae diversity.

Multi-integrated approach for unraveling small open reading frames potentially associated with secondary metabolism in Streptomyces

IMPORTANCE

Due to their small size and special chemical features, small open reading frame (smORF)-encoding peptides (SEPs) are often neglected. However, they may play critical roles in regulating gene expression, enzyme activity, and metabolite production. Studies on bacterial microproteins have mainly focused on pathogenic bacteria, which are importance to systematically investigate SEPs in streptomycetes and are rich sources of bioactive secondary metabolites. Our study is the first to perform a global identification of smORFs in streptomycetes. We established a peptidogenomic workflow for non-model microbial strains and identified multiple novel smORFs that are potentially linked to secondary metabolism in streptomycetes. Our multi-integrated approach in this study is meaningful to improve the quality and quantity of the detected smORFs. Ultimately, the workflow we established could be extended to other organisms and would benefit the genome mining of microproteins with critical functions for regulation and engineering useful microorganisms.

Comparative analyses of Linderniaceae plastomes, with implications for its phylogeny and evolution

Introduction

The recently established Linderniaceae, separated from the traditionally defined Scrophulariaceae, is a taxonomically complicated family. Although previous phylogenetic studies based on a few short DNA markers have made great contributions to the taxonomy of Linderniaceae, limited sampling and low resolution of the phylogenetic tree have failed to resolve controversies between some generic circumscriptions. The plastid genome exhibits a powerful ability to solve phylogenetic relationships ranging from shallow to deep taxonomic levels. To date, no plastid phylogenomic studies have been carried out in Linderniaceae.

Methods

In this study, we newly sequenced 26 plastid genomes of Linderniaceae, including eight genera and 25 species, to explore the phylogenetic relationships and genome evolution of the family through plastid phylogenomic and comparative genomic analyses.

Results

The plastid genome size of Linderniaceae ranged from 152,386 bp to 154,402 bp, exhibiting a typical quartile structure. All plastomes encoded 114 unique genes, comprising 80 protein-coding genes, 30 tRNA genes, and four rRNA genes. The inverted repeat regions were more conserved compared with the single-copy regions. A total of 1803 microsatellites and 1909 long sequence repeats were identified, and five hypervariable regions (petN-psbM, rps16-trnQ, rpl32-trnL, rpl32, and ycf1) were screened out. Most protein-coding genes were relatively conserved, with only the ycf2 gene found under positive selection in a few species. Phylogenomic analyses confirmed that Linderniaceae was a distinctive lineage and revealed that the presently circumscribed Vandellia and Torenia were non-monophyletic.

Discussion

Comparative analyses showed the Linderniaceae plastomes were highly conservative in terms of structure, gene order, and gene content. Combining morphological and molecular evidence, we supported the newly established Yamazakia separating from Vandellia and the monotypic Picria as a separate genus. These findings provide further evidence to recognize the phylogenetic relationships among Linderniaceae and new insights into the evolution of the plastid genomes.

Analysis of complete chloroplast genome sequences and insight into the phylogenetic relationships of Ferula L

Background

Ferula L. is one of the largest and most taxonomically complicated genera as well as being an important medicinal plant resource in the family Apiaceae. To investigate the plastome features and phylogenetic relationships of Ferula and its neighboring genera Soranthus Ledeb., Schumannia Kuntze., and Talassia Korovin, we sequenced 14 complete plastomes of 12 species. 

Results

The size of the 14 complete chloroplast genomes ranged from 165,607 to 167,013 base pairs (bp) encoding 132 distinct genes (87 protein-coding, 37 tRNA, and 8 rRNA genes), and showed a typical quadripartite structure with a pair of inverted repeats (IR) regions. Based on comparative analysis, we found that the 14 plastomes were similar in codon usage, repeat sequence, simple sequence repeats (SSRs), and IR borders, and had significant collinearity. Based on our phylogenetic analyses, Soranthus, Schumannia, and Talassia should be considered synonymous with Ferula. Six highly divergent regions (rps16/trnQ-UUG, trnS-UGA/psbZ, psbH/petB, ycf1/ndhF, rpl32, and ycf1) were also detected, which may represent potential molecular markers, and combined with selective pressure analysis, the weak positive selection gene ccsA may be a discriminating DNA barcode for Ferula species.

Conclusion

Plastids contain abundant informative sites for resolving phylogenetic relationships. Combined with previous studies, we suggest that there is still much room for improvement in the classification of Ferula. Overall, our study provides new insights into the plastome evolution, phylogeny, and taxonomy of this genus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08868-z.