Fine dissection of limber pine resistance to Cronartium ribicola using targeted sequencing of the NLR family

Genome-Wide Association Study Reveals Polygenic Architecture for Limber Pine Quantitative Disease Resistance to White Pine Blister Rust

Development of durable resistance effective against a broad range of pathotypes is crucial for restoration of pathogen-damaged ecosystems. This study dissected the complex genetic architecture for limber pine quantitative disease resistance (QDR) to Cronartium ribicola using a genome-wide association study. Eighteen-month-old seedlings were inoculated for resistance screening under controlled conditions. Disease development was quantitatively assessed for QDR-related traits over 4 years postinoculation. To reveal the genomic architecture contributing to QDR-related traits, a set of genes related to disease resistance with genome-wide distribution was selected for targeted sequencing for genotyping of single-nucleotide polymorphisms (SNPs). The genome-wide association study revealed a set of SNPs significantly associated with quantitative traits for limber pine QDR to white pine blister rust, including number of needle spots and stem cankers, as well as survival 4 years postinoculation. The peaks of marker-trait associations displayed a polygenic pattern, with genomic regions as potential resistant quantitative trait loci, distributed over 10 of the 12 linkage groups (LGs) of Pinus. None of them was linked to the Cr4-controlled major gene resistance previously mapped on LG08. Both normal canker and bole infection were mapped on LG05, and the associated SNPs explained their phenotypic variance up to 52%, tagging a major resistant quantitative trait locus. Candidate genes containing phenotypically associated SNPs encoded putative nucleotide-binding site leucine-rich repeat proteins, leucine-rich repeat-receptor-like kinase, cytochrome P450 superfamily protein, heat shock cognate protein 70, glutamate receptor, RNA-binding family protein, and unknown protein. The confirmation of resistant quantitative trait loci broadens the genetic pool of limber pine resistance germplasm for resistance breeding.

Maritime Pine Rootstock Genotype Modulates Gene Expression Associated with Stress Tolerance in Grafted Stems

Climate change-induced hazards, such as drought, threaten forest resilience, particularly in vulnerable regions such as the Mediterranean Basin. Maritime pine (Pinus pinaster Aiton), a model species in Western Europe, plays a crucial role in the Mediterranean forest due to its genetic diversity and ecological plasticity. This study characterizes transcriptional profiles of scion and rootstock stems of four P. pinaster graft combinations grown under well-watered conditions. Our grafting scheme combined drought-sensitive and drought-tolerant genotypes for scions (GAL1056: drought-sensitive scion; and Oria6: drought-tolerant scion) and rootstocks (R1S: drought-sensitive rootstock; and R18T: drought-tolerant rootstock). Transcriptomic analysis revealed expression patterns shaped by genotype provenance and graft combination. The accumulation of differentially expressed genes (DEGs) encoding proteins, involved in defense mechanisms and pathogen recognition, was higher in drought-sensitive scion stems and also increased when grafted onto drought-sensitive rootstocks. DEGs involved in drought tolerance mechanisms were identified in drought-tolerant genotypes as well as in drought-sensitive scions grafted onto drought-tolerant rootstocks, suggesting their establishment prior to drought. These mechanisms were associated with ABA metabolism and signaling. They were also involved in the activation of the ROS-scavenging pathways, which included the regulation of flavonoid and terpenoid metabolisms. Our results reveal DEGs potentially associated with the conifer response to drought and point out differences in drought tolerance strategies. These findings suggest genetic trade-offs between pine growth and defense, which could be relevant in selecting more drought-tolerant Pinus pinaster trees.

Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome

Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

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.

A Genome Sequence for the Threatened Whitebark Pine

Whitebark pine (WBP, Pinus albicaulis ) is a white pine of subalpine regions in western contiguous US 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 megametophyte 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 Gbp of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gbp). Approximately 87.2% (24.0 Gbp) of total sequence was placed on the twelve 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 three subclasses of NLRs (TNL, CNL, RNL). Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo assembled transcriptomes.

Beyond the genomes of Fulvia fulva (syn. Cladosporium fulvum) and Dothistroma septosporum: New insights into how these fungal pathogens interact with their host plants

Fulvia fulva and Dothistroma septosporum are closely related apoplastic pathogens with similar lifestyles but different hosts: F. fulva is a pathogen of tomato, whilst D. septosporum is a pathogen of pine trees. In 2012, the first genome sequences of these pathogens were published, with F. fulva and D. septosporum having highly fragmented and near-complete assemblies, respectively. Since then, significant advances have been made in unravelling their genome architectures. For instance, the genome of F. fulva has now been assembled into 14 chromosomes, 13 of which have synteny with the 14 chromosomes of D. septosporum, suggesting these pathogens are even more closely related than originally thought. Considerable advances have also been made in the identification and functional characterization of virulence factors (e.g., effector proteins and secondary metabolites) from these pathogens, thereby providing new insights into how they promote host colonization or activate plant defence responses. For example, it has now been established that effector proteins from both F. fulva and D. septosporum interact with cell-surface immune receptors and co-receptors to activate the plant immune system. Progress has also been made in understanding how F. fulva and D. septosporum have evolved with their host plants, whilst intensive research into pandemics of Dothistroma needle blight in the Northern Hemisphere has shed light on the origins, migration, and genetic diversity of the global D. septosporum population. In this review, we specifically summarize advances made in our understanding of the F. fulva-tomato and D. septosporum-pine pathosystems over the last 10 years.

Transcriptomic Reprogramming and Genetic Variations Contribute to Western Hemlock Defense and Resistance Against Annosus Root and Butt Rot Disease

Western hemlock (Tsuga heterophylla) is highly susceptible to Annosus root and butt rot disease, caused by Heterobasidion occidentale across its native range in western North America. Understanding molecular mechanisms of tree defense and dissecting genetic components underlying disease resistance will facilitate forest breeding and disease control management. The aim of this study was to profile host transcriptome reprogramming in response to pathogen infection using RNA-seq analysis. Inoculated seedlings were clearly grouped into three types: quantitative resistant (QR), susceptible (Sus), and un-infected (Uif), based on profiles of H. occidentale genes expressed in host tissues. Following de novo assembly of a western hemlock reference transcriptome with more than 33,000 expressed genes, the defensive transcriptome reprogramming was characterized and a set of differentially expressed genes (DEGs) were identified with gene ontology (GO) annotation. The QR seedlings showed controlled and coordinated molecular defenses against biotic stressors with enhanced biosynthesis of terpenoids, cinnamic acids, and other secondary metabolites. The Sus seedlings showed defense responses to abiotic stimuli with a few biological processes enhanced (such as DNA replication and cell wall organization), while others were suppressed (such as killing of cells of other organism). Furthermore, non-synonymous single nucleotide polymorphisms (ns-SNPs) of the defense- and resistance-related genes were characterized with high genetic variability. Both phylogenetic analysis and principal coordinate analysis (PCoA) revealed distinct evolutionary distances among the samples. The QR and Sus seedlings were well separated and grouped into different phylogenetic clades. This study provides initial insight into molecular defense and genetic components of western hemlock resistance against the Annosus root and butt rot disease. Identification of a large number of genes and their DNA variations with annotated functions in plant resistance and defense promotes the development of genomics-based breeding strategies for improved western hemlock resistance to H. occidentale.

Comparative Association Mapping Reveals Conservation of Major Gene Resistance to White Pine Blister Rust in Southwestern White Pine (Pinus strobiformis) and Limber Pine (P. flexilis)

All native North American white pines are highly susceptible to white pine blister rust (WPBR) caused by Cronartium ribicola. Understanding genomic diversity and molecular mechanisms underlying genetic resistance to WPBR remains one of the great challenges in improvement of white pines. To compare major gene resistance (MGR) present in two species, southwestern white pine (Pinus strobiformis) Cr3 and limber pine (P. flexilis) Cr4, we performed association analyses of Cr3-controlled resistant traits using single nucleotide polymorphism (SNP) assays designed with Cr4-linked polymorphic genes. We found that ∼70% of P. flexilis SNPs were transferable to P. strobiformis. Furthermore, several Cr4-linked SNPs were significantly associated with the Cr3-controlled traits in P. strobiformis families. The most significantly associated SNP (M326511_1126R) almost colocalized with Cr4 on the Pinus consensus linkage group 8, suggesting that Cr3 and Cr4 might be the same R locus, or have localizations very close to each other in the syntenic region of the P. strobiformis and P. flexilis genomes. M326511_1126R was identified as a nonsynonymous SNP, causing amino acid change (Val₃₇₆Ile) in a putative pectin acetylesterase, with coding sequences identical between the two species. Moreover, top Cr3-associated SNPs were further developed as TaqMan genotyping assays, suggesting their usefulness as marker-assisted selection (MAS) tools to distinguish genotypes between quantitative resistance and MGR. This work demonstrates the successful transferability of SNP markers between two closely related white pine species in the hybrid zone, and the possibility for deployment of MAS tools to facilitate long-term WPBR management in P. strobiformis breeding and conservation.