Adaptation Insights from Comparative Transcriptome Analysis of Two Opisthopappus Species in the Taihang Mountains

doi:10.21203/rs.3.rs-1023687/v1

PDF

Research Article

Adaptation Insights from Comparative Transcriptome Analysis of Two Opisthopappus Species in the Taihang Mountains

https://doi.org/10.21203/rs.3.rs-1023687/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Opisthopappus is a major wild source of Asteraceae with good cold and drought resistance. Two species of this genus (Opisthopappus taihangensis and Opisthopappus longilobus) have been employed as model systems to address the evolutionary history of perennial herb biomes in the Taihang Mountains of China. However, further studies on the adaptive divergence processes of these two species are currently impeded by the lack of genomic resources. To elucidate the molecular mechanisms involved, a comparative analysis of these two species was conducted. Among the identified transcription factors, the bHLH members were most prevalent, which exhibited significantly different expression levels in the terpenoid metabolic pathway. O. longilobus revealed a higher expression than did O. taihangensis in terms of terpenes biosynthesis and metabolism, particularly regarding monoterpenoids and diterpenoids. Analyses of the positive selection genes (PSGs) identified from O. taihangensis and O. longilobus, 1203 genes were found that related to adaptative divergence, which were under rapid evolution and/or have signs of positive selection. Different PSG expressions occurred primarily in the mitochondrial electron transport, starch degradation, secondary metabolism, as well as nucleotide synthesis and S-metabolism pathway processes. Two PSGs were obviously differentially expressed in terpenes biosynthesis that might result in the fragrances divergence between O. longilobus and O. taihangensis, which would provide insights as to how the two species adapted to different environments, characterized by sub-humid warm temperate and temperate continental monsoon climates. The comparative analysis for these two species of Opisthopappus not only revealed how the divergence occurred from molecular perspective, but also provided novel insights into how differential adaptations occurred in Taihang Mountains.

Plant Molecular Biology and Genetics

Opisthopappus

comparative transcriptome

differentiation

adaptation

Adaptation is a topic of fundamental interest to researchers in ecology, evolution, and conservation biology [1–4]. When different populations experience heterogeneous environments, natural selection can drive phenotypic divergence while modulating underlying genomic architectures, which induces populations to adapt to their habitats [3, 5]. Adaptative divergence has been well verified as a major mechanism that initiates evolutionary diversification and speciation [6, 7]. There are myriad variations in the combination of phenotypic traits and environmental factors [8]. Variation created by mutation, the raw material for evolutionary change, is translated into phenotypes by flux through metabolic pathways and by the topography and dynamics of molecular networks [9]. Thus, empirical adaptation tests are of principal importance as they explore the balance between the evolutionary processes that shape populations (e.g., strength of selection in promoting local adaptation) relative to non-adaptive or even maladaptive factors such as gene flow, recombination, mutation, and genetic drift [10, 11, 12].

With the development of high-throughput sequencing, researchers have been provided with exciting opportunities to conduct large-scale comparisons at genomic, exomic, or transcriptomic levels [13, 14, 15, 4]. This has significantly increased the capacity to obtain insights into the molecular mechanisms that drive the evolution and differentiation of organisms [16, 17]. Transcriptome sequencing technology is a fast and cost-effective approach that can obtain all of the RNA data of an organism [18]. It has increasingly been employed to identify the functional elements in the genomes, metabolic pathways, and differentially expressed genes (DEGs) of plants [19, 20, 21], particularly those that lack a sequenced reference genome, or non-model species [22].

As growing numbers of species have been included in transcriptome sequencing projects, comparative transcriptomics has received greater attention from researchers [16]. Moreover, comparative transcriptome analysis has become a powerful tool for the investigation of plant growth, development, and physiology at the transcriptional and metabolic levels [23], which provides opportunities to estimate transcriptome-wide divergence and identify loci under selection [16]. Especially, comparative RNA-sequencing studies between intimately related species can, in principle, not only provide additional genomic resources but also offer data related to the processes of speciation or adaptive evolution [15]. Even more, they can assess the functional differences between orthologous genes following divergence [24].

The Opisthopappus Shih genus belonging to the family Asteraceae, which includes two closely related species (Opisthopappus taihangensis Shih and Opisthopappus longilobus (Ling) Shih) [25]. These two species generally grow within the cliff cracks and rock gaps of the Taihang Mountains in China. They present similar morphological characteristics; the two diagnostic attributes used to distinguish them including the style of pinnatisect leaves and the presence or absence of bracteal leaves [26, 27]. For O. longilobus, most stems and leaves are pinnatifid with one pair of bracteal leaves; O. taihangensis has bipinnatifid stems and leaves with no bracteal leaves [27]. Based on our observations in a common-garden, the stems of O. longilobus are more brittle than those of O. taihangensis, and O. longilobus individuals are relatively smaller than those of O. taihangensis. Further, O. taihangensis more easily survives transplantation into a common garden than does O. longilobus. Meanwhile, both O. taihangensis and O. longilobus can exude unique fragrances. Using HPLC (High Performance Liquid Chromatography) and FTIR (Fourier Transform Infrared Spectrometery), it was found that the chemical compositions and contents of chlorogenic acid, apigenin, and quercetin were different between the two species (data not published).

During the long-term evolution of Opisthopappus, the dramatic geological changes and significant topographical barriers of Taihang Mountains could limit the large-scale population expansion and gradually introduce heterogeneity within different populations of O. taihangensis and O. longilobus [28]. O. taihangensis and O. longilobus likely diverged during the early Miocene under a strengthening East Asian monsoon [29]. Subsequently, the dramatic climatic transformation and complex topography of the Taihang Mountains further promoted variations between the two species and their populations [29]. At the present, O. taihangensis and O. longilobus occupy different niches due to difference in growth habits [30]. Through nuclear and organelle markers, our previous genetic and phylogeographic survey indicated that significant genetic variations were found between the two species [30-34). Utilizing transcriptome data, it was revealed that there were 3410 differentially expressed genes among two species, which primarily involved lipid, carbohydrate, and amino acid metabolism, as well as xenobiotic biodegradation and metabolism [28, 34]. Thus, according to the current geographical distribution of O. taihangensis and O. longilobus, we presumed that adaptive divergence processes might have played an additional role in maintaining the differentiation of the two Opisthopappus species that inhabit different regions of the Taihang Mountains [29].

Whether driven by morphological, genetic, or metabolic aspects, differentiation proceeded between O. taihangensis and O. longilobus. Nevertheless, the molecular mechanisms underlying these phenomena remains unclear, where a paucity of genetic resources has made further studies on the neutral and adaptive divergence processes of Opisthopappus a challenging task.

For this study, the transcriptome sequencing through a series of detailed comparative analyses were performed to highlight the characterization of baseline transcriptomic differences between O. taihangensis and O. longilobus. Further, the analyses of divergent gene expression and major metabolic pathway patterns for these two species were carried to provide useful data regarding their adaptation to specific habitats and enable the exploration of adaptive evolutionary signatures between the two species (i.e., types of genes with putatively accelerated nonsynonymous divergence, and those that showed signs of positive selection and their potential correlations with adaptability). It is anticipated that the data generated in this study will enhance our knowledge of the regulatory networks and pathways involved in assessing the roles of evolutionary factors (e.g., natural selection) on the adaptation, and speciation of Opisthopappus.

Plant materials

O. taihangensis distributes mostly across Henan and Shanxi Provinces, which are located along the southern Taihang Mountains and belong to a sub-humid warm temperate climate; while O. longilobus generally grows in Hebei and Shanxi Provinces, which reside along the northern Taihang Mountains with a temperate continental monsoon climate.

To elucidate the molecular mechanisms of differentiation proceeded between O. taihangensis and O. longilobus, the perennial plants of the two species were sampled from Shanxi Province where the two have adjacent distribution as closely as possible (Study protocol comply with relevant institutional, national, and international guidelines and legislation, and collection of O. taihangensis and O. longilobus have obtained the permission from the College of Life Sciences of Shanxi Normal University). During the same growth period of O. taihangensis and O. longilobus at July in 2015, fresh mature leaves from the same branch and position of each species were sampled. Each individual from the same population was collected from different locations at least 10 m apart. Then mixed leaves of five individuals were collected as a sample. Three samples (biological replicates) were tested for each population, including taihangensis1, taihangensis2 and taihangensis3 for O. taihangensis and longilobus1, longilobus2 and longilobus3 for O. longilobus. The samples were frozen in liquid nitrogen (immediately stored at -80°C) and then sent to Shanghai Personal Biotechnology Co., Ltd. The data were generated in accordance with a transcriptome analysis program displayed in our previous studies [28].

Numerous studies have shown that transcription factors (TFs) are involved in the physiological processes and adaptation of plants, which can regulate gene expression during all developmental and growth stages. Positive selection genes (PSGs) are those that operate under natural selection, which reflect the great adaptation of species to factors that are intimately linked to their survival under specific environmental conditions. Thus, TFs and PSGs were analyzed to explore the variations in the development and adaptation of the two Opisthopappus species for this study.

TF identification and expression analysis

PlantTFDB (http://planttfdb.cbi.pku.edu.cn/index.php) is a plant transcription factor database that includes the sequences of 58 plant transcription factor families from 165 plant species [35 Jin et al. 2017]. Clean reads were assembled de novo using Trinity software, and the sequences of unigenes were Blastx aligned to the transcription factor database PlantTFDB, where the best of these with E values of less than 1e-5, were screened as the annotation data of the unigenes. Candidates that contained DNA binding domains were recognized by GO annotation for final TF identification.

An in-depth transcriptome analysis would definitely facilitate an improved perception of the species under study [15, 21, 36]. An advanced bioinformatics tool (MapMan] could be used to analyze and compare transcriptional responses between phylogenetically species, and especially useful at the level of functional categories in cross-species comparisons [37]. This analysis allows us to explore gene categories from large data sets to extract meaningful information. In this study, MapMan [38, 39] was employed to comprehensively interpret transcriptome data and visualize the functionalities of associated genes. Through the use of MapMan, it was evident that significant variations between two species were conspicuous with respect to differentially expressed genes [40]. At the same time, we employed metabolic overview installed in the MapMan tool to identify the primary and secondary metabolic pathways related to the genes associated with transcription factors (TFs) [41]. MapMan metabolism encompasses photosynthesis pathways, carbohydrate metabolism, N-dependent pathways (e.g., amino acid metabolism), as well as cell wall, lipid, and secondary metabolism. Meanwhile, the diﬀerentially expressed TFs and transporters were isolated by Fisher’s Exact Test (P-value < 0.05) and enrichment fold ≥1.5 compared with the whole genome background.

PSG identification and expression analysis

Orthologous contigs and estimation of substitution rates

The assembled sequences were initially employed to predict the ORF (Open Read Frame) regions. Subsequently, the ORFs were blasted using OrthoMCL (http:// orthomcl.org/orthomcl/) software. Orthologous pairs with identities of < 60% were excluded, and only 1:1 orthologous pair in both lineages were retained. Next, based on the BLAST results, Markov clustering was used with the Markov Cluster algorithm (MCL) method. Aligned sequences showing a 90% identity were defined as pairs of putative orthologues. The best-hit sequences of each cluster were then used in subsequent analyses.

The KAKS_CALCULATOR software [42, 43] was used to calculate nonsynonymous (dN) and synonymous (dS) substitution rates of the single-copy orthologous genes, and the dN/dS ratios of each putative orthologous pair using the YN [44] algorithm in PAML software. The dN/dS ratio is an estimate of natural selection acting on the genes. The estimated values and its proportion were assigned to one out of five selection scenarios: (ⅰ) dN/dS = 0–0.6 (strong negative selection), (ⅱ) dN/dS = 0.61–0.9 (weak negative selection), (ⅲ) dN/dS = 0.91–1.5 (neutral selection), (ⅳ) dN/dS = 1.51–5 (weak positive selection), and (ⅴ) dN/dS > 5 (strong positive selection) [45]. Alignments with dS values of < 0.01 and dS or dN > 2 were discarded. Very high dN/dS (>10) were also removed as they indicated bias.

According to the standard errors computed from YN00 for a Student’s t-test, the statistical significance of the differences between dN and dS was analyzed [44 Yang & Nielsen, 2000]. The pair of othologous with a P-value < 0.05 and dN/dS > 1.5 was designated as preliminary candidate positive selection. A P-value < 0.01 and dN/dS > 5 were used as threshold for candidate genes with stronger positive selection.

GO classification analysis

The identified positive selection genes were processed to retrieve associated Gene Ontology (GO) terms describing biological processes, molecular functions, and cellular components using TBtools software. Expression levels for all of the positive selection genes in replicates were assessed by mapping high quality (HQ) filtered reads using BOWTIE2 [46]. Mapped reads were further normalized using the Fragments Per Kilobase Per Millions (FPKM) method. The significance of the difference in GO term abundance between the two datasets was tested using the Fisher’s exact test of the GOSSIP package [47] implemented in BLAST2GO [48].

Pathway enrichment analysis

To analyze the PSG pathways between the two species, we employed two annotation methods (KEGG and MapMan) to compare the expressional differentiation between O. taihangensis and O. longilobus. KEGG (http://www.kegg.jp/) [49] pathway annotation was performed using the BLASTX algorithm with E-values < 1.00E^-5 and enabled ‘sensitive mode’ against the KEGG (Kyoto Encyclopedia of Genes and Genomes) gene peptide database. Subsequently, the corresponding KO identifiers, EC numbers and pathway categories were parsed using customized Perl scripts. For the KEGG function classification, unigenes were annotated into five categories of KEGG metabolic pathways (cellular processes, environmental information processing, organismal systems, and metabolic and genetic information processing).

For MapMan analysis [41], transcripts expressed in both species were annotated with the TAIR database (Arabidopsis homologs). A MapMan BIN file, with a hierarchical ontology system for basil genes, was prepared using Mercator [50, 51], by comparing transcripts against already-classified proteins. BIN level information and gene identifiers were also derived from the same database.

Transcription factors (TFs)

Of the 483 transcript unigenes, 406 were transcription factors that belonged to 41 transcription factors families, including bHLH, bZIP, etc. Among them, the bHLH family was the most abundant with 67 transcription factors, whereas the second most abundant family was FAR1 with 50 transcription factors, followed by the NAC (40 transcription factors), M-type (24 transcription factors), MYB and MYB related (24 transcription factors), and WRKY (22 transcription factors) families (Table 1).

While, families having 10-19 transcription factors included Nin-like (19 factors), B3 (17 factors), EIL (13 factors), C3H (11 factors), G2-like (10 factors), and Trihelix (10 factors). Other families had less than 10 transcription factors, such as BES1 (8 factors), SBP (7 factors), and TCP (4 factors). Seven families (ZF-HD, RAV, NF-YA, HRT-like, GATA, E2F/DP, and CAMTA) had only one factor. Among the 406 TFs, 309 were shared by O. taihangensis and O. longilobus, while 46 were exclusive to O. longilobus and 13 were exclusive to O. taihangensis.

To systematically assign functions to the 406 TF genes, we employed the diverse overview tools installed in MapMan. The metabolism overview divided the 406 TF genes as: seven for amino acid metabolism, four for secondary metabolism, three for lipid metabolism, one for cell wall metabolism, one for major carbohydrate metabolism, one for photosystems, and one for fermentation (Figure1).

Significant expression differentiation between O. taihangensis and O. longilobus occurred in the terpenoid metabolic pathway, in which the transcription factors of O. longilobus were obviously more active than those of O. taihangensis, which had a relatively lower expression level. Furthermore, for the amino acid metabolism pathways, the TFs were observed to be more numerously expressed in O. taihangensis, whereas in other pathways the TFs were almost equally expressed between the two species (Figure 1).

Positive selection genes (PSGs)

Based on the predicted ORF, we identified 59753 pairs of putative orthologous contigs between O. taihangensis and O. longilobus using ORTHOMCL. The one-to-one, reciprocal best method for elucidating orthologous proteins generated 50587 putative orthologous pairs. After filtering the gene pairs annotated with different protein products in the Swiss-Prot database, 38986 pairs of putative orthologs were finally identified and used in downstream analyses.

Of these, 11598 pairs had only synonymous or non-synonymous substitutions, and 38986 pairs had both types of substitutions, for which the dN/dS ratios were calculated. Among the 38986 pairs of putative orthologs, 1596 putative othologous pairs had dN/dS ratios = 0.91–1.5, which indicated ongoing neutral evolution. And 36187 had dN/dS ratios between 0.0001 and 0.9 under negative selection that might have experienced relaxed purifying selection and/or unfixed mutations without altering the encoded amino acid sequence during the speciation period.

There were five othologous pairs with dN/dS ratios >10 and 1203 with dN/dS ratios > 1.51 (P < 0.05). According the frequency distribution plots of all dN/dS ratios of two species, we took a more appropriate threshold of 1.5 for the dN/dS ratio as an indicator of positive selection [22, 16, 45], thus 1203 pairs with dN/dS values >1.51 were identified as undergoing positive selection (P < 0.05) (Figure 2A). Among the positive selection genes, 63 were shared by the two species, and 1125 were exclusive to O. taihangensis and O. longilobus, respectively (Figure 2A). All positively selection genes were shown in a boxplot distribution with a median of 2.0308 (Figure 2B). These results suggested divergent evolutionary processes in O. taihangensis and O. longilobus. Meanwhile, we clustered the above 1203 putative orthologous pairs into three main GO categories (biological process, cellular component, and molecular function). Within the biological process category, the term ‘menaquinone metabolic process’ and ‘menaquinone biosynthetic process’ were the most dominant. Within the cellular component category, ‘nuclear origin of replication recognition complex and ‘origin recognition complex’ represented the major subcategories. Within the molecular function category, the main functional subcategories were ‘2-succinyl-5-enolpyruvyl-6-hydroxy-3-cyclohexene-1’ and ‘metal ion transmembrane transporter activity’ (Figure 3A).

Furthermore, we used KEGG pathways to investigate the biological functions of positive selection genes. Overall, the PSGs were assigned to 22 KEGG pathways. The most represented pathways were “Nicotinate and nicotinamide metabolism”, “Zeatin biosynthesis (Metabolism of terpenoids and polyketides)”, “Other glycan degradation”, “Histidine metabolism”, and “Sphingolipid metabolism” (Figure 3B). For a comprehensive assessment of the variations/similarities in the transcriptomes of both species, PSG transcripts were plotted using the MapMan tool, which were separated by bins based on their functional ontology (Figure 4). The functional pathway classification of 1203 PSGs mapped to the “overview of metabolism”, with 4 to lipid metabolism, 12 to secondary metabolism, 7 to amino acid metabolism, 5 to mitochondrial electron transport, 4 to nucleotide metabolism, 3 to major carbohydrate metabolism, one to the tricarboxylic acid (TCA) cycle, 2 redox ascorbate and glutathione ascorbate, 1 to oxidative PP.6-phosphogluconolactonase, one to C1-metabolism, one to cell wall metabolism, one to photosystems and one to tetrapyrrole synthesis.

Although major differences were not observed between O. longilobus and O. taihangensis, and most transcripts were shared by the two species, there were still a number of PSGs that differed. For example, the lipids pathway, major carbohydrate metabolism pathways, nucleotide metabolism pathways, tetrapyrrole metabolism pathways, mitochondrial electron transport pathways, presented higher expression levels in O. longilobus than in O. taihangensis. On the contrary, valine, leucine, and isoleucine biosynthesis in the amino acid pathway and C-1 metabolism pathways were obviously more active in O. taihangensis than O. longilobus (Figure 4A).

We simultaneously annotated the mitochondrial electron transport / ATP synthesis pathway, nucleotide metabolism pathway and sucrose-starch pathway. An overview of mitochondrial electron transport / ATP synthesis pathway showed that the various pathways presented mainly on NADH-DH complex I and cytochrome c oxidase of the mitochondrial electron transport chains complex I and IV (Figure 4B). An overview of major carbohydrate metabolism pathways suggested that the PSG expression levels involved in beta-amylase for O. longilobus was significantly higher than those for O. taihangensis (Figure 4C).

Secondary metabolites related to TFs and PSGs

Secondary metabolites play major roles in the adaptation of plants to the environment and in overcoming stress conditions, while also contributing to the specific odors, tastes, and colors of plants. The genes related to secondary metabolism are pertinent for the responses of plants to various environmental stimuli. O. taihangensis and O. longilobus secrete chemicals and exude aromas to protect against biological and abiotic stresses [52]. These chemicals are primarily secondary metabolites, including polyphenols, alkaloids, and terpenoids.

Based on the above results, either TFs or PSGs, or both, can enrich the metabolic process pathway; thus, we further annotated the secondary metabolite pathway of the metabolic process. The analysis of pathways enriched by MapMan revealed that terpenoid biosynthesis was simultaneously enriched at PSGs (Figure 5A) and TFs (Figure 5B), and only the terpene biosynthesis pathway exhibited higher expression in O. longilobus than O. taihangensis. Except for the terpene biosynthesis pathway, the PSGs were primarily identified on the flavanones biosynthesis and xanthophyll biosynthesis pathways (Figure 5A). Furthermore, there were significantly different expression levels between the two species in these pathways.

Among the PSGs, terpenoids biosynthesis were enriched on mono-/sesquiterpene-/diterpene synthase where the expression levels were higher in O. longilobus than O. taihangensis (Figure 6A). Except for the monoterpene and diterpene biosynthesis pathways, the genes were mainly enriched on non-mevalonate pathway secondary metabolism pathways, especially ISPD (4-diphosphocytidyl-2-Cmethyl-D-erythritol synthase) and ISPE (4-diphosphocytidyl-2-C-methyl-D-erythritol kinase), where the expression levels were still higher in O. longilobus than O. taihangensis (Figure 6A). These involved in trinity_dn129788_c0_g1 gene and trinity_dn155665_c0_g2 gene.

The terpenoid biosynthesis pathway was further analyzed by MapMan tools. The results indicated that the enriched pathway was carotenoid biosynthesis pathway with higher expression levels in O. longilobus than O. taihangensis for TFs. Moreover, this differentiation primarily occurred along the monoterpene and diterpene biosynthesis pathways. For these two pathways, the involved TFs all belonged to the bHLH family (Figure 6B).

Although the differentiation and variations between the Opisthopappus species were performed using diverse molecular markers [28, 30, 32-34], the availability of data on the molecular mechanisms of evolution and adaptation remained generally limited. To reveal the molecular mechanisms of the two Opisthopappus species associated with their surroundings, comparative transcriptome analysis was performed, which provided some clues for an improved elucidation of major differences under natural selection.

Adaptive characteristics of TFs

Transcription factors (TFs) are significant key elements of various physiological and biochemical pathways in higher plants. The activation and repression of TFs can primarily regulate the expression of genes related to growth, stress responses of plants, and other developmental processes [53, 40]. The activities of TFs often depend on the developmental stage, exogenous stimuli, and/or the presence of co-regulatory proteins. Various TFs are employed for the regulation of the expression of different genes. Thus, TF-based gene expression regulation allows plants to respond to changes in their environments [54, 40].

Based on transcriptome data, we initially obtained 406 TFs (Table 1). The bHLH, FAR1, NAC, MYB, and WRKY family members were highest among the obtained factors. These factors were found to be significantly expressed for lipid, secondary, amino acid, and nucleotide metabolism, which constituted the main biological processes. The responses of plants to various environments are pertinent for successful growth and reproduction [40]. The proactive responses to environmental/developmental cues were explicitly depicted in O. taihangensis, as many TFs were involved in the amino acid metabolism pathways (Figure 1).

In particular, for our study, O. longilobus had a higher expression level in the terpenoid metabolic pathway than did O. taihangensis (Figure 5). Terpenoids are one of the most common compounds in the secondary metabolites of plants, which have numerous physiological and ecological functions, such as attracting pollinated insects, regulating growth and development, resisting environmental stress, and participating in the defense of pests [55]. Meanwhile terpenoids are important components of the fragrances of flowers. For the two species of Opisthopappus under study, these secondary metabolites are primarily emitted by their leaves and flowers.

According to the Flora Reipublicae Popularis Sinicae (FRPS), the surfaces of O. longilobus leaves are smooth and glabrous, while O. taihangensis is sparsely pubescent on both leaf surfaces, which include glandular trichomes and non-glandular trichomes [52]. The terpenoids produced in O. taihangensis can be released through glandular trichomes give rise to the unique aroma for communication and defense [52]. For O. longilobus, owing to the lack of glandular trichomes on the surfaces of its leaves, terpenoid biosynthesis genes are more highly expressed to compensate for this deficit through the use of alternate emission strategies.

Basic helix-loop-helix (bHLH) transcription factors, which are one of the largest families and a large superfamily in plants, play relevant roles in a variety of developmental and evolutionarily conserved processes, including cell-fate specification, tissue differentiation, and secondary metabolites [56, 57, 24]. For example, the bHLH factor VvMYC1 can regulate anthocyanin and/or proanthocyanidin (PA) synthesis of the flavonoid pathway [58]. In the terpenoid biosynthesis pathway of Opisthopappus, the involved TFs are all bHLH family members (Figure 5). Several bHLH TFs have been found in Medicago truncatula [59], Chenopodium quinoa [60], Panax notoginseng [61], and Panax ginseng [62]. Meanwhile, many bHLH genes respond to various forms of stress such as drought, salt, and cold stresses [59, 53, 63]. Other TFs (e.g., NAC, MYB, and WRKY) have been investigated for their capacities to improve tolerance and resistance in many plants [64-68], while participating in abiotic and biotic stress responses [69 , 70,17]. Both bHLH and MYB can modulate secondary metabolism biosynthesis in plants that are subject to environmental stress. This suggested that TFs (particularly bHLH family members) might play critical roles in the responses of the two Opisthopappus species to heterogeneous environments.

Adaptability of PSGs

To describe genome-wide levels of coding sequence evolution and estimate the effects of natural selection on lineage divergence, we calculated the dN/dS ratios for O. longilobus and O. taihangensis orthologs. This ratio has frequently been employed as an indicator of frequency and selection modes under which protein-coding genes evolved [16]. In this study, the average dN/dS ratios across the pairs of the Opisthopappus genus were much lower than 0.9, which suggested that purifying selection had a general influence on the evolution of most protein-coding regions (ORF) of the two species, as has been observed in other plants [71, 16]. Actually, the genes that were most under the influence of purifying selection primarily contained structural or “housekeeping” genes, for example tyrosyl-tRNA synthetase. As these genes are involved in processes that are crucial for organisms, purifying selection serves to eliminate deleterious, nonsynonymous mutations [72, 73, 16]. Among the 38986 orthologs shared between lineages, 1203 gene pairs significantly exhibited dN/dS ratios >1.5 (P < 0.05) (Figure 2A, 2B). These genes were involved in several biological functions (e.g., metabolic processes) and may constitute candidates that are under the effect of positive selection; thus, potentially associated with species divergence.

Adaptive divergence at molecular level may be reflected by an increased rate of non-synonymous changes within genes involved in adaptation [74]. Analyses of the PSGs identified from O. taihangensis and O. longilobus yielded significantly different landscapes of biological processes, cellular component and molecular function (Figure 3). O. taihangensis and O. longilobus would experience divergent adaptation during their evolution, and that of genes related to adaptation was under rapid evolution and/or have signs of positive selection. These findings implicated that Opisthopappus were experienced the ongoing accelerated evolution under different environment. Thus, more than 1000 genes under positive selection between the two Opisthopappus species might have played a role in shaping the divergence of this genus, given that O. longilobus and O. taihangensis are exposed to sub-humid warm temperate, and temperate continental monsoon climates, respectively.

Subsequently, KEGG enrichment and MapMan analyses suggested that significant GO terms, mainly related to metabolic regulation (Figure 3B), were important among the Opisthopappus species (Figure 4). Secondary metabolism was also significantly represented for O. longilobus and O. taihangensis in our PSG analysis (Figure 5).

Beta-amylase, a member of family 14 of glycosyl hydrolases [75], can hydrolyzes the α-1,4-glucosidic linkages in starch, removing successive maltose units from the non-reducing ends of the chains [76]. Starch is the main form of carbon storage in higher plants and it accumulates in different organs, which can be consumed by the cell in which they are produced, transported to nonphotosynthetic sink tissues, or stored for later use. According to studies using transgenic and mutant plants, beta-amylase seems to be important for normal degradation of the transient starch accumulated in chloroplasts [76]. Meanwhile, beta-amylases possibly facilitating the rapid starch degradation under heat and drought [77]. In this study, we found that the PSG expression levels involved in beta-amylase for O. longilobus was significantly higher than those for O. taihangensis (Figure 4A). Starch degradation as a major response in conditions prohibiting photosynthesis or prolonged drought [77], therefore, O. longilobus seems better adapted to weak light or drought conditions than O. taihangensis.

Beta-amylase can convert starch into maltose and glucose as the predominant form in which sugars are translocated in plants presented in Figure 4C. Consistent with this notion, O. longilobus exhibited higher levels of cellular glucose and fructose than did O. taihangensis. Collectively, these data signified that the higher expression levels of beta-amylase genes in O. longilobus might accelerate its cell growth processes (in contrast to O. taihangensis), which have important roles in the evolution of the two species. Ultimately, glycan degradation resulted in the accumulation of sucrose, which is critical for the protection of plants against xenobiotic and oxidative stresses [13], and which might facilitate the enhanced survival of O. longilobus for relatively more drought habitats over than O. taihangensis.

The electron transport chain (ETC, respiratory chain) is a series of protein complexes that transfer electrons from electron donors to electron acceptors via redox reactions (both reduction and oxidation occurring simultaneously) and couples this electron transfer with the transfer of protons (H+ ions) across a membrane. The electron transport chain is built up of peptides, enzymes, and other molecules. In this study, the expression levels of relative PSGs involved in complex I (NADH ubiquinone oxireductase) and complex IV (cytochrome c oxidase, EC 1.9.3.1) for O. longilobus was higher than those for O. taihangensis (Figure 4B). The ETC sustains the major mitochondrial function of ATP generation, in relation to the metabolic dynamics of tricarboxylic acids (TCAs), acetyl-CoA, ADP, oxidized (NAD+) or reduced (NADH) β-nicotinamide adenine dinucleotide, oxidized (FAD) or reduced (FADH2) flavin adenine dinucleotide [78]. Thus, we concluded that the changes in ETC can affect nicotinate and nicotinamide metabolism, glycan degradation metabolism and transmembrane transport activity, which was supported by the KEGG and MapMan results (Figures 3 and 4). The various expression levels of relative PSGs between these two species were also regarded as different adaption signals to habitats.

Additionally, the responses to sulfur between O. longilobus and O. taihangensis were significantly different, which suggested that the two species possessed different stress tolerances for S in response to heterogeneous environments.

Adaptive function of Terpenoids

Secondary metabolism produces large populations of specialized molecules that are required for plants to survive in their environments and are essential for communicating with other organisms in a mutualistic (e.g., to attract beneficial organisms such as pollinators), or antagonistic (e.g., to combat herbivores and pathogens) manner. Under baseline or non-stress conditions, it is expected that mutualistic metabolites, or those required for normal physiological processes, are expressed [40]. This scenario, particularly in terms of terpenoid biosynthesis, was unambiguously observed in the present study as a difference between the two species TFs and PSGs (Figure 5).

The terpenoid pathway is intricately regulated by endogenous and environmental factors that enable spatially and temporally controlled metabolite production [79-81]. Terpenes are the major components of plant terpenoids, which comprise the largest and the most diverse class of natural products as a homologous series of molecules as isoprene polymers. They often play important roles in plant defenses by attracting the enemies or predators of herbivores and repelling herbivorous insect feeders [82].

Isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP) are the main intermediate compounds in the terpenoid pathway for terpene biosynthesis [83, 20]. Their precursor molecules (C5) are generated via the process of isoprenogenesis [84, 85]. In plants, isoprenogenesis occurs through two discrete biosynthetic pathways: the mevalonic acid (MVA) pathway in the cytosol, and the 2-C-methyl-D-erythritol 4-phosphate/1-deoxyD-xylulose 5-phosphate (MEP/DOXP) pathway within plastids. In the MEP pathways, the C5 units are catalyzed step by step via a series of enzymes to synthesize IPP. The 4-diphosphocytidyl-2C methyl-D-erythritol synthase (ISPD/CMS) and 4-diphosphocytidyl-2C-methyl-D-erythritol kinase (ISPE/CMK) catalyze the synthesis of 4-diphosphocytidyl-2C-methyl-D-erythritol (CDP-ME) and 4-diphosphocytidyl-2C-methyl-D-erythritol 2-phosphate (CDP-ME2P), respectively. CDP-ME and CDP-ME2P are the intermediate compounds of IPP biosynthesis.

IPP or DMAPP may be linked through head to tail condensation reactions to generate terpenes of different classes (e.g., mono, di, and triterpenes), which are catalyzed by geranyl diphosphate synthase (GPPS) to synthesize geranyl pyrophosphate (GPP). The GPP results in farnesyl pyrophosphate (FPP) and geranylgeranyl diphosphate (GGPP) via farnesyl pyrophosphate synthase (FPPS) and geranylgeranyl diphosphate synthase (GGPPS) enzymes, respectively. GPP and GGPP are substrates for monoterpene and diterpene biosynthesis in the terpene pathway [86 -88].

Monoterpenes and diterpenes can impart the distinct flavors and aromas of plants [89, 20]. These terpenes are produced in plants as secondary bioactive metabolites, often for ecological adjustment and protection from microbial pathogens, fungi, pests, and predation [90]. An overview of terpenoid metabolism of the two species revealed that the putative PSGs (trinity_dn129788_c0_g1 and trinity_dn155665_c0_g2) displayed different levels of expression for IPP biosynthesis, particularly in the synthesis of CDP-ME and CDP-ME2P (Figure 6).

Meanwhile, both the PSGs and TFs involved in the downstream steps of the MEP pathway exhibited a relatively higher expression in contrast to the two species under study, which displayed different expression levels for monoterpene and diterpene biosynthesis (Figure 6). These results indicated variable adaptability in the responses of the two species to different environments. During the evolutionary process, the genes involved in terpenoid biosynthesis pathway (in particular, monoterpene and diterpene biosynthesis), begun to diverge from their expression levels under different long-term environment stresses. Simultaneously, several bHLH TFs involved this metabolic process to further regulate the expression of these positive selection genes. Ultimately, the adaptive phenotypic characteristics occurred in O. longilobus and O. taihangensis, such as the presence of trichomes or not, and different fragrances.

This work presents a comparative genome-wide transcriptome analysis between two Opisthopappus species, toward the identification of biological processes and functional genes that facilitate differentiation and adaption of these two species during evolutionary processes. Under natural selection, O. taihangensis gradually modified its gene expression model to adapt to its surroundings, which differed from O. longilobus. The screened several PSGs and TFs involved in the terpenoid biosynthesis of monoterpenoids and diterpenoids and proposed adaptive signal responses to heterogeneous environments. The enhancement or reduction of the expression level of several genes between O. longilobus and O. taihangensis, it provided a comprehensive characterization of the expression profiles during the developmental stages of O. longilobus and O. taihangensis. These results will not only shed light on how differentiations between the two Opisthopappus species occur, but also open the door to an increased understanding of how plants thrive in mountains environments. Additionally, our research lays a theoretical foundation for the further investigation of the molecular mechanisms of specific characteristics that have important academic and application value for other Asteraceae species.

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Funding

This study was supported by the National Natural Science Foundation of China (31970358), the Research Project Supported by Shanxi Scholarship Council of China (2020-090) and the Graduate Education Innovation Project of Shanxi (2020SY323).

Conflict of interest

No conflict of interest exits in the submission of this manuscript, and the manuscript is approved by all authors for publication.

Author contributions

Chen N, Genlou Sun and Yiling Wang designed experiments; Chen N and Yiling Wang carried out experiments, analyzed experimental results and data, developed analysis tools and wrote the main manuscript text. Hao Zhang, En Zang, Zhixia Liu, Yafei Lan, Weili Hao, Shan He, Xing Fan, Genlou Sun, Yiling Wang assisted with writing the manuscript. All authors reviewed the manuscript.

Acknowledgements

Data availability statement

Data deposited via the Genebank: accession SUB4002238: PRJNA471103.

Ferchaud, A.L., Hansen, M.M. (2016). The impact of selection, gene flow and demographic history on heterogeneous genomic divergence: threespine sticklebacks in divergent environments. Mol Ecol, 25(1), 238-259.
Svensson, E.I., Goedert, D., Gómez-Llano, M.A., Spagopoulou, F., Nava-Bolaños, A., Booksmythe, I. (2018). Sex differences in local adaptation: what can we learn from reciprocal transplant experiments? Philos T R Soc B, 373, 20170420.
Li, Y.L., Xue, D.X., Zhang, B.D., Liu, J.X. (2019). Environmental selection in the catadromous roughskin sculpin Trachidermus fasciatus. Genome Biol Evol, 11(7), 1751-1764.
Enbodya, E.D., Petterssona, M.E., Sprehna, C.G., Palmb, S., Wickströmb, H., Andersson. L. (2021). Ecological adaptation in European eels is based on phenotypic plasticity. P Natl Acad Sci USA, 118: 4e2022620118.
Zhang, X., Rayner, J.G., Blaxter, M., Bailey, N.W. (2021). Rapid parallel adaptation despite gene flow in silent crickets. Nat Commun, 12, 50.
Delph, L.F. (2017). The study of local adaptation: a thriving field of research. J Hered, 109,1-2.
Ranade, S.S., García-Gil, M.R. (2021). Molecular signatures of local adaptation to light in Norway spruce. Planta, 253(2), 53.
Honjo, M.N., Kudoh, H. (2019). Arabidopsis halleri: a perennial model system for studying population differentiation and local adaptation. AoB Plants, 11, plz076.
Olson-Manning, C.F., Wagner, M.R., Mitchell-Olds, T. (2012). Adaptive evolution: evaluating empirical support for theoretical predictions[J]. Nature Reviews Genetics, 13(12):867.
Kawecki, T.J., Ebert, D. (2004). Conceptual issues in local adaptation. Ecol Lett, 7, 1225-1241.
Hereford, J. (2009). A quantitative survey of local adaptation and fitness trade-offs. Am Nat, 173, 579-588.
Whitlock, M.C. (2015). Modern approaches to local adaptation. Am Nat, 186, S1-S4.
Yan, L., Liu, Z., Xu, H., Zhang, X., Zhao, A., Liang, F., … Ni, Z. (2018). Transcriptome analysis reveals potential mechanisms for different grain size between natural and resynthesized allohexaploid wheats with near-identical AABB genomes. BMC Plant Biol, 18, 28.
Kim, J., Park, S.J., Lee, I.H., Chu, H., Penfold, C.A., Kim, J.H., … Lim, P.O. (2018). Comparative transcriptome analysis in Arabidopsis ein2/ore3 and ahk3/ore12 mutants during dark-induced leaf senescence. J Exp Bot, 69, 3023-3036.
Li, J., Ye, C., Chang, C. (2020). Comparative transcriptomics analysis revealing flower trichome development during flower development in two Lonicera japonica Thunb. cultivars using RNA-seq. BMC Plant Biol 20, 341.
Wang, Y., Jiang, W., Ye, W., Fu, C., Gitzendanner, M.A. Soltis, P.S., … Qiu, Y. (2018). Evolutionary insights from comparative transcriptome and transcriptome-wide coalescence analyses in Tetrastigma hemsleyanum. BMC Plant Biol, 18, 208.
Ma, Z., Lin, S. (2019). Transcriptomic revelation of phenolic compounds involved in aluminum toxicity responses in roots of Cunninghamia lanceolata (Lamb.) Hook. Genes, 10, 835.
Zhang, H.B., Xia, E.H., Huang, H., Jiang, J.J., Liu, B.Y., Gao, L.Z. (2015). De novo transcriptome assembly of the wild relative of tea tree (Camellia taliensis) and comparative analysis with tea transcriptome identified putative genes associated with tea quality and stress response. BMC Genomics, 16, 298.
Onda, Y., Mochida, K., Yoshida, T., Sakurai, T., Seymour, R.S., Umekawa, Y., … Ito, K. (2015). Transcriptome analysis of thermogenic Arum concinnatum reveals the molecular components of floral scent production. Scie Rep, 5, 8753.
Torre, S., Tattini, M., Brunetti, C., Guidi, L., Gori, A., Marzano, C,. … Landi, M, (2016). De novo assembly and comparative transcriptome analyses of red and green morphs of sweet basil grown in full sunlight. PLoS One, 11(8), e0160370.
Huang, G., Zeng, Y., Wei, L., Yao, Y., Dai, J., Liu, G., … Gui, Z. (2020). Comparative transcriptome analysis of mulberry reveals anthocyanin biosynthesis mechanisms in black (Morus atropurpurea Roxb.) and white (Morus alba L.) fruit genotypes. BMC Plant Biol, 20, 279.
Zhang, L., Yan, H.F., Wu, W., Yu, H., Ge, X.J. (2013). Comparative transcriptome analysis and marker development of two closely related Primrose species (Primula poissonii and Primula wilsonii). BMC Genomics, 14, 39.
Bonasio, R., Zhang, G., Ye, C., Mutti, N.S., Fang, X., Qin, N., … Liebig, J. (2010). Genomic comparison of the ants Camponotus floridanus and Harpegnathos saltator. Science, 329, 1068-1071.
Bao, L., Gao, H., Zheng, Z., Zhao, X., Zhang, M., Jiao, F., . . . Qian, Y. (2020). Integrated transcriptomic and un-targeted metabolomics analysis reveals mulberry fruit (Morus atropurpurea) in response to sclerotiniose pathogen Ciboria shiraiana infection. Int J Mol Sci, 21, 1789.
Shih, C. (1979). Opisthopappus shih-a new genus of compositae from China. Acta Phytotaxon Sin, 32, 110-112
Wang, F.Z., Tang, J., Chen, X.Q., Liang, S.J., Dai, L.K. (1978). Liliaceae. Flora of China Editorial Committee of China Academy of Science. Beijing: Science Press 73-74.
Ding, B.Z., Wang, S.Y. (1997). Flora of Henan. Zhengzhou: Henan Science and Technology Press.
Chai, M., Wang, S., He, J., Chen, W., Fan, Z., Li, J., … Wang, Y. (2018). De novo assembly and transcriptome characterization of Opisthopappus (Asteraceae) for population differentiation and adaption. Front Genet, 9, 371.
Ye, H., Wang, Z., Hou, H.M., Wu, J.H., Gao, Y., Han, W., … Wang, Y.L. (2021). Localized environmental heterogeneity drives the population differentiation of two endangered and endemic Opisthopappus Shih species. BMC Ecol Evo, 21, 56.
Wang, Y.L., Yan, G.Q. (2014). Molecular phylogeography and population genetic structure of O. longilobus and O. taihangensis (Opisthopappus) on the Taihang Mountains. PLoS One, 9, e104773.
Wang, Y.L. (2013). Chloroplast microsatellite diversity of Opisthopappus shih (Asteraceae) endemic to China. Plant Syst Evol, 299, 1849-1858.
Wang, Y.L., Yan, G.Q. (2013). Genetic diversity and population structure of Opisthopappus longilobus and Opisthopappus taihangensis (Asteraceae) in China determined using sequence related amplified polymorphism markers. Biochem Syst Ecol, 49, 115-124.
Wang, Y.L., Zhang, C.Q., Lin, L.L., Yuan, L.H. (2015). ITS sequence analysis of Opisthopappus taihangensis and Opisthopappus longilobus. Acta Hortic Sin, 42, 178-186.
Chai, M., Ye, H., Wang, Z., Zhou, Y., Wu, J., Gao, Y., … Wang, Y. (2020). Genetic divergence and relationship among Opisthopappus species identified by development of EST-SSR markers. Front Genet, 11, 177.
Jin, J., Tian, F., Yang, D.C., Meng, Y.Q., Kong, L., Lou, J., … Gao, G. (2017). PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res, 45 (D1), D1040-1045,
Bathke, J., Konzer, A., Remes, B., McIntosh, M., Klug, G. (2019). Comparative analyses of the variation of the transcriptome and proteome of Rhodobacter sphaeroides throughout growth. BMC Genomics, 20, 358.
Usadel, B., Poree, F., Nagel, A., et al. (2010). A guide to using MapMan to visualize and compare Omics data in plants: a case study in the crop species, Maize[J]. Plant Cell & Environment, 32(9):1211-1229.
Jung, K.H., An, G. (2012). Application of MapMan and rice net drives systematic analyses of the early heat stress transcriptome in rice seedlings. J Plant Biol, 55(6), 436-449.
Chandran, A.K., Lee, G.S., Yoo, Y.H., Yoon, U.H., Ahn, B.O., Yun, D.W., … Jung, K.H. (2016). Functional classification of rice flanking sequence tagged genes using MapMan terms and global understanding on metabolic and regulatory pathways affected by dxr mutant having defects in light response. Rice (NY), 9, 17.
Rathinam, M., Mishra, P., Vasudevan, M., Budhwar, R., Mahato, A., Prabha, A.L., … Sreevathsa, R. (2019). Comparative transcriptome analysis of pigeonpea, Cajanus cajan (L.) and one of its wild relatives Cajanus platycarpus (Benth.) Maesen. PLoS One, 14 (7), e0218731.
Thimm, O., Bläesing, O., Gibon, Y., Nagel, A., Meyer, S., Krüger, P., … Stitt, M. (2004). MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J, 37, 914-939.
Zhang, Z., Li, J., Zhao, X.Q., Wang, J., Wong, G.K., Yu, J.M. (2006). Ka/Ks_Calculator: calculating Ka and Ks through model selection and model averaging. Genom, Proteom Bioinf, 4, 259-263.
Wang, D., Zhang, Y., Zhang, Z., Zhu, J., Yu, J. (2010). KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genom. Proteom. Bioinformics, 8, 77-80.
Yang, Z., Nielsen, R. (2000). Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol Biol Evol, 17, 32-43.
Minaya, M., A Díaz-Pérez, Mason-Gamer, R., Pimentel, M., & P Catalán. (2015). Evolution of the beta-amylase gene in the temperate grasses: Non-purifying selection, recombination, semiparalogy, homeology and phylogenetic signal [J]. Molecular Phylogenetics & Evolution. 91:68-85.
Langmead, B., Salzberg, S. (2012). Fast gapped-read alignment with Bowtie 2. Nat Methods, 9, 357-359.
Blüthgen, N., Brand, K., Cajavec, B., Swat, M., Herzel, H., Beule, D. (2005). Biological profiling of gene groups utilizing Gene Ontology. Genome Inform, 16,106-115.
Conesa, A., Gotz, S., Garcia-Gomez, J.M., Terol, J., Talon, M., Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics, 21, 3674-3676.
Kanehisa, M., Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res, 28, 27-30.
Hanley, M.E., Lamont, B., Fairbanks, M.M., Rafferty, C.M. (2007). Plant structural traits and their role in anti-herbivore defence. Perspect Plant Ecol, 8,157-178.
Lohse, M., Nagel, A., Herter, T., May, P., Schroda, M., Zrenner, R., … Usadel, B. (2014). Mercator: a fast and simple web server for genome scale functional annotation of plant sequence data. Plant Cell Environ, 37, 1250–1258.
Guo, Y., Zhang, T., Zhong, J., B, T., Xu, T., Zhang Q., … Sun, M. (2020). Identification of the volatile compounds and observation of the glandular trichomes in Opisthopappus taihangensis and four species of Chrysanthemum. Plants, 9, 855.
Niu, X., Guan, Y., Chen, S., Li, H. (2017). Genome-wide analysis of basic helix-loophelix (bHLH) transcription factors in Brachypodium distachyon. BMC Genomics, 18, 619.
Gu, H., Yang, Y., Xing, M., Yue, C., Wei, F., Zhang, Y., … Huang J. (2019). Physiological and transcriptome analyses of Opisthopappus taihangensis in response to drought stress. Cell Biosci, 9, 56.
Dudareva, N., Klempien, A., Muhlemann, J.K., Kaplan, I. (2013). Biosynthesis, function and metabolic engineering of plant volatile organic compounds. New Phytol, 198, 16-32.
Gyoja, F. (2017). Basic helix-loop-helix transcription factors in evolution: Roles in development of mesoderm and neural tissues. Genesis, 55(9), e23051.
Zheng, Y., Wang, P., Chen, X., Sun, Y., Yue, C., Ye, N. (2019). Transcriptome and metabolite profiling reveal novel insights into volatile heterosis in the tea plant (Camellia sinensis). Molecules, 24, 3380.
Hichri, I., Heppel, S.C., Pillet, J., Léon, G., Czemmel, S., Delrot, V.L., … Bogs, J. (2010). The basic helix-loop-helix transcription factor MYC1 is involved in the regulation of the flavonoid biosynthesis pathway in grapevine. Mol Plant, 3, 509-523.
Mertens, J., Pollier, J., Bossche, R.V., Lopez-Vidriero, I., Franco-Zorrilla, J.M., Goossens, A. (2016). The bHLH transcription factors TSAR1 and TSAR2 regulate triterpene saponin biosynthesis in Medicago truncatula. Plant Physiol, 170, 194-210.
Jarvis, D.E., Ho, Y.S., Lightfoot, D.J., Schmöckel, S.M., Li, B., Borm, T.J.A., … Tester, M. (2017). The genome of Chenopodium quinoa. Nature, 542, 307-312.
Zhang, X., Ge, F., Deng, B., Shah, T., Huang, Z., Liu, D., … Chen, C. (2017). Molecular cloning and characterization of PnbHLH1 transcription factor in Panax notoginseng. Molecules, 22, 1268.
Chu, Y., Xiao, S., Su, H., Liao, B., Zhang, J., Xu, J., … Chen, S. (2018). Genome-wide characterization and analysis of bHLH transcription factors in Panax ginseng. Acta Pharm Sin B, 8(4), 666-677.
Miao, L., Gao, Y., Zhao, K., Kong, L., Yu, S., Li, R., … Yu, X. (2020). Comparative analysis of basic helix–loop– helix gene family among Brassica oleracea, Brassica rapa, and Brassica napus. BMC Genomics, 21, 178.
Chen, L., Song, Y., Li, S., Zhang, L., Zou, C., Yu, D. (2012). The role of WRKY transcription factors in plant abiotic stresses. BBA-Gene Regul Mech, 1819, 120-128.
Liu, C., Mao, B., Ou, S., Wang, W., Liu, L., Wu, Y., … Wang, X. (2014). OsbZIP71, a bZIP transcription factor, confers salinity and drought tolerance in rice. Plant Mol Biol, 84, 19-36.
Liu, J., Osbourn, A., Ma, P. (2015). MYB transcription factors as regulators of phenylpropanoid metabolism in plants. Mol Plant, 8, 689-708.
Zhang, X., Zhang, B., Li, M.J., Yin, X.M., Huang, L.F., Cui, Y.C., … Xia, X. (2016). OsMSR15 encoding a rice C2H2-type zinc finger protein confers enhanced drought tolerance in transgenic Arabidopsis. J Plant Biol, 59, 271-281.
Yoo, Y.H., Chandran, N.A.K., Park, J.C., Gho, Y.S., Lee, S.W., An, G., … Jung, K.H. (2017). OsPhyB-mediating novel regulatory pathway for drought tolerance in rice root identified by a global RNA-Seq transcriptome analysis of rice genes in response to water deficiencies. Front Plant Sci, 8, 850.
Kumari, M., Taylor, G.J., Deyholos, M.K. (2008). Transcriptomic responses to aluminum stress in roots of Arabidopsis thaliana. Mol Genet Genomics, 279, 339-357.
Tang, N., Zhang, H., Li, X., Xiao, J., Xiong, L. (2012). Constitutive activation of transcription factor OsbZIP46 improves drought tolerance in rice. Plant Physiol, 158, 1755-1768.
Buschiazzo, E., Ritland, C., Bohlmann, J., Ritland, K. (2012). Slow but not low: genomic comparisons reveal slower evolutionary rate and higher dN/dS in conifers compared to angiosperms. BMC Evol Biol, 12, 8.
Hughes, A.L., Packer, B., Welch, R., Bergen, A.W., Chanock, S.J., Yeager, M. (2003). Widespread purifying selection at polymorphic sites in human proteincoding loci. P Natl Acad Sci USA, 100,15754-15757.
Renaut, S., Nolte, A.W., Bernatchez, L. (2010). Mining transcriptome sequences towards identifying adaptive single nucleotide polymorphisms in lake whitefish species pairs (Coregonus spp. Salmonidae). Mol Ecol, 19, 115-131.
Qiu Q, Zhang G, Ma T, Qian W, Wang J, Ye Z, Cao C, Hu Q, Kim J, Larkin DM. (2012). The yak genome and adaptation to life at high altitude. Nature Genetics. 44(8):946–9.
Henrissat, B. (1991). A classification of glycosyl hydrolases based on amino-acid sequence similarities. Biochem. J., 280, pp. 309-316
Lloyd, J.R., Kossmann, J., Ritte, G. (2005). Leaf starch degradation comes out of the shadows. Trends Plant Sci., 10 (3), pp. 130-137
Tsamir‐Rimon, M., Ben‐Dor, S., Feldmesser, E., Oppenhimer‐Shaanan, Y., David‐Schwartz, R., Samach, A., & Klein, T. (2021). Rapid starch degradation in the wood of olive trees under heat and drought is permitted by three stress‐specific beta amylases. New Phytologist, 229(3), 1398-1414.
Bohovych, I., Khalimonchuk, O. (2016). Sending out an SOS: mitochondria as a signaling hub. Frontiers in Cell and Developmental Biology 4: 109.
Nieuwenhuizen, N.J., Chen, X., Wang, M.Y., Matich, A.J., Perez, R.L., Allan, A.C., … Atkinson, R.G. (2015). Natural variation in monoterpene synthesis in kiwifruit: transcriptional regulation of terpene synthases by NAC and ETHYLENEINSENSITIVE3-like transcription factors. Plant Physiol, 167,1243-1258.
Wang, Q., Reddy, V.A., Panicker, D., Mao, H.Z., Kumar, N., Rajan, C., … Sarojam, R. (2016). Metabolic engineering of terpene biosynthesis in plants using a trichomespecific transcription factor MsYABBY5 from spearmint (Mentha spicata). Plant Biotechnol J, 14, 1619-1632.
Costantini, L., Kappel, C.D., Trenti, M., Battilana, J., Emanuelli, F., Sordo, M., … Grando, M.S. (2017). Drawing links from transcriptome to metabolites: the evolution of aroma in the ripening berry of Moscato Bianco (Vitis vinifera L.). Front Plant Sci, 8, 780.
Niu, L., Pan, L., Zeng, W., Lu, Z., Cui, G., Fan, M., … Li, G. (2018). Dynamic transcriptomes of resistant and susceptible peach lines after infestation by green peach aphids (Myzus persicae Sülzer) reveal defence responses controlled by the Rm3 locus. BMC Genomics, 19, 846.
Chang, W., Song, H., Liu, H., Liu, P. (2013). Current development in isoprenoid precursor biosynthesis and regulation. Curr Opin Chem Biol, 17, 571-579.
Karami, A., Khorshidi, R., Jowkar, A. (2015). Diurnal changes in the volatile oil compounds of rose-scented geranium (Pelargonium graveolens). Anal Chem Lett, 5, 103-108.
Ravindra, N.S., Kulkarni, R. (2015). Essential oil yield and quality in rose-scented geranium: Variation among clones and plant parts. Sci Hortic, 184, 31-35.
Tholl, D. (2006). Terpene synthases and the regulation, diversity and biological roles of terpene metabolism. Curr Opin Plant Biol, 9, 297-304.
Cheng, A.X., Lou, Y.G., Mao, Y.B., Lu, S., Wang, L.J., Chen, X.Y. (2007). Plant terpenoids: biosynthesis and ecological functions. J Integr Plant Biol, 49,179-186.
Narnoliya, L.K., Kaushal, G., Singh, S.P., Sangwan, R.S. (2017). De novo transcriptome analysis of rosescented geranium provides insights into the metabolic specificity of terpene and tartaric acid biosynthesis. BMC Genomics, 18, 74.
Gang, D.R. (2001). An investigation of the storage and biosynthesis of phenylpropenes in sweet basil. Plant Physiol, 125, 539-555.
Nazzaro, F., Fratianni, F., De-Martino, L., Coppola, R., De-Feo, V. (2013). Effect of essential oils on pathogenic bacteria. Pharmaceuticals, 6,1451-1474.

Table 1 Transcription factors (TFs) identified in this study.

TFs	Unigene counts	TFs	Unigene counts	TFs	Unigene counts	TFs	Unigene counts
bHLH	67	G2-like	10	C2H2	4	HD-ZIP	2
FAR1	50	Trihelix	10	CO-like	4	LBD	2
NAC	40	bZIP	10	Dof	4	STAT	2
M-type	24	MYB	9	S1Fa-like	4	CAMTA	1
WRKY	22	BES1	8	TCP	4	E2F/DP	1
Nin-like	19	GRAS	8	ERF	3	GATA	1
B3	17	GeBP	7	HB-other	3	HRT-like	1
MYB-related	15	SBP	7	NF-YB	3	NY-YA	1
EIL	13	YABBY	7	NF-YC	3	RAV	1
C3H	11	HSF	5	ARF	2	ZF-HD	1

No competing interests reported.

PDF

Editorial decision: Major revision
04 Jan, 2022
Reviews received at journal
23 Dec, 2021
Reviewers agreed at journal
13 Dec, 2021
Reviewers invited by journal
13 Dec, 2021
Editor assigned by journal
13 Dec, 2021
Editor invited by journal
27 Nov, 2021
Submission checks completed at journal
27 Nov, 2021
First submitted to journal
27 Oct, 2021

You are reading this latest preprint version

Adaptation Insights from Comparative Transcriptome Analysis of Two Opisthopappus Species in the Taihang Mountains

Status:

Version 1

Abstract

Figures

Introduction

Materials And Methods

Results

Discussion

Conclusion

Declarations

References

Tables

Additional Declarations

Status:

Version 1