Divergent diversity patterns of NBS and LRR domains of resistance gene analogs in wild emmer wheat populations

Identification and characterization of genome-wide resistance gene analogs (RGAs) of durian (Durio zibethinus L.)

BACKGROUND

Durian (Durio zibethinus L.) is a tropical fruit crop which is popular in Southeast Asia but recently gaining popularity in other parts of the world. In this study, we analyzed the resistance gene analogs (RGAs) of durian through mining of the currently available reference genome of its 'Musang King' cultivar (PRJNA400310).

RESULTS

A total of 2586 RGAs were identified in the durian genome consisting of 47 nucleotide binding site proteins (NBS), 158 NBS-leucine rich repeat proteins (NL), 400 coiled-coil NBS-LRR (CNL), 72 toll/interleukin-1 receptor NBS-LRR (TNL), 54 coiled-coil NBS (CN), 10 toll/interleukin-1 receptor NBS (TN), 19 toll/interleukin-1 receptor with unknown domain (TX), 246 receptor-like proteins (RLP), 1,377 receptor-like kinases (RLK), 185 TM-CC, and 18 other NBS-containing proteins with other domains. These RGAs were functionally annotated and characterized via gene ontology (GO) analysis. Among the RGAs with the highest copies in durian genome include the putative disease resistance RPP13-like protein 1, disease resistance protein At4g27190, disease resistance protein RPS6, Probable disease resistance protein At4g27220, and putative disease resistance protein RGA3, while 35 RGAs were found to be novel. Phylogenetic analyses revealed that the genome-wide RGAs were broadly clustered into four major clades based on their domain classification.

CONCLUSION

To our knowledge, this is the most comprehensive analysis of durian RGAs which provides a valuable resource for genetic, agronomic, and other biological research of this important tropical fruit crop.

Subsets of NLR genes show differential signatures of adaptation during colonization of new habitats

Nucleotide binding site, leucine-rich repeat receptors (NLRs) are canonical resistance (R) genes in plants, fungi and animals, functioning as central (helper) and peripheral (sensor) genes in a signalling network. We investigate NLR evolution during the colonization of novel habitats in a model tomato species, Solanum chilense. We used R-gene enrichment sequencing to obtain polymorphism data at NLRs of 140 plants sampled across 14 populations covering the whole species range. We inferred the past demographic history of habitat colonization by resequencing whole genomes from three S. chilense plants from three key populations and performing approximate Bayesian computation using data from the 14 populations. Using these parameters, we simulated the genetic differentiation statistics distribution expected under neutral NLR evolution and identified small subsets of outlier NLRs exhibiting signatures of selection across populations. NLRs under selection between habitats are more often helper genes, whereas those showing signatures of adaptation in single populations are more often sensor-NLRs. Thus, centrality in the NLR network does not constrain NLR evolvability, and new mutations in central genes in the network are key for R-gene adaptation during colonization of different habitats.

Diversity and evolution of Rp1 rust resistance genes in four maize lines

This manuscript provides genome-level analysis of disease resistance genes in four maize lines, including studies of haplotype and resistance gene number as well as selection and recombination. The Rp1 locus of maize is a complex resistance gene (R-gene) cluster that confers race-specific resistance to Puccinia sorghi, the causal agent of common leaf rust. Rp1 NB-LRR disease resistance genes were isolated from two Rp1 haplotypes (HRp1-B and HRp1-M) and two maize inbred lines (B73 and H95). Sixty-one Rp1 genes were isolated from Rp1-B, Rp1-M, B73 and H95 with a PCR-based approach. The four maize lines carried from 12 to 19 Rp1 genes. From 4 to 9 of the identified Rp1 genes were transcribed in the four maize lines. The Rp1 gene nucleotide diversity was higher in HRp1-B and HRp1-M than in B73 and H95. Phylogenic analysis of 69 Rp1 genes revealed that the Rp1 genes maintained in HRp1-B, HRp1-M and H95 are evolving independently of each other, while Rp1 genes in B73 and HRp1-D appear more like each other than they do genes in the other lines. The results also revealed that the analysed Rp1 R-genes were under positive selection in HRp1-M and B73. Intragenic recombination was detected in Rp1 genes maintained in the four maize lines. This demonstrates that a genetic process that has the potential to generate new resistance genes with new specificities is active at the Rp1 locus in the four analysed maize lines and that the new resistance genes may act against newly arising pathogen races that become prevalent in the pathogen population.

Three-dimensional modeling and diversity analysis reveals distinct AVR recognition sites and evolutionary pathways in wild and domesticated wheat Pm3 R genes

The Pm3 gene confers resistance against wheat powdery mildew. Studies of Pm3 diversity have shown that Pm3 alleles isolated from southern populations of wild emmer wheat located in Lebanon, Jordan, Israel, and Syria are more diverse and more distant from bread wheat alleles than alleles from the northern wild wheat populations located in Turkey, Iran, and Iraq. Therefore, southern populations from Israel were studied extensively to reveal novel Pm3 alleles that are absent from the cultivated gene pool. Candidate Pm3 genes were isolated via a polymerase chain reaction cloning approach. Known and newly identified Pm3 genes were subjected to variation analysis and polymorphic amino acid residues were superimposed on a three-dimensional (3D) model of PM3. The region of highest interspecies diversity between Triticum aestivum and T. dicoccoides lies in leucine-rich repeats (LRR) 19 to 24, whereas most intraspecies diversity in T. aestivum is located in LRR 25 to 28. Interestingly, these two regions are separated by one large LRR whose propensity for flexibility facilitates the conformation of the PM3 LRR domain into two differently structured models. The combination of evolutionary and protein 3D structure analysis revealed that Pm3 genes in wild and domesticated wheat show different evolutionary histories which might have been triggered through different interactions with the powdery mildew pathogen.

Isolation and diversity analysis of resistance gene homologues from switchgrass

Resistance gene homologs (RGHs) were isolated from the switchgrass variety Alamo by a combination of polymerase chain reaction and expressed sequence tag (EST) database mining. Fifty-eight RGHs were isolated by polymerase chain reaction and 295 RGHs were identified in 424,545 switchgrass ESTs. Four nucleotide binding site--leucine-rich repeat RGHs were selected to investigate RGH haplotypic diversity in seven switchgrass varieties chosen for their representation of a broad range of the switchgrass germplasm. Lowland and upland ecotypes were found to be less similar, even from nearby populations, than were more distant populations with similar growth environments. Most (83.5%) of the variability in these four RGHs was found to be attributable to the within-population component. The difference in nucleotide diversity between and within populations was observed to be small, whereas this diversity is maintained to similar degrees at both population and ecotype levels. The results also revealed that the analyzed RGHs were under positive selection in the studied switchgrass accessions. Intragenic recombination was detected in switchgrass RGHs, thereby demonstrating an active genetic process that has the potential to generate new resistance genes with new specificities that might act against newly-arising pathogen races.

Genetic variation within and between winter wheat genotypes from Turkey, Kazakhstan, and Europe as determined by nucleotide-binding-site profiling

The genetic diversity within wheat breeding programs across Turkey and Kazakhstan was compared with a selection of European cultivars that represented the genetic diversity across eight European countries and six decades of wheat breeding. To focus the measure of genetic diversity on that relevant to disease-resistant phenotypes, nucleotide-binding-site (NBS) profiling was used to detect polymorphisms associated with the NBS motifs found within the NBS--leucine-rich repeat (LRR) class of resistance (R) genes. Cereal-specific NBS primers, designed specifically to the conserved NBS motifs found within cereal R-genes, provided distinct NBS profiles. Although the genetic diversity associated with NBS motifs was only slightly higher within the Eastern wheat genotypes, the NBS profiles produced by Eastern and European wheat lines differed considerably. Structure analysis divided the wheat genotypes into four groups, which compared well with the origin of the wheat genotypes. The highest levels of genetic diversity were seen for the wheat genotypes from the Genetic Resource Collection held in Ankara, Turkey, as wheat genotypes within breeding programs were genetically more similar. The wheat genotypes from Kazakhstan were the most similar to the European cultivars, reflecting the significant number of eastern European cultivars used in the breeding program in Kazakhstan. In general, the NBS profiles suggested that NBS-LRR R-gene usage in winter wheat breeding in Turkey and Kazakhstan differed from that deployed in European cultivars.

Rapid linkage disequilibrium decay in the Lr10 gene in wild emmer wheat (Triticum dicoccoides) populations

INTRODUCTION

Recombination is a key evolutionary factor enhancing diversity. However, the effect of recombination on diversity in inbreeding species is expected to be low. To estimate this effect, recombination and diversity patterns of Lr10 gene were studied in natural populations of the inbreeder species, wild emmer wheat (Triticum dicoccoides). Wild emmer wheat is the progenitor of most cultivated wheats and it harbors rich genetic resources for disease resistance. Lr10 is a leaf rust resistance gene encoding three domains: a coiled-coil, nucleotide-binding site, and leucine-rich repeat (CC-NBS-LRR).

RESULTS

Lr10 was sequenced from 58 accessions representing 12 diverse habitats in Israel. Diversity analysis revealed a high rate of synonymous and non-synonymous substitutions (d (S) = 0.029, d (N) = 0.018, respectively) in the NBS-LRR domains. Moreover, in contrast to other resistance genes, in Lr10 the CC domain was more diverse than the NBS-LRR domains (d (S) = 0.069 vs. 0.029, d (N) = 0.094 vs. 0.018) and was subjected to positive selection in some of the populations. Seventeen recombination events were detected between haplotypes, especially in the CC domain. Linkage disequilibrium (LD) analysis has shown a rapid decay from r (2) = 0.5 to r (2) = 0.1 within a 2-kb span.

CONCLUSION

These results suggest that recombination is a diversifying force for the R-gene, Lr10, in the selfing species T. dicoccoides. This is the first report of a short-range LD decay in wild emmer wheat.

Pathogen race determines the type of resistance response in the stripe rust-Triticum dicoccoides pathosystem

Wild relatives of crop plants may serve as a promising source for screening for new disease resistance genes that can be utilized in breeding programs. Triticum dicoccoides, the wild progenitor of most cultivated wheats, was shown to harbor many resistance genes against the major diseases attacking cultivated wheat. Stripe rust is a devastating fungal disease that attacks wheat in many regions of the world. New races of Puccinia striiformis Westend. f. sp. tritici, the causative agent of stripe rust, have overcome most of the known Yr resistance genes in wheat. Therefore, there is a need to search for new resistance genes in the T. dicoccoides gene pool. A set of 120 T. dicoccoides accessions, collected from 13 populations representing different habitats in Israel and vicinity, was tested for resistance to three prevalent stripe rust races (38E134, 6E16 and 6E0). Of these 120 accessions, 14, 8 and 12% were resistant to races 38E134, 6E16 and 6E0, respectively, while 57, 2 and 4% were moderately resistant to these races, respectively. A unique resistance was found in the population of Mt Hermon where >80% of the accessions showed resistance to all races. Distribution of infection types (ITs) of race 38E134 showed a normal distribution that can fit a quantitative pattern of response, while the distributions of ITs of races 6E16 and 6E0 had excess of extreme values and therefore showing a qualitative pattern of response. anova testing the main factor effects and interaction showed significant effects of population, race and their interaction on IT. Significant positive correlations were obtained between the resistance to races 6E16 and 6E0 and humidity variables of the collections sites, while resistance to race 38E134 was positively correlated with temperature variables. These results show that the pathogen race can determine the type of resistance response, qualitative or quantitative, in the stripe rust-T. dicoccoides pathosystem. The obtained results also reveal that the distribution of resistance to different pathogen races can be affected by different climatic factors.