Atlas of tissue-specific and tissue-preferential gene expression in ecologically and economically significant conifer Pinus sylvestris

Comparative transcriptomics of seed nourishing tissues: uncovering conserved and divergent pathways in seed plants

The evolutionary and ecological success of spermatophytes is intrinsically linked to the seed habit, which provides a protective environment for the initial development of the new generation. This environment includes an ephemeral nourishing tissue that supports embryo growth. In gymnosperms this tissue originates from the asexual proliferation of the maternal megagametophyte, while in angiosperms it is a product of fertilization, and is called the endosperm. The emergence of these nourishing tissues is of profound evolutionary value, and they are also food staples for most of the world's population. Here, using Orthofinder to infer orthologue genes among newly generated and previously published datasets, we provide a comparative transcriptomic analysis of seed nourishing tissues from species of several angiosperm clades, including those of early diverging lineages, as well as of one gymnosperm. Our results show that, although the structure and composition of seed nourishing tissues has seen significant divergence along evolution, there are signatures that are conserved throughout the phylogeny. Conversely, we identified processes that are specific to species within the clades studied, and thus illustrate their functional divergence. With this, we aimed to provide a foundation for future studies on the evolutionary history of seed nourishing structures, as well as a resource for gene discovery in future functional studies.

Strong Purifying Selection in Haploid Tissue-Specific Genes of Scots Pine Supports the Masking Theory

The masking theory states that genes expressed in a haploid stage will be under more efficient selection. In contrast, selection will be less efficient in genes expressed in a diploid stage, where the fitness effects of recessive deleterious or beneficial mutations can be hidden from selection in heterozygous form. This difference can influence several evolutionary processes such as the maintenance of genetic variation, adaptation rate, and genetic load. Masking theory expectations have been confirmed in single-cell haploid and diploid organisms. However, in multicellular organisms, such as plants, the effects of haploid selection are not clear-cut. In plants, the great majority of studies indicating haploid selection have been carried out using male haploid tissues in angiosperms. Hence, evidence in these systems is confounded with the effects of sexual selection and intraspecific competition. Evidence from other plant groups is scarce, and results show no support for the masking theory. Here, we have used a gymnosperm Scots pine megagametophyte, a maternally derived seed haploid tissue, and four diploid tissues to test the strength of purifying selection on a set of genes with tissue-specific expression. By using targeted resequencing data of those genes, we obtained estimates of genetic diversity, the site frequency spectrum of 0-fold and 4-fold sites, and inferred the distribution of fitness effects of new mutations in haploid and diploid tissue-specific genes. Our results show that purifying selection is stronger for tissue-specific genes expressed in the haploid megagametophyte tissue and that this signal of strong selection is not an artifact driven by high expression levels.

SMRT and Illumina sequencing provide insights into mechanisms of lignin and terpenoids biosynthesis in Pinus massoniana Lamb

Wood and oleoresin are important industrial raw materials with high economic value; however, their molecular formation and biosynthesis mechanisms in different tissues of Pinus massoniana remain unexplored. Therefore, we used single-molecule real-time sequencing technology (SMRT) and Illumina RNA sequencing to establish a transcriptome dataset and explore the expression pattern of genes related to secondary metabolites involved in wood formation and oleoresin biosynthesis in six different P. massoniana tissues. In total, 63.58 Gb of polymerase reads were obtained, including 41,407 isoforms with an average length of 1822 bp. We identified 3939 and 8785 isoforms and 161 and 481 transcription factors with tissue expression specificity and in the reproductive and vegetative organs, respectively. Eighty isoforms were annotated as cellulose synthases and 224 isoforms involved in lignin biosynthesis were enriched. Additionally, we identified 217 isoforms involved in the terpenoid biosynthesis pathway, with needles having the most tissue-specific genes for terpenoid biosynthesis. Some isoforms related to lignin biosynthesis were highly expressed in the xylem, according to the results of transcriptome sequencing and real-time quantitative reverse-transcription polymerase chain reaction. Our research confirmed the advantages of SMRT sequencing and provided valuable information for the transcriptional annotation of P. massoniana, which will be beneficial for producing better raw wood and oleoresin materials.

Taming the massive genome of Scots pine with PiSy50k, a new genotyping array for conifer research

Pinus sylvestris (Scots pine) is the most widespread coniferous tree in the boreal forests of Eurasia, with major economic and ecological importance. However, its large and repetitive genome presents a challenge for conducting genome-wide analyses such as association studies, genetic mapping and genomic selection. We present a new 50K single-nucleotide polymorphism (SNP) genotyping array for Scots pine research, breeding and other applications. To select the SNP set, we first genotyped 480 Scots pine samples on a 407 540 SNP screening array and identified 47 712 high-quality SNPs for the final array (called 'PiSy50k'). Here, we provide details of the design and testing, as well as allele frequency estimates from the discovery panel, functional annotation, tissue-specific expression patterns and expression level information for the SNPs or corresponding genes, when available. We validated the performance of the PiSy50k array using samples from Finland and Scotland. Overall, 39 678 (83.2%) SNPs showed low error rates (mean = 0.9%). Relatedness estimates based on array genotypes were consistent with the expected pedigrees, and the level of Mendelian error was negligible. In addition, array genotypes successfully discriminate between Scots pine populations of Finnish and Scottish origins. The PiSy50k SNP array will be a valuable tool for a wide variety of future genetic studies and forestry applications.

Genomics and adaptation in forest ecosystems

Rapid human-induced environmental changes like climate warming represent a challenge for forest ecosystems. Due to their biological complexity and the long generation time of their keystone tree species, genetic adaptation in these ecosystems might not be fast enough to keep track with conditions changing at such a fast pace. The study of adaptation to environmental change and its genetic mechanisms is therefore key for ensuring a sustainable support and management of forests. The 4-day conference of the European Research Group EvolTree (https://www.evoltree.eu) on the topic of "Genomics and Adaptation in Forest Ecosystems" brought together over 130 scientists to present and discuss the latest developments and findings in forest evolutionary research. Genomic studies in forest trees have long been hampered by the lack of high-quality genomics resources and affordable genotyping methods. This has dramatically changed in the last few years; the conference impressively showed how such tools are now being applied to study past demography, adaptation and interactions with associated organisms. Moreover, genomic studies are now finally also entering the world of conservation and forest management, for example by measuring the value or cost of interspecific hybridization and introgression, assessing the vulnerability of species and populations to future change, or accurately delineating evolutionary significant units. The newly launched conference series of EvolTree will hopefully play a key role in the exchange and synthesis of such important investigations.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11295-022-01542-1.

Seed-derived defensins from Scots pine: structural and functional features

The recombinant PsDef5.1 defensin inhibits the growth of phytopathogenic fungi, Gram-positive and Gram-negative bacteria, and human pathogen Candida albicans. Expression of seed-derived Scots pine defensins is tissue-specific and developmentally regulated. Plant defensins are ubiquitous antimicrobial peptides that possess a broad spectrum of activities and multi-functionality. The genes for these antimicrobial proteins form a multigenic family in the plant genome and are expressed in every organ. Most of the known defensins have been isolated from seeds of various monocot and dicot species, but seed-derived defensins have not yet been characterized in gymnosperms. This study presents the isolation of two new 249 bp cDNA sequences from Scots pine seeds with 97.9% nucleotide homology named PsDef5.1 and PsDef5.2. Their deduced amino acid sequences have typical plant defensin features, including an endoplasmic reticulum signal sequence of 31 amino acids (aa), followed by a characteristic defensin domain of 51 aa. To elucidate the functional activity of new defensins, we expressed the mature form of PsDef5.1 in a prokaryotic system. The purified recombinant peptide exhibited activity against the phytopathogenic fungi and Gram-negative and Gram-positive bacteria with the IC₅₀ of 5-18 µM. Moreover, it inhibited the growth of the human pathogen Candida albicans with the IC₅₀ of 6.0 µM. Expression analysis showed that transcripts of PsDef5.1-2 genes were present in immature and mature pine seeds and different parts of seedlings at the early stage of germination. In addition, unlike the PsDef5.2, the PsDef5.1 gene was expressed in the reproductive organs. Our findings indicate that novel defensins are promising candidates for transgenic application and the development of new antimicrobial drugs.