Haplotype-resolved genome assembly and resequencing analysis provide insights into genome evolution and allelic imbalance in Pinus densiflora.

Haplotype-level allelic characterization facilitates research on the functional, evolutionary and breeding-related features of extremely large and complex plant genomes. We report a 21.7-Gb chromosome-level haplotype-resolved assembly in Pinus densiflora. We found genome rearrangements involving translocations and inversions between chromosomes 1 and 3 of Pinus species and a proliferation of specific long terminal repeat (LTR) retrotransposons (LTR-RTs) in P. densiflora. Evolutionary analyses illustrated that tandem and LTR-RT-mediated duplications led to an increment of transcription factor (TF) genes in P. densiflora. The haplotype sequence comparison showed allelic imbalances, including presence-absence variations of genes (PAV genes) and their functional contributions to flowering and abiotic stress-related traits in P. densiflora. Allele-aware resequencing analysis revealed PAV gene diversity across P. densiflora accessions. Our study provides insights into key mechanisms underlying the evolution of genome structure, LTR-RTs and TFs within the Pinus lineage as well as allelic imbalances and diversity across P. densiflora.

An Improved Chromosome-scale Genome Assembly and Population Genetics resource for Populus tremula.

Aspen (Populus tremula L.) is a keystone species and a model system for forest tree genomics. We present an updated resource comprising a chromosome-scale assembly, population genetics and genomics data. Using the resource, we explore the genetic basis of natural variation in leaf size and shape, traits with complex genetic architecture. We generated the genome assembly using long-read sequencing, optical and high-density genetic maps. We conducted whole-genome resequencing of the Umeå Aspen (UmAsp) collection. Using the assembly and re-sequencing data from the UmAsp, Swedish Aspen (SwAsp) and Scottish Aspen (ScotAsp) collections we performed genome-wide association analyses (GWAS) using Single Nucleotide Polymorphisms (SNPs) for 26 leaf physiognomy phenotypes. We conducted Assay of Transposase Accessible Chromatin sequencing (ATAC-Seq), identified genomic regions of accessible chromatin, and subset SNPs to these regions, improving the GWAS detection rate. We identified candidate long non-coding RNAs in leaf samples, quantified their expression in an updated co-expression network, and used this to explore the functions of candidate genes identified from the GWAS. A GWAS found SNP associations for seven traits. The associated SNPs were in or near genes annotated with developmental functions, which represent candidates for further study. Of particular interest was a ~177-kbp region harbouring associations with several leaf phenotypes in ScotAsp. We have incorporated the assembly, population genetics, genomics, and GWAS data into the PlantGenIE.org web resource, including updating existing genomics data to the new genome version, to enable easy exploration and visualisation. We provide all raw and processed data to facilitate reuse in future studies.

Multivariate Analysis among Marker Compounds, Environmental Factors, and Fruit Quality of Schisandra chinensis at Different Locations in South Korea.

This study aimed to investigate the correlation among the contents of marker compounds, growth characteristics, and environmental factors of Schisandra chinensis fruits across South Korea. The fruits were collected from 36 cultivation sites in 28 regions across the country. We investigated nine growth characteristics, twelve soil physicochemical properties, eight meteorological data, and three marker compounds in this study. We optimized and validated an optimized method for quantifying marker compounds using UPLC and performed correlation analysis among the contents of marker compounds, growth characteristics, and environmental factors. The UPLC-UV method for analyzing marker compounds was validated by measuring linearity, LOD, LOQ, precision, and accuracy. The marker compounds were negatively correlated with the fruit size and sugar contents, and growth characteristics were negatively correlated with some physicochemical properties of the soil. The results of this study can be used as basic data for the standard cultural practices and quality control of S. chinensis fruits.

Seasonal Developing Xylem Transcriptome Analysis of Pinus densiflora Unveils Novel Insights for Compression Wood Formation.

Wood is the most important renewable resource not only for numerous practical utilizations but also for mitigating the global climate crisis by sequestering atmospheric carbon dioxide. The compressed wood (CW) of gymnosperms, such as conifers, plays a pivotal role in determining the structure of the tree through the reorientation of stems displaced by environmental forces and is characterized by a high content of lignin. Despite extensive studies on many genes involved in wood formation, the molecular mechanisms underlying seasonal and, particularly, CW formation remain unclear. This study examined the seasonal dynamics of two wood tissue types in Pinus densiflora: CW and opposite wood (OW). RNA sequencing of developing xylem for two consecutive years revealed comprehensive transcriptome changes and unique differences in CW and OW across seasons. During growth periods, such as spring and summer, we identified 2255 transcripts with differential expression in CW, with an upregulation in lignin biosynthesis genes and significant downregulation in stress response genes. Notably, among the laccases critical for monolignol polymerization, PdeLAC17 was found to be specifically expressed in CW, suggesting its vital role in CW formation. PdeERF4, an ERF transcription factor preferentially expressed in CW, seems to regulate PdeLAC17 activity. This research provides an initial insight into the transcriptional regulation of seasonal CW development in P. densiflora, forming a foundation for future studies to enhance our comprehension of wood formation in gymnosperms.

Chromosome-level genome assembly of the Asian aspen Populus davidiana Dode.

The genome of Populus davidiana, a keystone aspen species, has been sequenced to improve our understanding of the evolutionary and functional genomics of the Populus genus. The Hi-C scaffolding genome assembly resulted in a 408.1 Mb genome with 19 pseudochromosomes. The BUSCO assessment revealed that 98.3% of the genome matched the embryophytes dataset. A total of 31,862 protein-coding sequences were predicted, of which 31,619 were functionally annotated. The assembled genome was composed of 44.9% transposable elements. These findings provide new knowledge about the characteristics of the P. davidiana genome and will facilitate comparative genomics and evolutionary research on the genus Populus.

Draft Genome Sequence for the Symbiotic Pine Mushroom Tricholoma matsutake.

We report the high-quality genome sequence of Tricholoma matsutake strain 2001, which was isolated from a mushroom fruiting body in South Korea. The genome has 80 contigs, a size of 162.6 Mb, and an N50 value of 5,103,859 bp and will provide insight into the symbiotic association between T. matsutake and Pinus densiflora.

Comparative functional analysis of PdeNAC2 and AtVND6 in the tracheary element formation.

Tracheary elements (i.e. vessel elements and tracheids) are highly specialized, non-living cells present in the water-conducting xylem tissue. In angiosperms, proteins in the VASCULAR-RELATED NAC-DOMAIN (VND) subgroup of the NAC (NAM, ATAF1,2, and CUC2) transcription factor family (e.g. AtVND6) are required for the differentiation of vessel elements through transcriptional regulation of genes responsible for secondary cell wall formation and programmed cell death. Gymnosperms, however, produce only tracheids, the mechanism of which remains elusive. Here, we report functional characteristics of PdeNAC2, a VND homolog in Pinus densiflora, as a key regulator of tracheid formation. Interestingly, our molecular genetic analyses show that PdeNAC2 can induce the formation of vessel element-like cells in angiosperm plants, demonstrated by transgenic overexpression of either native or NAC domain-swapped synthetic genes of PdeNAC2 and AtVND6 in both Arabidopsis and hybrid poplar. Subsequently, genome-wide identification of direct target (DT) genes of PdeNAC2 and AtVND6 revealed 138 and 174 genes as putative DTs, respectively, but only 17 genes were identified as common DTs. Further analyses have found that PdeNAC2 does not control some AtVND6-dependent vessel differentiation genes in angiosperm plants, such as AtVRLK1, LBD15/30 and pit-forming Rho-like GTPases from plant (ROP) signaling genes. Collectively, our results suggest that different target gene repertoires of PdeNAC2 and AtVND6 may contribute to the evolution of tracheary elements.

Identification, In Silico Characterization, and Differential Expression Profiles of Carotenoid, Xanthophyll, Apocarotenoid Biosynthetic Pathways Genes, and Analysis of Carotenoid and Xanthophyll Accumulation in Heracleum moellendorffii Hance.

Heracleum moellendorffii Hance is a non-woody forest plant widely used in China, Korea, and Japan because of its various therapeutic properties. However, the genetic details of the carotenoid pathway (CP), xanthophyll pathway (XP), and apocarotenoid pathway (AP) genes have not been studied. Thus, the CP, XP, and AP genes of H. moellendorffii were detected and analyzed. A total of fifteen genes were identified, of which eight, four, and three belonged to CP, XP, and AP, respectively. All identified genes possessed full open reading frames. Phylogenetic characterization of the identified gene sequences showed the highest similarity with other higher plants. Multiple alignments and 3D dimensional structures showed several diverse conserved motifs, such as the carotene-binding motif, dinucleotide-binding motif, and aspartate or glutamate residues. The results of real-time PCR showed that the CP, XP, and AP genes were highly expressed in leaves, followed by the stems and roots. In total, eight different individual carotenoids were identified using HPLC analysis. The highest individual and total carotenoid content were achieved in the leaves, followed by the stems and roots. This study will provide more information on the gene structure of the CP, XP, and AP genes, which may help to increase the accumulation of carotenoids in H. moellendorffii through genetic engineering. These results could be helpful for further molecular and functional studies of CP, XP, and AP genes.

Chromosome-level genome assembly of the fully mycoheterotrophic orchid Gastrodia elata.

Gastrodia elata, an obligate mycoheterotrophic orchid, requires complete carbon and mineral nutrient supplementation from mycorrhizal fungi during its entire life cycle. Although full mycoheterotrophy occurs most often in the Orchidaceae family, no chromosome-level reference genome from this group has been assembled to date. Here, we report a high-quality chromosome-level genome assembly of G. elata, using Illumina and PacBio sequencing methods with Hi-C technique. The assembled genome size was found to be 1045 Mb, with an N50 of 50.6 Mb and 488 scaffolds. A total of 935 complete (64.9%) matches to the 1440 embryophyte Benchmarking Universal Single-Copy Orthologs were identified in this genome assembly. Hi-C scaffolding of the assembled genome resulted in 18 pseudochromosomes, 1008 Mb in size and containing 96.5% of the scaffolds. A total of 18,844 protein-coding sequences (CDSs) were predicted in the G. elata genome, of which 15,619 CDSs (82.89%) were functionally annotated. In addition, 74.92% of the assembled genome was found to be composed of transposable elements. Phylogenetic analysis indicated a significant contraction of genes involved in various biosynthetic processes and cellular components and an expansion of genes for novel metabolic processes and mycorrhizal association. This result suggests an evolutionary adaptation of G. elata to a mycoheterotrophic lifestyle. In summary, the genomic resources generated in this study will provide a valuable reference genome for investigating the molecular mechanisms of G. elata biological functions. Furthermore, the complete G. elata genome will greatly improve our understanding of the genetics of Orchidaceae and its mycoheterotrophic evolution.

PtrMYB120 functions as a positive regulator of both anthocyanin and lignin biosynthetic pathway in a hybrid poplar.

Both anthocyanins and lignins are essential secondary metabolites in plant growth and development. Their biosynthesis is metabolically interconnected and diverges in the central metabolite 4-coumaroyl CoA of the phenylpropanoid pathway. Considerable progress has been made in understanding transcriptional regulation of genes involved in lignin and anthocyanin synthesis pathways, but the concerted regulation of these pathways is not yet fully understood. Here, we functionally characterized PtrMYB120, a R2R3-MYB transcription factor from Populus trichocarpa. Overexpression of PtrMYB120 in a hybrid poplar (i.e., 35S::PtrMYB120) was associated with increased anthocyanin (i.e., cyanidin 3-O-glucoside) accumulation and upregulation of anthocyanin biosynthetic genes. However, transgenic poplars with dominant suppression of PtrMYB120 function achieved by fusing the ERF-associated amphiphilic repression motif to PtrMYB120 (i.e., 35S::PtrMYB120-SRDX) had a dramatic decrease in not only anthocyanin but also Klason lignin content with downregulation of both anthocyanin and lignin biosynthetic genes. Indeed, 35S::PtrMYB120-SRDX poplars had irregularly shaped xylem vessels with reduced S-lignin content in stems, which was proportionally related to the level of the introduced PtrMYB120-SRDX gene. Furthermore, protoplast-based transcriptional activation assay using the PtrMYB120-GR system suggested that PtrMYB120 directly regulates genes involved in both anthocyanin and lignin biosynthesis, including chalcone synthase and ferulate-5 hydroxylase. Interestingly, the saccharification efficiency of line #6 of 35S::PtrMYB120-SRDX poplars, which had slightly reduced lignin content with a normal growth phenotype, was dramatically enhanced (>45%) by NaOH treatment. Taken together, our results suggest that PtrMYB120 functions as a positive regulator of both anthocyanin and lignin biosynthetic pathways and can be targeted to enhance saccharification efficiency in woody perennials.

Wood transcriptome analysis of Pinus densiflora identifies genes critical for secondary cell wall formation and NAC transcription factors involved in tracheid formation.

Although conifers have significant ecological and economic value, information on transcriptional regulation of wood formation in conifers is still limited. Here, to gain insight into secondary cell wall (SCW) biosynthesis and tracheid formation in conifers, we performed wood tissue-specific transcriptome analyses of Pinus densiflora (Korean red pine) using RNA sequencing. In addition, to obtain full-length transcriptome information, PacBio single molecule real-time iso-sequencing was carried out using RNAs from 28 tissues of P. densiflora. Subsequent comparative tissue-specific transcriptome analysis successfully pinpointed critical genes encoding key proteins involved in biosynthesis of the major secondary wall components (cellulose, galactoglucomannan, xylan and lignin). Furthermore, we predicted a total of 62 NAC (NAM, ATAF1/2 and CUC2) family transcription factor members and identified seven PdeNAC genes preferentially expressed in developing xylem tissues in P. densiflora. Protoplast-based transcriptional activation analysis found that four PdeNAC genes, homologous to VND, NST and SND/ANAC075, upregulated GUS activity driven by an SCW-specific cellulose synthase promoter. Consistently, transient overexpression of the four PdeNACs induced xylem vessel cell-like SCW deposition in both tobacco (Nicotiana benthamiana) and Arabidopsis leaves. Taken together, our data provide a foundation for further research to unravel transcriptional regulation of wood formation in conifers, especially SCW formation and tracheid differentiation.

Overexpression of a Poplar RING-H2 Zinc Finger, Ptxerico, Confers Enhanced Drought Tolerance via Reduced Water Loss and Ion Leakage in Populus.

Drought stress is one of the major environmental problems in the growth of crops and woody perennials, but it is getting worse due to the global climate crisis. XERICO, a RING (Really Interesting New Gene) zinc-finger E3 ubiquitin ligase, has been shown to be a positive regulator of drought tolerance in plants through the control of abscisic acid (ABA) homeostasis. We characterized a poplar (Populus trichocarpa) RING protein family and identified the closest homolog of XERICO called PtXERICO. Expression of PtXERICO is induced by both salt and drought stress, and by ABA treatment in poplars. Overexpression of PtXERICO in Arabidopsis confers salt and ABA hypersensitivity in young seedlings, and enhances drought tolerance by decreasing transpirational water loss. Consistently, transgenic hybrid poplars overexpressing PtXERICO demonstrate enhanced drought tolerance with reduced transpirational water loss and ion leakage. Subsequent upregulation of genes involved in the ABA homeostasis and drought response was confirmed in both transgenic Arabidopsis and poplars. Taken together, our results suggest that PtXERICO will serve as a focal point to improve drought tolerance of woody perennials.

The complete chloroplast genome of Korean Gastrodia elata Blume.

The completed chloroplast genome of Gastrodia elata Blume (G. elata) from Korea was determined in this study. The cpDNA is 35,230 bp in length and lacked the large and small single copy (LSC and SSC) regions, due to the lost inverted repeat (IR). The overall AT content is 73.30%, and the cpDNA contains 20 protein-coding genes, 5 tRNA genes, and 3 rRNA genes. Remarkably, the Korean G. elata cp genome was 74 bp smaller than that of the Chinese G. elata. It revealed substantial sequence variants 495 SNPs and 75 InDels between the two G. elata genomes.

Overexpression of Populus transcription factor PtrTALE12 increases axillary shoot development by regulating WUSCHEL expression.

The TALE (Three Amino acid Loop Extension) transcription factor family has been shown to control meristem formation and organogenesis in plants. To understand the functional roles of the TALE family in woody perennials, each of the TALE members of Populus trichocarpa was overexpressed in Arabidopsis as a proxy. Among them, the overexpression of PtrTALE12 (i.e., 35S::PtrTALE12) resulted in a dramatic increase of axillary shoot development with early flowering. Interestingly, expression of WUSCHEL (WUS), a central regulator of both apical and axillary meristem formation, was significantly increased in the 35S::PtrTALE12 Arabidopsis plants. Conversely, WUS expression was downregulated in 35S::PtrTALE12-SRDX (short transcriptional repressor domain) plants. Further analysis found that PtrTALE12, expressed preferentially in meristem tissues, directly regulates WUS expression in transient activation assays using Arabidopsis leaf protoplast. Yeast two-hybrid assays showed that PtrTALE12 interacts with SHOOT MERISTEMLESS (STM); however, the interaction does not affect the WUS expression. In addition, expression of both CIRCADIAN CLOCK ASSOCIATED1 (CCA1) and LATE ELONGATED HYPOCOTYL (LHY) genes was suppressed accordingly for early flowering 35S::PtrTALE12 Arabidopsis. Indeed, transgenic poplars overexpressing PtrTALE12 as well as Arabidopsis plants overexpressing AtBLH11, a close homolog of PtrTALE12, phenocopied the 35S::PtrTALE12 Arabidopsis (i.e., increased axillary shoot development). Taken together, our results suggest that PtrTALE12 functions as a positive regulator of axillary shoot formation in both Arabidopsis and poplar.

Evidence for widespread selection in shaping the genomic landscape during speciation of Populus.

Increasing our understanding of how evolutionary processes drive the genomic landscape of variation is fundamental to a better understanding of the genomic consequences of speciation. However, genome-wide patterns of within- and between- species variation have not been fully investigated in most forest tree species despite their global ecological and economic importance. Here, we use whole-genome resequencing data from four Populus species spanning the speciation continuum to reconstruct their demographic histories and investigate patterns of diversity and divergence within and between species. Using Populus trichocarpa as an outgroup species, we further infer the genealogical relationships and estimate the extent of ancient introgression among the three aspen species (Populus tremula, Populus davidiana and Populus tremuloides) throughout the genome. Our results show substantial variation in these patterns along the genomes with this variation being strongly predicted by local recombination rates and the density of functional elements. This implies that the interaction between recurrent selection and intrinsic genomic features has dramatically sculpted the genomic landscape over long periods of time. In addition, our findings provide evidence that, apart from background selection, recent positive selection and long-term balancing selection have also been crucial components in shaping patterns of genome-wide variation during the speciation process.

Identification of transcriptome-wide, nut weight-associated SNPs in Castanea crenata.

Nut weight is one of the most important traits that can affect a chestnut grower's returns. Due to the long juvenile phase of chestnut trees, the selection of desired characteristics at early developmental stages represents a major challenge for chestnut breeding. In this study, we identified single nucleotide polymorphisms (SNPs) in transcriptomic regions, which were significantly associated with nut weight in chestnuts (Castanea crenata), using a genome-wide association study (GWAS). RNA-sequencing (RNA-seq) data were generated from large and small nut-bearing trees, using an Illumina HiSeq. 2000 system, and 3,271,142 SNPs were identified. A total of 21 putative SNPs were significantly associated with chestnut weight (false discovery rate [FDR] < 10-5), based on further analyses. We also applied five machine learning (ML) algorithms, support vector machine (SVM), C5.0, k-nearest neighbour (k-NN), partial least squares (PLS), and random forest (RF), using the 21 SNPs to predict the nut weights of a second population. The average accuracy of the ML algorithms for the prediction of chestnut weights was greater than 68%. Taken together, we suggest that these SNPs have the potential to be used during marker-assisted selection to facilitate the breeding of large chestnut-bearing varieties.

Wood Transcriptome Profiling Identifies Critical Pathway Genes of Secondary Wall Biosynthesis and Novel Regulators for Vascular Cambium Development in Populus.

Wood, the most abundant biomass on Earth, is composed of secondary xylem differentiated from vascular cambium. However, the underlying molecular mechanisms of wood formation remain largely unclear. To gain insight into wood formation, we performed a series of wood-forming tissue-specific transcriptome analyses from a hybrid poplar (Populus alba × P. glandulosa, clone BH) using RNA-seq. Together with shoot apex and leaf tissue, cambium and xylem tissues were isolated from vertical stem segments representing a gradient of secondary growth developmental stages (i.e., immature, intermediate, and mature stem). In a comparative transcriptome analysis of the 'developing xylem' and 'leaf' tissue, we could identify critical players catalyzing each biosynthetic step of secondary wall components (e.g., cellulose, xylan, and lignin). Several candidate genes involved in the initiation of vascular cambium formation were found via a co-expression network analysis using abundantly expressed genes in the 'intermediate stem-derived cambium' tissue. We found that transgenic Arabidopsis plants overexpressing the PtrHAM4-1, a GRAS family transcription factor, resulted in a significant increase of vascular cambium development. This phenotype was successfully reproduced in the transgenic poplars overexpressing the PtrHAM4-1. Taken together, our results may serve as a springboard for further research to unravel the molecular mechanism of wood formation, one of the most important biological processes on this planet.

The complete chloroplast genome of Castanea crenata Sieb. & Zucc.

The complete chloroplast genome sequence of Castanea crenata was sequenced and assembled using PacBio Sequel data. The cpDNA was 160,787 bp in length, containing a pair of inverted repeats (IRs) of 25,654 bp each separated by a large and small single copy (LSC and SSC) regions of 90,645 bp and 18,836 bp, respectively. The cpDNA contained 102 genes, including 65 protein-coding genes, 8 ribosomal RNA genes and 37 transfer RNA genes. Phylogenetic analysis indicated that C. crenata was closest to C. pumila var. pumila, which is known as a typical variety of American chinquapin or dwarf chestnut.

Wood forming tissue-specific bicistronic expression of PdGA20ox1 and PtrMYB221 improves both the quality and quantity of woody biomass production in a hybrid poplar.

With the exponential growth of the human population and industrial developments, research on renewable energy resources is required to alleviate environmental and economic impacts caused by the consumption of fossil fuels. In this study, we present a synthetic biological application of a wood forming tissue-specific bicistronic gene expression system to improve both the quantity and quality of woody biomass to minimize undesirable growth penalties. Our transgenic poplars, designed to express both PdGA20ox1 (a GA20-oxidase from Pinus densiflora producing bioactive gibberellin, GA) and PtrMYB221 (a MYB transcription factor negatively regulating lignin biosynthesis) under the developing xylem (DX) tissue-specific promoter (i.e., DX15::PdGA20ox1-2A-PtrMYB221 poplar), resulted in a 2-fold increase in biomass quantity compared to wild-type (WT), without undesirable growth defects. A similar phenotype was observed in transgenic Arabidopsis plants harboring the same gene constructs. These phenotypic consequences were further verified in the field experiments. Importantly, our transgenic poplars exhibited an improved quality of biomass with reduced lignin content (~16.0 wt%) but increased holocellulose content (~6.6 wt%). Furthermore, the saccharification efficiency of our transgenic poplar increased significantly by up to 8%. Our results demonstrate that the controlled production of both GA and a secondary wall modifying regulator in the same spatio-temporal manner can be utilized as an efficient biotechnological tool for producing the desired multi-purpose woody biomass.

Functional and evolutionary genomic inferences in Populus through genome and population sequencing of American and European aspen.

The Populus genus is one of the major plant model systems, but genomic resources have thus far primarily been available for poplar species, and primarily Populus trichocarpa (Torr. & Gray), which was the first tree with a whole-genome assembly. To further advance evolutionary and functional genomic analyses in Populus, we produced genome assemblies and population genetics resources of two aspen species, Populus tremula L. and Populus tremuloides Michx. The two aspen species have distributions spanning the Northern Hemisphere, where they are keystone species supporting a wide variety of dependent communities and produce a diverse array of secondary metabolites. Our analyses show that the two aspens share a similar genome structure and a highly conserved gene content with P. trichocarpa but display substantially higher levels of heterozygosity. Based on population resequencing data, we observed widespread positive and negative selection acting on both coding and noncoding regions. Furthermore, patterns of genetic diversity and molecular evolution in aspen are influenced by a number of features, such as expression level, coexpression network connectivity, and regulatory variation. To maximize the community utility of these resources, we have integrated all presented data within the PopGenIE web resource (PopGenIE.org).