Transcriptome sequencing of Pinus kesiya var. langbianensis and comparative analysis in the Pinus phylogeny

Comparative transcriptomics reveals potential regulators of climate adaptation in Santalum album L. (Indian Sandalwood)

Santalum album L. (Indian Sandalwood), a valued tree species known for its fragrant heartwood and essential oil is facing increasing threat due to severe anthropogenic pressures compounded by climate change which has resulted in depletion of its adaptive gene pool. The present study investigates the transcriptome-level responses of nine sandalwood genotypes sourced from diverse climatic zones to identify adaptive genes in the species. Comparative transcriptomics predicted 727, 1141 and 479 differentially expressed transcripts (DETs) across wet vs. dry; monsoon vs. dry and wet vs. monsoon conditions, respectively, and majority of DETs were up-regulated in samples sourced from high rainfall areas. Transcripts including heat shock proteins, Zinc finger binding protein, ribosomal proteins, transcription factors and protein kinase were identified as probable regulators of climate adaptation in S. album. The expression changes of eight selected transcripts were further validated by real-time quantitative PCR. Protein-protein interaction analysis revealed key hub transcripts involved in climate response, while alternative splicing events in transcripts such as SURP and G-patch domain-containing protein 1-like protein, G-type lectin S-receptor-like serine/threonine protein kinase B120, Tetraspanin-3 and ARM repeat superfamily protein indicated the probable role of alternate splicing in increasing the transcript diversity during adaptation. This study presents the first insight into the molecular mechanisms of climate adaptation in the species and can form the basis for specific interventions such as selective breeding, genetic manipulation, and habitat management for conservation and long-term survival of sandalwood.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-025-04218-4.

Full-length transcriptome characterization and comparative analysis of Gleditsia sinensis

As an economically important tree, Gleditsia sinensis Lam. is widely planted. A lack of background genetic information on G. sinensis hinders molecular breeding. Based on PacBio single-molecule real-time (SMRT) sequencing and analysis of G. sinensis, a total of 95,183 non-redundant transcript sequences were obtained, of which 93,668 contained complete open reading frames (ORFs), 2,858 were long non-coding RNAs (LncRNAs) and 18,855 alternative splicing (AS) events were identified. Genes orthologous to different Gleditsia species pairs were identified, stress-related genes had been positively selected during the evolution. AGA, AGG, and CCA were identified as the universal optimal codon in the genus of Gleditsia. EIF5A was selected as a suitable fluorescent quantitative reference gene. 315 Cytochrome P450 monooxygenases (CYP450s) and 147 uridine diphosphate (UDP)-glycosyltransferases (UGTs) were recognized through the PacBio SMRT transcriptome. Randomized selection of GsIAA14 for cloning verified the reliability of the PacBio SMRT transcriptome assembly sequence. In conclusion, the research data lay the foundation for further analysis of the evolutionary mechanism and molecular breeding of Gleditsia.

Identification and Analysis of the CBF Gene Family in Three Species of Acer under Cold Stress

The C-Repeat Binding Factor (CBF) gene family has been identified and characterized in multiple plant species, and it plays a crucial role in responding to low temperatures. Presently, only a few studies on tree species demonstrate the mechanisms and potential functions of CBFs associated with cold resistance, while our study is a novel report on the multi-aspect differences of CBFs among three tree species, compared to previous studies. In this study, genome-wide identification and analysis of the CBF gene family in Acer truncatum, Acer pseudosieboldianum, and Acer yangbiense were performed. The results revealed that 16 CBF genes (five ApseCBFs, four AcyanCBFs, and seven AtruCBFs) were unevenly distributed across the chromosomes, and most CBF genes were mapped on chromosome 2 (Chr2) and chromosome 11 (Chr11). The analysis of phylogenetic relationships, gene structure, and conserved motif showed that 16 CBF genes could be clustered into three subgroups; they all contained Motif 1 and Motif 5, and most of them only spanned one exon. The cis-acting elements analysis showed that some CBF genes might be involved in hormone and abiotic stress responsiveness. In addition, CBF genes exhibited tissue expression specificity. High expressions of ApseCBF1, ApseCBF3, AtruCBF1, AtruCBF4, AtruCBF6, AtruCBF7, and ApseCBF3, ApseCBF4, ApseCBF5 were detected on exposure to low temperature for 3 h and 24 h. Low expressions of AtruCBF2, AtruCBF6, AtruCBF7 were detected under cold stress for 24 h, and AtruCBF3 and AtruCBF5 were always down-regulated under cold conditions. Taken together, comprehensive analysis will enhance our understanding of the potential functions of the CBF genes on cold resistance, thereby providing a reference for the introduction of Acer species in our country.

Comprehensive collection of genes and comparative analysis of full-length transcriptome sequences from Japanese larch (Larix kaempferi) and Kuril larch (Larix gmelinii var. japonica)

BACKGROUND

Japanese larch (Larix kaempferi) is an economically important deciduous conifer species that grows in cool-temperate forests and is endemic to Japan. Kuril larch (L. gmelinii var. japonica) is a variety of Dahurian larch that is naturally distributed in the Kuril Islands and Sakhalin. The hybrid larch (L. gmelinii var. japonica × L. kaempferi) exhibits heterosis, which manifests as rapid juvenile growth and high resistance to vole grazing. Since these superior characteristics have been valued by forestry managers, the hybrid larch is one of the most important plantation species in Hokkaido. To accelerate molecular breeding in these species, we collected and compared full-length cDNA isoforms (Iso-Seq) and RNA-Seq short-read, and merged them to construct candidate gene as reference for both Larix species. To validate the results, candidate protein-coding genes (ORFs) related to some flowering signal-related genes ​were screened from the reference sequences, and the phylogenetic relationship with closely related species was elucidated.

RESULTS

Using the isoform sequencing of PacBio RS ll and the de novo assembly of RNA-Seq short-read sequences, we identified 50,690 and 38,684 ORFs in Japanese larch and Kuril larch, respectively. BUSCO completeness values were 90.5% and 92.1% in the Japanese and Kuril larches, respectively. After comparing the collected ORFs from the two larch species, a total of 19,813 clusters, comprising 22,571 Japanese larch ORFs and 22,667 Kuril larch ORFs, were contained in the intersection of the Venn diagram. In addition, we screened several ORFs related to flowering signals (SUPPRESSER OF OVEREXPRESSION OF CO1: SOC1, LEAFY: LFY, FLOWERING Locus T: FT, CONSTANCE: CO) from both reference sequences, and very similar found in other species.

CONCLUSIONS

The collected ORFs will be useful as reference sequences for molecular breeding of Japanese and Kuril larches, and also for clarifying the evolution of the conifer genome and investigating functional genomics.

Comparative Transcriptomics and Metabolites Analysis of Two Closely Related Euphorbia Species Reveal Environmental Adaptation Mechanism and Active Ingredients Difference

Roots of Euphorbia fischeriana and Euphorbia ebracteolata are recorded as the source plant of traditional Chinese medicine "Langdu," containing active ingredients with anticancer and anti-AIDS activity. However, the two species have specific patterns in the graphic distribution. Compared with E. ehracteolata, E. fischeriana distributes in higher latitude and lower temperature areas and might have experienced cold stress adaptation. To reveal the molecular mechanism of environmental adaptation, RNA-seq was performed toward the roots, stems, and leaves of E. fischeriana and E. ehracteolata. A total of 6,830 pairs of putative orthologs between the two species were identified. Estimations of non-synonymous or synonymous substitution rate ratios for these orthologs indicated that 533 of the pairs may be under positive selection (Ka/Ks > 0.5). Functional enrichment analysis revealed that significant proportions of the orthologs were in the TCA cycle, fructose and mannose metabolism, starch and sucrose metabolism, fatty acid biosynthesis, and terpenoid biosynthesis providing insights into how the two closely related Euphorbia species adapted differentially to extreme environments. Consistent with the transcriptome, a higher content of soluble sugars and proline was obtained in E. fischeriana, reflecting the adaptation of plants to different environments. Additionally, 5 primary or secondary metabolites were screened as the biomarkers to distinguish the two species. Determination of 4 diterpenoids was established and performed, showing jolkinolide B as a representative component in E. fischeriana, whereas ingenol endemic to E. ebracteolate. To better study population genetics, EST-SSR markers were generated and tested in 9 species of Euphorbia. A total of 33 of the 68 pairs were screened out for producing clear fragments in at least four species, which will furthermore facilitate the studies on the genetic improvement and phylogenetics of this rapidly adapting taxon. In this study, transcriptome and metabolome analyses revealed the evolution of genes related to cold stress tolerance, biosynthesis of TCA cycle, soluble sugars, fatty acids, and amino acids, consistent with the molecular strategy that genotypes adapting to environment. The key active ingredients of the two species were quantitatively analyzed to reveal the difference in pharmacodynamic substance basis and molecular mechanism, providing insights into rational crude drug use.

Applications of machine learning in pine nuts classification

Pine nuts are not only the important agent of pine reproduction and afforestation, but also the commonly consumed nut with high nutritive values. However, it is difficult to distinguish among pine nuts due to the morphological similarity among species. Therefore, it is important to improve the quality of pine nuts and solve the adulteration problem quickly and non-destructively. In this study, seven pine nuts (Pinus bungeana, Pinus yunnanensis, Pinus thunbergii, Pinus armandii, Pinus massoniana, Pinus elliottii and Pinus taiwanensis) were used as study species. 210 near-infrared (NIR) spectra were collected from the seven species of pine nuts, five machine learning methods (Decision Tree (DT), Random Forest (RF), Multilayer Perceptron (MLP), Support Vector Machine (SVM) and Naive Bayes (NB)) were used to identify species of pine nuts. 303 images were used to collect morphological data to construct a classification model based on five convolutional neural network (CNN) models (VGG16, VGG19, Xception, InceptionV3 and ResNet50). The experimental results of NIR spectroscopy show the best classification model is MLP and the accuracy is closed to 0.99. Another experimental result of images shows the best classification model is InceptionV3 and the accuracy is closed to 0.964. Four important range of wavebands, 951–957 nm, 1,147–1,154 nm, 1,907–1,927 nm, 2,227–2,254 nm, were found to be highly related to the classification of pine nuts. This study shows that machine learning is effective for the classification of pine nuts, providing solutions and scientific methods for rapid, non-destructive and accurate classification of different species of pine nuts.

Morphological and Comparative Transcriptome Analysis of Three Species of Five-Needle Pines: Insights Into Phenotypic Evolution and Phylogeny

Pinus koraiensis, Pinus sibirica, and Pinus pumila are the major five-needle pines in northeast China, with substantial economic and ecological values. The phenotypic variation, environmental adaptability and evolutionary relationships of these three five-needle pines remain largely undecided. It is therefore important to study their genetic differentiation and evolutionary history. To obtain more genetic information, the needle transcriptomes of the three five-needle pines were sequenced and assembled. To explore the relationship of sequence information and adaptation to a high mountain environment, data on needle morphological traits [needle length (NL), needle width (NW), needle thickness (NT), and fascicle width (FW)] and 19 climatic variables describing the patterns and intensity of temperature and precipitation at six natural populations were recorded. Geographic coordinates of altitude, latitude, and longitude were also obtained. The needle morphological data was combined with transcriptome information, location, and climate data, for a comparative analysis of the three five-needle pines. We found significant differences for needle traits among the populations of the three five-needle pine species. Transcriptome analysis showed that the phenotypic variation and environmental adaptation of the needles of P. koraiensis, P. sibirica, and P. pumila were related to photosynthesis, respiration, and metabolites. Analysis of orthologs from 11 Pinus species indicated a closer genetic relationship between P. koraiensis and P. sibirica compared to P. pumila. Our study lays a foundation for genetic improvement of these five-needle pines and provides insights into the adaptation and evolution of Pinus species.

Inferring the Phylogenetic Positions of Two Fig Wasp Subfamilies of Epichrysomallinae and Sycophaginae Using Transcriptomes and Mitochondrial Data

Fig wasps are a group of insects (Hymenoptera: Chalcidoidea) that live in the compact syconia of fig trees (Moraceae: Ficus). Accurate classification and phylogenetic results are very important for studies of fig wasps, but the taxonomic statuses of some fig wasps, especially the non-pollinating subfamilies are difficult to determine, such as Epichrysomallinae and Sycophaginae. To resolve the taxonomic statuses of Epichrysomallinae and Sycophaginae, we obtained transcriptomes and mitochondrial genome (mitogenome) data for four species of fig wasps. These newly added data were combined with the data of 13 wasps (data on 11 fig wasp species were from our laboratory and two wasp species were download from NCBI). Based on the transcriptome and genome data, we obtained 145 single-copy orthologous (SCO) genes in 17 wasp species, and based on mitogenome data, we obtained 13 mitochondrial protein-coding genes (PCGs) for each of the 17 wasp species. Ultimately, we used 145 SCO genes, 13 mitochondrial PCGs and combined SCO genes and mitochondrial genes data to reconstruct the phylogenies of fig wasps using both maximum likelihood (ML) and Bayesian inference (BI) analyses. Our results suggest that both Epichrysomallinae and Sycophaginae are more closely related to Agaonidae with a high statistical support.

Pharmaceutical resource discovery from traditional medicinal plants: Pharmacophylogeny and pharmacophylogenomics

The worldwide botanical and medicinal culture diversity are astonishing and constitute a Pierian spring for innovative drug R&D. Here, the latest awareness and the perspectives of pharmacophylogeny and pharmacophylogenomics, as well as their expanding utility in botanical drug R&D, are systematically summarized and highlighted. Chemotaxonomy is based on the fact that closely related plants contain the same or similar chemical profiles. Correspondingly, it is better to combine morphological characters, DNA markers and chemical markers in the inference of medicinal plant phylogeny. Medicinal plants within the same phylogenetic groups may have the same or similar therapeutic effects, thus forming the core of pharmacophylogeny. Here we systematically review and comment on the versatile applications of pharmacophylogeny in (1) looking for domestic resources of imported drugs, (2) expanding medicinal plant resources, (3) quality control, identification and expansion of herbal medicines, (4) predicting the chemical constituents or active ingredients of herbal medicine and assisting in the identification and determination of chemical constituents, (5) the search for new drugs sorting out, and (6) summarizing and improving herbal medicine experiences, etc. Such studies should be enhanced within the context of deeper investigations of molecular biology and genomics of traditional medicinal plants, phytometabolites and metabolomics, and ethnomedicine-based pharmacological activity, thus enabling the sustainable conservation and utilization of traditional medicinal resources.

The complete chloroplast genome sequence of Pinus kesiya var. langbianensis

Pinus kesiya var. langbianensis is an important oleoresin resource tree species in Yunnan province, China. In this study, the complete chloroplast genome (cpDNA) sequence of P. kesiya var. langbianensis was determined from Illumina HiSeq pair-end sequencing data. The cpDNA is 119,780 bp in length, contains a large single copy region (LSC) of 65,863 bp and a small single copy region (SSC) of 53,117 bp, which were separated by a pair of inverted repeat (IR) regions of 400 bp. The genome contains 113 genes, including 73 protein-coding genes, 4 ribosomal RNA genes, and 36 transfer RNA genes. The overall GC content of the whole genome is 38.5%, and the corresponding values of the LSC, SSC, and IR regions are 37.9, 39.4, and 35.5%, respectively. Further phylogenomic analysis showed that P. kesiya var. langbianensis clustered together with Pinus kesiya, Pinus densata, and Pinus taiwanensis.

Comparative genomics and transcriptomics analysis reveals evolution patterns of selection in the Salix phylogeny

BACKGROUND

Willows are widely distributed in the northern hemisphere and have good adaptability to different living environment. The increasing of genome and transcriptome data provides a chance for comparative analysis to study the evolution patterns with the different origin and geographical distributions in the Salix phylogeny.

RESULTS

Transcript sequences of 10 Salicaceae species were downloaded from public databases. All pairwise of orthologues were identified by comparative analysis in these species, from which we constructed a phylogenetic tree and estimated the rate of diverse. Divergence times were estimated in the 10 Salicaceae using comparative transcriptomic analysis. All of the fast-evolving positive selection sequences were identified, and some cold-, drought-, light-, universal-, and heat- resistance genes were discovered.

CONCLUSIONS

The divergence time of subgenus Vetrix and Salix was about 17.6-16.0 Mya during the period of Middle Miocene Climate Transition (21-14 Mya). Subgenus Vetrix diverged to migratory and resident groups when the climate changed to the cool and dry trend by 14 Mya. Cold- and light- stress genes were involved in positive selection among the resident Vetrix, and which would help them to adapt the cooling stage. Universal- stress genes exhibited positive selection among the migratory group and subgenus Salix. These data are useful for comprehending the adaptive evolution and speciation in the Salix lineage.