Forest genomics: Advancing climate adaptation, forest health, productivity, and conservation

The Complex Genetic Legacy of Hybridization and Introgression between the Rare Ocotea loxensis van der Werff and the Widespread O. infrafoveolata van der Werff (Lauraceae)

Hybridization and introgression are complex evolutionary mechanisms that can increase species diversity and lead to speciation, but may also lead to species extinction. In this study, we tested the presence and genetic consequences of hybridization between the rare and Ecuadorian endemic O. loxensis van der Werff and the widespread O. infrafoveolata van der Werff (Lauraceae). Phenotypically, some trees are difficult to identify, and we expect that some might in fact be cryptic hybrids. Thus, we developed nuclear microsatellites to assess the existence of hybrids, as well as the patterns of genetic diversity and population structure in allopatric and sympatric populations. The results revealed high levels of genetic diversity, even in the rare O. loxensis, being usually significantly higher in sympatric than in allopatric populations. The Bayesian assignment of individuals into different genetic classes revealed a complex scenario with different hybrid generations occurring in all sympatric populations, but also in allopatric ones. The absence of some backcrossed hybrids suggests the existence of asymmetric gene flow, and that some hybrids might be more fitted than others might. The existence of current and past interspecific gene flow also explains the blurring of species boundaries in these species and could be linked to the high rates of species found in Ocotea.

Evolutionary divergence of CXE gene family in green plants unveils that PtoCXEs overexpression reduces fungal colonization in transgenic Populus

Plant enzymes significantly contribute to the rapidly diversified metabolic repertoire since the colonization of land by plants. Carboxylesterase is just one of the ubiquitous, multifunctional and ancient enzymes that has particularly diversified during plant evolution. This study provided a status on the carboxylesterase landscape within Viridiplantae. A total of 784 carboxylesterases were identified from the genome of 31 plant species representing nine major lineages of sequenced Viridiplantae and divided into five clades based on phylogenetic analysis. Clade I carboxylesterase genes may be of bacterial origin and then expanded and diversified during plant evolution. Clade II was first gained in the ancestor of bryophytes after colonization of land by plants, Clade III and Clade IV in ferns which were considered the most advanced seedless vascular plants, while Clade V was gained in seed plants. To date, the functions of carboxylesterase genes in woody plants remain unclear. In this study, 51 carboxylesterase genes were identified from the genome of Populus trichocarpa and further divided into eight classes. Tandem and segmental duplication events both contributed to the expansion of carboxylesterase genes in Populus. Although carboxylesterase genes were proven to enhance resistance to pathogens in many herbaceous species, relevant researches on forest trees are still needed. In this study, pathogen incubation assays showed that overexpressing of six Class VI carboxylesterases in Populus tomentosa, to a greater or lesser degree, reduced colonization of detached leaves by fungus Cytospora chrysosperma. A significant difference was also found in functional divergence patterns for genes derived from different gene duplication events. Functional differentiation of duplicated carboxylesterase genes in Populus was proved for the first time by in vivo physiological analysis. The identification of the potentially anti-fungal PtoCXE06 gene also laid a theoretical foundation for promoting the genetic improvement of disease-resistance traits in forest trees.

Planting long-lived trees in a warming climate: Theory shows the importance of stage-dependent climatic tolerance

Climate change poses a particular threat to long-lived trees, which may not adapt or migrate fast enough to keep up with rising temperatures. Assisted gene flow could facilitate adaptation of populations to future climates by using managed translocation of seeds from a warmer location (provenance) within the current range of a species. Finding the provenance that will perform best in terms of survival or growth is complicated by a trade-off. Because trees face a rapidly changing climate during their long lives, the alleles that confer optimal performance may vary across their lifespan. For instance, trees from warmer provenances could be well adapted as adults but suffer from colder temperatures while juvenile. Here we use a stage-structured model, using both analytical predictions and numerical simulations, to determine which provenance would maximize the survival of a cohort of long-lived trees in a changing climate. We parameterize our simulations using empirically estimated demographic transition matrices for 20 long-lived tree species. Unable to find reliable quantitative estimates of how climatic tolerance changes across stages in these same species, we varied this parameter to study its effect. Both our mathematical model and simulations predict that the best provenance depends strongly on how fast the climate changes and also how climatic tolerance varies across the lifespan of a tree. We thus call for increased empirical efforts to measure how climate tolerance changes over life in long-lived species, as our model suggests that it should strongly influence the best provenance for assisted gene flow.

Genome sequences and population genomics reveal climatic adaptation and genomic divergence between two closely related sweetgum species

Understanding the genetic basis of population divergence and adaptation is an important goal in population genetics and evolutionary biology. However, the relative roles of demographic history, gene flow, and/or selective regime in driving genomic divergence, climatic adaptation, and speciation in non-model tree species are not yet fully understood. To address this issue, we generated whole-genome resequencing data of Liquidambar formosana and L. acalycina, which are broadly sympatric but altitudinally segregated in the Tertiary relict forests of subtropical China. We integrated genomic and environmental data to investigate the demographic history, genomic divergence, and climatic adaptation of these two sister species. We inferred a scenario of allopatric species divergence during the late Miocene, followed by secondary contact during the Holocene. We identified multiple genomic islands of elevated divergence that mainly evolved through divergence hitchhiking and recombination rate variation, likely fostered by long-term refugial isolation and recent differential introgression in low-recombination genomic regions. We also found some candidate genes with divergent selection signatures potentially involved in climatic adaptation and reproductive isolation. Our results contribute to a better understanding of how late Tertiary/Quaternary climatic change influenced speciation, genomic divergence, climatic adaptation, and introgressive hybridization in East Asia's Tertiary relict flora. In addition, they should facilitate future evolutionary, conservation genomics, and molecular breeding studies in Liquidambar, a genus of important medicinal and ornamental values.

Sustainable Development versus Extractivist Deforestation in Tropical, Subtropical, and Boreal Forest Ecosystems: Repercussions and Controversies about the Mother Tree and the Mycorrhizal Network Hypothesis

This research reviews the phenomenon of extractive deforestation as a possible trigger for cascade reactions that could affect part of the forest ecosystem and its biodiversity (surface, aerial, and underground) in tropical, subtropical, and boreal forests. The controversy and disparities in criteria generated in the international scientific community around the hypothesis of a possible link between "mother trees" and mycorrhizal networks in coopetition for nutrients, nitrogen, and carbon are analyzed. The objective is to promote awareness to generate more scientific knowledge about the eventual impacts of forest extraction. Public policies are emphasized as crucial mediators for balanced sustainable development. Currently, the effects of extractive deforestation on forest ecosystems are poorly understood, which requires caution and forest protection. Continued research to increase our knowledge in molecular biology is advocated to understand the adaptation of biological organisms to the new conditions of the ecosystem both in the face of extractive deforestation and reforestation. The environmental impacts of extractive deforestation, such as the loss of biodiversity, soil degradation, altered water cycles, and the contribution of climate change, remain largely unknown. Long-term and high-quality research is essential to ensure forest sustainability and the preservation of biodiversity for future generations.

Under Stress: Searching for Genes Involved in the Response of Abies pinsapo Boiss to Climate Change

Currently, Mediterranean forests are experiencing the deleterious effects of global warming, which mainly include increased temperatures and decreased precipitation in the region. Relict Abies pinsapo fir forests, endemic in the southern Iberian Peninsula, are especially sensitive to these recent environmental disturbances, and identifying the genes involved in the response of this endangered tree species to climate-driven stresses is of paramount importance for mitigating their effects. Genomic resources for A. pinsapo allow for the analysis of candidate genes reacting to warming and aridity in their natural habitats. Several members of the complex gene families encoding late embryogenesis abundant proteins (LEAs) and heat shock proteins (HSPs) have been found to exhibit differential expression patterns between wet and dry seasons when samples from distinct geographical locations and dissimilar exposures to the effects of climate change were analyzed. The observed changes were more perceptible in the roots of trees, particularly in declining forests distributed at lower altitudes in the more vulnerable mountains. These findings align with previous studies and lay the groundwork for further research on the molecular level. Molecular and genomic approaches offer valuable insights for mitigating climate stress and safeguarding this endangered conifer.

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.

Comparative transcriptomics reveals divergence in pathogen response gene families amongst 20 forest tree species

Forest trees provide critical ecosystem services for humanity that are under threat due to ongoing global change. Measuring and characterizing genetic diversity are key to understanding adaptive potential and developing strategies to mitigate negative consequences arising from climate change. In the area of forest genetic diversity, genetic divergence caused by large-scale changes at the chromosomal level has been largely understudied. In this study, we used the RNA-seq data of 20 co-occurring forest trees species from genera including Acer, Alnus, Amelanchier, Betula, Cornus, Corylus, Dirca, Fraxinus, Ostrya, Populus, Prunus, Quercus, Ribes, Tilia, and Ulmus sampled from Upper Peninsula of Michigan. These data were used to infer the origin and maintenance of gene family variation, species divergence time, as well as gene family expansion and contraction. We identified a signal of common whole genome duplication events shared by core eudicots. We also found rapid evolution, namely fast expansion or fast contraction of gene families, in plant-pathogen interaction genes amongst the studied diploid species. Finally, the results lay the foundation for further research on the genetic diversity and adaptive capacity of forest trees, which will inform forest management and conservation policies.

Genomic signals of local adaptation in Picea crassifolia

BACKGROUND

Global climate change poses a grave threat to biodiversity and underscores the importance of identifying the genes and corresponding environmental factors involved in the adaptation of tree species for the purposes of conservation and forestry. This holds particularly true for spruce species, given their pivotal role as key constituents of the montane, boreal, and sub-alpine forests in the Northern Hemisphere.

RESULTS

Here, we used transcriptomes, species occurrence records, and environmental data to investigate the spatial genetic distribution of and the climate-associated genetic variation in Picea crassifolia. Our comprehensive analysis employing ADMIXTURE, principal component analysis (PCA) and phylogenetic methodologies showed that the species has a complex population structure with obvious differentiation among populations in different regions. Concurrently, our investigations into isolation by distance (IBD), isolation by environment (IBE), and niche differentiation among populations collectively suggests that local adaptations are driven by environmental heterogeneity. By integrating population genomics and environmental data using redundancy analysis (RDA), we identified a set of climate-associated single-nucleotide polymorphisms (SNPs) and showed that environmental isolation had a more significant impact than geographic isolation in promoting genetic differentiation. We also found that the candidate genes associated with altitude, temperature seasonality (Bio4) and precipitation in the wettest month (Bio13) may be useful for forest tree breeding.

CONCLUSIONS

Our findings deepen our understanding of how species respond to climate change and highlight the importance of integrating genomic and environmental data in untangling local adaptations.

Genome-wide association studies: an intuitive solution for SNP identification and gene mapping in trees

Analysis of natural diversity in wild/cultivated plants can be used to understand the genetic basis for plant breeding programs. Recent advancements in DNA sequencing have expanded the possibilities for genetically altering essential features. There have been several recently disclosed statistical genetic methods for discovering the genes impacting target qualities. One of these useful methods is the genome-wide association study (GWAS), which effectively identifies candidate genes for a variety of plant properties by examining the relationship between a molecular marker (such as SNP) and a target trait. Conventional QTL mapping with highly structured populations has major limitations. The limited number of recombination events results in poor resolution for quantitative traits. Only two alleles at any given locus can be studied simultaneously. Conventional mapping approach fails to work in perennial plants and vegetatively propagated crops. These limitations are sidestepped by association mapping or GWAS. The flexibility of GWAS comes from the fact that the individuals being examined need not be linked to one another, allowing for the use of all meiotic and recombination events to increase resolution. Phenotyping, genotyping, population structure analysis, kinship analysis, and marker-trait association analysis are the fundamental phases of GWAS. With the rapid development of sequencing technologies and computational methods, GWAS is becoming a potent tool for identifying the natural variations that underlie complex characteristics in crops. The use of high-throughput sequencing technologies along with genotyping approaches like genotyping-by-sequencing (GBS) and restriction site associated DNA (RAD) sequencing may be highly useful in fast-forward mapping approach like GWAS. Breeders may use GWAS to quickly unravel the genomes through QTL and association mapping by taking advantage of natural variances. The drawbacks of conventional linkage mapping can be successfully overcome with the use of high-resolution mapping and the inclusion of multiple alleles in GWAS.

Leveraging explainable machine learning models to assess forest health: A case study in Hainan, China

Global forest area has declined over the past few years, forest quality has declined, and ecological and environmental events have increased with climate change and human activity. In the context of ecological civilization, forest health issues have received unprecedented attention. By improving forest health, forests can better perform their ecosystem service functions and promote green development. This study was carried out in the WuZhi Shan area of Hainan Tropical Rainforest National Park. We employed a decision tree algorithm, a machine learning technique, for our modeling due to its high accuracy and interpretability. The objective weighted method using criteria of importance through intercriteria correlation (CRITIC) was used to determine forest health classes based on survey and experimental data from 132 forest samples. The results showed that species diversity is the most important metric to measure forest health. An interpretable decision tree machine learning model was proposed to incorporate forest health indicators, providing up to 90% accuracy in the classification of forest health conditions. The model demonstrated a high degree of effectiveness, achieving an average precision of 90%, a recall of 67%, and an F1 score of 70.2% in predicting forest health. The interpretable decision tree classification results showed that breast height diameter is the most important variable in classifying the health status of both primary and secondary forests. This study highlights the importance of using interpretable machine learning methods for the decision-making process. Our work contributes to the scientific underpinnings of sustainable forest development and effective conservation planning.

Applying molecular and genetic methods to trees and their fungal communities

Forests provide invaluable economic, ecological, and social services. At the same time, they are exposed to several threats, such as fragmentation, changing climatic conditions, or increasingly destructive pests and pathogens. Trees, the inherent species of forests, cannot be viewed as isolated organisms. Manifold (micro)organisms are associated with trees playing a pivotal role in forest ecosystems. Of these organisms, fungi may have the greatest impact on the life of trees. A multitude of molecular and genetic methods are now available to investigate tree species and their associated organisms. Due to their smaller genome sizes compared to tree species, whole genomes of different fungi are routinely compared. Such studies have only recently started in forest tree species. Here, we summarize the application of molecular and genetic methods in forest conservation genetics, tree breeding, and association genetics as well as for the investigation of fungal communities and their interrelated ecological functions. These techniques provide valuable insights into the molecular basis of adaptive traits, the impacts of forest management, and changing environmental conditions on tree species and fungal communities and can enhance tree-breeding cycles due to reduced time for field testing. It becomes clear that there are multifaceted interactions among microbial species as well as between these organisms and trees. We demonstrate the versatility of the different approaches based on case studies on trees and fungi. KEY POINTS: • Current knowledge of genetic methods applied to forest trees and associated fungi. • Genomic methods are essential in conservation, breeding, management, and research. • Important role of phytobiomes for trees and their ecosystems.

Small-scale genetic structure of coral populations in Palau based on whole mitochondrial genomes: Implications for future coral resilience

The ability of local populations to adapt to future climate conditions is facilitated by a balance between short range dispersal allowing local buildup of adaptively beneficial alleles, and longer dispersal moving these alleles throughout the species range. Reef building corals have relatively low dispersal larvae, but most population genetic studies show differentiation only over 100s of km. Here, we report full mitochondrial genome sequences from 284 tabletop corals (Acropora hyacinthus) from 39 patch reefs in Palau, and show two signals of genetic structure across reef scales from 1 to 55 km. First, divergent mitochondrial DNA haplotypes exist in different proportions from reef to reef, causing Phi_ST values of 0.02 (p = 0.02). Second, closely related sequences of mitochondrial Haplogroups are more likely to be co-located on the same reefs than expected by chance alone. We also compared these sequences to prior data on 155 colonies from American Samoa. In these comparisons, many Haplogroups in Palau were disproportionately represented or absent in American Samoa, and inter-regional Phi_ST = 0.259. However, we saw three instances of identical mitochondrial genomes between locations. Together, these data sets suggest two features of coral dispersal revealed by occurrence patterns in highly similar mitochondrial genomes. First, the Palau-American Samoa data suggest that long distance dispersal in corals is rare, as expected, but that it is common enough to deliver identical mitochondrial genomes across the Pacific. Second, higher than expected co-occurrence of Haplogroups on the same Palau reefs suggests greater retention of coral larvae on local reefs than predicted by many current oceanographic models of larval movement. Increased attention to local scales of coral genetic structure, dispersal, and selection may help increase the accuracy of models of future adaptation of corals and of assisted migration as a reef resilience intervention.

How ancient forest fragmentation and riparian connectivity generate high levels of genetic diversity in a microendemic Malagasy tree

Understanding landscape changes is central to predicting evolutionary trajectories and defining conservation practices. While human-driven deforestation is intense throughout Madagascar, exceptions in areas such as the Loky-Manambato region (north) raise questions regarding the causes and age of forest fragmentation. The Loky-Manambato region also harbours a rich and endemic flora, whose evolutionary origin remains poorly understood. We assessed the genetic diversity of an endangered microendemic Malagasy olive species (Noronhia spinifolia Hong-Wa) to better understand the vegetation dynamics in the Loky-Manambato region and its influence on past evolutionary processes. We characterized 72 individuals sampled across eight forests through nuclear and mitochondrial restriction-associated DNA sequencing data and chloroplast microsatellites. Combined population and landscape genetics analyses indicate that N. spinifolia diversity is largely explained by the current forest cover, highlighting a long-standing habitat mosaic in the region. This sustains a major and long-term role of riparian corridors in maintaining connectivity across these antique mosaic habitats, calling for the study of organismal interactions that promote gene flow.

Improving Coconut Using Modern Breeding Technologies: Challenges and Opportunities

Coconut (Cocos nucifera L.) is a perennial palm with a wide range of distribution across tropical islands and coastlines. Multitude use of coconut by nature is important in the socio-economic fabric framework among rural smallholders in producing countries. It is a major source of income for 30 million farmers, while 60 million households rely on the coconut industry directly as farm workers and indirectly through the distribution, marketing, and processing of coconut and coconut-based products. Stagnant production, inadequate planting materials, the effects of climate change, as well as pests and diseases are among the key issues that need to be urgently addressed in the global coconut industry. Biotechnology has revolutionized conventional breeding approaches in creating genetic variation for trait improvement in a shorter period of time. In this review, we highlighted the challenges of current breeding strategies and the potential of biotechnological approaches, such as genomic-assisted breeding, next-generation sequencing (NGS)-based genotyping and genome editing tools in improving the coconut. Also, combining these technologies with high-throughput phenotyping approaches and speed breeding could speed up the rate of genetic gain in coconut breeding to solve problems that have been plaguing the industry for decades.

Genomic insights into local adaptation and future climate-induced vulnerability of a keystone forest tree in East Asia

Rapid global climate change is posing a substantial threat to biodiversity. The assessment of population vulnerability and adaptive capacity under climate change is crucial for informing conservation and mitigation strategies. Here we generate a chromosome-scale genome assembly and re-sequence genomes of 230 individuals collected from 24 populations for Populus koreana, a pioneer and keystone tree species in temperate forests of East Asia. We integrate population genomics and environmental variables to reveal a set of climate-associated single-nucleotide polymorphisms, insertion/deletions and structural variations, especially numerous adaptive non-coding variants distributed across the genome. We incorporate these variants into an environmental modeling scheme to predict a highly spatiotemporal shift of this species in response to future climate change. We further identify the most vulnerable populations that need conservation priority and many candidate genes and variants that may be useful for forest tree breeding with special aims. Our findings highlight the importance of integrating genomic and environmental data to predict adaptive capacity of a key forest to rapid climate change in the future.

Assisted migration is plausible for a boreal tree species under climate change: A quantitative and population genetics study of trembling aspen (Populus tremuloides Michx.) in western Canada

A novel method was tested for improving tree breeding strategies that integrate quantitative and population genetics based on range‐wide reciprocal transplant experiments. Five reciprocal common garden tests of Populus tremuloides were investigated including 6450 trees across western Canada focusing on adaptation traits and growth. Both genetic parameters and home‐site transplant models were evaluated. We found a genetic trade‐off between growth and early spring leaf flush and late fall senescence. Coefficients of phenotypic variation (CVp) of cell lysis (CL), a measure of freezing injury, shrank from 0.28 to 0.10 during acclimation in the fall, and the CVp slope versus the freezing temperature was significantly different from zero (R² = 0.33, p = .02). There was more between‐population genetic variation in fall phenology than in spring leaf phenology. We suggest that P. tremuloides demonstrated a discrepancy between the ecological optimum and the physiological optimum minimum winter temperature. The sub‐optimal growing condition of P. tremuloides is potentially caused by the warmer ecological optimum than the physiological optimum. Assisted migration and breeding of fast growers to reforest cooler plantation sites can improve productivity. Transferring the study populations to less than 4°C of extreme minimum temperature appears safe for reforestation aligning with the historical recolonization direction of the species. This is equivalent to a 5–10° latitudinal northward movement. Fall frost hardiness is an effective criterion for family selection in the range tested in this study.

Using Genomic Selection to Develop Performance-Based Restoration Plant Materials

Effective native plant materials are critical to restoring the structure and function of extensively modified ecosystems, such as the sagebrush steppe of North America’s Intermountain West. The reestablishment of native bunchgrasses, e.g., bluebunch wheatgrass (Pseudoroegneria spicata [Pursh] À. Löve), is the first step for recovery from invasive species and frequent wildfire and towards greater ecosystem resiliency. Effective native plant material exhibits functional traits that confer ecological fitness, phenotypic plasticity that enables adaptation to the local environment, and genetic variation that facilitates rapid evolution to local conditions, i.e., local adaptation. Here we illustrate a multi-disciplinary approach based on genomic selection to develop plant materials that address environmental issues that constrain local populations in altered ecosystems. Based on DNA sequence, genomic selection allows rapid screening of large numbers of seedlings, even for traits expressed only in more mature plants. Plants are genotyped and phenotyped in a training population to develop a genome model for the desired phenotype. Populations with modified phenotypes can be used to identify plant syndromes and test basic hypotheses regarding relationships of traits to adaptation and to one another. The effectiveness of genomic selection in crop and livestock breeding suggests this approach has tremendous potential for improving restoration outcomes for species such as bluebunch wheatgrass.

Combining Spectral, Spatial-Contextual, and Structural Information in Multispectral UAV Data for Spruce Crown Delineation

Precise delineation of individual tree crowns is critical for accurate forest biophysical parameter estimation, species classification, and ecosystem modelling. Multispectral optical remote sensors mounted on low-flying unmanned aerial vehicles (UAVs) can rapidly collect very-high-resolution (VHR) photogrammetric optical data that contain the spectral, spatial, and structural information of trees. State-of-the-art tree crown delineation approaches rely mostly on spectral information and underexploit the spatial-contextual and structural information in VHR photogrammetric multispectral data, resulting in crown delineation errors. Here, we propose the spectral, spatial-contextual, and structural information-based individual tree crown delineation (S3-ITD) method, which accurately delineates individual spruce crowns by minimizing the undesirable effects due to intracrown spectral variance and nonuniform illumination/shadowing in VHR multispectral data. We evaluate the performance of the S3-ITD crown delineation method over a white spruce forest in Quebec, Canada. The highest mean intersection over union (IoU) index of 0.83, and the lowest mean crown-area difference (CAD) of 0.14 m2, proves the superior crown delineation performance of the S3-ITD method over state-of-the-art methods. The reduction in error by 2.4 cm and 1.0 cm for the allometrically derived diameter at breast height (DBH) estimates compared with those from the WS-ITD and the BF-ITD approaches, respectively, demonstrates the effectiveness of the S3-ITD method to accurately estimate biophysical parameters.

Climate Adaptation, Drought Susceptibility, and Genomic-Informed Predictions of Future Climate Refugia for the Australian Forest Tree Eucalyptus globulus

Understanding the capacity of forest tree species to adapt to climate change is of increasing importance for managing forest genetic resources. Through a genomics approach, we modelled spatial variation in climate adaptation within the Australian temperate forest tree Eucalyptus globulus, identified putative climate drivers of this genomic variation, and predicted locations of future climate refugia and populations at-risk of future maladaptation. Using 812,158 SNPs across 130 individuals from 30 populations (i.e., localities) spanning the species’ natural range, a gradientForest algorithm found 1177 SNPs associated with locality variation in home-site climate (climate-SNPs), putatively linking them to climate adaptation. Very few climate-SNPs were associated with population-level variation in drought susceptibility, signalling the multi-faceted nature and complexity of climate adaptation. Redundancy analysis (RDA) showed 24% of the climate-SNP variation could be explained by annual precipitation, isothermality, and maximum temperature of the warmest month. Spatial predictions of the RDA climate vectors associated with climate-SNPs allowed mapping of genomically informed climate selective surfaces across the species’ range under contemporary and projected future climates. These surfaces suggest over 50% of the current distribution of E. globulus will be outside the modelled adaptive range by 2070 and at risk of climate maladaptation. Such surfaces present a new integrated approach for natural resource managers to capture adaptive genetic variation and plan translocations in the face of climate change.

Pedigree reconstruction and spatial analysis for genetic testing and selection in a Larix kaempferi (Lamb.) Carrière plantation

BACKGROUND

Larix kaempferi is one of the major timber species in Northeast Asia. Demand for the reforestation of the species is rising in South Korea due to an increase in large timber production and utilization. However, progeny trials for the species have not been explored, making it challenging to foster advanced generations of tree improvement. In the present study, genetic testing and selection for diameter growth were conducted using pedigree reconstruction and phenotypic spatial distribution analysis in a plantation of L. kaempferi. The aim of the present study was to select the superior larch individuals using the pedigree reconstruction and phenotypic spatial distribution to substitute progeny trials. The plantation of seed orchard crops was established in 1990 and one-hundred and eighty-eight trees were selected as the study material. Genetic variation was investigated first to validate its adequacy as breeding material. Genetic testing was carried out using a model considering pedigree information and spatial autoregression of the phenotypes.

RESULTS

The expected heterozygosity of the mother trees and offspring were 0.672 and 0.681 presenting the corresponding level of genetic variation between two groups. The pedigree reconstruction using maternity analysis assigned one to six progenies to ninety-two candidate mothers. The accuracy of genetic testing was exceedingly increased with the animal model considering AR1 ⊗ AR1 structure compared to the animal model only. The estimated genetic variance of the former was 9.086 whereas that of the latter was 4.9E-5 for DBH. The predicted breeding values of the offspring for DBH were ranged from -5.937 cm to 5.655 cm and the estimated heritability of diameter growth was 0.344.

CONCLUSIONS

The genetic testing approach based on pedigree reconstruction and phenotypic spatial distribution analysis was considered a useful analytical scheme that could replace or supplement progeny trials.

Global evidence on the asymmetric response of gross primary productivity to interannual precipitation changes

Understanding gross primary productivity (GPP) response to precipitation (PPT) changes is essential for predicting land carbon uptake under increasing PPT variability and extremes. Previous studies found that ecosystem GPP may have an asymmetric response to PPT changes, leading to the inconsistency of GPP gains in wet years compared to GPP declines in dry years. However, it is unclear how the asymmetric responses vary among vegetation types and under different PPT variabilities. This study evaluated the global patterns of asymmetries of GPP response to different PPT changes using two state-of-science global GPP datasets. The result shows that under mild PPT changes (|ΔPPT| ≤ 25%), grasslands, savannas, shrublands, and tundra show positive asymmetric responses (i.e., larger GPP gains in wet years than GPP losses in dry years), while other vegetation types show negative asymmetric responses (i.e., larger GPP losses in dry years than GPP gains in wet years). Conversely, all vegetation types show negative GPP asymmetric responses to moderate (25% < |ΔPPT| ≤ 50%) and extreme (|ΔPPT| > 50%) PPT changes. Thus, we propose a new non-linear asymmetric GPP-PPT model that incorporates three modes with regards to vegetation types. Meanwhile, we found that the spatial patterns of asymmetry were mainly driven by PPT amount and variability. Stronger and negative asymmetries were found in areas with smaller PPT amount and variability, while positive asymmetries were found in areas with higher PPT variability. These findings promote our understanding of carbon dynamics under increased PPT variability and extremes and provide new insights for land models to better predict future carbon uptake and its feedback to climate change.

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.

Difficult climate-adaptive decisions in forests as complex social-ecological systems

Climate change threatens the social, ecological, and economic benefits enjoyed by forest-dependent communities worldwide. Climate-adaptive forest management strategies such as genomics-based assisted migration (AM) may help protect many of these threatened benefits. However, such novel technological interventions in complex social-ecological systems will generate new risks, benefits, and uncertainties that interact with diverse forest values and preexisting risks. Using data from 16 focus groups in British Columbia, Canada, we show that different stakeholders (forestry professionals, environmental nongovernmental organizations, local government officials, and members of local business communities) emphasize different kinds of risks and uncertainties in judging the appropriateness of AM. We show the difficulty of climate-adaptive decisions in complex social-ecological systems in which both climate change and adaptation will have widespread and cascading impacts on diverse nonclimate values. Overarching judgments about AM as an adaptation strategy, which may appear simple when elicited in surveys or questionnaires, require that participants make complex trade-offs among multiple domains of uncertain and unknown risks. Overall, the highest-priority forest management objective for most stakeholders is the health and integrity of the forest ecosystem from which all other important forest values derive. The factor perceived as riskiest is our lack of knowledge of how forest ecosystems work, which hinders stakeholders in their assessment of AM's acceptability. These results are further evidence of the inherent risk in privileging natural science above other forms of knowledge at the science-policy interface. When decisions are framed as technical, the normative and ethical considerations that define our fundamental goals are made invisible.

Meta-Analysis as a Tool to Identify Candidate Genes Involved in the Fagus sylvatica L. Abiotic Stress Response

In this study, we aimed to evaluate whether candidate genes for abiotic stresses in Fagus sylvatica L. are also candidate genes for herbaceous plants, with the purpose of better defining the abiotic stress response model of F. sylvatica. Therefore, a meta-analysis was performed on published papers related to abiotic stress. Firstly, we carried out a systematic review regarding the activity of 24 candidate genes selected for F. sylvatica under abiotic stress reported in 503 articles. After choosing the inclusion criteria, 73 articles out of 503, regarding 12 candidate genes, were included in this analysis. We performed an exploratory meta-analysis based on the random-effect model and the combined effect-size approach (Cohen’s d). The results obtained through Forest and Funnel plots indicate that the candidate genes for F. sylvatica are considered to be candidate genes in other herbaceous species. These results allowed us to set up models of plants’ response to abiotic stresses implementing the stress models in forest species. The results of this study will serve to bridge knowledge gaps regarding the pathways of response to abiotic stresses in trees based on the meta-analysis. The study approach used could be extended to observe larger gene databases and different species.

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.

Demographic history and natural selection shape patterns of deleterious mutation load and barriers to introgression across Populus genome

Hybridization and resulting introgression are important processes shaping the tree of life and appear to be far more common than previously thought. However, how the genome evolution was shaped by various genetic and evolutionary forces after hybridization remains unresolved. Here we used whole-genome resequencing data of 227 individuals from multiple widespread Populus species to characterize their contemporary patterns of hybridization and to quantify genomic signatures of past introgression. We observe a high frequency of contemporary hybridization and confirm that multiple previously ambiguous species are in fact F₁ hybrids. Seven species were identified, which experienced different demographic histories that resulted in strikingly varied efficacy of selection and burdens of deleterious mutations. Frequent past introgression has been found to be a pervasive feature throughout the speciation of these Populus species. The retained introgressed regions, more generally, tend to contain reduced genetic load and to be located in regions of high recombination. We also find that in pairs of species with substantial differences in effective population size, introgressed regions are inferred to have undergone selective sweeps at greater than expected frequencies in the species with lower effective population size, suggesting that introgression likely have higher potential to provide beneficial variation for species with small populations. Our results, therefore, illustrate that demography and recombination have interplayed with both positive and negative selection in determining the genomic evolution after hybridization.

Landscape Genomics in Tree Conservation Under a Changing Environment

Understanding the genetic basis of how species respond to changing environments is essential to the conservation of species. However, the molecular mechanisms of adaptation remain largely unknown for long-lived tree species which always have large population sizes, long generation time, and extensive gene flow. Recent advances in landscape genomics can reveal the signals of adaptive selection linking genetic variations and landscape characteristics and therefore have created novel insights into tree conservation strategies. In this review article, we first summarized the methods of landscape genomics used in tree conservation and elucidated the advantages and disadvantages of these methods. We then highlighted the newly developed method "Risk of Non-adaptedness," which can predict the genetic offset or genomic vulnerability of species via allele frequency change under multiple scenarios of climate change. Finally, we provided prospects concerning how our introduced approaches of landscape genomics can assist policymaking and improve the existing conservation strategies for tree species under the ongoing global changes.

Evolutionary Genetics of Crop-Wild Complexes

Since Darwin's time, the role of crop wild relatives (CWR), landraces, and cultivated genepools in shaping plant diversity and boosting food resources has been a major question [...].

Integrative Pre-Breeding for Biotic Resistance in Forest Trees

Climate change is unleashing novel biotic antagonistic interactions for forest trees that may jeopardize populations' persistence. Therefore, this review article envisions highlighting major opportunities from ecological evolutionary genomics to assist the identification, conservation, and breeding of biotic resistance in forest tree species. Specifically, we first discuss how assessing the genomic architecture of biotic stress resistance enables us to recognize a more polygenic nature for a trait typically regarded Mendelian, an expectation from the Fisherian runaway pathogen-host concerted arms-race evolutionary model. Secondly, we outline innovative pipelines to capture and harness natural tree pre-adaptations to biotic stresses by merging tools from the ecology, phylo-geography, and omnigenetics fields within a predictive breeding platform. Promoting integrative ecological genomic studies promises a better understanding of antagonistic co-evolutionary interactions, as well as more efficient breeding utilization of resistant phenotypes.

Induced Genetic Variations in Fruit Trees Using New Breeding Tools: Food Security and Climate Resilience

Fruit trees provide essential nutrients to humans by contributing to major agricultural outputs and economic growth globally. However, major constraints to sustainable agricultural productivity are the uncontrolled proliferation of the population, and biotic and abiotic stresses. Tree mutation breeding has been substantially improved using different physical and chemical mutagens. Nonetheless, tree plant breeding has certain crucial bottlenecks including a long life cycle, ploidy level, occurrence of sequence polymorphisms, nature of parthenocarpic fruit development and linkage. Genetic engineering of trees has focused on boosting quality traits such as productivity, wood quality, and resistance to biotic and abiotic stresses. Recent technological advances in genome editing provide a unique opportunity for the genetic improvement of woody plants. This review examines application of the CRISPR-Cas system to reduce disease susceptibility, alter plant architecture, enhance fruit quality, and improve yields. Examples are discussed of the contemporary CRISPR-Cas system to engineer easily scorable PDS genes, modify lignin, and to alter the flowering onset, fertility, tree architecture and certain biotic stresses.

Pests, diseases, and aridity have shaped the genome of Corymbia citriodora

Corymbia citriodora is a member of the predominantly Southern Hemisphere Myrtaceae family, which includes the eucalypts (Eucalyptus, Corymbia and Angophora; ~800 species). Corymbia is grown for timber, pulp and paper, and essential oils in Australia, South Africa, Asia, and Brazil, maintaining a high-growth rate under marginal conditions due to drought, poor-quality soil, and biotic stresses. To dissect the genetic basis of these desirable traits, we sequenced and assembled the 408 Mb genome of Corymbia citriodora, anchored into eleven chromosomes. Comparative analysis with Eucalyptus grandis reveals high synteny, although the two diverged approximately 60 million years ago and have different genome sizes (408 vs 641 Mb), with few large intra-chromosomal rearrangements. C. citriodora shares an ancient whole-genome duplication event with E. grandis but has undergone tandem gene family expansions related to terpene biosynthesis, innate pathogen resistance, and leaf wax formation, enabling their successful adaptation to biotic/abiotic stresses and arid conditions of the Australian continent.

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.

Genetic mixing for population management: From genetic rescue to provenancing

Animal and plant species around the world are being challenged by the deleterious effects of inbreeding, loss of genetic diversity, and maladaptation due to widespread habitat destruction and rapid climate change. In many cases, interventions will likely be needed to safeguard populations and species and to maintain functioning ecosystems. Strategies aimed at initiating, reinstating, or enhancing patterns of gene flow via the deliberate movement of genotypes around the environment are generating growing interest with broad applications in conservation and environmental management. These diverse strategies go by various names ranging from genetic or evolutionary rescue to provenancing and genetic resurrection. Our aim here is to provide some clarification around terminology and to how these strategies are connected and linked to underlying genetic processes. We draw on case studies from the literature and outline mechanisms that underlie how the various strategies aim to increase species fitness and impact the wider community. We argue that understanding mechanisms leading to species decline and community impact is a key to successful implementation of these strategies. We emphasize the need to consider the nature of source and recipient populations, as well as associated risks and trade-offs for the various strategies. This overview highlights where strategies are likely to have potential at population, species, and ecosystem scales, but also where they should probably not be attempted depending on the overall aims of the intervention. We advocate an approach where short- and long-term strategies are integrated into a decision framework that also considers nongenetic aspects of management.

Consistent community genetic effects in the context of strong environmental and temporal variation in Eucalyptus

Provenance translocations of tree species are promoted in forestry, conservation, and restoration in response to global climate change. While this option is driven by adaptive considerations, less is known of the effects translocations can have on dependent communities. We investigated the relative importance and consistency of extended genetic effects in Eucalyptus using two species-E. globulus and E. pauciflora. In E. globulus, the dependent arthropod and pathogen canopy communities were quantified based on the abundance of 49 symptoms from 722 progeny from 13 geographic sub-races across 2 common gardens. For E. pauciflora, 6 symptoms were quantified over 2 years from 238 progeny from 16 provenances across 2 common gardens. Genetic effects significantly influenced communities in both species. However, site and year effects outweighed genetic effects with site explaining approximately 3 times the variation in community traits in E. globulus and site and year explaining approximately 6 times the variation in E. pauciflora. While the genetic effect interaction terms were significant in some community traits, broad trends in community traits associated with variation in home-site latitude for E. globulus and home-site altitude for E. pauciflora were evident. These broad-scale trends were consistent with patterns of adaptive differentiation within each species, suggesting there may be extended consequences of local adaptation. While small in comparison to site and year, the consistency of genetic effects highlights the importance of provenance choice in tree species, such as Eucalyptus, as adaptive divergence among provenances may have significant long-term effects on biotic communities.

Inheritance of Rootstock Effects in Avocado (Persea americana Mill.) cv. Hass

Grafting is typically utilized to merge adapted seedling rootstocks with highly productive clonal scions. This process implies the interaction of multiple genomes to produce a unique tree phenotype. However, the interconnection of both genotypes obscures individual contributions to phenotypic variation (rootstock-mediated heritability), hampering tree breeding. Therefore, our goal was to quantify the inheritance of seedling rootstock effects on scion traits using avocado (Persea americana Mill.) cv. Hass as a model fruit tree. We characterized 240 diverse rootstocks from 8 avocado cv. Hass orchards with similar management in three regions of the province of Antioquia, northwest Andes of Colombia, using 13 microsatellite markers simple sequence repeats (SSRs). Parallel to this, we recorded 20 phenotypic traits (including morphological, biomass/reproductive, and fruit yield and quality traits) in the scions for 3 years (2015-2017). Relatedness among rootstocks was inferred through the genetic markers and inputted in a "genetic prediction" model to calculate narrow-sense heritabilities (h 2) on scion traits. We used three different randomization tests to highlight traits with consistently significant heritability estimates. This strategy allowed us to capture five traits with significant heritability values that ranged from 0.33 to 0.45 and model fits (r) that oscillated between 0.58 and 0.73 across orchards. The results showed significance in the rootstock effects for four complex harvest and quality traits (i.e., total number of fruits, number of fruits with exportation quality, and number of fruits discarded because of low weight or thrips damage), whereas the only morphological trait that had a significant heritability value was overall trunk height (an emergent property of the rootstock-scion interaction). These findings suggest the inheritance of rootstock effects, beyond root phenotype, on a surprisingly wide spectrum of scion traits in "Hass" avocado. They also reinforce the utility of polymorphic SSRs for relatedness reconstruction and genetic prediction of complex traits. This research is, up to date, the most cohesive evidence of narrow-sense inheritance of rootstock effects in a tropical fruit tree crop. Ultimately, our work highlights the importance of considering the rootstock-scion interaction to broaden the genetic basis of fruit tree breeding programs while enhancing our understanding of the consequences of grafting.

Modern Strategies to Assess and Breed Forest Tree Adaptation to Changing Climate

Studying the genetics of adaptation to new environments in ecologically and industrially important tree species is currently a major research line in the fields of plant science and genetic improvement for tolerance to abiotic stress. Specifically, exploring the genomic basis of local adaptation is imperative for assessing the conditions under which trees will successfully adapt in situ to global climate change. However, this knowledge has scarcely been used in conservation and forest tree improvement because woody perennials face major research limitations such as their outcrossing reproductive systems, long juvenile phase, and huge genome sizes. Therefore, in this review we discuss predictive genomic approaches that promise increasing adaptive selection accuracy and shortening generation intervals. They may also assist the detection of novel allelic variants from tree germplasm, and disclose the genomic potential of adaptation to different environments. For instance, natural populations of tree species invite using tools from the population genomics field to study the signatures of local adaptation. Conventional genetic markers and whole genome sequencing both help identifying genes and markers that diverge between local populations more than expected under neutrality, and that exhibit unique signatures of diversity indicative of "selective sweeps." Ultimately, these efforts inform the conservation and breeding status capable of pivoting forest health, ecosystem services, and sustainable production. Key long-term perspectives include understanding how trees' phylogeographic history may affect the adaptive relevant genetic variation available for adaptation to environmental change. Encouraging "big data" approaches (machine learning-ML) capable of comprehensively merging heterogeneous genomic and ecological datasets is becoming imperative, too.

Predicting Thermal Adaptation by Looking Into Populations' Genomic Past

Molecular evolution offers an insightful theory to interpret the genomic consequences of thermal adaptation to previous events of climate change beyond range shifts. However, disentangling often mixed footprints of selective and demographic processes from those due to lineage sorting, recombination rate variation, and genomic constrains is not trivial. Therefore, here we condense current and historical population genomic tools to study thermal adaptation and outline key developments (genomic prediction, machine learning) that might assist their utilization for improving forecasts of populations' responses to thermal variation. We start by summarizing how recent thermal-driven selective and demographic responses can be inferred by coalescent methods and in turn how quantitative genetic theory offers suitable multi-trait predictions over a few generations via the breeder's equation. We later assume that enough generations have passed as to display genomic signatures of divergent selection to thermal variation and describe how these footprints can be reconstructed using genome-wide association and selection scans or, alternatively, may be used for forward prediction over multiple generations under an infinitesimal genomic prediction model. Finally, we move deeper in time to comprehend the genomic consequences of thermal shifts at an evolutionary time scale by relying on phylogeographic approaches that allow for reticulate evolution and ecological parapatric speciation, and end by envisioning the potential of modern machine learning techniques to better inform long-term predictions. We conclude that foreseeing future thermal adaptive responses requires bridging the multiple spatial scales of historical and predictive environmental change research under modern cohesive approaches such as genomic prediction and machine learning frameworks.

High-throughput drone-based remote sensing reliably tracks phenology in thousands of conifer seedlings

Phenology is an important indicator of environmental variation and climate change impacts on tree responses. In conifers, monitoring phenology of photosynthesis through remote sensing has been unreliable, because needle foliage varies little throughout the year. This is challenging for modelling ecosystem carbon uptake and monitoring phenology for enhanced breeding (genomic selection) and forest health. Here, we demonstrate that drone-based carotenoid-sensitive spectral indices, such as the Chl/carotenoid index (CCI), can be used to track phenology in conifers by taking advantage of the close relationship between seasonally changing carotenoid levels and the variation of photosynthetic activity. Physiological ground measurements, including photosynthetic pigments and maximum quantum yield of Chl fluorescence, indicated that CCI tracked the variation of photosynthetic activity better than other vegetation indices for 30 white spruce seedlings measured over 1 yr. A machine-learning approach, using CCI derived from drone-based multispectral imagery, was used to model phenology of photosynthesis for the entire pedigree population (6000 seedlings). This high-throughput drone-based phenotyping approach is suitable for studying climate change impacts and environmental variation on the physiological status of thousands of field-grown conifers at unprecedented speed and scale.