Transcriptomic response is more sensitive to water deficit in shoots than roots of Vitis riparia (Michx.)

The physiology of drought stress in two grapevine cultivars: Photosynthesis, antioxidant system, and osmotic regulation responses

Drought stress impedes viticultural plant growth and development by modifying various metabolic pathways. However, the regulatory network response underlying drought stress is not yet clear. In this study, the leaves and roots of "Shine Muscat" ("SM," Vitis labruscana × Vitis vinifera) and "Thompson Seedless" ("TS," V. vinifera L. cv.) were subjected to drought stress to study the regulatory network used by drought stress. Morphophysiological results showed that the malondialdehyde content after 28 days of drought stress increased more significantly in "TS" than "SM." Furthermore, the multiomics analysis studies showed that a total of 3036-6714 differentially expressed genes and 379-385 differentially abundant metabolites were identified in "SM" and "TS" grapevine cultivars under drought stress. Furthermore, the retained intron was the major form of differential alternative splicing event under drought stress. The photosynthesis pathway, antioxidant system, plant hormone signal transduction, and osmotic adjustment were the primary response systems in the two grapevine cultivars under drought stress. We have identified GRIK1, RFS2, and LKR/SDH as the hub genes in the coexpression network of drought stress. In addition, the difference in the accumulation of pheophorbide-a reveals different drought resistance mechanisms in the two grapevine cultivars. Our study explained the difference in drought response between cultivars and tissues and identified drought stress-responsive genes, which provides reference data for further understanding the regulatory network of drought tolerance in grapevine.

Network Biology Analyses and Dynamic Modeling of Gene Regulatory Networks under Drought Stress Reveal Major Transcriptional Regulators in Arabidopsis

Drought is one of the most serious abiotic stressors in the environment, restricting agricultural production by reducing plant growth, development, and productivity. To investigate such a complex and multifaceted stressor and its effects on plants, a systems biology-based approach is necessitated, entailing the generation of co-expression networks, identification of high-priority transcription factors (TFs), dynamic mathematical modeling, and computational simulations. Here, we studied a high-resolution drought transcriptome of Arabidopsis. We identified distinct temporal transcriptional signatures and demonstrated the involvement of specific biological pathways. Generation of a large-scale co-expression network followed by network centrality analyses identified 117 TFs that possess critical properties of hubs, bottlenecks, and high clustering coefficient nodes. Dynamic transcriptional regulatory modeling of integrated TF targets and transcriptome datasets uncovered major transcriptional events during the course of drought stress. Mathematical transcriptional simulations allowed us to ascertain the activation status of major TFs, as well as the transcriptional intensity and amplitude of their target genes. Finally, we validated our predictions by providing experimental evidence of gene expression under drought stress for a set of four TFs and their major target genes using qRT-PCR. Taken together, we provided a systems-level perspective on the dynamic transcriptional regulation during drought stress in Arabidopsis and uncovered numerous novel TFs that could potentially be used in future genetic crop engineering programs.

Improving abiotic stress tolerance of forage grasses - prospects of using genome editing

Due to an increase in the consumption of food, feed, and fuel and to meet global food security needs for the rapidly growing human population, there is a necessity to obtain high-yielding crops that can adapt to future climate changes. Currently, the main feed source used for ruminant livestock production is forage grasses. In temperate climate zones, perennial grasses grown for feed are widely distributed and tend to suffer under unfavorable environmental conditions. Genome editing has been shown to be an effective tool for the development of abiotic stress-resistant plants. The highly versatile CRISPR-Cas system enables increasingly complex modifications in genomes while maintaining precision and low off-target frequency mutations. In this review, we provide an overview of forage grass species that have been subjected to genome editing. We offer a perspective view on the generation of plants resilient to abiotic stresses. Due to the broad factors contributing to these stresses the review focuses on drought, salt, heat, and cold stresses. The application of new genomic techniques (e.g., CRISPR-Cas) allows addressing several challenges caused by climate change and abiotic stresses for developing forage grass cultivars with improved adaptation to the future climatic conditions. Genome editing will contribute towards developing safe and sustainable food systems.

Genetic analysis of grapevine root system architecture and loci associated gene networks

Own-rooted grapevines and grapevine rootstocks are vegetatively propagated from cuttings and have an adventitious root system. Unraveling the genetic underpinnings of the adventitious root system architecture (RSA) is important for improving own-rooted and grafted grapevine sustainability for a changing climate. Grapevine RSA genetic analysis was conducted in an Vitis sp. 'VRS-F2' population. Nine root morphology, three total root system morphology, and two biomass traits that contribute to root anchorage and water and nutrient uptake were phenotyped. Quantitative trait loci (QTL) analysis was performed using a high density integrated GBS and rhAmpSeq genetic map. Thirty-one QTL were detected for eleven of the RSA traits (surface area, root volume, total root length, fresh weight, number of tips, forks or links, longest root and average root diameter, link length, and link surface area) revealing many small effects. Several QTL were colocated on chromosomes 1, 9, 13, 18, and 19. QTL with identical peak positions on chromosomes 1 or 13 were enriched for AP2-EREBP, AS2, C2C2-CO, HMG, and MYB transcription factors, and QTL on chromosomes 9 or 13 were enriched for the ALFIN-LIKE transcription factor and regulation of autophagy pathways. QTL modeling for individual root traits identified eight models explaining 13.2 to 31.8% of the phenotypic variation. 'Seyval blanc' was the grandparent contributing to the allele models that included a greater surface area, total root length, and branching (number of forks and links) traits promoting a greater root density. In contrast, V. riparia 'Manitoba 37' contributed the allele for greater average branch length (link length) and diameter, promoting a less dense elongated root system with thicker roots. LATERAL ORGAN BOUNDARY DOMAIN (LBD or AS2/LOB) and the PROTODERMAL FACTOR (PFD2 and ANL2) were identified as important candidate genes in the enriched pathways underlying the hotspots for grapevine adventitious RSA. The combined QTL hotspot and trait modeling identified transcription factors, cell cycle and circadian rhythm genes with a known role in root cell and epidermal layer differentiation, lateral root development and cortex thickness. These genes are candidates for tailoring grapevine root system texture, density and length in breeding programs.

Combined analysis of transcriptome and metabolome reveals the molecular mechanism and candidate genes of Haloxylon drought tolerance

Haloxylon ammodendron and Haloxylon persicum, as typical desert plants, show strong drought tolerance and environmental adaptability. They are ideal model plants for studying the molecular mechanisms of drought tolerance. Transcriptomic and metabolomic analyses were performed to reveal the response mechanisms of H. ammodendron and H. persicum to a drought environment at the levels of transcription and physiological metabolism. The results showed that the morphological structures of H. ammodendron and H. persicum showed adaptability to drought stress. Under drought conditions, the peroxidase activity, abscisic acid content, auxin content, and gibberellin content of H. ammodendron increased, while the contents of proline and malondialdehyde decreased. The amino acid content of H. persicum was increased, while the contents of proline, malondialdehyde, auxin, and gibberellin were decreased. Under drought conditions, 12,233 and 17,953 differentially expressed genes (DEGs) were identified in H. ammodendron and H. persicum , respectively, including members of multiple transcription factor families such as FAR1, AP2/ERF, C2H2, bHLH, MYB, C2C2, and WRKY that were significantly up-regulated under drought stress. In the positive ion mode, 296 and 452 differential metabolites (DEMs) were identified in H. ammodendron and H. persicum, respectively; in the negative ion mode, 252 and 354 DEMs were identified, primarily in carbohydrate and lipid metabolism. A combined transcriptome and metabolome analysis showed that drought stress promoted the glycolysis/gluconeogenesis pathways of H. ammodendron and H. persicum and increased the expression of amino acid synthesis pathways, consistent with the physiological results. In addition, transcriptome and metabolome were jointly used to analyze the expression changes of the genes/metabolites of H. ammodendron and H. persicum that were associated with drought tolerance but were regulated differently in the two plants. This study identified drought-tolerance genes and metabolites in H. ammodendron and H. persicum and has provided new ideas for studying the drought stress response of Haloxylon.

Wild grapevines as rootstock regulate the oxidative defense system of in vitro grafted scion varieties under drought stress

The narrow genetic base of modern cultivars is becoming a key bottleneck for crop improvement and the use of wild relatives is an appropriate approach to improve the genetic diversity of crops to manage the sustainable production under different abiotic and biotic constraints. In Pakistan, wild germplasm of grapevine viz Dakh, Toran, and Zarishk belong to Vitis vinifera subsp. sylvestris and Fatati belong to Vitis vinifera subsp. sativa is naturally present in humid and sub-humid areas of mountainous and sub-mountainous regions and showed varying level of tolerance against drought stress but have not been evaluated as rootstock. In this study, different tolerant behavior of wild grapevines as rootstock in grafted scion varieties were explored under different levels of PEG-6000 mediated drought stress i.e., -4.00, -6.00, and -8.00 bars. In response to drought stress, wild grapevines evoked several non-enzymatic and enzymatic activities. Among non-enzymatic activities, total chlorophyll contents of commercial varieties were sustained at higher level when grafted on wild grapevines Dakh and Fatati which subsequently reduced the damage of cell membrane via MDA. Whereas, to cope the membranous damage due to excessive cellular generation of ROS, wild grapevines triggered the enhanced activities of SOD to dismutase the free oxygen radicals into H₂O₂, then CAT enzyme convert the H₂O₂ into water molecules. Higher accumulation of ROS in commercial scion varieties were also coped by wild grapevines Dakh and Fatati through the upregulation of POD and APX enzymes activities. Based on these enzymatic and non-enzymatic indices, biplot and cluster analysis classified the wild grapevines as rootstock into three distinct categories comprises on relatively tolerant i.e., Dakh (Vitis vinifera subsp. sylvestris) and Fatati (Vitis vinifera subsp. sativa), moderate tolerant i.e., Toran (Vitis vinifera subsp. sylvestris) and relatively susceptible category i.e., Zarishk (Vitis vinifera subsp. sylvestris).

Characterization of the gene expression profile response to drought stress in Haloxylon using PacBio single-molecule real-time and Illumina sequencing

Haloxylon ammodendron and Haloxylon persicum are important drought-tolerant plants in northwest China. The whole-genome sequencing of H. ammodendron and H. persicum grown in their natural environment is incomplete, and their transcriptional regulatory network in response to drought environment remains unclear. To reveal the transcriptional responses of H. ammodendron and H. persicum to an arid environment, we performed single-molecule real-time (SMRT) and Illumina RNA sequencing. In total, 20,246,576 and 908,053 subreads and 435,938 and 210,334 circular consensus sequencing (CCS) reads were identified by SMRT sequencing of H. ammodendron and H. persicum, and 15,238 and 10,135 unigenes, respectively, were successfully obtained. In addition, 9,794 and 7,330 simple sequence repeats (SSRs) and 838 and 71 long non-coding RNAs were identified. In an arid environment, the growth of H. ammodendron was restricted; plant height decreased significantly; basal and branch diameters became thinner and hydrogen peroxide (H₂O₂) content and peroxidase (POD) activity were increased. Under dry and wet conditions, 11,803 and 15,217 differentially expressed genes (DEGs) were identified in H. ammodendron and H. persicum, respectively. There were 319 and 415 DEGs in the signal transduction pathways related to drought stress signal perception and transmission, including the Ca²⁺ signal pathway, the ABA signal pathway, and the MAPK signal cascade. In addition, 217 transcription factors (TFs) and 398 TFs of H. ammodendron and H. persicum were differentially expressed, including FAR1, MYB, and AP2/ERF. Bioinformatic analysis showed that under drought stress, the expression patterns of genes related to active oxygen [reactive oxygen species (ROS)] scavenging, functional proteins, lignin biosynthesis, and glucose metabolism pathways were altered. Thisis the first full-length transcriptome report concerning the responses of H. ammodendron and H. persicum to drought stress. The results provide a foundation for further study of the adaptation to drought stress. The full-length transcriptome can be used in genetic engineering research.

Comparative Transcriptomic Analysis of Root and Leaf Transcript Profiles Reveals the Coordinated Mechanisms in Response to Salinity Stress in Common Vetch

Owing to its strong environmental suitability to adverse abiotic stress conditions, common vetch (Vicia sativa) is grown worldwide for both forage and green manure purposes and is an important protein source for human consumption and livestock feed. The germination of common vetch seeds and growth of seedlings are severely affected by salinity stress, and the response of common vetch to salinity stress at the molecular level is still poorly understood. In this study, we report the first comparative transcriptomic analysis of the leaves and roots of common vetch under salinity stress. A total of 6361 differentially expressed genes were identified in leaves and roots. In the roots, the stress response was dominated by genes involved in peroxidase activity. However, the genes in leaves focused mainly on Ca²⁺ transport. Overexpression of six salinity-inducible transcription factors in yeast further confirmed their biological functions in the salinity stress response. Our study provides the most comprehensive transcriptomic analysis of common vetch leaf and root responses to salinity stress. Our findings broaden the knowledge of the common and distinct intrinsic molecular mechanisms within the leaves and roots of common vetch and could help to develop common vetch cultivars with high salinity tolerance.

Transcriptome analysis and phenotyping of walnut seedling roots under nitrogen stresses

Nitrogen is an essential core element in walnut seedling growth and development. However, nitrogen starvation and excessive nitrogen stress can cause stunted growth and development of walnut seedlings, and environmental pollution is also of concern. Therefore, it is necessary to study the mechanism of walnut seedling resistance to nitrogen stress. In this study, morphological and physiological observations and transcriptome sequencing of walnut seedlings under nitrogen starvation and excess nitrogen stress were performed. The results showed that walnut seedlings under nitrogen starvation and excess stress could adapt to the changes in the nitrogen environment by changing the coordination of their root morphology and physiological indexes. Based on an analysis of transcriptome data, 4911 differential genes (DEGs) were obtained (2180 were upregulated and 2731 were downregulated) in a comparison of nitrogen starvation and control groups. A total of 9497 DEGs (5091 upregulated and 4406 downregulated) were obtained in the comparison between the nitrogen overdose and control groups. When these DEGs were analysed, the differential genes in both groups were found to be significantly enriched in the plant’s circadian pathway. Therefore, we selected the circadian rhythm as the focus for further analysis. We made some discoveries by analysing the gene co-expression network of nitrogen metabolism, circadian rhythm, and hormone signal transduction. (a) Nitrite nitrogen (NO₂⁻) or Glu may act as a nitrogen signal to the circadian clock. (b) Nitrogen signalling may be input into the circadian clock by regulating changes in the abundance of the CRY1 gene. (c) After the nitrogen signal enters the circadian clock, the expression of the LHY gene is upregulated, which causes a phase shift in the circadian clock. (d) The RVE protein may send information about the circadian clock’s response to nitrogen stress back to the nitrogen metabolic pathway via the hormone transduction pathway. In conclusion, various metabolic pathways in the roots of walnut seedlings coordinated with one another to resist the ill effects of nitrogen stress on the root cells, and these coordination relationships were regulated by the circadian clock. This study is expected to provide valuable information on the circadian clock regulation of plant resistance to nitrogen stress.

Characterization of the Gene Expression Profile Response to Drought Stress in Populus ussuriensis Using PacBio SMRT and Illumina Sequencing

In this study, we characterized the gene expression profile in the roots of Populus ussuriensis at 0, 6, 12, 24, 48 and 120 h after the start of polyethylene glycol (PEG)-induced drought stress using PacBio single-molecule real-time sequencing (SMRT-seq) and Illumina RNA sequencing. Compared to the control, 2244 differentially expressed genes (DEGs) were identified, and many of these DEGs were associated with the signal transduction, antioxidant system, ion accumulation and drought-inducing proteins. Changes in certain physiological and biochemical indexes, such as antioxidant activity and the contents of Ca²⁺, proline, and total soluble sugars, were further confirmed in P. ussuriensis roots. Furthermore, most of the differentially expressed transcription factors were members of the AP2/ERF, C2H2, MYB, NAC, C2C2 and WRKY families. Additionally, based on PacBio SMRT-seq results, 5955 long non-coding RNAs and 700 alternative splicing events were identified. Our results provide a global view of the gene expression profile that contributes to drought resistance in P. ussuriensis and meaningful information for genetic engineering research in the future.

Grapevine aquaporins: Diversity, cellular functions, and ecophysiological perspectives

High-scored premium wines are typically produced under moderate drought stress, suggesting that the water status of grapevine is crucial for wine quality. Aquaporins greatly influence the plant water status by facilitating water diffusion across the plasma membrane in a tightly regulated manner. They adjust the hydraulic conductance of the plasma membrane rapidly and reversibly, which is essential in specific physiological events, including adaptation to soil water scarcity. The comprehension of the sophisticated plant-water relations at the molecular level are thus important to optimize agricultural practices or to assist plant breeding programs. This review explores the recent progresses in understanding the water transport in grapevine at the cellular level through aquaporins and its regulation. Important aspects, including aquaporin structure, diversity, cellular localization, transport properties, and regulation at the cellular and whole plant level are addressed. An ecophysiological perspective about the roles of grapevine aquaporins in plant response to drought stress is also provided.

Comparative transcriptome analyses between cultivated and wild grapes reveal conservation of expressed genes but extensive rewiring of co-expression networks

The transcriptomes of wild and cultivated grapes consists of similar expressed genes but distinct wiring of co-expressed modules associated with environmental conditions. Grapevine is an important fruit crop worldwide, with high economic value and widespread distribution. Commercial production is based on Vitis vinifera, and, to a lesser extent, on hybrids with American grapes, such as V. labrusca. Wild grape relatives are important sources of resistance against biotic and abiotic factors; however, their global gene expression patterns remain poorly characterized. We associated genome-wide transcript profiling to phenotypic analyses to investigate the responses of cultivated and wild vines to vineyard conditions. The expressed genes in the Vitis reference transcriptome are largely shared by wild grapes, V. labrusca hybrids and vinifera cultivars. In contrast, significant differential regulation between wild and vinifera genotypes represents 80% of gene expression variation, regardless of the environment. In wild grapes, genes associated to regulatory processes are downregulated, whereas those involved in metabolic pathways are upregulated, in comparison to vinifera. Photosynthesis-related ontologies are overrepresented in the induced genes, in agreement with higher contents of chlorophyll in wild grapes. Co-regulated gene network analyses provide evidence of more complex transcriptome organization in vinifera. In wild grapes, genes involved in signaling pathways of stress-related hormones are overrepresented in modules associated with the environment. Consensus network analyses revealed high preservation within co-regulated gene modules between cultivated and wild grapes, but divergent relationships among the expression clusters. In conclusion, the distinct phenotypes of wild and cultivated grapes are underlain by differences in gene expression, but also by distinct higher-order organization of the transcriptome and contrasting association of co-expressed gene clusters with the environment.

Gaining Insight into Exclusive and Common Transcriptomic Features Linked to Drought and Salinity Responses across Fruit Tree Crops

The present study aimed at identifying and mapping key genes expressed in root tissues involved in drought and salinity tolerance/resistance conserved among different fruit tree species. Twenty-six RNA-Seq samples were analyzed from six published studies in five plant species (Olea europaea, Vitis riparia Michx, Prunus mahaleb, Prunus persica, Phoenix dactylifera). This meta-analysis used a bioinformatic pipeline identifying 750 genes that were commonly modulated in three salinity studies and 683 genes that were commonly regulated among three drought studies, implying their conserved role in resistance/tolerance/response to these environmental stresses. A comparison was done on the genes that were in common among both salinity and drought resulted in 82 genes, of which 39 were commonly regulated with the same trend of expression (23 were upregulated and 16 were downregulated). Gene set enrichment and pathway analysis pointed out that pathways encoding regulation of defense response, drug transmembrane transport, and metal ion binding are general key molecular responses to these two abiotic stress responses. Furthermore, hormonal molecular crosstalk plays an essential role in the fine-tuning of plant responses to drought and salinity. Drought and salinity induced a different molecular "hormonal fingerprint". Dehydration stress specifically enhanced multiple genes responsive to abscisic acid, gibberellin, brassinosteroids, and the ethylene-activated signaling pathway. Salt stress mostly repressed genes encoding for key enzymes in signaling proteins in auxin-, gibberellin-(gibberellin 2 oxidase 8), and abscisic acid-related pathways (aldehyde oxidase 4, abscisic acid-responsive element-binding protein 3). Abiotic stress-related genes were mapped into the chromosome to identify molecular markers usable for the improvement of these complex quantitative traits. This meta-analysis identified genes that serve as potential targets to develop cultivars with enhanced drought and salinity resistance and/or tolerance across different fruit tree crops in a biotechnological sustainable way.

Plant growth under suboptimal water conditions: early responses and methods to study them

Drought stress forms a major environmental constraint during the life cycle of plants, often decreasing plant yield and in extreme cases threatening survival. The molecular and physiological responses induced by drought have been the topic of extensive research during the past decades. Because soil-based approaches to studying drought responses are often challenging due to low throughput and insufficient control of the conditions, osmotic stress assays in plates were developed to mimic drought. Addition of compounds such as polyethylene glycol, mannitol, sorbitol, or NaCl to controlled growth media has become increasingly popular since it offers the advantage of accurate control of stress level and onset. These osmotic stress assays enabled the discovery of very early stress responses, occurring within seconds or minutes following osmotic stress exposure. In this review, we construct a detailed timeline of early responses to osmotic stress, with a focus on how they initiate plant growth arrest. We further discuss the specific responses triggered by different types and severities of osmotic stress. Finally, we compare short-term plant responses under osmotic stress versus in-soil drought and discuss the advantages, disadvantages, and future of these plate-based proxies for drought.

Drought tolerance of the grapevine, Vitis champinii cv. Ramsey, is associated with higher photosynthesis and greater transcriptomic responsiveness of abscisic acid biosynthesis and signaling

BACKGROUND

Grapevine is an economically important crop for which yield and berry quality is strongly affected by climate change. Large variations in drought tolerance exist across Vitis species. Some of these species are used as rootstock to enhance abiotic and biotic stress tolerance. In this study, we investigated the physiological and transcriptomic responses to water deficit of four different genotypes that differ in drought tolerance: Ramsey (Vitis champinii), Riparia Gloire (Vitis riparia), Cabernet Sauvignon (Vitis vinifera), and SC2 (Vitis vinifera x Vitis girdiana).

RESULTS

Ramsey was particularly more drought tolerant than the other three genotypes. Ramsey maintained a higher stomatal conductance and photosynthesis at equivalent levels of moderate water deficit. We identified specific and common transcriptomic responses shared among the four different Vitis species using RNA sequencing analysis. A weighted gene co-expression analysis identified a water deficit core gene set with the ABA biosynthesis and signaling genes, NCED3, RD29B and ABI1 as potential hub genes. The transcript abundance of many abscisic acid metabolism and signaling genes was strongly increased by water deficit along with genes associated with lipid metabolism, galactinol synthases and MIP family proteins. This response occurred at smaller water deficits in Ramsey and with higher transcript abundance than the other genotypes. A number of aquaporin genes displayed differential and unique responses to water deficit in Ramsey leaves. Genes involved in cysteine biosynthesis and metabolism were constitutively higher in the roots of Ramsey; thus, linking the gene expression of a known factor that influences ABA biosynthesis to this genotype's increased NCED3 transcript abundance.

CONCLUSION

The drought tolerant Ramsey maintained higher photosynthesis at equivalent water deficit than the three other grapevine genotypes. Ramsey was more responsive to water deficit; its transcriptome responded at smaller water deficits, whereas the other genotypes did not respond until more severe water deficits were reached. There was a common core gene network responding to water deficit for all genotypes that included ABA metabolism and signaling. The gene clusters and sub-networks identified in this work represent interesting gene lists to explore and to better understand drought tolerance molecular mechanisms.

Water Deficit Transcriptomic Responses Differ in the Invasive Tamarix chinensis and T. ramosissima Established in the Southern and Northern United States

Tamarix spp. (saltcedar) were introduced from Asia to the southern United States as windbreak and ornamental plants and have spread into natural areas. This study determined differential gene expression responses to water deficit (WD) in seedlings of T. chinensis and T. ramosissima from established invasive stands in New Mexico and Montana, respectively. A reference de novo transcriptome was developed using RNA sequences from WD and well-watered samples. Blast2GO analysis of the resulting 271,872 transcripts yielded 89,389 homologs. The reference Tamarix (Tamaricaceae, Carophyllales order) transcriptome showed homology with 14,247 predicted genes of the Beta vulgaris subsp. vulgaris (Amaranthaceae, Carophyllales order) genome assembly. T. ramosissima took longer to show water stress symptoms than T. chinensis. There were 2068 and 669 differentially expressed genes (DEG) in T. chinensis and T. ramosissima, respectively; 332 were DEG in common between the two species. Network analysis showed large biological process networks of similar gene content for each of the species under water deficit. Two distinct molecular function gene ontology networks (binding and transcription factor-related) encompassing multiple up-regulated transcription factors (MYB, NAC, and WRKY) and a cellular components network containing many down-regulated photosynthesis-related genes were identified in T. chinensis, in contrast to one small molecular function network in T. ramosissima.