Full text

Turn on search term navigation

Introduction

Western white pine (Pinus monticola Dougl.) is an ecologically dominant, and highly important for forestry, species in the western North American ecosystems ([20]). In 1980, two P. monticola provenances from Jackson, Oregon were introduced in Greece in comparative species field trials of the Aristole University of Thessaloniki, Laboratory of Forest Genetics Experimental Centers located in Pertouli and Taxiarchis ([26], [27]). The aim of this P. monticola introduction was twofold: (a) the use of the species in plantation forestry for wood production, and (b) the evaluation of resistance to the white pine blister rust (Cronartium ribicola J.C. Fisch.). Trial results indicate large variation within provenance, growth comparable to that of local Abies borisii regis Mattf., but reduced survival ([28], [27]).

Cronartium ribicola is among the most destructive white pine fungal diseases which has caused severe mortality in North American white pines and especially in P. monticola, which is currently considered a species under threat ([21]). P. monticola species exhibits genetic variation in disease resistance to C. ribicola involving mainly a dominant disease resistance gene (Cr2 - [12], [19]). Disease resistance and maintenance of genetic diversity in production populations are very important aspects for P. monticola breeding ([4]). The development through breeding of trees resistant to fungi and the maintenance of a variety of genetic resources are of pivotal significance towards the substantial establishment of P. monticola reforestation and plantation forestry programs ([18]).

Classical breeding approaches, such as interspecific hybridization with the indigenous to Greece and resistant to C. ribicola Balkan white pine Pinus peuce Grisb., are time-consuming without warranties of specific results in a definite timeframe. Engagement of innovative molecular breeding approaches allows breeders to establish sustainable resistance more accurately and at a faster pace ([14]). PCR approaches implementing mainly degenerate primers were extensively applied in order to isolate members of the NBS family (nucleotide-binding site) / LRRs (Leucine-Rich Repeats) - containing genes, the most highly expanded group of genes linked directly to resistance (R) genes mediated functions in plants ([5]). This protein family is further classified into two major subfamilies based on structural characteristics of their putative amino-terminal signaling domain ([13]). One subfamily, TIR-NBS-LRR, contains an amino-terminal signaling domain that depicts homology with Toll and interleukin- 1 receptor (TIR), whilst the other subfamily, CC-NBS-LRR, instead of the TIR domain, carries a putative coiled-coil structure (CC) one ([13]). Both these subfamilies have been reported to mediate resistance to a wide range of phytopathenic fungi ([34]).

Resistance gene analogs (RGAs) are a large class of potential R-genes, being an important resource for the development of molecular markers, hence enhancing numerous disease resistance breeding efforts ([1], [30]). Moreover, these RGAs commonly contain tandem repeats such as LRRs which give rise to novel ligand-binding specificities under positive selection ([11]) via a number of different modes, such as point mutations, variations in repeat numbers and gene duplications ([8], [24], [37], [39]). Positive selection is defined as the rapid fixation of beneficial non-synonymous mutations and is an important evolutionary driving force for a number of paralogous genes towards their functional diversification, which might be selectively favored under biotic pressures in plants ([6]). Furthermore, the location of positively selected amino acid residues is an important aspect for obtaining novel gene functions ([25]). Previous studies have reported that solvent-exposed regions of the LRRs repeats are subjected to positive selection; such evidence has been interpreted as a sign of the involvement of these regions in pathogen recognition ([29], [39], [11]). Overall, the evolution patterns of plant NBS-LRR R genes and of their RGAs counterparts is likely a complex process ([40], [31]).

The TIR-NBS-LRR RGAs from P. monticola were previously investigated for their diversity and expression, in order to gain insights regarding their potential usage in disease resistance breeding programs ([17]). In this study, we address the question whether signatures of positive selection are systematically employed as an evolutionary driven force acting on TIR-NBS-LRR-encoding RGAs genes at P. monticola. Thus, we gain insights into the evolutionary profiles of the P. monticola TIR-NBS-LRR-containing RGAs genes, hypothesizing that those successive episodes of positive selection might contribute in the acquisition of novel pathogens recognition repertoires. Overall, our findings enrich and support the hypothesis that these P. monticola RGAs could be employed as biomarkers for both, monitoring the ongoing and, more effectively applying future P. monticola disease resistance breeding efforts.

Materials and methods

P. monticola TIR-NBS-LRR RGAs mining and phylogenetic analysis

All 33 P. monticola RGAs partial proteins previously identified by Liu & Ekramoddoullah ([17]) were downloaded from the GenBank by keyword searches based on their NB-ARC domain profile (PF00931.19/ CL0023). Using these amino acid sequence IDs, the respective ORF nucleotide sequences were also retrieved. All amino acid sequences had to fulfill the following criteria in older to be evaluated for evolutionary signatures acting upon them: (a) existence of valid structures by means of protein lengths up to 90 residues, positions of NB-ARC domain motifs, blast(p) hits against known plant resistance genes; and (b) verification of functional domains using the predictions against the Pfam database (⇒ http:/­/­pfam.­sanger.­ac.­uk/­search#tabview= tab1) and the InterProScan 5 software ([9]).

The phylogenetic relationships among these P. monticola RGAs amino acid sequences were revealed by aligning them using the MUSCLE software ([7]) and performing a tree reconstruction using the RAxML software ([32]) inferring a gamma model of rate heterogeneity. All these analyses were performed on the Geneious R7 platform ([10]).

P. monticola TIR-NBS-LRR RGAs evolutionary analysis

The YN00, CODEML and CODEMLSITES programs were employed from the PAML version 4.8 software package ([36]) in order to assess evidence of positive selection signatures acting in these P. monticola TIR-NBS-LRR RGAs. These tests compare the substitution rates of non-synonymous (dN) to synonymous (dS) mutations, and estimate the derived omega (ω) values, both among the nucleotide codon sequences and along the branch lineages. The rates of dN/dS nucleotide substitutions per site across all possible pairwise sequence comparisons were estimated using the counting approach ([35]) implemented in the YN00 program. The total dN, dS, and ω values for all pairwise sequence comparisons among these RGAs are shown in Tab. S1 (Supplementary material).

An NJ tree at amino acid level and the respective nucleotide alignment were both used in the CODEML and CODEMLSITES analyses according to Lynn et al. ([22], [23]); details of the exact methods for these analyses are found in Zambounis et al. ([39]). Briefly, the amino acid sequence alignment was conducted by employing the MUSCLE software ([7]), whilst the construction of the respective phylogenetic tree was performed using the MEGA 5 software ([33]). A Poisson correction was employed with a gamma distribution modeling among sites and rates variations ([39]). Finally, the Baseml software (PAML version 4.8 - [36]) was implemented for estimating the alpha shape parameter using the REV model.

CODEML was tested in P. monticola TIR-NBS-LRR RGAs dataset for estimating the variable selective pressures acting among lineages in the phylogeny, whilst CODEMLSITES was employed to test for site-specific codon substitution models allowing the identification of selective signatures across amino acid residues (Tab. 1). In both programs log-likelihood calculations were performed for each model and comparisons between each other were implemented by likelihood ratio tests (LRTs). The empirical Bayesian method further allowed the detection of codon residues subjected to positive selection by higher posterior probabilities. All data sets as well as the outputs of CODEML and CODEMLSITES programs are available upon request.

Tab. 1 - Overview of CODEMLSITES method statistics and distribution of positively selected amino acid sites across TIR-NBS-LRR RGAs in P. monticola.

Model P Parameters l dn/ds Positively Selected Sites Statistical significance
M0: one ratio 2 ω = 0.395 -2855.76 = ω - -
M1: neutral 3 p0 = 0.72, ω0 = 0.22 -2821.47 0.44 - -
p1 = 0.27, ω1 = 1
M2: selection 5 p0 = 0.72, ω0 = 0.22 -2821.47 0.44 - P

Results

P. monticola TIR-NBS-LRR RGAs structure and phylogenetic analysis

The 33 P. monticola TIR-NBS-LRR RGAs studied were found to share Pfam hits with the NB-ARC domain (PF00931.19) and the InterProScan ones with LRR-containing proteins (PTHR23155). The RGAs were corresponding to regions between conserved kinase-2 and P-loop domains. A RAxML phylogenetical analysis at the amino acid level was conducted in order to explore their expansion (Fig. 1). The backbone topology of this tree was properly resolved with 63 branches (Fig. 1), whilst the overall pairwise identity among the amino acid sequences was 59.6%, indicative of a high degree of divergence.

Fig. 1 - RAxML phylogenetic reconstruction of the amino acid sequences of P. monticola TIR-NBS-LRR RGAs. In various colors the six clades which were being assigned whilst, asterisks (*) indicate tree branches under positive selection ([omega] values higher than 1) with posterior probabilities using the CODEML method.

Discussion

Conifer genomes are in general large with DNA content reaching more than 26.500 Mbp per haploid genome ([2] and references therein). Taking also into consideration their long generation times, classical breeding, especially for disease resistance is a very slow process. On the contrary, molecular breeding approaches may provide promising and effective strategies for both R genes and RGAs counterpart identification based on their highly conserved motifs ([20]). Previously, it was postulated that there is a highly diverse group of TIR-NBS-LRR genes in P. monticola ([17]) along with the existence of numerous CC-NBS-LRR-containing proteins ([20]). Similar abundance and diversity of these genes was also observed in other plant genomes ([40], [3], [31]). Hitherto, TIR-NBS-LRR RGAs in P. monticola exhibited high variation at amino acid levels, whilst almost half of them were found to be pseudogenes representing a rather complex TIR-NBS-LRR subfamily ([17]).

In our study, the analysis of 33 P. monticola TIR-NBS-LRR RGAs showed a high level of diversity that may contribute to their evolutionary fitness, and the putatively functional character of their transcription activity ([17]). These RGAs are a potential source for the development of molecular markers, especially regarding selection of resistance to white pine blister rust. Furthermore, RGAs were found to share conserved structures rendering them ideal candidates for applying evolutionary bioinformatics approaches.

The CODEML analysis, showed recent episodes of positive selection acting extensively in clades I, IV and V, which in turn were overlapping similar more ancient events. Similar positive selective pressure acting on plant RGAs has also been reported in other surveys ([15], [34], [30], [11]). In principle, this reservoir of genetic variation could drive the evolution of new R-gene specificities ([30], [31]). To our knowledge this is the first systematically report of evolutionary pressures acting in P. monticola RGAs. Positively selected RGAs branches were not only at terminal positions, implying rather recent positive selective signatures, but also at ancestral positions on the phylogeny. This finding is in accordance with previous reports which showed that numerous selective episodes occur even from the birth of a paralogous locus ([37], [39]).

Signs of positively selected residues were found to be distributed rather evenly across the RGAs amino acid sequences. Positive selection on NBS-containing gene families in plants and especially among their LRR domains, were also reported in various plant hosts ([3], [16], [39], [38], [30], [11]). These findings support the view that selection for durable disease resistance may be a crucial component of plant breeding programs. Additionally, comparative analyses of Rosaceae trees regarding their R genes revealed that solvent-exposed residues of the LRRs domains are hyper-variable with intensive positive selective pressures acting on them ([30]). Moreover, such evidence is consistent with host-pathogen co-evolution processes leading to the acquisition of novel resistance functional specificities ([39], [30]). These findings would facilitate the development of informative RGA-derived molecular markers for identifying fungal resistant P. monticola and P. monticola × P. peuce genotypes in breeding programs using segregating progenies. They may also prove to be of importance for P. peuce itself, as despite the fact that Balkan white pine is a resistant species, future selection pressures on C. ribicola in light of climatic change may alter the virulence of the fungus towards its host.

Conclusions

The P. monticola RGAs gene homologues encoding for R gene members of the TIR-NBS-LRR subfamily were found to be subjected to intensive positive selection, acting widely in three clades across their phylogeny, both on terminal and ancestral lineages. The majority of positively selected residues sites were found to be localized widely across these RGAs sequences, putatively affecting the structures of their ligand-binding domains. Practical applications in the development of disease resistant plants through breeding applications may also be considered, facilitating for instance, the development of informative RGAs-derived sequence characterized amplified region (SCAR) markers, based on RGAs sequences with extensive positive selection acting upon them. Such markers might be particularly useful for identifying fungal resistant genotypes in routine marker-assisted selection (MAS) breeding programs.

List of Abbreviations

The following abbreviations have been used throughout the paper:

1. CC: N-terminal coiled-coil domain;

2. LRR: Leucine Rich Repeat;

3. NBS: nucleotide-binding site domain;

4. PAML: Phylogenetic Analysis by Maximum Likelihood;

5. RaxML: Randomized Axelerated Maximum Likelihood;

6. RGA: Resistance Gene Analog;

7. TIR: Toll/interleukin-1 receptor domain;

8. WWP: Western white pine.

Acknowledgments

The co-authors have not a conflict of interest to declare.

(1) Ameline-Torregrosa C, Wang B-B, O’Bleness MS, Deshpande S, Zhu H, Roe B, Young ND, Cannon SB (2008). Identification and characterization of nucleotide-binding site-leucine-rich repeat genes in the model plant Medicago truncatula. Plant Physiology 146 (1): 5-21.

(2) Aravanopoulos FA (2014). Genomics of trees. In: “Tree Biotechnology” (Ramawat KG, Merillon JM, Ahuja MR eds). CRC Press, Boca Raton, USA, pp. 514-557.

(3) Chen Q, Han Z, Jiang H, Tian D, Yang S (2010). Strong positive selection drives rapid diversification of R-genes in Arabidopsis relatives. Journal of Molecular Evolution 70 (2): 137-148.

(4) Dangl JL, Horvath DM, Staskawicz BJ (2013). Pivoting the plant immune system from dissection to deployment. Science 341 (6147): 746-751.

(5) Debener T, Byrne DH (2014). Disease resistance breeding in rose: current status and potential of biotechnological tools. Plant Science 228: 107-117.

(6) Delph LF, Kelly JK (2014). On the importance of balancing selection in plants. New Phytologist 201 (1): 45-56.

(7) Edgar RC (2004). MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 5: 1-19.

(8) Friedman AR, Baker BJ (2007). The evolution of resistance genes in multi-protein plant resistance systems. Current Opinion in Genetics and Development 17 (6): 493-499.

(9) Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, McWilliam H, Maslen J, Mitchell A, Nuka G, Pesseat S, Quinn AF, Sangrador-Vegas A, Scheremetjew M, Yong SY, Lopez R, Hunter S (2014). InterProScan 5: genome-scale protein function classification. Bioinformatics 30 (9): 1236-1240.

(10) Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, Thierer T, Ashton B, Mentjies P, Drummond A (2012). Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28 (12): 1647-1649.

(11) Khan AM, Khan AA, Azhar MT, Amrao L, Cheema HM (2016). Comparative analysis of resistance gene analogues encoding NBS-LRR domains in cotton. Journal of the Science of Food and Agriculture 96 (2): 530-538.

(12) Kinloch BB, Sniezko RA, Barnes GD, Greathouse TE (1999). A major gene for resistance to white pine blister rust in western white pine from the western Cascade Range. Phytopathology 89: 861-867.

(13) Kim J, Lim CJ, Lee BW, Choi JP, Oh SK, Ahmad R, Kwon SY, Ahn J, Hur CG (2012). A genome-wide comparison of NB-LRR type of resistance gene analogs (RGA) in the plant kingdom. Molecules and Cells 33 (4): 385-392.

(14) Lalli DA, Decroocq V, Blenda AV, Schurdi-Levraud V, Garay L, Le Gall O, Damsteegt V, Reighard GL, Abbott AG (2005). Identification and mapping of resistance gene analogs (RGAs) in Prunus: a resistance map for Prunus. Theoretical and Applied Genetics 111 (8): 1504-1513.

(15) Li J, Xu Y, Fei S, Li L (2006). Isolation, characterization and evolutionary analysis of resistance gene analogs in annual ryegrass, perennial ryegrass and their hybrid. Physiologia Plantarum 126 (4): 627-638.

(16) Li J, Ding J, Zhang W, Zhang Y, Tang P, Chen J-Q, Tian D, Yang S (2010). Unique evolutionary pattern of numbers of Gramineous NBS-LRR genes. Molecular Genetics and Genomics 283 (5): 427-438.

(17) Liu JJ, Ekramoddoullah AK (2003). Isolation, genetic variation and expression of TIR-NBS-LRR resistance gene analogs from western white pine (Pinus monticola Dougl. ex. D. Don.). Molecular Genetics and Genomics 270: 432-441.

(18) Liu JJ, Hunt R, Ekramoddoullah AK (2004). Recent insights into genetic resistance of western white pine to white pine blister rust. Recent Research Developments in Biotechnology and Bioengineering 6: 65-76. Online | Gscholar

(19) Liu JJ, Ekramoddoullah AK, Hunt RS, Zamani A (2006). Identification and characterization of random amplified polymorphic DNA markers linked to a major gene (Cr2) for resistance to Cronartium ribicola in Pinus monticola. Phytopathology 4: 395-399.

(20) Liu JJ, Ekramoddoullah AK (2007). The CC-NBS-LRR subfamily in Pinus monticola: targeted identification, gene expression, and genetic linkage with resistance to Cronartium ribicola. Phytopathology 97 (6): 728-736.

(21) Loehman RA, Clark JA, Keane RE (2011). Modeling effects of climate change and fire management on western white pine (Pinus monticola) in the northern Rocky Mountains, USA. Forests 2 (4): 832-860.

(22) Lynn DJ, Higgs R, Gaines S, Tierney J, James T, Lloyd AT, Fares MA, Mulcahy G, Farrelly C (2004a). Bioinformatic discovery and initial characterisation of nine novel antimicrobial peptide genes in the chicken. Immunogenetics 56 (3): 170-177.

(23) Lynn DJ, Lloyd AT, Fares MA, O’Farrelly C (2004b). Evidence of positively selected sites in mammalian α-defensins. Molecular Biology and Evolution 21 (5): 819-827.

(24) Meyers BC, Kozik A, Griego A, Kuang H, Michelmore RW (2003). Genome-wide analysis of NBS-LRR-encoding genes in Arabidopsis. Plant Cell 15 (4): 809-834.

(25) Mondragón-Palomino M, Meyers BC, Michelmore RW, Gaut BS (2002). Patterns of positive selection in the complete NBS-LRR gene family of Arabidopsis thaliana. Genome Research 12 (9): 1305-1315.

(26) Panetsos KP, Moulalis D, Mitsopoulos D, Skaltsoyiannes A (1983). Results of exotic and native forest species trials at the Aristotle University forests. Laboratory of Forest Genetics and Tree Breeding, Aristotle University, Publ. No. 6, Thessaloniki, Greece, pp. 30. [in Greek with English summary] Gscholar

(27) Panetsos KP, Alizoti PG, Moulalis D, Skaltsoyiannes A, Aravanopoulos FA (1996). Results of exotic forest species testing in Pertouli. In: Proceedings of the “6th Pan-Hellenic Conference of the Hellenic Scientific Society of Plant Genetics and Breeding”. Florina (Greece) 2-4 Oct 1996. Hellenic Scientific Society for Plant Genetics and Breeding , Thessaloniki, Greece, pp. 310-316. [in Greek with English summary] Gscholar

(28) Papadima A, Sotikopoulos N (1995). Results of exotic and native forest species trials at the Aristotle University forests. BSc Thesis, Faculty of Agriculture, Forestry and Natural Environment, Aristotle University of Thessaloniki, Aristotle University Publ., Thessaloniki, Greece, pp. 73. [in Greek with English summary] Gscholar

(29) Parniske M, Hammond-Kosack KE, Golstein C, Thomas CM, Jones DA, Harrison K, Wulff BBH, Jones JDG (1997). Novel disease resistance specificities result from sequence exchange between tandemly repeated genes at the Cf-4/9 locus of tomato. Cell 91 (6): 821-832.

(30) Perazzolli M, Malacarne G, Baldo A, Righetti L, Bailey A, Fontana P, Velasco R, Malnoy M (2014). Characterization of resistance gene analogues (RGAs) in apple (Malus x domestica Borkh.) and their evolutionary history of the Rosaceae family. PloS One 9 (2): e83844.

(31) Sekhwal MK, Li P, Lam I, Wang X, Cloutier S, You FM (2015). Disease resistance gene analogs (RGAs) in plants. International Journal of Molecular Sciences 16 (8): 19248-19290.

(32) Stamatakis A (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30 (9): 1312-1313.

(33) Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Molecular Biology and Evolution 28 (10): 2731-2739.

(34) Wan H, Zhao Z, Malik AA, Qian C, Chen J (2010). Identification and characterization of potential NBS-encoding resistance genes and induction kinetics of a putative candidate gene associated with downy mildew resistance in Cucumis. BMC Plant Biology 10 (1): 186.

(35) Yang Z, Nielsen R (2000). Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Molecular Biology and Evolution 17 (1): 32-43.

(36) Yang Z (2007). PAML 4: phylogenetic analysis by maximum likelihood. Molecular Biology and Evolution 24 (8): 1586-1591.

(37) Yang S, Zhang X, Yue J-X, Tian D, Chen J-Q (2008). Recent duplications dominate NBS-encoding gene expansion in two woody species. Molecular Genetics and Genomics 280 (3): 187-198.

(38) Yang S, Li J, Zhang X, Zhang Q, Huang J, Chen J-Q, Hartl DL, Tian D (2013). Rapidly evolving R genes in diverse grass species confer resistance to rice blast disease. Proceedings of the National Academy of Sciences USA 110 (46): 18572-18577.

(39) Zambounis A, Elias M, Sterck L, Maumus F, Gachon CMM (2012). Highly dynamic exon shuffling in candidate pathogen receptors… What if brown algae were capable of adaptive immunity? Molecular Biology and Evolution 29 (4): 1263-1276.

(40) Zhou T, Wang Y, Chen JQ, Araki H, Jing Z, Jiang K, Shen J, Tian D (2004). Genome-wide identification of NBS genes in japonica rice reveals significant expansion of divergent non-TIR NBS-LRR genes. Molecular Genetics and Genomics 271 (4): 402-415.

Word count: 3500

Show less

© 2017. This work is licensed under https://creativecommons.org/licenses/by-nc/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.

Abstract

Western white pine (Pinus monticola Dougl.) is an important forest tree species, which is intensively plagued by the fungus Cronartium ribicola. Resistance gene analogs (RGAs) are the most highly abundant class of potential resistance (R) genes sharing greatly conserved domains and structures. Hence RGAs are crucial components for disease resistance breeding programs on P. monticola serving as useful functional markers. A total of 33 P. monticola RGAs gene homologues were mined from GenBank, encoding for R gene members of the TIR-NBS-LRR subfamily. The existence of positive selection acting upon RGAs was determined using a series of maximum likelihood analyses. Robust evidence of positive selection was showed to be acting widely in three clades across RGA gene phylogeny, both on terminal and ancestral lineages. Furthermore, our analysis revealed that the majority of positively selected residues sites are localized widely across these RGAs sequences, putatively affecting the structures of their ligand-binding domains and offering novel specificities. These results may find immediate application in ongoing disease resistance breeding programs.

Details

Title
Adaptive response of Pinus monticola driven by positive selection upon resistance gene analogs (RGAs) of the TIR-NBS-LRR subfamily
Author
Zambounis, Antonios; Tsaftaris, Athanasios; Avramidou, Evangelia; Papadima, Aglaea; Barbas, Evaggelos; Aravanopoulos, Filippos A; Madesis, Panagiotis; Arzimanoglou, Iordanis
Pages
237-241
Section
Research Articles
Publication year
2017
Publication date
2017
Publisher
The Italian Society of Silviculture and Forest Ecology (SISEF)
ISSN
19717458
Source type
Scholarly Journal
Language of publication
English
DOI
https://doi.org/10.3832/ifor2050-009
ProQuest document ID
2661965661
Copyright
© 2017. This work is licensed under https://creativecommons.org/licenses/by-nc/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.