Chromosomal-scale genomes of two Rosa species provide insights into genome evolution and ascorbate accumulation.

Rosa roxburghii and Rosa sterilis, two species belonging to the Rosaceae family, are widespread in the southwest of China. These species have gained recognition for their remarkable abundance of ascorbate in their fresh fruits, making them an ideal vitamin C resource. In this study, we generated two high-quality chromosome-scale genome assemblies for R. roxburghii and R. sterilis, with genome sizes of 504 and 981.2 Mb, respectively. Notably, we present a haplotype-resolved, chromosome-scale assembly for diploid R. sterilis. Our results indicated that R. sterilis originated from the hybridization of R. roxburghii and R. longicuspis. Genome analysis revealed the absence of recent whole-genome duplications in both species and identified a series of duplicated genes that possibly contributing to the accumulation of flavonoids. We identified two genes in the ascorbate synthesis pathway, GGP and GalLDH, that show signs of positive selection, along with high expression levels of GDP-d-mannose 3', 5'-epimerase (GME) and GDP-l-galactose phosphorylase (GGP) during fruit development. Furthermore, through co-expression network analysis, we identified key hub genes (MYB5 and bZIP) that likely regulate genes in the ascorbate synthesis pathway, promoting ascorbate biosynthesis. Additionally, we observed the expansion of terpene synthase genes in these two species and tissue expression patterns, suggesting their involvement in terpenoid biosynthesis. Our research provides valuable insights into genome evolution and the molecular basis of the high concentration of ascorbate in these two Rosa species.

Chloroplast genome sequence of triploid Toxicodendron vernicifluum and comparative analyses with other lacquer chloroplast genomes.

BACKGROUND: Toxicodendron vernicifluum, belonging to the family Anacardiaceae, is an important commercial arbor species, which can provide us with the raw lacquer, an excellent adhesive and painting material used to make lacquer ware. Compared with diploid, triploid lacquer tree has a higher yield of raw lacquer and stronger resistance to stress. Triploid T. vernicifluum was a newly discovered natural triploid lacquer tree. However, the taxonomy of triploid T. vernicifluum has remained uncertain. Here, we sequenced and analyzed the complete chloroplast (cp) genome of triploid T. vernicifluum and compared it with related species of Toxicodendron genus based on chloroplast genome and SSR markers. RESULTS: The plastome of triploid T. vernicifluum is 158,221 bp in length, including a pair of inverted repeats (IRs) of 26,462 bp, separated by a large single-copy region of 86,951 bp and a small single-copy region of 18,346 bp. In total, 132 genes including 87 protein-coding genes, 37 tRNA genes and 8 rRNA genes were identified in the triploid T. vernicifluum. Among these, 16 genes were duplicated in the IR regions, 14 genes contain one intron, while three genes contain two introns. After nucleotide substitutions, seven small inversions were analyzed in the chloroplast genomes, eight hotspot regions were found, which could be useful molecular genetic markers for future population genetics. Phylogenetic analyses showed that triploid T. vernicifluum was a sister to T. vernicifluum cv. Dahongpao and T. vernicifluum cv. Hongpigaobachi. Moreover, phylogenetic clustering based on the SSR markers showed that all the samples of triploid T. vernicifluum, T. vernicifluum cv. Dahongpao and T. vernicifluum cv. Hongpigaobachi in one group, while the samples of T. vernicifluum and T. succedaneum in another group, which is consistent with the cp genome and morphological analysis. CONCLUSIONS: The current genomic datasets provide pivotal genetic resources to determine the phylogenetic relationships, variety identification, breeding and resource exploitation, and future genetic diversity-related studies of T. vernicifluum.

Comparative Transcriptomic and Metabolomic Analyses of Differences in Trunk Spiral Grain in Pinus yunnanensis.

Having a spiral grain is considered to be one of the most important wood properties influencing wood quality. Here, transcriptome profiles and metabolome data were analyzed in the straight grain and twist grain of Pinus yunnanensis. A total of 6644 differential expression genes were found between the straight type and the twist type. A total of 126 differentially accumulated metabolites were detected. There were 24 common differential pathways identified from the transcriptome and metabolome, and these pathways were mainly annotated in ABC transporters, arginine and proline metabolism, flavonoid biosynthesis, isoquinoline alkaloid biosynthesis, linoleic acid metabolism, phenylpropanoid, tryptophan metabolism, etc. A weighted gene coexpression network analysis showed that the lightblue4 module was significantly correlated with 2'-deoxyuridine and that transcription factors (basic leucine zipper (bZIP), homeodomain leucine zipper (HD-ZIP), basic helix-loop-helix (bHLH), p-coumarate 3-hydroxylase (C3H), and N-acetylcysteine (NAC)) play important roles in regulating 2'-deoxyuridine, which may be involved in the formation of spiral grains. Meanwhile, the signal transduction of hormones may be related to spiral grain, as previously reported. ARF7 and MKK4_5, as indoleacetic acid (IAA)- and ethylene (ET)-related receptors, may explain the contribution of plant hormones in spiral grain. This study provided useful information on spiral grain in P. yunnanensis by transcriptome and metabolome analyses and could lay the foundation for future molecular breeding.

Comparison of the diversity and structure of the rhizosphere microbial community between the straight and twisted trunk types of Pinus yunnanensis.

Background: Pinus yunnanensis is a major silvicultural species in Southwest China. Currently, large areas of twisted-trunk Pinus yunnanensis stands severely restrict its productivity. Different categories of rhizosphere microbes evolve alongside plants and environments and play an important role in the growth and ecological fitness of their host plant. However, the diversity and structure of the rhizosphere microbial communities between P. yunnanensis with two different trunk types-straight and twisted-remain unclear. Methods: We collected the rhizosphere soil of 5 trees with the straight and 5 trees with the twisted trunk type in each of three sites in Yunnan province. We assessed and compared the diversity and structure of the rhizosphere microbial communities between P. yunnanensis with two different trunk types by Illumina sequencing of 16S rRNA genes and internal transcribed spacer (ITS) regions. Results: The available phosphorus in soil differed significantly between P. yunnanensis with straight and twisted trunks. Available potassium had a significant effect on fungi. Chloroflexi dominated the rhizosphere soils of the straight trunk type, while Proteobacteria was predominant in the rhizosphere soils of the twisted trunk type. Trunk types significantly explained 6.79% of the variance in bacterial communities. Conclusion: This study revealed the composition and diversity of bacterial and fungal groups in the rhizosphere soil of P. yunnanensis with straight and twisted trunk types, providing proper microbial information for different plant phenotypes.

Genome-wide identification of polar auxin transporter gene families reveals a possible new polar auxin flow in inverted cuttings of Populus yunnanensis.

Inverted cuttings of Populus yunnanensis are characterized by enlarged stems and dwarfed new shoots, and phytohormones play a crucial role in the response to inversion. The polar auxin transport (PAT) system is distinct from the transport systems of other hormones and is controlled by three major transporter gene families: pin-formed (PIN), auxin-resistant/like aux (AUX/LAX) and ATP-binding cassette transporters of the B class (ABCB). Here, we identified these three families in P. trichocarpa, P. euphratica and P. yunnanensis through a genome-wide analysis. The Populus PIN, AUX/LAX and ABCB gene families comprised 15, 8 and 31 members, respectively. Most PAT genes in Populus and Arabidopsis were identified as clear sister pairs, and some had unique motifs. Transcriptome profiling revealed that the expression of most PAT genes was unrelated to cutting inversion and that only several genes showed altered expression when cuttings were inverted. The auxin content difference at positions was opposite in upright and inverted cutting bodies during rooting, which obeyed the original plant polarity. However, during plant growth, the two direction types exhibited similar auxin movements in the cutting bodies, and the opposite auxin changes were observed in new shoots. Four PAT genes with a positive response to cutting inversion, PyuPIN10, PyuPIN11, PyuLAX6 and PyuABCB27, showed diverse expression patterns between upright and inverted cuttings during rooting and plant growth. Furthermore, PAT gene expression retained its polarity, which differs from the results found for auxin flow during plant growth. The inconformity indicated that a new downward auxin flow in addition to the old upward flow might be established during the growth of inverted cuttings. Some highly polar PAT genes were involved in the maintenance of original auxin polarity, which might cause the enlarged stems of inverted cuttings. This work lays a foundation for understanding the roles of auxin transport in plant responses to inversion.

Comparative analysis of the complete chloroplast genomes of seven Populus species: Insights into alternative female parents of Populus tomentosa.

Populus tomentosa, of section Populus, is distributed mainly in northern China. This species has high resistance to many diseases and insects, and it plays key roles in shelterbelts and urban afforestation in northern China. It has long been suspected to be a hybrid, but its parents remain unknown. In the present study, we report four newly sequenced complete cp genomes from section Populus and comparative genomic analyses of these new sequences and three published cp genome sequences. The seven cp genomes ranged from 155,853 bp (P. tremula var. davidiana) to 156,746 bp (P. adenopoda) in length, and their gene orders, gene numbers and GC contents were similar. We analyzed SNPs, indels, SSRs and repeats among the seven cp genomes, and eight small inversions were detected in the ndhC-trnV, rbcL-accD, petA-psbJ, trnW-trnP, rpl16-rps3, trnL-ycf15, ycf15-trnL, and ndhF-trnL intergenic regions. Furthermore, seven divergent regions (trnH-psbA, matK, psbM-psbD, ndhC-trnV, ycf1, ndhF-ccsA and ccsA-ndhD) were found in more highly variable regions. The phylogenetic tree reveals that P. tomentosa is closely related to P. alba and P. alba var. pyramidalis. Hence, P. alba was involved in the formation of P. tomentosa.

Plastome Sequences Help to Resolve Deep-Level Relationships of Populus in the Family Salicaceae.

Populus, a core genus of Salicaceae, plays a significant ecological role as a source of pioneer species in boreal forests. However, interspecific hybridization and high levels of morphological variation among poplars have resulted in great difficulty in classifying species for systematic and comparative evolutionary studies. Here, we present phylogenetic analyses of 24 newly sequenced Populus plastomes and 36 plastomes from GenBank, which represent seven genera of Salicaceae, in combination with a matrix of eighteen morphological characters of 40 Populus taxa to reconstruct highly supported relationships of genus Populus. Relationships among the 60 taxa of Salicaceae strongly supported two monophyletic genera: Populus and Salix. Chosenia was nested within the genus Salix, and five clades within Populus were divided. Clade I included the three taxa P. euphratica, P. pruinosa, and P. ilicifolia. Clade II contained thirteen taxa [P. adenopoda, P. alba, P. bolleana, P. davidiana, P. hopeiensis, P. nigra, P. qiongdaoensis, P. rotundifolia, P. rotundifolia var. duclouxiana, P. tremula, P. tremula × alba, P. tomentosa, and P. tomentosa (NC)]. Clade III included the ten taxa P. haoana, P. kangdingensis, P. lasiocarpa, P. pseudoglauca, P. qamdoensis, P. schneideri, P. simonii, P. szechuanica, P. szechuanica var. tibetica, and P. yunnanensis. Clade IV included P. cathayana, P. gonggaensis, P. koreana, P. laurifolia, P. trinervis, P. wilsonii, and P. xiangchengensis. The last clade comprised P. angustifolia, P. balsamifera, P. deltoides, P. deltoides × nigra, P. fremontii, P. mexicana, and P. trichocarpa. This phylogeny is also supported by morphological traits, including bark smoothness, bud size, petiole shape, leaf inflorescence, male anther length and male anther tip.

Bark tissue transcriptome analyses of inverted Populus yunnanensis cuttings reveal the crucial role of plant hormones in response to inversion.

Inverted cuttings of Populus yunnanensis exhibit an interesting growth response to inversion. This response is characterized by enlargement of the stem above the shoot site, while the upright stem shows obvious outward growth below the shoot site. In this study, we examined transcriptome changes in bark tissue at four positions on upright and inverted cuttings of P. yunnanensis: position B, the upper portion of the stem; position C, the lower portion of the stem; position D, the bottom of new growth; and position E, the top of new growth. The results revealed major transcriptomic changes in the stem, especially at position B, but little alteration was observed in the bark tissue of the new shoot. The differentially expressed genes (DEGs) were mainly assigned to four pathways: plant hormone signal transduction, plant-pathogen interaction, mitogen-activated protein kinase (MAPK) signaling pathway-plant, and adenosine triphosphate-binding cassette (ABC) transporters. Most of these DEGs were involved in at least two pathways. The levels of many hormones, such as auxin (IAA), cytokinin (CTK), gibberellins (GAs), ethylene (ET), and brassinosteroids (BRs), underwent large changes in the inverted cuttings. A coexpression network showed that the top 20 hub unigenes at position B in the upright and inverted cutting groups were associated mainly with the BR and ET signaling pathways, respectively. Furthermore, brassinosteroid insensitive 1-associated receptor kinase 1 (BAK1) in the BR pathway and both ethylene response (ETR) and constitutive triple response 1 (CTR1) in the ET pathway were important hubs that interfaced with multiple pathways.

Physiological Analysis and Transcriptome Profiling of Inverted Cuttings of Populus yunnanensis Reveal That Cell Wall Metabolism Plays a Crucial Role in Responding to Inversion.

Inverted cuttings of Populus yunnanensis remain alive by rooting from the original morphological apex and sprouting from the base, but the lateral branches exhibit less vigorous growth than those of the upright plant. In this study, we examined the changes in hormone contents, oxidase activities, and transcriptome profiles between upright and inverted cuttings of P. yunnanensis. The results showed that the indole-3-acetic acid (IAA) and gibberellic acid (GA3) contents were significantly lower in inverted cuttings than in upright cuttings only in the late growth period (September and October), while the abscisic acid (ABA) level was always similar between the two direction types. The biosynthesis of these hormones was surprisingly unrelated to the inversion of P. yunnanensis during the vegetative growth stage (July and August). Increased levels of peroxidases (PODs) encoded by 13 differentially expressed genes (DEGs) served as lignification promoters that protected plants against oxidative stress. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis showed that most DEGs (107) were related to carbohydrate metabolism. Furthermore, altered activities of uridine diphosphate (UDP)-sugar pyrophosphorylase (USP, 15 DEGs) for nucleotide sugars, pectin methylesterase (PME, 7 DEGs) for pectin, and POD (13 DEGs) for lignin were important factors in the response of the trees to inversion, and these enzymes are all involved cell wall metabolism.

Analyzing the phylogeny of poplars based on molecular data.

Methods for constructing trees using DNA sequences, known as molecular phylogenetics, have been applied to analyses of phylogenetic origin, evolutionary relatedness and taxonomic classification. Combining data sequenced in this study and downloaded from GenBank, we sampled 112 (chloroplast data) / 122 (ITS data) specimens belonging to 49 (chloroplast data) / 46 (ITS data) poplar species or hybrids from six (chloroplast data) / five sections (ITS data). Maximum parsimony and Bayesian inference were used to analyze phylogenetic relationships within the genus Populus based on eight chloroplast combinations and ITS regions. The results suggested that Bayesian inference might be more suitable for the phylogenetic reconstruction of Populus. All Populus species could be divided into two clades: clade 1, including subclades 1 and 2, and clade 2, including subclades 3 and 4. Species within clade 1, involving five sections except for Leuce, clustered coinciding with their two specific main geographical distribution areas: China (subclade 1) and North America (subclade 2). Clustering in subclade 3, section Leuce was confirmed to be of monophyletic origin and independent evolution. Its two subsections, namely Albidae and Trepidae, could be separated by chloroplast data but had frequent gene flow based on ITS data. Phylogeny analysis based on chloroplast data demonstrated once more that section Aigeiros was paraphyletic and further showed that the P. deltoides lineage is restricted in subclade 2 and that P. nigra lineage, located in subclade 3, originated from a hybrid of which an Albidae ancestor species was the material parent. Similarly, section Tacamahaca was found to be paraphyletic and had two lineages: a clade 1 lineage, such as P. cathayana, and a clade 2 lineage, such as P. simonii. Section Leucoides was paraphyletic and closely linked to section Tacamahaca. Their section boundaries were not conclusively delimitated by sequencing information.