An improved assembly of the loblolly pine mega-genome using long-read single-molecule sequencing

Physiological and molecular response to drought in somatic plants from Pinus radiata embryonal masses induced at high temperatures

Drought and heat are among the major abiotic stresses in forest trees and are directly related with the consequences of climatic change. Many responses to abiotic stresses in plants have been associated with plant memory but mechanisms underlying this phenomenon remain unclear. Somatic embryogenesis, which is considered one of the most important methods for large-scale vegetative propagation of plants, is also used for stress induction and study the mechanisms involved in adaptation to abiotic stress. Specifically, heat stress during initiation stage of somatic embryogenesis has shown to have an impact in differential expression of stress related genes in pines. Modifications caused by a previous stress could eventually influence the stress tolerance of somatic plants years later. In this study we analysed the response to drought in 2-year-old radiata pine somatic plants, derived from embryonal masses initiated at 60 °C, at physiological, transcriptomic and amino acid accumulation level. Our results showed a more pronounce response to drought in plants coming from 60 °C treatment, which presented lower values in several physiological parameters as well as higher proline and tyrosine levels. Additionally, the transcriptomic response to drought was stronger in heat primed plants compared to control plants, suggesting a memory acquired two years before.

Optimising Exome Captures in Species With Large Genomes Using Species-Specific Repetitive DNA Blocker

Large and highly repetitive genomes are common. However, research interests usually lie within the non-repetitive parts of the genome, as they are more likely functional, and can be used to answer questions related to adaptation, selection and evolutionary history. Exome capture is a cost-effective method for providing sequencing data from protein-coding parts of the genes. C0t-1 DNA blockers consist of repetitive DNA and are used in exome captures to prevent the hybridisation of repetitive DNA sequences to capture baits or bait-bound genomic DNA. Universal blockers target repetitive regions shared by many species, while species-specific c0t-1 DNA is prepared from the DNA of the studied species, thus perfectly matching the repetitive DNA contents of the species. So far, the use of species-specific c0t-1 DNA has been limited to a few model species. Here, we evaluated the performance of blocker treatments in exome captures of Pinus sylvestris, a widely distributed conifer species with a large (> 20 Gbp) and highly repetitive genome. We compared treatment with a commercial universal blocker to treatments with species-specific c0t-1 (30,000 and 60,000 ng). Species-specific c0t-1 captured more unique exons than the initial set of targets leading to increased SNP discovery and reduced sequencing of tandem repeats compared to the universal blocker. Based on our results, we recommend optimising exome captures using at least 60,000 ng of species-specific c0t-1 DNA. It is relatively easy and fast to prepare and can also be used with existing bait set designs.

gymnotoa-db: a database and application to optimize functional annotation in gymnosperms

Gymnosperms are a clade of non-flowering plants that include about 1000 living species. Due to their complex genomes and lack of genomic resources, functional annotation in genomics and transcriptomics on gymnosperms suffers from limitations. Here we present gymnotoa-db, which is a novel, publicly accessible relational database designed to facilitate functional annotation in gymnosperms. This database stores non-redundant records of gymnosperm proteins, encompassing taxonomic and functional information. The complementary software, gymnotoa-app, enables users to download gymnotoa-db and execute a comprehensive functional annotation pipeline for high-throughput sequencing-derived DNA or cDNA sequences. gymnotoa-app's user-friendly interface and efficient algorithms streamline the functional annotation process, making it an invaluable tool for researchers studying gymnosperms. We compared gymnotoa-app's performance against other annotation tools utilizing disparate reference databases. Our results demonstrate gymnotoa-app's superior ability to accurately annotate gymnosperm transcripts, recovering a greater number of transcripts and unique, non-redundant Gene Ontology terms. gymnotoa-db's distinctive features include comprehensive coverage with a non-redundant dataset of gymnosperm protein sequences, robust functional information that integrates data from multiple ontology systems, including GO, KEGG, EC, and MetaCYC, while keeping the taxonomic context, including Arabidopsis homologs. Database URL: https://blogs.upm.es/gymnotoa-db/2024/09/19/gymnotoa-app/.

Development of a 101.6K liquid-phased probe for GWAS and genomic selection in pine wilt disease-resistance breeding in Masson pine

Masson pine (Pinus massoniana Lamb.), indigenous to southern China, faces serious threats from pine wilt disease (PWD). Several natural genotypes have survived PWD outbreaks. Conducting genetic breeding with these resistant genotypes holds promise for enhancing resistance to PWD in Masson pine at its source. We conducted a genome-wide association study (GWAS) and genomic selection (GS) on 1013 Masson pine seedlings from 72 half-sib families to advance disease-resistance breeding. A set of efficient 101.6K liquid-phased probes was developed for single-nucleotide polymorphisms (SNPs) genotyping through target sequencing. PWD inoculation experiments were then performed to obtain phenotypic data for these populations. Our analysis reveals that the targeted sequencing data successfully divided the experimental population into three subpopulations consistent with the provenance, verifying the reliability of the liquid-phased probe. A total of 548 SNPs were considerably associated with disease-resistance traits using four GWAS algorithms. Among them, 283 were located on or linked to 169 genes, including common plant disease resistance-related protein families such as NBS-LRR and AP2/ERF. The DNNGP (deep neural network-based method for genomic prediction) model demonstrated superior performance in GS, achieving a maximum predictive accuracy of 0.71. The accuracy of disease resistance predictions reached 90% for the top 20% of the testing population ordered by resistance genomic estimated breeding value. This study establishes a foundational framework for advancing research on disease-resistant genes in P. massoniana and offers preliminary evidence supporting the feasibility of utilizing GS for the early identification of disease-resistant individuals.

Epigenetic memory of temperature sensed during somatic embryo maturation in 2-yr-old maritime pine trees

Embryogenesis is a brief but potentially critical phase in the life cycle of a tree for adaptive phenotypic plasticity. Using somatic embryogenesis in maritime pine (Pinus pinaster Ait.), we found that temperature during the maturation phase affects embryo development and postembryonic tree growth for up to 3 yr. We examined whether this somatic stress memory could stem from temperature- and/or development-induced changes in DNA methylation. For this, we developed a 200 mb custom sequence capture bisulfite analysis of genes and promoters to identify differentially methylated cytosines (DMCs) between temperature treatments (18, 23, and 28 °C) and developmental stages (immature and cotyledonary embryos, shoot apical meristem of 2-yr-old plants) and investigate if these differences can be mitotically transmitted from embryonic to postembryonic development (epigenetic memory). We revealed a high prevalence of temperature-induced DMCs in genes (8% to 14%) compared to promoters (<1%) in all 3 cytosine contexts. Developmental DMCs showed a comparable pattern but only in the CG context and with a strong trend toward hypomethylation, particularly in the promoters. A high percentage of DMCs induced by developmental transitions were found memorized in genes (up to 45%-50%) and promoters (up to 90%). By contrast, temperature-induced memory was lower and confined to genes after both embryonic (up to 14%) and postembryonic development (up to 8%). Using stringent criteria, we identified 10 genes involved in defense responses and adaptation, embryo development, and chromatin regulation that are candidates for the establishment of a persistent epigenetic memory of temperature sensed during embryo maturation in maritime pine. Here, we provide evidence that DNA methylation marks established during the embryonic phase are transmitted to the postembryonic plant development phase.

Beneficial mutualistic fungus Suillus luteus provided excellent buffering insurance in Scots pine defense responses under pathogen challenge at transcriptome level

BACKGROUND

Mutualistic mycorrhiza fungi that live in symbiosis with plants facilitates nutrient and water acquisition, improving tree growth and performance. In this study, we evaluated the potential of mutualistic fungal inoculation to improve the growth and disease resistance of Scots pine (Pinus sylvestris L.) against the forest pathogen Heterobasidion annosum.

RESULTS

In co-inoculation experiment, Scots pine seedlings were pre-inoculated with mutualistic beneficial fungus (Suillus luteus) prior to H. annosum infection. The result revealed that inoculation with beneficial fungus promoted plant root growth. Transcriptome analyses revealed that co-inoculated plants and plants inoculated with beneficial fungus shared some similarities in defense gene responses. However, pathogen infection alone had unique sets of genes encoding pathogenesis-related (PR) proteins, phenylpropanoid pathway/lignin biosynthesis, flavonoid biosynthesis, chalcone/stilbene biosynthesis, ethylene signaling pathway, JA signaling pathway, cell remodeling and growth, transporters, and fungal recognition. On the other hand, beneficial fungus inoculation repressed the expression of PR proteins, and other defense-related genes such as laccases, chalcone/stilbene synthases, terpene synthases, cytochrome P450s. The co-inoculated plants did not equally enhance the induction of PR genes, chalcone/stilbene biosynthesis, however genes related to cell wall growth, water and nutrient transporters, phenylpropanoid/lignin biosynthesis/flavonoid biosynthesis, and hormone signaling were induced.

CONCLUSION

S. luteus promoted mutualistic interaction by suppressing plant defense responses. Pre-inoculation of Scots pine seedlings with beneficial fungus S. luteus prior to pathogen challenge promoted primary root growth, as well as had a balancing buffering role in plant defense responses and cell growth at transcriptome level.

Unveiling Key Genes and Unique Transcription Factors Involved in Secondary Cell Wall Formation in Pinus taeda

Pinus taeda is a key timber species, and extensive research has been conducted on its wood formation. However, a comprehensive investigation into the biosynthetic pathways of lignin, cellulose, and hemicellulose in P. taeda is lacking, resulting in an incomplete understanding of secondary cell wall (SCW) formation in this species. In this study, we systematically analyzed transcriptomic data from previously published sources and constructed detailed pathways for lignin, cellulose, and hemicellulose biosynthesis. We identified 188 lignin-related genes and 78 genes associated with cellulose and hemicellulose biosynthesis. An RT-qPCR highlighted 15 key lignin biosynthesis genes and 13 crucial genes for cellulose and hemicellulose biosynthesis. A STEM analysis showed that most essential enzyme-coding genes clustered into Profile 14, suggesting their significant role in SCW formation. Additionally, we identified seven NAC and six MYB transcription factors (TFs) from atypical evolutionary clades, with distinct expression patterns from those of the previously characterized NAC and MYB genes, indicating potentially unique functions in SCW formation. This research provides the first comprehensive overview of lignin, cellulose, and hemicellulose biosynthetic genes in P. taeda and underscores the importance of non-canonical NAC and MYB TFs, laying a genetic foundation for future studies on SCW regulatory mechanisms.

When less is more: sketching with minimizers in genomics

The exponential increase in sequencing data calls for conceptual and computational advances to extract useful biological insights. One such advance, minimizers, allows for reducing the quantity of data handled while maintaining some of its key properties. We provide a basic introduction to minimizers, cover recent methodological developments, and review the diverse applications of minimizers to analyze genomic data, including de novo genome assembly, metagenomics, read alignment, read correction, and pangenomes. We also touch on alternative data sketching techniques including universal hitting sets, syncmers, or strobemers. Minimizers and their alternatives have rapidly become indispensable tools for handling vast amounts of data.

Experimental Induction of Extreme Indented Growth Rings (Hazel Wood) in Pinus halepensis Miller by Wide and Long Parallel Bark and Vascular Cambium Woundings

Indented growth rings were found long ago to be experimentally induced in Pinus halepensis Miller by thin parallel axial scratching of the bark up to the vascular cambium with a sharp blade. Here, we show that when the bark and vascular cambium of P. halepensis are wounded by wide and long parallel axial wounds ("windows") rather than by thin scratches, the induced indented growth rings become dramatically more indented. All ten trees that were wounded by long parallel "windows" responded with very strong growth (especially in the first two years) that resulted in the formation of very conspicuous, extremely indented growth rings in the wood formed in between the long and wide woundings. This is true for both the trunks that were wounded all around their circumference and those that were wounded only in part of their circumference. We also suggest further lines of research.

Genome-wide identification and characterization of the AP2/ERF gene family in loblolly pine (Pinus taeda L.)

The loblolly pine (Pinus taeda L.) is one of the most profitable forest species worldwide owing to its quick growth, high wood yields, and strong adaptability. The AP2/ERF gene family plays a widespread role in the physiological processes of plant defense responses and the biosynthesis of metabolites. Nevertheless, there are no reports on this gene family in loblolly pine (P. taeda). In this study, a total of 303 members of the AP2/ERF gene family were identified. Through multiple sequence alignment and phylogenetic analysis, they were classified into four subfamilies, including AP2 (34), RAV (17), ERF (251), and Soloist (1). An analysis of the conservation domains, conserved motifs, and gene structure revealed that every PtAP2/ERF transcription factor (TF) had at least one AP2 domain. While evolutionary conservation was displayed within the same subfamilies, the distribution of conserved domains, conserved motifs, and gene architectures varied between subfamilies. Cis-element analysis revealed abundant light-responsive elements, phytohormone-responsive elements, and stress-responsive elements in the promoter of the PtAP2/ERF genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of potential target genes showed that the AP2/ERF gene family might play a critical role in plant growth and development, the response to environmental stresses, and metabolite biosynthesis. Utilizing quantitative real-time PCR (qRT-PCR), we examined the expression patterns of 10 randomly selected genes from Group IX after 6 h of treatments with mechanical injury, ethephon (Eth), and methyl jasmonate (MeJA). The AP2/ERF gene family in the loblolly pine was systematically analyzed for the first time in this study, offering a theoretical basis for exploring the functions and applications of AP2/ERF genes.

A genome sequence for the threatened whitebark pine

Whitebark pine (WBP, Pinus albicaulis) is a white pine of subalpine regions in the Western contiguous United States and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola) and additional threats from mountain pine beetle (Dendroctonus ponderosae), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short reads of haploid megagametophyte tissue and Oxford Nanopore long reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gb of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gb). Approximately 87.2% (24.0 Gb) of total sequence was placed on the 12 WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the 3 subclasses of NLRs. Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo-assembled transcriptomes.

Genetic transformation in conifers: current status and future prospects

Genetic transformation has been a cornerstone in plant molecular biology research and molecular design breeding, facilitating innovative approaches for the genetic improvement of trees with long breeding cycles. Despite the profound ecological and economic significance of conifers in global forestry, the application of genetic transformation in this group has been fraught with challenges. Nevertheless, genetic transformation has achieved notable advances in certain conifer species, while these advances are confined to specific genotypes, they offer valuable insights for technological breakthroughs in other species. This review offers an in-depth examination of the progress achieved in the genetic transformation of conifers. This discussion encompasses various factors, including expression vector construction, gene-delivery methods, and regeneration systems. Additionally, the hurdles encountered in the pursuit of a universal model for conifer transformation are discussed, along with the proposal of potential strategies for future developments. This comprehensive overview seeks to stimulate further research and innovation in this crucial field of forest biotechnology.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

Identification and high-throughput genotyping of single nucleotide polymorphism markers in a non-model conifer (Abies nordmanniana (Steven) Spach)

Single nucleotide polymorphism (SNP) markers are powerful tools for investigating population structures, linkage analysis, and genome-wide association studies, as well as for breeding and population management. The availability of SNP markers has been limited to the most commercially important timber species, primarily due to the cost of genome sequencing required for SNP discovery. In this study, a combination of reference-based and reference-free approaches were used to identify SNPs in Nordmann fir (Abies nordmanniana), a species previously lacking genomic sequence information. Using a combination of a genome assembly of the closely related Silver fir (Abies alba) species and a de novo assembly of low-copy regions of the Nordmann fir genome, we identified a high density of reliable SNPs. Reference-based approaches identified two million SNPs in common between the Silver fir genome and low-copy regions of Nordmann fir. A combination of one reference-free and two reference-based approaches identified 250 shared SNPs. A subset of 200 SNPs were used to genotype 342 individuals and thereby tested and validated in the context of identity analysis and/or clone identification. The tested SNPs successfully identified all ramets per clone and five mislabeled individuals via identity and genomic relatedness analysis. The identified SNPs will be used in ad hoc breeding of Nordmann fir in Denmark.

Identifying genetic variation associated with environmental gradients and drought-tolerance phenotypes in ponderosa pine

As climate changes, understanding the genetic basis of local adaptation in plants becomes an ever more pressing issue. Combining genotype-environment association (GEA) with genotype-phenotype association (GPA) analysis has an exciting potential to uncover the genetic basis of environmental responses. We use these approaches to identify genetic variants linked to local adaptation to drought in Pinus ponderosa. Over 4 million Single Nucleotide Polymorphisms (SNPs) were identified using 223 individuals from across the Sierra Nevada of California. 927,740 (22.3%) SNPs were retained after filtering for proximity to genes and used in our association analyses. We found 1374 associated with five major climate variables, with the largest number (1151) associated with April 1st snowpack. We also conducted a greenhouse study with various drought-tolerance traits measured in first-year seedlings of a subset of the genotyped trees grown in the greenhouse. 796 SNPs were associated with control-condition trait values, while 1149 were associated with responsiveness of these traits to drought. While no individual SNPs were associated with both the environmental variables and the measured traits, several annotated genes were associated with both, particularly those involved in cell wall formation, biotic and abiotic stress responses, and ubiquitination. However, the functions of many of the associated genes have not yet been determined due to the lack of gene annotation information for conifers. Future studies are needed to assess the developmental roles and ecological significance of these unknown genes.

The genome of the Wollemi pine, a critically endangered "living fossil" unchanged since the Cretaceous, reveals extensive ancient transposon activity

We present the genome of the living fossil, Wollemia nobilis, a southern hemisphere conifer morphologically unchanged since the Cretaceous. Presumed extinct until rediscovery in 1994, the Wollemi pine is critically endangered with less than 60 wild adults threatened by intensifying bushfires in the Blue Mountains of Australia. The 12 Gb genome is among the most contiguous large plant genomes assembled, with extremely low heterozygosity and unusual abundance of DNA transposons. Reduced representation and genome re-sequencing of individuals confirms a relictual population since the last major glacial/drying period in Australia, 120 ky BP. Small RNA and methylome sequencing reveal conservation of ancient silencing mechanisms despite the presence of thousands of active and abundant transposons, including some transferred horizontally to conifers from arthropods in the Jurassic. A retrotransposon burst 8-6 my BP coincided with population decline, possibly as an adaptation enhancing epigenetic diversity. Wollemia, like other conifers, is susceptible to Phytophthora, and a suite of defense genes, similar to those in loblolly pine, are targeted for silencing by sRNAs in leaves. The genome provides insight into the earliest seed plants, while enabling conservation efforts.

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.

A long-read and short-read transcriptomics approach provides the first high-quality reference transcriptome and genome annotation for Pseudotsuga menziesii (Douglas-fir)

Douglas-fir (Pseudotsuga menziesii) is native to western North America. It grows in a wide range of environmental conditions and is an important timber tree. Although there are several studies on the gene expression responses of Douglas-fir to abiotic cues, the absence of high-quality transcriptome and genome data is a barrier to further investigation. Like for most conifers, the available transcriptome and genome reference dataset for Douglas-fir remains fragmented and requires refinement. We aimed to generate a highly accurate, and complete reference transcriptome and genome annotation. We deep-sequenced the transcriptome of Douglas-fir needles from seedlings that were grown under nonstress control conditions or a combination of heat and drought stress conditions using long-read (LR) and short-read (SR) sequencing platforms. We used 2 computational approaches, namely de novo and genome-guided LR transcriptome assembly. Using the LR de novo assembly, we identified 1.3X more high-quality transcripts, 1.85X more "complete" genes, and 2.7X more functionally annotated genes compared to the genome-guided assembly approach. We predicted 666 long noncoding RNAs and 12,778 unique protein-coding transcripts including 2,016 putative transcription factors. We leveraged the LR de novo assembled transcriptome with paired-end SR and a published single-end SR transcriptome to generate an improved genome annotation. This was conducted with BRAKER2 and refined based on functional annotation, repetitive content, and transcriptome alignment. This high-quality genome annotation has 51,419 unique gene models derived from 322,631 initial predictions. Overall, our informatics approach provides a new reference Douglas-fir transcriptome assembly and genome annotation with considerably improved completeness and functional annotation.

Targeted and untargeted metabolomics reveals deep analysis of drought stress responses in needles and roots of Pinus taeda seedlings

Drought stress is one of major environmental stresses affecting plant growth and yield. Although Pinus taeda trees are planted in rainy southern China, local drought sometime occurs and can last several months, further affecting their growth and resin production. In this study, P. taeda seedlings were treated with long-term drought (42 d), and then targeted and untargeted metabolomics analysis were carried out to evaluate drought tolerance of P. taeda. Targeted metabolomics analysis showed that levels of some sugars, phytohormones, and amino acids significantly increased in the roots and needles of water-stressed (WS) P. taeda seedlings, compared with well-watered (WW) pine seedlings. These metabolites included sucrose in pine roots, the phytohormones abscisic acid and sacylic acid in pine needles, the phytohormone gibberellin (GA4) and the two amino acids, glycine and asparagine, in WS pine roots. Compared with WW pine seedlings, the neurotransmitter acetylcholine significantly increased in needles of WS pine seedlings, but significantly reduced in their roots. The neurotransmitters L-glutamine and hydroxytyramine significantly increased in roots and needles of WS pine seedlings, respectively, compared with WW pine seedlings, but the neurotransmitter noradrenaline significantly reduced in needles of WS pine seedlings. Levels of some unsaturated fatty acids significantly reduced in roots or needles of WS pine seedlings, compared with WW pine seedlings, such as linoleic acid, oleic acid, myristelaidic acid, myristoleic acid in WS pine roots, and palmitelaidic acid, erucic acid, and alpha-linolenic acid in WS pine needles. However, three saturated fatty acids significantly increased in WS pine seedlings, i.e., dodecanoic acid in WS pine needles, tricosanoic acid and heptadecanoic acid in WS pine roots. Untargeted metabolomics analysis showed that levels of some metabolites increased in WS pine seedlings, especially sugars, long-chain lipids, flavonoids, and terpenoids. A few of specific metabolites increased greatly, such as androsin, piceatanol, and panaxatriol in roots and needles of WS pine seedlings. Comparing with WW pine seedlings, it was found that the most enriched pathways in WS pine needles included flavone and flavonol biosynthesis, ABC transporters, diterpenoid biosynthesis, plant hormone signal transduction, and flavonoid biosynthesis; in WS pine roots, the most enriched pathways included tryptophan metabolism, caffeine metabolism, sesquiterpenoid and triterpenoid biosynthesis, plant hormone signal transduction, biosynthesis of phenylalanine, tyrosine, and tryptophan. Under long-term drought stress, P. taeda seedlings showed their own metabolomics characteristics, and some new metabolites and biosynthesis pathways were found, providing a guideline for breeding drought-tolerant cultivars of P. taeda.

Genetic diversity and structure of the 4th cycle breeding population of Chinese fir (Cunninghamia lanceolata (lamb.) hook)

Studying population genetic structure and diversity is crucial for the marker-assisted selection and breeding of coniferous tree species. In this study, using RAD-seq technology, we developed 343,644 high-quality single nucleotide polymorphism (SNP) markers to resolve the genetic diversity and population genetic structure of 233 Chinese fir selected individuals from the 4^th cycle breeding program, representing different breeding generations and provenances. The genetic diversity of the 4^th cycle breeding population was high with nucleotide diversity (P_i ) of 0.003, and H_o and H_e of 0.215 and 0.233, respectively, indicating that the breeding population has a broad genetic base. The genetic differentiation level between the different breeding generations and different provenances was low (F_st < 0.05), with population structure analysis results dividing the 233 individuals into four subgroups. Each subgroup has a mixed branch with interpenetration and weak population structure, which might be related to breeding rather than provenance, with aggregation from the same source only being in the local branches. Our results provide a reference for further research on the marker-assisted selective breeding of Chinese fir and other coniferous trees.

A draft Diabrotica virgifera virgifera genome: insights into control and host plant adaption by a major maize pest insect

BACKGROUND

Adaptations by arthropod pests to host plant defenses of crops determine their impacts on agricultural production. The larval host range of western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is restricted to maize and a few grasses. Resistance of D. v. virgifera to crop rotation practices and multiple insecticides contributes to its status as the most damaging pest of cultivated maize in North America and Europe. The extent to which adaptations by this pest contributes to host plant specialization remains unknown.

RESULTS

A 2.42 Gb draft D. v. virgifera genome, Dvir_v2.0, was assembled from short shotgun reads and scaffolded using long-insert mate-pair, transcriptome and linked read data. K-mer analysis predicted a repeat content of ≥ 61.5%. Ortholog assignments for Dvir_2.0 RefSeq models predict a greater number of species-specific gene duplications, including expansions in ATP binding cassette transporter and chemosensory gene families, than in other Coleoptera. A majority of annotated D. v. virgifera cytochrome P450s belong to CYP4, 6, and 9 clades. A total of 5,404 transcripts were differentially-expressed between D. v. virgifera larvae fed maize roots compared to alternative host (Miscanthus), a marginal host (Panicum virgatum), a poor host (Sorghum bicolor) and starvation treatments; Among differentially-expressed transcripts, 1,908 were shared across treatments and the least number were between Miscanthus compared to maize. Differentially-expressed transcripts were enriched for putative spliceosome, proteosome, and intracellular transport functions. General stress pathway functions were unique and enriched among up-regulated transcripts in marginal host, poor host, and starvation responses compared to responses on primary (maize) and alternate hosts.

CONCLUSIONS

Manual annotation of D. v. virgifera Dvir_2.0 RefSeq models predicted expansion of paralogs with gene families putatively involved in insecticide resistance and chemosensory perception. Our study also suggests that adaptations of D. v. virgifera larvae to feeding on an alternate host plant invoke fewer transcriptional changes compared to marginal or poor hosts. The shared up-regulation of stress response pathways between marginal host and poor host, and starvation treatments may reflect nutrient deprivation. This study provides insight into transcriptomic responses of larval feeding on different host plants and resources for genomic research on this economically significant pest of maize.

Proteomics research in forest trees: A 2012-2022 update

This review is a compilation of proteomic studies on forest tree species published in the last decade (2012-2022), mostly focused on the most investigated species, including Eucalyptus, Pinus, and Quercus. Improvements in equipment, platforms, and methods in addition to the increasing availability of genomic data have favored the biological knowledge of these species at the molecular, organismal, and community levels. Integration of proteomics with physiological, biochemical and other large-scale omics in the direction of the Systems Biology, will provide a comprehensive understanding of different biological processes, from growth and development to responses to biotic and abiotic stresses. As main issue we envisage that proteomics in long-living plants will thrive light on the plant responses and resilience to global climate change, contributing to climate mitigation strategies and molecular breeding programs. Proteomics not only will provide a molecular knowledge of the mechanisms of resilience to either biotic or abiotic stresses, but also will allow the identification on key gene products and its interaction. Proteomics research has also a translational character being applied to the characterization of the variability and biodiversity, as well as to wood and non-wood derived products, traceability, allergen and bioactive peptides identification, among others. Even thought, the full potential of proteomics is far from being fully exploited in forest tree research, with PTMs and interactomics being reserved to plant model systems. The most outstanding achievements in forest tree proteomics in the last decade as well as prospects are discussed.

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

The long road to bloom in conifers

More than 600 species of conifers (phylum Pinophyta) serve as the backbone of the Earth's terrestrial plant community and play key roles in global carbon and water cycles. Although coniferous forests account for a large fraction of global wood production, their productivity relies largely on the use of genetically improved seeds. However, acquisition of such seeds requires recurrent selection and testing of genetically superior parent trees, eventually followed by the establishment of a seed orchard to produce the improved seeds. The breeding cycle for obtaining the next generation of genetically improved seeds can be significantly lengthened when a target species has a long juvenile period. Therefore, development of methods for diminishing the juvenile phase is a cost-effective strategy for shortening breeding cycle in conifers. The molecular regulatory programs associated with the reproductive transition and annual reproductive cycle of conifers are modulated by environmental cues and endogenous developmental signals. Mounting evidence indicates that an increase in global average temperature seriously threatens plant productivity, but how conifers respond to the ever-changing natural environment has yet to be fully characterized. With the breakthrough of assembling and annotating the giant genome of conifers, identification of key components in the regulatory cascades that control the vegetative to reproductive transition is imminent. However, comparison of the signaling pathways that control the reproductive transition in conifers and the floral transition in Arabidopsis has revealed many differences. Therefore, a more complete understanding of the underlying regulatory mechanisms that control the conifer reproductive transition is of paramount importance. Here, we review our current understanding of the molecular basis for reproductive regulation, highlight recent discoveries, and review new approaches for molecular research on conifers.

Distinct composition and amplification dynamics of transposable elements in sacred lotus (Nelumbo nucifera Gaertn.)

Sacred lotus (Nelumbo nucifera Gaertn.) is a basal eudicot plant with a unique lifestyle, physiological features, and evolutionary characteristics. Here we report the unique profile of transposable elements (TEs) in the genome, using a manually curated repeat library. TEs account for 59% of the genome, and hAT (Ac/Ds) elements alone represent 8%, more than in any other known plant genome. About 18% of the lotus genome is comprised of Copia LTR retrotransposons, and over 25% of them are associated with non-canonical termini (non-TGCA). Such high abundance of non-canonical LTR retrotransposons has not been reported for any other organism. TEs are very abundant in genic regions, with retrotransposons enriched in introns and DNA transposons primarily in flanking regions of genes. The recent insertion of TEs in introns has led to significant intron size expansion, with a total of 200 Mb in the 28 455 genes. This is accompanied by declining TE activity in intergenic regions, suggesting distinct control efficacy of TE amplification in different genomic compartments. Despite the prevalence of TEs in genic regions, some genes are associated with fewer TEs, such as those involved in fruit ripening and stress responses. Other genes are enriched with TEs, and genes in epigenetic pathways are the most associated with TEs in introns, indicating a dynamic interaction between TEs and the host surveillance machinery. The dramatic differential abundance of TEs with genes involved in different biological processes as well as the variation of target preference of different TEs suggests the composition and activity of TEs influence the path of evolution.

The genome of homosporous maidenhair fern sheds light on the euphyllophyte evolution and defences

Euphyllophytes encompass almost all extant plants, including two sister clades, ferns and seed plants. Decoding genomes of ferns is the key to deep insight into the origin of euphyllophytes and the evolution of seed plants. Here we report a chromosome-level genome assembly of Adiantum capillus-veneris L., a model homosporous fern. This fern genome comprises 30 pseudochromosomes with a size of 4.8-gigabase and a contig N50 length of 16.22 Mb. Gene co-expression network analysis uncovered that homospore development in ferns has relatively high genetic similarities with that of the pollen in seed plants. Analysing fern defence response expands understanding of evolution and diversity in endogenous bioactive jasmonates in plants. Moreover, comparing fern genomes with those of other land plants reveals changes in gene families important for the evolutionary novelties within the euphyllophyte clade. These results lay a foundation for studies on fern genome evolution and function, as well as the origin and evolution of euphyllophytes.

Conifer Biotechnology: An Overview

The peculiar characteristics of conifers determine the difficulty of their study and their great importance from various points of view. However, their study faces numerous important scientific, methodological, cultural, economic, social, and legal challenges. This paper presents an approach to several of those challenges and proposes a multidisciplinary scientific perspective that leads to a holistic understanding of conifers from the perspective of the latest technical, computer, and scientific advances. This review highlights the deep connection that all scientific contributions to conifers can have in each other as fully interrelated communicating vessels.

The Larix kaempferi genome reveals new insights into wood properties

Here, through single-molecule real-time sequencing, we present a high-quality genome sequence of the Japanese larch (Larix kaempferi), a conifer species with great value for wood production and ecological afforestation. The assembled genome is 10.97 Gb in size, harboring 45,828 protein-coding genes. Of the genome, 66.8% consists of repeat sequences, of which long terminal repeat retrotransposons are dominant and make up 69.86%. We find that tandem duplications have been responsible for the expansion of genes involved in transcriptional regulation and stress responses, unveiling their crucial roles in adaptive evolution. Population transcriptome analysis reveals that lignin content in L. kaempferi is mainly determined by the process of monolignol polymerization. The expression values of six genes (LkCOMT7, LkCOMT8, LkLAC23, LkLAC102, LkPRX148, and LkPRX166) have significantly positive correlations with lignin content. These results indicated that the increased expression of these six genes might be responsible for the high lignin content of the larches' wood. Overall, this study provides new genome resources for investigating the evolution and biological function of conifer trees, and also offers new insights into wood properties of larches.

Evaluating the accuracy of variant calling methods using the frequency of parent-offspring genotype mismatch

The use of next‐generation sequencing (NGS) data sets has increased dramatically over the last decade, but there have been few systematic analyses quantifying the accuracy of the commonly used variant caller programs. Here we used a familial design consisting of diploid tissue from a single lodgepole pine (Pinus contorta) parent and the maternally derived haploid tissue from 106 full‐sibling offspring, where mismatches could only arise due to mutation or bioinformatic error. Given the rarity of mutation, we used the rate of mismatches between parent and offspring genotype calls to infer the single nucleotide polymorphism (SNP) genotyping error rates of FreeBayes, HaplotypeCaller, SAMtools, UnifiedGenotyper, and VarScan. With baseline filtering HaplotypeCaller and UnifiedGenotyper yielded more SNPs and higher error rates by one to two orders of magnitude, whereas FreeBayes, SAMtools and VarScan yielded lower numbers of SNPs and more modest error rates. To facilitate comparison between variant callers we standardized each SNP set to the same number of SNPs using additional filtering, where UnifiedGenotyper consistently produced the smallest proportion of genotype errors, followed by HaplotypeCaller, VarScan, SAMtools, and FreeBayes. Additionally, we found that error rates were minimized for SNPs called by more than one variant caller. Finally, we evaluated the performance of various commonly used filtering metrics on SNP calling. Our analysis provides a quantitative assessment of the accuracy of five widely used variant calling programs and offers valuable insights into both the choice of variant caller program and the choice of filtering metrics, especially for researchers using non‐model study systems.

Genomic Prediction of Complex Traits in Perennial Plants: A Case for Forest Trees

This chapter provides an overview of the genomic selection progress in long-lived forest tree species. Factors affecting the prediction accuracy in genomic prediction are assessed with examples from empirical studies. Infrastructure and resources required for the implementation of genomic selection are evaluated. Some general guidelines are provided for the successful application of genomic selection in forest tree breeding programs.

MicroRNA-mediated post-transcriptional regulation of Pinus pinaster response and resistance to pinewood nematode

Pine wilt disease (PWD), caused by the parasitic nematode Bursaphelenchus xylophilus, or pinewood nematode (PWN), is a serious threat to pine forests in Europe. Pinus pinaster is highly susceptible to the disease and it is currently the most affected European pine species. In this work, we investigated the role of small RNAs (sRNAs) in regulating P. pinaster–PWN interaction in an early stage of infection. After performing an artificial PWN inoculation assay, we have identified 105 plant microRNAs (miRNAs) responsive to PWN. Based on their predicted targets, part of these miRNAs was associated with roles in jasmonate-response pathway, ROS detoxification, and terpenoid biosynthesis. Furthermore, by comparing resistant and susceptible plants, eight miRNAs with putative functions in plant defence and resistance to PWN have been identified. Finally, we explored the possibility of bidirectional trans-kingdom RNA silencing, identifying several P. pinaster genes putatively targeted by PWN miRNAs, which was supported by degradome analysis. Targets for P. pinaster miRNAs were also predicted in PWN, suggesting a role for trans-kingdom miRNA transfer and gene silencing both in PWN parasitism as in P. pinaster resistance to PWD. Our results provide new insights into previously unexplored roles of sRNA post-transcriptional regulation in P. pinaster response and resistance to PWN.

Identification and Tissue-Specific Expression Analysis of CYP720B Subfamily Genes in Slash Pine and Loblolly Pine

Diterpene resin acids (DRAs) are major components of pine oleoresin that can effectively resist the invasion of insects and pathogenic microorganisms. The subfamily of cytochrome P450s, CYP720B, catalyzes diterpene products into DRAs. Identifying CYP720B subfamily members and revealing the characteristics of tissue-specific expression would help understand diterpene-rich structures and diverse types. Slash pine and loblolly pine are important pines that provide oleoresin products. In this study, we identified CYP720B candidate genes based on the Pinus taeda V2.0 genome and full-length transcriptome of slash pine by PacBio. A total of 17 genes in slash pine and 19 in loblolly pine were identified and classified into four main clades by phylogenetic analysis. An analysis of cis-acting elements showed that CYP720B genes were closely related to adversity resistance. The gene expression of these candidates in different tissues was quantified by real-time quantitative PCR (RT–qPCR) analysis. Most of the genes showed relatively higher expression levels in roots and stems than in the other tissues, corresponding with the results of DRA component detection by gas chromatography–mass spectrometry (GC–MS), which indicated that stems and roots might be important tissues in oleoresin biosynthesis. These results provide a valuable resource for a better understanding of the biological role of individual CYP720Bs in slash pine and loblolly pine.

Genome-Wide Prediction of Transcription Start Sites in Conifers

The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.

Identifying and testing marker-trait associations for growth and phenology in three pine species: Implications for genomic prediction

In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 single nucleotide polymorphisms (SNPs) from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.

Achievements and Challenges of Genomics-Assisted Breeding in Forest Trees: From Marker-Assisted Selection to Genome Editing

Forest tree breeding efforts have focused mainly on improving traits of economic importance, selecting trees suited to new environments or generating trees that are more resilient to biotic and abiotic stressors. This review describes various methods of forest tree selection assisted by genomics and the main technological challenges and achievements in research at the genomic level. Due to the long rotation time of a forest plantation and the resulting long generation times necessary to complete a breeding cycle, the use of advanced techniques with traditional breeding have been necessary, allowing the use of more precise methods for determining the genetic architecture of traits of interest, such as genome-wide association studies (GWASs) and genomic selection (GS). In this sense, main factors that determine the accuracy of genomic prediction models are also addressed. In turn, the introduction of genome editing opens the door to new possibilities in forest trees and especially clustered regularly interspaced short palindromic repeats and CRISPR-associated protein 9 (CRISPR/Cas9). It is a highly efficient and effective genome editing technique that has been used to effectively implement targetable changes at specific places in the genome of a forest tree. In this sense, forest trees still lack a transformation method and an inefficient number of genotypes for CRISPR/Cas9. This challenge could be addressed with the use of the newly developing technique GRF-GIF with speed breeding.

A genome-wide SNP genotyping resource for tropical pine tree species

We performed gene and genome targeted SNP discovery towards the development of a genome-wide, multi-species genotyping array for tropical pines. Pooled RNA-seq data from shoots of seedlings from five tropical pine species was used to identify transcript-based SNPs resulting in 1.3 million candidate Affymetrix SNP probe sets. In addition, we used a custom 40K probe set to perform capture-seq in pooled DNA from 81 provenances representing the natural ranges of six tropical pine species in Mexico and Central America resulting in 563K candidate SNP probe sets. Altogether, 300K RNA-seq (72%) and 120K capture-seq (28%) derived SNP probe sets were tiled on a 420K screening array that was used to genotype 576 trees representing the 81 provenances and commercial breeding material. Based on the screening array results, 50K SNPs were selected for commercial SNP array production (Axiom 384 format, Thermo Fisher Scientific) including 20K polymorphic SNPs for P. patula, P. tecunumanii, P. oocarpa and P. caribaea, 15K for P. greggii and P. maximinoi, 13K for P. elliottii and 8K for P. pseudostrobus. We included 9.7K ancestry informative SNPs that will be valuable for species and hybrid discrimination. Of the 50K SNP markers, 25% are polymorphic in only one species, while 75% are shared by two or more species. The Pitro50K SNP chip will be useful for population genomics and molecular breeding in this group of pine species that, together with their hybrids, represent the majority of fast-growing tropical and subtropical pine plantations globally.

Sequencing and Chromosome-Scale Assembly of Plant Genomes, Brassica rapa as a Use Case

With the rise of long-read sequencers and long-range technologies, delivering high-quality plant genome assemblies is no longer reserved to large consortia. Not only sequencing techniques, but also computer algorithms have reached a point where the reconstruction of assemblies at the chromosome scale is now feasible at the laboratory scale. Current technologies, in particular long-range technologies, are numerous, and selecting the most promising one for the genome of interest is crucial to obtain optimal results. In this study, we resequenced the genome of the yellow sarson, Brassica rapa cv. Z1, using the Oxford Nanopore PromethION sequencer and assembled the sequenced data using current assemblers. To reconstruct complete chromosomes, we used and compared three long-range scaffolding techniques, optical mapping, Omni-C, and Pore-C sequencing libraries, commercialized by Bionano Genomics, Dovetail Genomics, and Oxford Nanopore Technologies, respectively, or a combination of the three, in order to evaluate the capability of each technology.

Evolutionary history of the mediterranean Pinus halepensis-brutia species complex using gene-resequencing and transcriptomic approaches

Complementary gene-resequencing and transcriptomic approaches reveal contrasted evolutionary histories in a species complex. Pinus halepensis and Pinus brutia are closely related species that can intercross, but occupy different geographical ranges and bioclimates. To study the evolution of this species complex and to provide genomic resources for further research, we produce and analyze two new complementary sets of genetic resources: (i) a set of 172 re-sequenced genomic target loci analyzed in 45 individuals, and (ii) a set of 11 transcriptome assemblies. These two datasets provide insights congruent with previous studies: P. brutia displays high level of genetic diversity and no genetic sub-structure, while P. halepensis shows three main genetic clusters, the western Mediterranean and North African clusters displaying much lower genetic diversity than the eastern Mediterranean cluster, the latter cluster having similar genetic diversity to P. brutia. In addition, these datasets provide new insights on the timing of the species-complex history: the two species would have split at the end of the tertiary, and the changing climatic conditions of the Mediterranean region at the end of the Tertiary-beginning of the Quaternary, together with the distinct species tolerance to harsh climatic conditions would have resulted in different geographic distributions, demographic histories and genetic patterns of the two pines. The multiple glacial-interglacial cycles during the Quaternary would have led to the expansion of P. brutia in the Middle East, while P. halepensis would have been through bottlenecks. The last glaciations, from 0.6 Mya on, would have affected further the Western genetic pool of P. halepensis.

Toward genomic selection in Pinus taeda: Integrating resources to support array design in a complex conifer genome

PREMISE

An informatics approach was used for the construction of an Axiom genotyping array from heterogeneous, high-throughput sequence data to assess the complex genome of loblolly pine (Pinus taeda).

METHODS

High-throughput sequence data, sourced from exome capture and whole genome reduced-representation approaches from 2698 trees across five sequence populations, were analyzed with the improved genome assembly and annotation for the loblolly pine. A variant detection, filtering, and probe design pipeline was developed to detect true variants across and within populations. From 8.27 million variants, a total of 642,275 were evaluated and 423,695 of those were screened across a range-wide population.

RESULTS

The final informatics and screening approach delivered an Axiom array representing 46,439 high-confidence variants to the forest tree breeding and genetics community. Based on the annotated reference genome, 34% were located in or directly upstream or downstream of genic regions.

DISCUSSION

The Pita50K array represents a genome-wide resource developed from sequence data for an economically important conifer, loblolly pine. It uniquely integrates independent projects that assessed trees sampled across the native range. The challenges associated with the large and repetitive genome are addressed in the development of this resource.

Breeding for Climate Change Resilience: A Case Study of Loblolly Pine (Pinus taeda L.) in North America

Earth's atmosphere is warming and the effects of climate change are becoming evident. A key observation is that both the average levels and the variability of temperature and precipitation are changing. Information and data from new technologies are developing in parallel to provide multidisciplinary opportunities to address and overcome the consequences of these changes in forest ecosystems. Changes in temperature and water availability impose multidimensional environmental constraints that trigger changes from the molecular to the forest stand level. These can represent a threat for the normal development of the tree from early seedling recruitment to adulthood both through direct mortality, and by increasing susceptibility to pathogens, insect attack, and fire damage. This review summarizes the strengths and shortcomings of previous work in the areas of genetic variation related to cold and drought stress in forest species with particular emphasis on loblolly pine (Pinus taeda L.), the most-planted tree species in North America. We describe and discuss the implementation of management and breeding strategies to increase resilience and adaptation, and discuss how new technologies in the areas of engineering and genomics are shaping the future of phenotype-genotype studies. Lessons learned from the study of species important in intensively-managed forest ecosystems may also prove to be of value in helping less-intensively managed forest ecosystems adapt to climate change, thereby increasing the sustainability and resilience of forestlands for the future.

Candidate Genes for the High-Altitude Adaptations of Two Mountain Pine Taxa

Mountain plants, challenged by vegetation time contractions and dynamic changes in environmental conditions, developed adaptations that help them to balance their growth, reproduction, survival, and regeneration. However, knowledge regarding the genetic basis of species adaptation to higher altitudes remain scarce for most plant species. Here, we attempted to identify such corresponding genomic regions of high evolutionary importance in two closely related European pines, Pinus mugo and P. uncinata, contrasting them with a reference lowland relative—P. sylvestris. We genotyped 438 samples at thousands of single nucleotide polymorphism (SNP) markers, tested their genetic differentiation and population structure followed by outlier detection and gene ontology annotations. Markers clearly differentiated the species and uncovered patterns of population structure in two of them. In P. uncinata three Pyrenean sites were grouped together, while two outlying populations constituted a separate cluster. In P. sylvestris, Spanish population appeared distinct from the remaining four European sites. Between mountain pines and the reference species, 35 candidate genes for altitude-dependent selection were identified, including such encoding proteins responsible for photosynthesis, photorespiration and cell redox homeostasis, regulation of transcription, and mRNA processing. In comparison between two mountain pines, 75 outlier SNPs were found in proteins involved mainly in the gene expression and metabolism.

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.

Transcript Dynamics in Wounded and Inoculated Scots Pine

Comparative transcriptome analysis provides a useful tool for the exploration of plant-pathogen interaction by allowing in-depth comparison of gene expression between unaffected, inoculated and wounded organisms. Here we present the results of comparative transcriptome analysis in genetically identical one-year-old Scots pine ramets after wounding and inoculation with Heterobasidion annosum. We identified 230 genes that were more than 2-fold upregulated in inoculated samples (compared to controls) and 116 downregulated genes. Comparison of inoculated samp les with wounded samples identified 32 differentially expressed genes (30 were upregulated after inoculation). Several of the genes upregulated after inoculation are involved in protection from oxidative stress, while genes involved in photosynthesis, water transport and drought stress tolerance were downregulated. An NRT3 family protein was the most upregulated transcript in response to both inoculation and wounding, while a U-box domain-containing protein gene was the most upregulated gene comparing inoculation to wounding. The observed transcriptome dynamics suggest involvement of auxin, ethylene, jasmonate, gibberellin and reactive oxygen species pathways and cell wall modification regulation in response to H. annosum infection. The results are compared to methyl jasmonate induced transcriptome dynamics.

Comparative Study of Pine Reference Genomes Reveals Transposable Element Interconnected Gene Networks

Sequencing the giga-genomes of several pine species has enabled comparative genomic analyses of these outcrossing tree species. Previous studies have revealed the wide distribution and extraordinary diversity of transposable elements (TEs) that occupy the large intergenic spaces in conifer genomes. In this study, we analyzed the distribution of TEs in gene regions of the assembled genomes of Pinus taeda and Pinus lambertiana using high-performance computing resources. The quality of draft genomes and the genome annotation have significant consequences for the investigation of TEs and these aspects are discussed. Several TE families frequently inserted into genes or their flanks were identified in both species’ genomes. Potentially important sequence motifs were identified in TEs that could bind additional regulatory factors, promoting gene network formation with faster or enhanced transcription initiation. Node genes that contain many TEs were observed in multiple potential transposable element-associated networks. This study demonstrated the increased accumulation of TEs in the introns of stress-responsive genes of pines and suggests the possibility of rewiring them into responsive networks and sub-networks interconnected with node genes containing multiple TEs. Many such regulatory influences could lead to the adaptive environmental response clines that are characteristic of naturally spread pine populations.

Molecular study of drought response in the Mediterranean conifer Pinus pinaster Ait.: Differential transcriptomic profiling reveals constitutive water deficit-independent drought tolerance mechanisms

Adaptation of long-living forest trees to respond to environmental changes is essential to secure their performance under adverse conditions. Water deficit is one of the most significant stress factors determining tree growth and survival. Maritime pine (Pinus pinaster Ait.), the main source of softwood in southwestern Europe, is subjected to recurrent drought periods which, according to climate change predictions for the years to come, will progressively increase in the Mediterranean region. The mechanisms regulating pine adaptive responses to environment are still largely unknown. The aim of this work was to go a step further in understanding the molecular mechanisms underlying maritime pine response to water stress and drought tolerance at the whole plant level. A global transcriptomic profiling of roots, stems, and needles was conducted to analyze the performance of siblings showing contrasted responses to water deficit from an ad hoc designed full-sib family. Although P. pinaster is considered a recalcitrant species for vegetative propagation in adult phase, the analysis was conducted using vegetatively propagated trees exposed to two treatments: well-watered and moderate water stress. The comparative analyses led us to identify organ-specific genes, constitutively expressed as well as differentially expressed when comparing control versus water stress conditions, in drought-sensitive and drought-tolerant genotypes. Different response strategies can point out, with tolerant individuals being pre-adapted for coping with drought by constitutively expressing stress-related genes that are detected only in latter stages on sensitive individuals subjected to drought.

Comparative Transcriptome Analysis of Pine Trees Treated with Resistance-Inducing Substances against the Nematode Bursaphelenchus xylophilus

The pinewood nematode (PWN) Bursaphelenchus xylophilus causes pine wilt disease, which results in substantial economic and environmental losses across pine forests worldwide. Although systemic acquired resistance (SAR) is effective in controlling PWN, the detailed mechanisms underlying the resistance to PWN are unclear. Here, we treated pine samples with two SAR elicitors, acibenzolar-S-methyl (ASM) and methyl salicylic acid (MeSA) and constructed an in vivo transcriptome of PWN-infected pines under SAR conditions. A total of 252 million clean reads were obtained and mapped onto the reference genome. Compared with untreated pines, 1091 and 1139 genes were differentially upregulated following the ASM and MeSA treatments, respectively. Among these, 650 genes showed co-expression patterns in response to both SAR elicitors. Analysis of these patterns indicated a functional linkage among photorespiration, peroxisome, and glycine metabolism, which may play a protective role against PWN infection-induced oxidative stress. Further, the biosynthesis of flavonoids, known to directly control parasitic nematodes, was commonly upregulated under SAR conditions. The ASM- and MeSA-specific expression patterns revealed functional branches for myricetin and quercetin production in flavonol biosynthesis. This study will enhance the understanding of the dynamic interactions between pine hosts and PWN under SAR conditions.

Development of a single nucleotide polymorphism array for population genomic studies in four European pine species

Pines are some of the most ecologically and economically important tree species in the world, and many have enormous natural distributions or have been extensively planted. However, a lack of rapid genotyping capability is hampering progress in understanding the molecular basis of genetic variation in these species. Here, we deliver an efficient tool for genotyping thousands of single nucleotide polymorphism (SNP) markers across the genome that can be applied to genetic studies in pines. Polymorphisms from resequenced candidate genes and transcriptome sequences of P. sylvestris, P. mugo, P. uncinata, P. uliginosa and P. radiata were used to design a 49,829 SNP array (Axiom_PineGAP, Thermo Fisher). Over a third (34.68%) of the unigenes identified from the P. sylvestris transcriptome were represented on the array, which was used to screen samples of four pine species. The conversion rate for the array on all samples was 42% (N = 20,795 SNPs) and was similar for SNPs sourced from resequenced candidate gene and transcriptome sequences. The broad representation of gene ontology terms by unigenes containing converted SNPs reflected their coverage across the full transcriptome. Over a quarter of successfully converted SNPs were polymorphic among all species, and the data were successful in discriminating among the species and some individual populations. The SNP array provides a valuable new tool to advance genetic studies in these species and demonstrates the effectiveness of the technology for rapid genotyping in species with large and complex genomes.

Two Independent Plastid accD Transfers to the Nuclear Genome of Gnetum and Other Insights on Acetyl-CoA Carboxylase Evolution in Gymnosperms

Acetyl-CoA carboxylase (ACCase) is the key regulator of fatty acid biosynthesis. In most plants, ACCase exists in two locations (cytosol and plastids) and in two forms (homomeric and heteromeric). Heteromeric ACCase comprises four subunits, three of them (ACCA-C) are nuclear encoded (nr) and the fourth (ACCD) is usually plastid encoded. Homomeric ACCase is encoded by a single nr-gene (ACC). We investigated the ACCase gene evolution in gymnosperms by examining the transcriptomes of newly sequenced Gnetum ula, combined with 75 transcriptomes and 110 plastomes of other gymnosperms. AccD-coding sequences are elongated through the insertion of repetitive DNA in four out of five cupressophyte families (except Sciadopityaceae) and were functionally transferred to the nucleus of gnetophytes and Sciadopitys. We discovered that, among the three genera of gnetophytes, only Gnetum has two copies of nr-accD. Furthermore, using protoplast transient expression assays, we experimentally verified that the nr-accD precursor proteins in Gnetum and Sciadopitys can be delivered to the plastids. Of the two nr-accD copies of Gnetum, one dually targets plastids and mitochondria, whereas the other potentially targets plastoglobuli. The distinct transit peptides, gene architectures, and flanking sequences between the two Gnetum accDs suggest that they have independent origins. Our findings are the first account of two distinctly targeted nr-accDs of any green plants and the most comprehensive analyses of ACCase evolution in gymnosperms to date.

Construction of genetic linkage map and identification of a novel major locus for resistance to pine wood nematode in Japanese black pine (Pinus thunbergii)

BACKGROUND

Pine wilt disease (PWD), which is caused by the pine wood nematode (PWN) Bursaphelenchus xylophilus, is currently the greatest threat to pine forests in Europe and East Asian countries including Japan. Constructing a detailed linkage map of DNA markers and identifying PWD resistance genes/loci lead to improved resistance in Pinus thunbergii, as well as other Pinus species that are also susceptible to PWD.

RESULTS

A total F₁ mapping population of 188 individuals derived from a cross between the PWD-resistant P. thunbergii varieties 'Tanabe 54' (resistant rank 2 to PWD) and 'Tosashimizu 63' (resistant rank 4 to PWD) was inoculated with PWN, and was evaluated for disease symptoms. To perform linkage analysis for PWN resistance, a set of three maps was constructed; two parental maps generated using the integrated two-way pseudo-testcross method, and a consensus map with population-type cross-pollination. The linkage map of 'Tanabe 54' consisted of 167 loci, and covered 14 linkage groups (LGs), with a total genetic distance of 1214.6 cM. The linkage map of 'Tosashimizu 63' consisted of 252 loci, and covered 14 LGs, with a total genetic distance of 1422.1 cM. The integrated consensus map comprised 12 LGs with the basic chromosome number of P. thunbergii, and a total genetic distance of 1403.6 cM. Results from quantitative trait loci (QTL) analysis using phenotype data and linkage maps indicated that PWN resistance is controlled by a single dominant allele, which was derived from the 'Tanabe 54' female parent. This major QTL was located on linkage group 3 and was designated PWD1 for PINE WILT DISEASE 1.

CONCLUSIONS

The PWD1 locus is a major resistance QTL located on the Pinus consensus LG03 that acts in a dominant manner to confer pine wood nematode resistance. Information from the present study will be useful for P. thunbergii breeding programs to improve resistance to PWD, and also to help identify susceptibility genes in Pinus species.