De Novo Sequencing and Assembly Analysis of Transcriptome in Pinus bungeana Zucc. ex Endl

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Mass spectrometry-based forest tree metabolomics

Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.

De novo transcriptomic analysis and identification of EST-SSR markers in Stephanandra incisa

Stephanandra incisa is a wild-type shrub with beautiful leaves and white flowers and is commonly used as a garden decoration accessory. However, the limited availability of genomic data of S. incisa has restricted its breeding process. Here, we identified EST-SSR markers using de novo transcriptome sequencing. In this study, a transcriptome database containing 35,251 unigenes, having an average length of 985 bp, was obtained from S. incisa. From these unigene sequences, we identified 5,555 EST-SSRs, with a distribution density of one SSR per 1.60 kb. Dinucleotides (52.96%) were the most detected SSRs, followed by trinucleotides (34.64%). From the EST-SSR loci, we randomly selected 100 sites for designing primer and used the DNA of 60 samples to verify the polymorphism. The average value of the effective number of alleles (Ne), Shannon's information index (I), and expective heterozygosity (He) was 1.969, 0.728, and 0.434, respectively. The polymorphism information content (PIC) value was in the range of 0.108 to 0.669, averaging 0.406, which represented a middle polymorphism level. Cluster analysis of S. incisa were also performed based on the obtained EST-SSR data in our work. As shown by structure analysis, 60 individuals could be classified into two groups. Thus, the identification of these novel EST-SSR markers provided valuable sequence information for analyzing the population structure, genetic diversity, and genetic resource assessment of S. incisa and other related species.

Transcriptome Analysis of Elm (Ulmus pumila) Fruit to Identify Phytonutrients Associated Genes and Pathways

Plant fruit is an important source of natural active phytonutrients that are profitable for human health. Elm (Ulmus pumila) fruit is considered as natural plant food in China that is rich in nutrients. In the present study, high-throughput RNA sequencing was performed in U. pumila edible fruits and leaves and 11,386 unigenes were filtered as dysregulated genes in fruit samples, including 5231 up- and 6155 downregulated genes. Hundreds of pathways were predicted to participate in seed development and phytonutrient biosynthesis in U. pumila by GO, MapMan, and KEGG enrichment analysis, including “seed maturation”, “glycine, serine, and threonine metabolism” and “phenylpropanoid biosynthesis”. ABA-mediated glucose response-related ethylene-activated signaling pathway (e.g., ABI4) were supposed to associate with elm fruit development; unsaturated fatty acids pathway (e.g., ACX2 and SAD) were predicted to participate in determination of fatty acid composition in elm fruit; flavonoid and coumarins biosynthesis (e.g., CYP98A3 and CCoAOMT1) were demonstrated to correlate with the bioactivity of elm fruits in human cancer and inflammation resistance. To provide more information about fruit developmental status, the qRT-PCR analysis for key genes of “phenylpropanoid biosynthesis” and “alpha-Linolenic acid metabolism” were conducted in samples of young fruits, ripe fruit, old fruit, and leaves. Two biosynthetic pathways for unsaturated fatty acid and Jasmonic acid (JA) were deduced to be involved in fruit development in U. pumila and the phenylpropanoid glycoside, syringin, was speculated to accumulate in the early development stages of elm fruit. Our transcriptome data supports molecular clues for seed development and biologically active substances in elm fruits.

Physiological and Transcriptome Analysis of a Yellow-Green Leaf Mutant in Birch (Betula platyphylla × B. Pendula)

Chlorophyll (Chl)-deficient mutants are ideal materials for the study of Chl biosynthesis, chloroplast development, and photosynthesis. Although the genes encoding key enzymes related to Chl biosynthesis have been well-characterized in herbaceous plants, rice (Oryza sativa L.), Arabidopsis (Arabidopsis thaliana), and maize (Zea mays L.), yellow-green leaf mutants have not yet been fully studied in tree species. In this work, we explored the molecular mechanism of the leaf color formation in a yellow-green leaf mutant (yl). We investigated the differentially expressed genes (DEGs) between yl and control plants (wild type birch (WT) and BpCCR1 overexpression line 11, (C11)) by transcriptome sequencing. Approximately 1163 genes (874 down-regulated and 289 up-regulated) and 930 genes (755 down-regulated and 175 up-regulated) were found to be differentially expressed in yl compared with WT and C11, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis for DEGs revealed that photosynthesis antenna proteins represent the most significant enriched pathway. The expressions of photosynthesis antenna proteins are crucial to the leaf color formation in yl. We also found that Chl accumulate, leaf anatomical structure, photosynthesis, and growth were affected in yl. Taken together, our results not only provide the difference of phenomenal, physiological, and gene expression characteristics in leaves between yl mutant and control plants, but also provide a new insight into the mutation underlying the chlorotic leaf phenotype in birch.

De Novo Transcriptomic Analysis and Development of EST–SSRs for Styrax japonicus

Styrax japonicus sieb. et Zucc. is widely distributed in China with ornamental and medicinal values. However, the transcriptome of S. japonicus has not yet been reported. In this study, we carried out the first transcriptome analysis of S. japonicus and developed a set of expressed sequence tag–simple sequence repeats (EST–SSRs). We obtained 338,570,222 clean reads in total, of which the mean GC content was 41.58%. In total, 136,071 unigenes were obtained having an average length of 611 bp and 71,226 unigenes were favorably annotated in the database. In total, we identified 55,977 potential EST–SSRs from 38,611 unigenes, of which there was 1 SSR per 6.73 kb. The di-nucleotide repeats (40.40%) were the most identified SSRs. One set of 60 primer pairs was randomly selected, and the amplified products in S. japonicus were validated; 28 primer pairs successfully produced clear amplicons. A total of 21 (35%) polymorphic genic SSR markers were identified between two populations. In total, 15 alleles were detected and the average number was 6. The average of observed heterozygosity and expected heterozygosity was 0.614 and 0.552, respectively. The polymorphism information content (PIC) value fluctuated between 0.074 and 0.855, with a mean value of 0.504, which was also the middle level. This study provides useful information for diversity studies and resource assessments of S. japonicus.