Transcriptional Profiling of Early Defense Response to White Pine Blister Rust Infection in Pinus albicaulis (Whitebark Pine).

Pathogen perception generates the activation of signal transduction cascades to host defense. White pine blister rust (WPBR) is caused by Cronartium ribicola J.C. Fisch and affects a number of species of Pinus. One of the most severely affected species is Pinus albicaulis Engelm (whitebark pine). WPBR resistance in the species is a polygenic and complex trait that requires an optimized immune response. We identified early responses in 2-year-old seedlings after four days of fungal inoculation and compared the underlying transcriptomic response with that of healthy non-inoculated individuals. A de novo transcriptome assembly was constructed with 56,796 high quality-annotations derived from the needles of susceptible and resistant individuals in a resistant half-sib family. Differential expression analysis identified 599 differentially expressed transcripts, from which 375 were upregulated and 224 were downregulated in the inoculated seedlings. These included components of the initial phase of active responses to abiotic factors and stress regulators, such as those involved in the first steps of flavonoid biosynthesis. Four days after the inoculation, infected individuals showed an overexpression of chitinases, reactive oxygen species (ROS) regulation signaling, and flavonoid intermediates. Our research sheds light on the first stage of infection and emergence of disease symptoms among whitebark pine seedlings. RNA sequencing (RNA-seq) data encoding hypersensitive response, cell wall modification, oxidative regulation signaling, programmed cell death, and plant innate immunity were differentially expressed during the defense response against C. ribicola.

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.

The Global Forest Health Crisis: A Public-Good Social Dilemma in Need of International Collective Action.

Society is confronted by interconnected threats to ecological sustainability. Among these is the devastation of forests by destructive non-native pathogens and insects introduced through global trade, leading to the loss of critical ecosystem services and a global forest health crisis. We argue that the forest health crisis is a public-good social dilemma and propose a response framework that incorporates principles of collective action. This framework enables scientists to better engage policymakers and empowers the public to advocate for proactive biosecurity and forest health management. Collective action in forest health features broadly inclusive stakeholder engagement to build trust and set goals; accountability for destructive pest introductions; pooled support for weakest-link partners; and inclusion of intrinsic and nonmarket values of forest ecosystems in risk assessment. We provide short-term and longer-term measures that incorporate the above principles to shift the societal and ecological forest health paradigm to a more resilient state.

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 (Val376Ile) 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.

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

BACKGROUND: Proteins with nucleotide binding site (NBS) and leucine-rich repeat (LRR) domains (NLR) make up one of most important resistance (R) families for plants to resist attacks from various pathogens and pests. The available transcriptomes of limber pine (Pinus flexilis) allow us to characterize NLR genes and related resistance gene analogs (RGAs) in host resistance against Cronartium ribicola, the causal fungal pathogen of white pine blister rust (WPBR) on five-needle pines throughout the world. We previously mapped a limber pine major gene locus (Cr4) that confers complete resistance to C. ribicola on the Pinus consensus linkage group 8 (LG-8). However, genetic distribution of NLR genes as well as their divergence between resistant and susceptible alleles are still unknown. RESULTS: To identify NLR genes at the Cr4 locus, the present study re-sequenced a total of 480 RGAs using targeted sequencing in a Cr4-segregated seed family. Following a call of single nucleotide polymorphisms (SNPs) and genetic mapping, a total of 541 SNPs from 155 genes were mapped across 12 LGs. Three putative NLR genes were newly mapped in the Cr4 region, including one that co-segregated with Cr4. The tight linkage of NLRs with Cr4-controlled phenotypes was further confirmed by bulked segregation analysis (BSA) using extreme-phenotype genome-wide association study (XP-GWAS) for significance test. Local tandem duplication in the Cr4 region was further supported by syntenic analysis using the sugar pine genome sequence. Significant gene divergences have been observed in the NLR family, revealing that diversifying selection pressures are relatively higher in local duplicated genes. Most genes showed similar expression patterns at low levels, but some were affected by genetic background related to disease resistance. Evidence from fine genetic dissection, evolutionary analysis, and expression profiling suggests that two NLR genes are the most promising candidates for Cr4 against WPBR. CONCLUSION: This study provides fundamental insights into genetic architecture of the Cr4 locus as well as a set of NLR variants for marker-assisted selection in limber pine breeding. Novel NLR genes were identified at the Cr4 locus and the Cr4 candidates will aid deployment of this R gene in combination with other major/minor genes in the limber pine breeding program.

Comparative Transcriptomics and RNA-Seq-Based Bulked Segregant Analysis Reveals Genomic Basis Underlying Cronartium ribicola vcr2 Virulence.

Breeding programs of five-needle pines have documented both major gene resistance (MGR) and quantitative disease resistance (QDR) to Cronartium ribicola (Cri), a non-native, invasive fungal pathogen causing white pine blister rust (WPBR). WPBR is one of the most deadly forest diseases in North America. However, Cri virulent pathotypes have evolved and can successfully infect and kill trees carrying resistance (R) genes, including vcr2 that overcomes MGR conferred by the western white pine (WWP, Pinus monticola) R gene (Cr2). In the absence of a reference genome, the present study generated a vcr2 reference transcriptome, consisting of about 20,000 transcripts with 1,014 being predicted to encode secreted proteins (SPs). Comparative profiling of transcriptomes and secretomes revealed vcr2 was significantly enriched for several gene ontology (GO) terms relating to oxidation-reduction processes and detoxification, suggesting that multiple molecular mechanisms contribute to pathogenicity of the vcr2 pathotype for its overcoming Cr2. RNA-seq-based bulked segregant analysis (BSR-Seq) revealed genome-wide DNA variations, including about 65,617 single nucleotide polymorphism (SNP) loci in 7,749 polymorphic genes shared by vcr2 and avirulent (Avcr2) pathotypes. An examination of the distribution of minor allele frequency (MAF) uncovered a high level of genomic divergence between vcr2 and Avcr2 pathotypes. By integration of extreme-phenotypic genome-wide association (XP-GWAS) analysis and allele frequency directional difference (AFDD) mapping, we identified a set of vcr2-associated SNPs within functional genes, involved in fungal virulence and other molecular functions. These included six SPs that were top candidate effectors with putative activities of reticuline oxidase, proteins with common in several fungal extracellular membrane (CFEM) domain or ferritin-like domain, polysaccharide lyase, rds1p-like stress responsive protein, and two Cri-specific proteins without annotation. Candidate effectors and vcr2-associated genes provide valuable resources for further deciphering molecular mechanisms of virulence and pathogenicity by functional analysis and the subsequent development of diagnostic tools for monitoring the virulence landscape in the WPBR pathosystems.

In-vitro anti-fungal assay and association analysis reveal a role for the Pinus monticola PR10 gene (PmPR10-3.1) in quantitative disease resistance to white pine blister rust.

Pathogenesis-related (PR) proteins play important roles in plant defense response. However, functional investigation of PR10 genes is still limited and their physiological roles have not been conclusively characterized in biological processes of conifer trees. Here, we identified multiple novel members in the western white pine (Pinus monticola) PmPR10 family by bioinformatic mining available transcriptomic data. Phylogenetic analysis of protein sequences revealed four PR10 and two PR10-like clusters with a high synteny across different species of five-needle pines. Of 10 PmPR10 genes, PmPR10-3.1 was selected and expressed in Escherichia coli. The purified recombinant protein exhibited inhibitory effects on spore hyphal growth of fungal pathogens Cronartium ribicola, Phoma exigua, and Phoma argillacea by in-vitro anti-fungal analysis. Genetic variation analysis detected a total of 21 single nucleotide polymorphisms (SNPs) within PmPR10-3.1 in a collection of P. monticola seed families. A nonsynonymous SNP (t178g) showed significant association with relative levels of quantitative disease resistance (QDR), explaining about 8.7% of phenotypic variation as the peak value across all SNPs. Our results provide valuable insight into the genetic architecture underlying P. monticola QDR and imply that PmPR10-3.1 may function as an important component in conifer basal immunity for non-specific resistance to a wide spectrum of pathogens.

Association Mapping and Development of Marker-Assisted Selection Tools for the Resistance to White Pine Blister Rust in the Alberta Limber Pine Populations.

Since its introduction to North America in the early 1900s, white pine blister rust (WPBR) caused by the fungal pathogen Cronartium ribicola has resulted in substantial economic losses and ecological damage to native North American five-needle pine species. The high susceptibility and mortality of these species, including limber pine (Pinus flexilis), creates an urgent need for the development and deployment of resistant germplasm to support recovery of impacted populations. Extensive screening for genetic resistance to WPBR has been underway for decades in some species but has only started recently in limber pine using seed families collected from wild parental trees in the USA and Canada. This study was conducted to characterize Alberta limber pine seed families for WPBR resistance and to develop reliable molecular tools for marker-assisted selection (MAS). Open-pollinated seed families were evaluated for host reaction following controlled infection using C. ribicola basidiospores. Phenotypic segregation for presence/absence of stem symptoms was observed in four seed families. The segregation ratios of these families were consistent with expression of major gene resistance (MGR) controlled by a dominant R locus. Based on linkage disequilibrium (LD)-based association mapping used to detect single nucleotide polymorphism (SNP) markers associated with MGR against C. ribicola, MGR in these seed families appears to be controlled by Cr4 or other R genes in very close proximity to Cr4. These associated SNPs were located in genes involved in multiple molecular mechanisms potentially underlying limber pine MGR to C. ribicola, including NBS-LRR genes for recognition of C. ribicola effectors, signaling components, and a large set of defense-responsive genes with potential functions in plant effector-triggered immunity (ETI). Interactions of associated loci were identified for MGR selection in trees with complex genetic backgrounds. SNPs with tight Cr4-linkage were further converted to TaqMan assays to confirm their effectiveness as MAS tools. This work demonstrates the successful translation and deployment of molecular genetic knowledge into specific MAS tools that can be easily applied in a selection or breeding program to efficiently screen MGR against WPBR in Alberta limber pine populations.

Genomic basis of white pine blister rust quantitative disease resistance and its relationship with qualitative resistance.

The genomic architecture and molecular mechanisms controlling variation in quantitative disease resistance loci are not well understood in plant species and have been barely studied in long-generation trees. Quantitative trait loci mapping and genome-wide association studies were combined to test a large single nucleotide polymorphism (SNP) set for association with quantitative and qualitative white pine blister rust resistance in sugar pine. In the absence of a chromosome-scale reference genome, a high-density consensus linkage map was generated to obtain locations for associated SNPs. Newly discovered associations for white pine blister rust quantitative disease resistance included 453 SNPs involved in wide biological functions, including genes associated with disease resistance and others involved in morphological and developmental processes. In addition, NBS-LRR pathogen recognition genes were found to be involved in quantitative disease resistance, suggesting these newly reported genes are qualitative genes with partial resistance, they are the result of defeated qualitative resistance due to avirulent races, or they have epistatic effects on qualitative disease resistance genes. This study is a step forward in our understanding of the complex genomic architecture of quantitative disease resistance in long-generation trees, and constitutes the first step towards marker-assisted disease resistance breeding in white pine species.

Characterization of Cronartium ribicola dsRNAs reveals novel members of the family Totiviridae and viral association with fungal virulence.

BACKGROUND: Mycoviruses were recently discovered in the white pine blister rust (WPBR) fungus Cronartium ribicola (J.C. Fisch.). Detection and characterization of their double stranded RNA (dsRNA) would facilitate understanding of pathogen virulence and disease pathogenesis in WPBR systems. METHODS: Full-length cDNAs were cloned from the dsRNAs purified from viral-infected C. ribicola, and their cDNA sequences were determined by DNA sequencing. Evolutionary relationships of the dsRNAs with related mycoviruses were determined by phylogenetic analysis. Dynamic distributions of the viral RNAs within samples of their fungal host C. ribicola were investigated by measurement of viral genome prevalence and viral gene expression. RESULTS: In this study we identified and characterized five novel dsRNAs from C. ribicola, designated as Cronartium ribicola totivirus 1-5 (CrTV1 to CrTV5). These dsRNA sequences encode capsid protein and RNA-dependent RNA polymerase with significant homologies to dsRNA viruses of the family Totiviridae. Phylogenetic analysis showed that the CrTVs were grouped into two distinct clades. CrTV2 through CrTV5 clustered within the genus Totivirus. CrTV1 along with a few un-assigned dsRNAs constituted a distinct phyletic clade that is genetically distant from presently known genera in the Totiviridae family, indicating that CrTV1 represents a novel genus in the Totiviridae family. The CrTVs were prevalent in fungal samples obtained from infected western white pine, whitebark pine, and limber pines. Viral RNAs were generally expressed at higher levels during in planta mycelium growth than in aeciospores and urediniospores. CrTV4 was significantly associated with C. ribicola virulent pathotype and specific C. ribicola host tree species, suggesting dsRNAs as potential tools for dissection of pathogenic mechanisms of C. ribicola and diagnosis of C. ribicola pathotypes. CONCLUSION: Phylogenetic and expression analyses of viruses in the WPBR pathogen, C. ribicola, have enchanced our understanding of virus diversity in the family Totiviridae, and provided a potential strategy to utilize pathotype-associated mycoviruses to control fungal forest diseases.

Are traded forest tree seeds a potential source of nonnative pests?

The international seed trade is considered relatively safe from a phytosanitary point of view and is therefore less regulated than trade in other plants for planting. However, the pests carried by traded seeds are not well known. We assessed insects and fungi in 58 traded seed lots of 11 gymnosperm and angiosperm tree species from North America, Europe, and Asia. Insects were detected by X-raying and molecular methods. The fungal community was characterized using high-throughput sequencing (HTS) and by growing fungi on non-selective agar. About 30% of the seed lots contained insect larvae. Gymnosperms contained mostly hymenopteran (Megastigmus spp.) and dipteran (Cecidomyiidae) larvae, while angiosperms contained lepidopteran (Cydia latiferreana) and coleopteran (Curculio spp.) larvae. HTS indicated the presence of fungi in all seed lots and fungi grew on non-selective agar from 96% of the seed lots. Fungal abundance and diversity were much higher than insect diversity and abundance, especially in angiosperm seeds. Almost 50% of all fungal exact sequence variants (ESVs) found in angiosperms were potential pathogens, in comparison with around 30% of potentially pathogenic ESVs found in gymnosperms. The results of this study indicate that seeds may pose a greater risk of pest introduction than previously believed or accounted for. A rapid risk assessment suggests that only a small number of species identified in this study is of phytosanitary concern. However, more research is needed to enable better risk assessment, especially to increase knowledge about the potential for transmission of fungi to seedlings and the host range and impact of identified species.

Limber pine (Pinus flexilis James) genetic map constructed by exome-seq provides insight into the evolution of disease resistance and a genomic resource for genomics-based breeding.

Limber pine (Pinus flexilis) is a keystone species of high-elevation forest ecosystems of western North America, but some parts of the geographic range have high infection and mortality from the non-native white pine blister rust caused byCronartium ribicola. Genetic maps can provide essential knowledge for understanding genetic disease resistance as well as local adaptation to changing climates. Exome-seq was performed to construct high-density genetic maps in two seed families. Composite maps positioned 9612 unigenes across 12 linkage groups (LGs). Syntenic analysis of genome structure revealed that the majority of orthologs were positional orthologous genes (POGs) with localization on homologousLGs among conifer species. Gene ontology (GO) enrichment analysis showed relatively fewer constraints forPOGs with putative roles in adaptation to environments and relatively more conservation forPOGs with roles in basic cell function and maintenance. The mapped genes included 639 nucleotide-binding site leucine-rich repeat genes (NBS-LRRs), 290 receptor-like protein kinase genes (RLKs), and 1014 genes with potential roles in the defense response and induced systemic resistance to attack by pathogens. Orthologous loci for resistance to rust pathogens were identified and were co-positioned with multiple members of theR gene family, revealing the evolutionary pressure acting upon them. This high-density genetic map provides a genomic resource and practical tool for breeding and genetic conservation programs, with applications in genome-wide association studies (GWASs), the characterization of functional genes underlying complex traits, and the sequencing and assembly of the full-length genomes of limber pine and relatedPinus species.

The role of hybridization during ecological divergence of southwestern white pine (Pinus strobiformis) and limber pine (P. flexilis).

Interactions between extrinsic factors, such as disruptive selection and intrinsic factors, such as genetic incompatibilities among loci, often contribute to the maintenance of species boundaries. The relative roles of these factors in the establishment of reproductive isolation can be examined using species pairs characterized by gene flow throughout their divergence history. We investigated the process of speciation and the maintenance of species boundaries between Pinus strobiformis and Pinus flexilis. Utilizing ecological niche modelling, demographic modelling and genomic cline analyses, we illustrated a divergence history with continuous gene flow. Our results supported an abundance of advanced generation hybrids and a lack of loci exhibiting steep transition in allele frequency across the hybrid zone. Additionally, we found evidence for climate-associated variation in the hybrid index and niche divergence between parental species and the hybrid zone. These results are consistent with extrinsic factors, such as climate, being an important isolating mechanism. A build-up of intrinsic incompatibilities and of coadapted gene complexes is also apparent, although these appear to be in the earliest stages of development. This supports previous work in coniferous species demonstrating the importance of extrinsic factors in facilitating speciation. Overall, our findings lend support to the hypothesis that varying strength and direction of selection pressures across the long lifespans of conifers, in combination with their other life history traits, delays the evolution of strong intrinsic incompatibilities.

The influence of genetics, defensive chemistry and the fungal microbiome on disease outcome in whitebark pine trees.

The invasive fungal pathogen Cronartium ribicola infects and kills whitebark pine (Pinus albicaulis) throughout western North America. Whitebark pine has been proposed for listing under the Endangered Species Act in the USA, and the loss of this species is predicted to have severe impacts on ecosystem composition and function in high-elevation forests. Numerous fungal endophytes live inside whitebark pine tissues and may influence the severity of C. ribicola infection, either directly by inhibition of pathogen growth or indirectly by the induction of chemical defensive pathways in the tree. Terpenes, a form of chemical defence in pine trees, can also influence disease. In this study, we characterized fungal endophyte communities in whitebark pine seedlings before and after experimental inoculation with C. ribicola, monitored disease progression and compared fungal community composition in susceptible vs. resistant seedlings in a common garden. We analysed the terpene composition of these same seedlings. Seed family identity or maternal genetics influenced both terpenes and endophyte communities. Terpene and endophyte composition correlated with disease severity, and terpene concentrations differed in resistant vs. susceptible seedlings. These results suggest that the resistance to C. ribicola observed in natural whitebark pine populations is caused by the combined effects of genetics, endophytes and terpenes within needle tissue, in which initial interactions between microbes and hosts take place. Tree genotype, terpene and microbiome combinations associated with healthy trees could help to predict or reduce disease severity and improve outcomes of future tree breeding programmes.

Profiling methyl jasmonate-responsive transcriptome for understanding induced systemic resistance in whitebark pine (Pinus albicaulis).

KEY MESSAGE: RNA-seq analysis on whitebark pine needles demonstrated that methyl jasmonate (MeJA)-triggered transcriptome re-programming substantially overlapped with defense responses against insects and fungal pathogens in Pinus species, increasing current knowledge regarding induced systemic resistance (ISR) to pathogens and pests in whitebark pine. Many whitebark pine populations are in steep decline due to high susceptibility to mountain pine beetle and the non-native white pine blister rust (WPBR). Resistance, including induced systemic resistance (ISR), is not well characterized in whitebark pine, narrowing the current options for increasing the success of restoration and breeding programs. Exogenous jasmonates are known to trigger ISR by activating the plant's immune system through regulation of gene expression to produce chemical defense compounds. This study reports profiles of whitebark pine needle transcriptomes, following methyl jasmonate (MeJA) treatment using RNA-seq. A MeJA-responsive transcriptome was de novo assembled and transcriptome profiling identified a set of differentially expressed genes (DEGs), revealing 1422 up- and 999 down-regulated transcripts with at least twofold change (FDR corrected p < 0.05) in needle tissues in response to MeJA application. GO analysis revealed that these DEGs have putative functions in plant defense signalling, transcription regulation, biosyntheses of secondary metabolites, and other biological processes. Lineage-specific expression of defense-related genes was characterized through comparison with MeJA signalling in model plants. In particular, MeJA-triggered transcriptome re-programming substantially overlapped with defense responses against WPBR and insects in related Pinus species, suggesting that MeJA may be used to improve whitebark pine resistance to pathogens/pests. Our study provides new insights into molecular mechanisms and metabolic pathways involved in whitebark pine ISR. DEGs identified in this study can be used as candidates to facilitate identification of genomic variation contributing to host resistance and aid in breeding selection of elite genotypes with better adaptive fitness to environmental stressors in this endangered tree species.

Saturated genic SNP mapping identified functional candidates and selection tools for the Pinus monticola Cr2 locus controlling resistance to white pine blister rust.

Molecular breeding incorporates efficient tools to increase rust resistance in five-needle pines. Susceptibility of native five-needle pines to white pine blister rust (WPBR), caused by the non-native invasive fungus Cronartium ribicola (J.C. Fisch.), has significantly reduced wild populations of these conifers in North America. Major resistance (R) genes against specific avirulent pathotypes have been found in several five-needle pine species. In this study, we screened genic SNP markers by comparative transcriptome and genetic association analyses and constructed saturated linkage maps for the western white pine (Pinus monticola) R locus (Cr2). Phenotypic segregation was measured by a hypersensitive reaction (HR)-like response on the needles and disease symptoms of cankered stems post inoculation by the C. ribicola avcr2 race. SNP genotypes were determined by HRM- and TaqMan-based SNP genotyping. Saturated maps of the Cr2-linkage group (LG) were constructed in three seed families using a total of 34 SNP markers within 21 unique genes. Cr2 was consistently flanked by contig_2142 (encoding a ruvb-like protein) and contig_3772 (encoding a delta-fatty acid desaturase) across the three seed families. Cr2 was anchored to the Pinus consensus LG-1, which differs from LGs where other R loci of Pinus species were mapped. GO annotation identified a set of NBS-LRR and other resistance-related genes as R candidates in the Cr2 region. Association of one nonsynonymous SNP locus of an NBS-LRR gene with Cr2-mediated phenotypes provides a valuable tool for marker-assisted selection (MAS), which will shorten the breeding cycle of resistance screening and aid in the restoration of WPBR-disturbed forest ecosystems.

Genetic Diversity and Population Structure of Whitebark Pine (Pinus albicaulis Engelm.) in Western North America.

Whitebark pine (WBP, Pinus albicaulis Engelm.) is an endangered conifer species due to heavy mortality from white pine blister rust (WPBR, caused by Cronartium ribicola) and mountain pine beetle (Dendroctonus ponderosae). Information about genetic diversity and population structure is of fundamental importance for its conservation and restoration. However, current knowledge on the genetic constitution and genomic variation is still limited for WBP. In this study, an integrated genomics approach was applied to characterize seed collections from WBP breeding programs in western North America. RNA-seq analysis was used for de novo assembly of the WBP needle transcriptome, which contains 97,447 protein-coding transcripts. Within the transcriptome, single nucleotide polymorphisms (SNPs) were discovered, and more than 22,000 of them were non-synonymous SNPs (ns-SNPs). Following the annotation of genes with ns-SNPs, 216 ns-SNPs within candidate genes with putative functions in disease resistance and plant defense were selected to design SNP arrays for high-throughput genotyping. Among these SNP loci, 71 were highly polymorphic, with sufficient variation to identify a unique genotype for each of the 371 individuals originating from British Columbia (Canada), Oregon and Washington (USA). A clear genetic differentiation was evident among seed families. Analyses of genetic spatial patterns revealed varying degrees of diversity and the existence of several genetic subgroups in the WBP breeding populations. Genetic components were associated with geographic variables and phenotypic rating of WPBR disease severity across landscapes, which may facilitate further identification of WBP genotypes and gene alleles contributing to local adaptation and quantitative resistance to WPBR. The WBP genomic resources developed here provide an invaluable tool for further studies and for exploitation and utilization of the genetic diversity preserved within this endangered conifer and other five-needle pines.

Genetic mapping of Pinus flexilis major gene (Cr4) for resistance to white pine blister rust using transcriptome-based SNP genotyping.

BACKGROUND: Linkage of DNA markers with phenotypic traits provides essential information to dissect clustered genes with potential phenotypic contributions in a target genome region. Pinus flexilis E. James (limber pine) is a keystone five-needle pine species in mountain-top ecosystems of North America. White pine blister rust (WPBR), caused by a non-native fungal pathogen Cronartium ribicola (J.C. Fisch.), has resulted in mortality in this conifer species and is still spreading through the distribution. The objective of this research was to develop P. flexilis transcriptome-wide single nucleotide polymorphism (SNP) markers using RNA-seq analysis for genetic mapping of the major gene (Cr4) that confers complete resistance to C. ribicola. RESULTS: Needle tissues of one resistant and two susceptible seedling families were subjected to RNA-seq analysis. In silico SNP markers were uncovered by mapping the RNA-seq reads back to the de novo assembled transcriptomes. A total of 110,573 in silico SNPs and 2,870 indels were identified with an average of 3.7 SNPs per Kb. These SNPs were distributed in 17,041 unigenes. Of these polymorphic P. flexilis unigenes, 6,584 were highly conserved as compared to the genome sequence of P. taeda L (loblolly pine). High-throughput genotyping arrays were designed and were used to search for Cr4-linked genic SNPs in megagametophyte populations of four maternal trees by haploid-segregation analysis. A total of 32 SNP markers in 25 genes were localized on the Cr4 linkage group (LG). Syntenic relationships of this Cr4-LG map with the model conifer species P. taeda anchored Cr4 on Pinus consensus LG8, indicating that R genes against C. ribicola have evolved independently in different five-needle pines. Functional genes close to Cr4 were annotated and their potential roles in Cr4-mediated resistance were further discussed. CONCLUSIONS: We demonstrated a very effective, low-cost strategy for developing a SNP genetic map of a phenotypic trait of interest. SNP discovery through transcriptome comparison was integrated with high-throughput genotyping of a small set of in silico SNPs. This strategy may be applied to mapping any trait in non-model plant species that have complex genomes. Whole transcriptome sequencing provides a powerful tool for SNP discovery in conifers and other species with complex genomes, for which sequencing and annotation of complex genomes is still challenging. The genic SNP map for the consensus Cr4-LG may help future molecular breeding efforts by enabling both Cr4 positional characterization and selection of this gene against WPBR.

Characterization of Five Novel Mitoviruses in the White Pine Blister Rust Fungus Cronartium ribicola.

The white pine blister rust (WPBR) fungus Cronartium ribicola (J.C. Fisch.) is an exotic invasive forest pathogen causing severe stem canker disease of native white pine trees (subgenus Strobus) in North America. The present study reports discovery of five novel mitoviruses in C. ribicola by deep RNA sequencing. The complete genome of each mitovirus was determined by rapid amplification of cDNA ends (RACE) and reverse transcriptase-polymerase chain reaction (RT-PCR). A single open reading frame (ORF) encoding a putative RNA-dependent RNA polymerase (RdRp) was detected in each of the viral genomes using mitochondrial genetic codes. Phylogenetic analysis indicated that the C. ribicola mitoviruses (CrMV1 to CrMV5) are new putative species of the genus Mitovirus. qRT-PCR and RNA-Seq analyses revealed that viral RNAs were significantly increased in fungal mycelia in cankered pine stems compared to expression during two different stages of spore development, suggesting that viral genome replication and transcription benefit from active growth of the host fungus. CrMVs were widespread with relatively high levels of minor allele frequency (MAF) in western North America. As the first report of mitoviruses in the Class Pucciniomycetes, this work allows further investigation of the dynamics of a viral community in the WPBR pathosystem, including potential impacts that may affect pathogenicity and virulence of the host fungus.

Transcriptome analysis of the white pine blister rust pathogen Cronartium ribicola: de novo assembly, expression profiling, and identification of candidate effectors.

BACKGROUND: The fungus Cronartium ribicola (Cri) is an economically and ecologically important forest pathogen that causes white pine blister rust (WPBR) disease on five-needle pines. To cause stem cankers and kill white pine trees the fungus elaborates a life cycle with five stages of spore development on five-needle pines and the alternate host Ribes plants. To increase our understanding of molecular WP-BR interactions, here we report genome-wide transcriptional profile analysis of C. ribicola using RNA-seq. RESULTS: cDNA libraries were constructed from aeciospore, urediniospore, and western white pine (Pinus monticola) tissues post Cri infection. Over 200 million RNA-seq 100-bp paired-end (PE) reads from rust fungal spores were de novo assembled and a reference transcriptome was generated with 17,880 transcripts that were expressed from 13,629 unigenes. A total of 734 unique proteins were predicted as a part of the Cri secretome from complete open reading frames (ORFs), and 41 % of them were Cronartium-specific. This study further identified a repertoire of candidate effectors and other pathogenicity determinants. Differentially expressed genes (DEGs) were identified to gain an understanding of molecular events important during the WPBR fungus life cycle by comparing Cri transcriptomes at different infection stages. Large-scale changes of in planta gene expression profiles were observed, revealing that multiple fungal biosynthetic pathways were enhanced during mycelium growth inside infected pine stem tissues. Conversely, many fungal genes that were up-regulated at the urediniospore stage appeared to be signalling components and transporters. The secreted fungal protein genes that were up-regulated in pine needle tissues during early infection were primarily associated with cell wall modifications, possibly to mask the rust pathogen from plant defenses. CONCLUSION: This comprehensive transcriptome profiling substantially improves our current understanding of molecular WP-BR interactions. The repertoire of candidate effectors and other putative pathogenicity determinants identified here are valuable for future functional analysis of Cri virulence and pathogenicity.