A genome-guided strategy for climate resilience in American chestnut restoration populations.

American chestnut (Castanea dentata) is a deciduous tree species of eastern North America that was decimated by the introduction of the chestnut blight fungus (Cryphonectria parasitica) in the early 20th century. Although millions of American chestnuts survive as root collar sprouts, these trees rarely reproduce. Thus, the species is considered functionally extinct. American chestnuts with improved blight resistance have been developed through interspecific hybridization followed by conspecific backcrossing, and by genetic engineering. Incorporating adaptive genomic diversity into these backcross families and transgenic lines is important for restoring the species across broad climatic gradients. To develop sampling recommendations for ex situ conservation of wild adaptive genetic variation, we coupled whole-genome resequencing of 384 stump sprouts with genotype-environment association analyses and found that the species range can be subdivided into three seed zones characterized by relatively homogeneous adaptive allele frequencies. We estimated that 21 to 29 trees per seed zone will need to be conserved to capture most extant adaptive diversity. We also resequenced the genomes of 269 backcross trees to understand the extent to which the breeding program has already captured wild adaptive diversity, and to estimate optimal reintroduction sites for specific families on the basis of their adaptive portfolio and future climate projections. Taken together, these results inform the development of an ex situ germplasm conservation and breeding plan to target blight-resistant breeding populations to specific environments and provides a blueprint for developing restoration plans for other imperiled tree species.

Validation of an Alternative Small Stem Assay for Blight Resistance in Chestnut Seedlings and Recommendations for Broader Use.

We evaluated an alternative small stem assay (AltSSA) for blight resistance in backcross hybrid chestnut trees (Castanea dentata/mollissima). Whereas standard small stem assays (SSAs) are done by inoculating small incisions in stems, in our AltSSA, 4- to 5-mm stems are cut off, and the exposed (living) stem tips are inoculated with discs of Cryphonectria parasitica inoculum and temporarily covered with plastic sleeves. Intended primarily for forward selection, this method was designed to be easy to implement, to consistently induce cankering, and to better enable seedling recovery via the development of lateral shoots from the lower stem. After 90+ days, cankers are evaluated and removed, and seedlings are prepared for out-planting. Previous results showed that AltSSAs performed at least as well as a common SSA method in distinguishing resistant and susceptible types. In this follow-up analysis of 35 lines of backcross seedlings studied in 2020 and 2021, we showed that mean orange zone canker length (OZCL) and a multifactor principal components analysis-based blight resistance index gave results consistent with predictions derived from two methods of blight resistance phenotyping and percentage of American chestnut ancestry of the parents of each line. As expected, based upon the apparent polygenic inheritance of blight resistance in backcross chestnut trees, mean OZCL of backcross families ranged from intermediate (F1 hybrid-level) to low (wild-type American chestnut-level). Consistent with prior results, canker production was near 100%, survivorship after out-planting was very high, and postinoculation stem dieback was not apparently related to the stem tip inoculations. Altogether, these results suggest that the AltSSA is a viable method for early detection of relative blight resistance in seedlings and may enable a reduction in the numbers of trees out-planted and placed under care for long-term evaluation and breeding. Thus, the AltSSA can prevent time, resources, and orchard space from being used on susceptible trees.

Frozen in time: Rangewide genomic diversity, structure, and demographic history of relict American chestnut populations.

American chestnut (Castanea dentata) was once the most economically and ecologically important hardwood species in the eastern United States. In the first half of the 20th century, an exotic fungal pathogen-Cryphonectria parasitica-decimated the species, killing billions of chestnut trees. Two approaches to developing blight-resistant American chestnut populations show promise, but both will require introduction of adaptive genomic diversity from wild germplasm to produce diverse, locally adapted restoration populations. Here we characterize population structure, demographic history, and genomic diversity in a range-wide sample of 384 wild American chestnuts to inform conservation and breeding with blight-resistant varieties. Population structure analyses suggest that the chestnut range can be roughly divided into northeast, central, and southwest populations. Within-population genomic diversity estimates revealed a clinal pattern with the highest diversity in the southwest, which likely reflects bottleneck events associated with Quaternary glaciation. Finally, we identified genomic regions under positive selection within each population, which suggests that defence against fungal pathogens is a common target of selection across all populations. Taken together, these results show that American chestnut underwent a postglacial expansion from the southern portion of its range leading to three extant genetic populations. These populations will serve as management units for breeding adaptive genetic variation into the blight-resistant tree populations for targeted reintroduction efforts.

Optimizing genomic selection for blight resistance in American chestnut backcross populations: A trade-off with American chestnut ancestry implies resistance is polygenic.

American chestnut was once a foundation species of eastern North American forests, but was rendered functionally extinct in the early 20th century by an exotic fungal blight (Cryphonectria parasitica). Over the past 30 years, the American Chestnut Foundation (TACF) has pursued backcross breeding to generate hybrids that combine the timber-type form of American chestnut with the blight resistance of Chinese chestnut based on a hypothesis of major gene resistance. To accelerate selection within two backcross populations that descended from two Chinese chestnuts, we developed genomic prediction models for five presence/absence blight phenotypes of 1,230 BC3F2 selection candidates and average canker severity of their BC3F3 progeny. We also genotyped pure Chinese and American chestnut reference panels to estimate the proportion of BC3F2 genomes inherited from parent species. We found that genomic prediction from a method that assumes an infinitesimal model of inheritance (HBLUP) has similar accuracy to a method that tends to perform well for traits controlled by major genes (Bayes C). Furthermore, the proportion of BC3F2 trees' genomes inherited from American chestnut was negatively correlated with the blight resistance of these trees and their progeny. On average, selected BC3F2 trees inherited 83% of their genome from American chestnut and have blight resistance that is intermediate between F1 hybrids and American chestnut. Results suggest polygenic inheritance of blight resistance. The blight resistance of restoration populations will be enhanced through recurrent selection, by advancing additional sources of resistance through fewer backcross generations, and by potentially by breeding with transgenic blight-tolerant trees.

Resistance to Phytophthora cinnamomi in American Chestnut (Castanea dentata) Backcross Populations that Descended from Two Chinese Chestnut (Castanea mollissima) Sources of Resistance.

Restoration of American chestnut (Castanea dentata) depends on combining resistance to both the chestnut blight fungus (Cryphonectria parasitica) and Phytophthora cinnamomi, which causes Phytophthora root rot, in a diverse population of C. dentata. Over a 14-year period (2004 to 2017), survival and root health of American chestnut backcross seedlings after inoculation with P. cinnamomi were compared among 28 BC3, 66 BC4, and 389 BC3F3 families that descended from two BC1 trees (Clapper and Graves) with different Chinese chestnut grandparents. The 5% most resistant Graves BC3F3 families survived P. cinnamomi infection at rates of 75 to 100% but had mean root health scores that were intermediate between resistant Chinese chestnut and susceptible American chestnut families. Within Graves BC3F3 families, seedling survival was greater than survival of Graves BC3 and BC4 families and was not genetically correlated with chestnut blight canker severity. Only low to intermediate resistance to P. cinnamomi was detected among backcross descendants from the Clapper tree. Results suggest that major-effect resistance alleles were inherited by descendants from the Graves tree, that intercrossing backcross trees enhances progeny resistance to P. cinnamomi, and that alleles for resistance to P. cinnamomi and C. parasitica are not linked. To combine resistance to both C. parasitica and P. cinnamomi, a diverse Graves backcross population will be screened for resistance to P. cinnamomi, survivors bred with trees selected for resistance to C. parasitica, and progeny selected for resistance to both pathogens will be intercrossed.

Hardwood Tree Genomics: Unlocking Woody Plant Biology.

Woody perennial angiosperms (i.e., hardwood trees) are polyphyletic in origin and occur in most angiosperm orders. Despite their independent origins, hardwoods have shared physiological, anatomical, and life history traits distinct from their herbaceous relatives. New high-throughput DNA sequencing platforms have provided access to numerous woody plant genomes beyond the early reference genomes of Populus and Eucalyptus, references that now include willow and oak, with pecan and chestnut soon to follow. Genomic studies within these diverse and undomesticated species have successfully linked genes to ecological, physiological, and developmental traits directly. Moreover, comparative genomic approaches are providing insights into speciation events while large-scale DNA resequencing of native collections is identifying population-level genetic diversity responsible for variation in key woody plant biology across and within species. Current research is focused on developing genomic prediction models for breeding, defining speciation and local adaptation, detecting and characterizing somatic mutations, revealing the mechanisms of gender determination and flowering, and application of systems biology approaches to model complex regulatory networks underlying quantitative traits. Emerging technologies such as single-molecule, long-read sequencing is being employed as additional woody plant species, and genotypes within species, are sequenced, thus enabling a comparative ("evo-devo") approach to understanding the unique biology of large woody plants. Resource availability, current genomic and genetic applications, new discoveries and predicted future developments are illustrated and discussed for poplar, eucalyptus, willow, oak, chestnut, and pecan.

Leaf cellulose density as the key determinant of inter- and intra-specific variation in leaf fracture toughness in a species-rich tropical forest.

Leaves as the main photosynthetic organ of plants must be well protected against various hazards to achieve their optimal lifespans. Yet, within-species variation and the material basis of leaf strength have been explored for very few species. Here, we present a large dataset of leaf fracture toughness from a species-rich humid tropical forest on Barro Colorado Island, Panama, reporting both among- and within-species variation in relation to light environment (sun-lit canopy versus shaded understorey) and ontogeny (seedlings versus adults). In this dataset encompassing 281 free-standing woody species and 428 species-light combinations, lamina fracture toughness varied ca 10 times. A central objective of our study was to identify generalizable patterns in the structural and material basis for interspecific variation in leaf lamina fracture toughness. The leaf lamina is a heterogeneous structure in which strong materials in cell walls, such as cellulose and lignin, contribute disproportionately to fracture toughness. We found significant increases in leaf fracture toughness from shade to sun and from seedling leaves to adult leaves. Both within and across species, leaf fracture toughness increased with total bulk density (dry biomass per unit volume) and cellulose mass concentration, but decreased with mass concentrations of lignin and hemicelluose. These bivariate relationships shift between light environments, but leaf cellulose density (cellulose mass per unit leaf volume) exhibits a common relationship with lamina fracture toughness between light environments and through ontogeny. Hence, leaf cellulose density is probably a universal predictor of leaf fracture toughness.

A Consensus Genetic Map for Pinus taeda and Pinus elliottii and Extent of Linkage Disequilibrium in Two Genotype-Phenotype Discovery Populations of Pinus taeda.

A consensus genetic map for Pinus taeda (loblolly pine) and Pinus elliottii (slash pine) was constructed by merging three previously published P. taeda maps with a map from a pseudo-backcross between P. elliottii and P. taeda. The consensus map positioned 3856 markers via genotyping of 1251 individuals from four pedigrees. It is the densest linkage map for a conifer to date. Average marker spacing was 0.6 cM and total map length was 2305 cM. Functional predictions of mapped genes were improved by aligning expressed sequence tags used for marker discovery to full-length P. taeda transcripts. Alignments to the P. taeda genome mapped 3305 scaffold sequences onto 12 linkage groups. The consensus genetic map was used to compare the genome-wide linkage disequilibrium in a population of distantly related P. taeda individuals (ADEPT2) used for association genetic studies and a multiple-family pedigree used for genomic selection (CCLONES). The prevalence and extent of LD was greater in CCLONES as compared to ADEPT2; however, extended LD with LGs or between LGs was rare in both populations. The average squared correlations, r(2), between SNP alleles less than 1 cM apart were less than 0.05 in both populations and r(2) did not decay substantially with genetic distance. The consensus map and analysis of linkage disequilibrium establish a foundation for comparative association mapping and genomic selection in P. taeda and P. elliottii.

Discovering candidate genes that regulate resin canal number in Pinus taeda stems by integrating genetic analysis across environments, ages, and populations.

Genetically improving constitutive resin canal development in Pinus stems may enhance the capacity to synthesize terpenes for bark beetle resistance, chemical feedstocks, and biofuels. To discover genes that potentially regulate axial resin canal number (RCN), single nucleotide polymorphisms (SNPs) in 4027 genes were tested for association with RCN in two growth rings and three environments in a complex pedigree of 520 Pinus taeda individuals (CCLONES). The map locations of associated genes were compared with RCN quantitative trait loci (QTLs) in a (P. taeda × Pinus elliottii) × P. elliottii pseudo-backcross of 345 full-sibs (BC1). Resin canal number was heritable (h(2) ˜ 0.12-0.21) and positively genetically correlated with xylem growth (rg ˜ 0.32-0.72) and oleoresin flow (rg ˜ 0.15-0.51). Sixteen well-supported candidate regulators of RCN were discovered in CCLONES, including genes associated across sites and ages, unidirectionally associated with oleoresin flow and xylem growth, and mapped to RCN QTLs in BC1. Breeding is predicted to increase RCN 11% in one generation and could be accelerated with genomic selection at accuracies of 0.45-0.52 across environments. There is significant genetic variation for RCN in loblolly pine, which can be exploited in breeding for elevated terpene content.

Association genetics of oleoresin flow in loblolly pine: discovering genes and predicting phenotype for improved resistance to bark beetles and bioenergy potential.

Rapidly enhancing oleoresin production in conifer stems through genomic selection and genetic engineering may increase resistance to bark beetles and terpenoid yield for liquid biofuels. We integrated association genetic and genomic prediction analyses of oleoresin flow (g 24 h(-1)) using 4854 single nucleotide polymorphisms (SNPs) in expressed genes within a pedigreed population of loblolly pine (Pinus taeda) that was clonally replicated at three sites in the southeastern United States. Additive genetic variation in oleoresin flow (h(2) ≈ 0.12-0.30) was strongly correlated between years in which precipitation varied (r(a) ≈ 0.95), while the genetic correlation between sites declined from 0.8 to 0.37 with increasing differences in soil and climate among sites. A total of 231 SNPs were significantly associated with oleoresin flow, of which 81% were specific to individual sites. SNPs in sequences similar to ethylene signaling proteins, ABC transporters, and diterpenoid hydroxylases were associated with oleoresin flow across sites. Despite this complex genetic architecture, we developed a genomic prediction model to accelerate breeding for enhanced oleoresin flow that is robust to environmental variation. Results imply that breeding could increase oleoresin flow 1.5- to 2.4-fold in one generation.

What makes a leaf tough? Patterns of correlated evolution between leaf toughness traits and demographic rates among 197 shade-tolerant woody species in a neotropical forest.

Slow-growing juveniles of shade-tolerant plant species are predicted to have tough leaves because of the high cost of leaf replacement in shade relative to potential carbon gain. We assessed the degree of correlated evolution among eight traits associated with leaf toughness and the relationships of those traits with the growth and mortality rates of 197 tree and shrub species from the understory of the 50-ha forest dynamics plot on Barro Colorado Island, Panama. Path analysis with phylogenetically independent contrasts revealed that leaves attained material toughness (resistance to fracture per unit fracture area) through increases in tissue density, percent cellulose per unit dry mass, and vein fracture toughness. Lamina density and cellulose content evolved independently and thus represent different paths to material toughness. Structural toughness (resistance to fracture per unit fracture length) depended on material toughness and lamina thickness. Mortality rates of individuals 1-10 cm in stem diameter were negatively correlated with material toughness and lamina density but were independent of structural toughness and cell wall fiber contents. Leaf toughness traits were uncorrelated with relative growth rates. Results imply that material toughness enhances resistance to natural enemies, which increases survival and offsets the biomass allocation cost of producing tough leaves in the shaded understory.