Pervasive purifying selection characterizes the evolution of I2 homologs

Evolutionary and Structural Analysis of PP16 in Viridiplantae

Members of the phloem protein 16 (PP16) gene family are induced by elicitors in rice and the corresponding proteins from cucurbits, which display RNA binding and intercellular transport activities, are accumulated in phloem sap. These proteins facilitate the movement of protein complexes through the phloem translocation flow and may be involved in the response to water deficit, among other functions. However, there is scant information regarding their function in other plants, including the identification of paralog genes in non-vascular plants and chlorophytes. In the present work, an evolutionary and structural analysis of the PP16 family in green plants (Viridiplantae) was carried out. Data mining in different databases indicated that PP16 likely originated from a larger gene present in an ancestral lineage that gave rise to chlorophytes and multicellular plants. This gene encodes a protein related to synaptotagmin, which is involved in vesicular transport in animal systems, although other members of this family play a role in lipid turnover in endomembranes and organelles. These proteins contain a membrane-binding C2 domain shared with PP16 proteins in vascular plants. In silico analysis of the predicted structure of the PP16 protein family identified several β-sheets, one α-helix, and intrinsically disordered regions. PP16 may have been originally involved in vesicular trafficking and/or membrane maintenance but specialized in long-distance signaling during the emergence of the plant vascular system.

Host-pathogen coevolution promotes the evolution of general, broad-spectrum resistance and reduces foreign pathogen spillover risk

Genetic variation for disease resistance within host populations can strongly impact the spread of endemic pathogens. In plants, recent work has shown that within-population variation in resistance can also affect the transmission of foreign spillover pathogens if that resistance is general. However, most hosts also possess specific resistance mechanisms that provide strong defenses against coevolved endemic pathogens. Here we use a modeling approach to ask how antagonistic coevolution between hosts and their endemic pathogen at the specific resistance locus can affect the frequency of general resistance, and therefore a host's vulnerability to foreign pathogens. We develop a two-locus model with variable recombination that incorporates both general (resistance to all pathogens) and specific (resistance to endemic pathogens only). We find that introducing coevolution into our model greatly expands the regions where general resistance can evolve, decreasing the risk of foreign pathogen invasion. Furthermore, coevolution greatly expands which conditions maintain polymorphisms at both resistance loci, thereby driving greater genetic diversity within host populations. This genetic diversity often leads to positive correlations between host resistance to foreign and endemic pathogens, similar to those observed in natural populations. However, if resistance loci become linked, the resistance correlations can shift to negative. If we include a third, linkage modifying locus into our model, we find that selection often favors complete linkage. Our model demonstrates how coevolutionary dynamics with an endemic pathogen can mold the resistance structure of host populations in ways that affect its susceptibility to foreign pathogen spillovers, and that the nature of these outcomes depends on resistance costs, as well as the degree of linkage between resistance genes.

Diversity, structure and function of the coiled-coil domains of plant NLR immune receptors

Plant nucleotide-binding, leucine-rich repeat receptors (NLRs) perceive pathogen avirulence effectors and activate defense responses. Nucleotide-binding, leucine-rich repeat receptors are classified into coiled-coil (CC)-containing and Toll/interleukin-1 receptor (TIR)-containing NLRs. Recent advances suggest that NLR CC domains often function in signaling activation, especially for induction of cell death. In this review, we outline our current understanding of NLR CC domains, including their diversity/classification and structure, their roles in cell death induction, disease resistance, and interaction with other proteins. Furthermore, we provide possible directions for future work.

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.

Non-host Plant Resistance against Phytophthora capsici Is Mediated in Part by Members of the I2 R Gene Family in Nicotiana spp

The identification of host genes associated with resistance to Phytophthora capsici is crucial to developing strategies of control against this oomycete pathogen. Since there are few sources of resistance to P. capsici in crop plants, non-host plants represent a promising source of resistance genes as well as excellent models to study P. capsici - plant interactions. We have previously shown that non-host resistance to P. capsici in Nicotiana spp. is mediated by the recognition of a specific P. capsici effector protein, PcAvr3a1 in a manner that suggests the involvement of a cognate disease resistance (R) genes. Here, we have used virus-induced gene silencing (VIGS) and transgenic tobacco plants expressing dsRNA in Nicotiana spp. to identify candidate R genes that mediate non-host resistance to P. capsici. Silencing of members of the I2 multigene family in the partially resistant plant N. edwardsonii and in the resistant N. tabacum resulted in compromised resistance to P. capsici. VIGS of two other components required for R gene-mediated resistance, EDS1 and SGT1, also enhanced susceptibility to P. capsici in N. edwardsonii, as well as in the susceptible plants N. benthamiana and N. clevelandii. The silencing of I2 family members in N. tabacum also compromised the recognition of PcAvr3a1. These results indicate that in this case, non-host resistance is mediated by the same components normally associated with race-specific resistance.

Whole genome sequencing of a banana wild relative Musa itinerans provides insights into lineage-specific diversification of the Musa genus

Crop wild relatives are valuable resources for future genetic improvement. Here, we report the de novo genome assembly of Musa itinerans, a disease-resistant wild banana relative in subtropical China. The assembled genome size was 462.1 Mb, covering 75.2% of the genome (615.2Mb) and containing 32, 456 predicted protein-coding genes. Since the approximate divergence around 5.8 million years ago, the genomes of Musa itinerans and Musa acuminata have shown conserved collinearity. Gene family expansions and contractions enrichment analysis revealed that some pathways were associated with phenotypic or physiological innovations. These include a transition from wood to herbaceous in the ancestral Musaceae, intensification of cold and drought tolerances, and reduced diseases resistance genes for subtropical marginally distributed Musa species. Prevalent purifying selection and transposed duplications were found to facilitate the diversification of NBS-encoding gene families for two Musa species. The population genome history analysis of M. itinerans revealed that the fluctuated population sizes were caused by the Pleistocene climate oscillations, and that the formation of Qiongzhou Strait might facilitate the population downsizing on the isolated Hainan Island about 10.3 Kya. The qualified assembly of the M. itinerans genome provides deep insights into the lineage-specific diversification and also valuable resources for future banana breeding.

Comparative transcriptomics of Central Asian Vitis vinifera accessions reveals distinct defense strategies against powdery mildew

Grape powdery mildew (PM), caused by the biotrophic ascomycete Erysiphe necator, is a devastating fungal disease that affects most Vitis vinifera cultivars. We have previously identified a panel of V. vinifera accessions from Central Asia with partial resistance to PM that possess a Ren1-like local haplotype. In this study, we show that in addition to the typical Ren1-associated late post-penetration resistance, these accessions display a range of different levels of disease development suggesting that alternative alleles or additional genes contribute to determining the outcome of the interaction with the pathogen. To identify potential Ren1-dependent transcriptional responses and functions associated with the different levels of resistance, we sequenced and analyzed the transcriptomes of these Central Asian accessions at two time points of PM infection. Transcriptomes were compared to identify constitutive differences and PM-inducible responses that may underlie their disease resistant phenotype. Responses to E. necator in all resistant accessions were characterized by an early up-regulation of 13 genes, most encoding putative defense functions, and a late down-regulation of 32 genes, enriched in transcriptional regulators and protein kinases. Potential Ren1-dependent responses included a hotspot of co-regulated genes on chromosome 18. We also identified 81 genes whose expression levels and dynamics correlated with the phenotypic differences between the most resistant accessions 'Karadzhandahal', DVIT3351.27, and O34-16 and the other genotypes. This study provides a first exploration of the functions associated with varying levels of partial resistance to PM in V. vinifera accessions that can be exploited as sources of genetic resistance in grape breeding programs.

Solanum resistance genes against Phytophthora infestans and their corresponding avirulence genes

Resistance genes against Phytophthora infestans (Rpi genes), the most important potato pathogen, are still highly valued in the breeding of Solanum spp. for enhanced resistance. The Rpi genes hitherto explored are localized most often in clusters, which are similar between the diverse Solanum genomes. Their distribution is not independent of late maturity traits. This review provides a summary of the most recent important revelations on the genomic position and cloning of Rpi genes, and the structure, associations, mode of action and activity spectrum of Rpi and corresponding avirulence (Avr) proteins. Practical implications for research into and application of Rpi genes are deduced and combined with an outlook on approaches to address remaining issues and interesting questions. It is evident that the potential of Rpi genes has not been exploited fully.

Identification of immunity related genes to study the Physalis peruviana--Fusarium oxysporum pathosystem

The Cape gooseberry (Physalisperuviana L) is an Andean exotic fruit with high nutritional value and appealing medicinal properties. However, its cultivation faces important phytosanitary problems mainly due to pathogens like Fusarium oxysporum, Cercosporaphysalidis and Alternaria spp. Here we used the Cape gooseberry foliar transcriptome to search for proteins that encode conserved domains related to plant immunity including: NBS (Nucleotide Binding Site), CC (Coiled-Coil), TIR (Toll/Interleukin-1 Receptor). We identified 74 immunity related gene candidates in P. peruviana which have the typical resistance gene (R-gene) architecture, 17 Receptor like kinase (RLKs) candidates related to PAMP-Triggered Immunity (PTI), eight (TIR-NBS-LRR, or TNL) and nine (CC–NBS-LRR, or CNL) candidates related to Effector-Triggered Immunity (ETI) genes among others. These candidate genes were categorized by molecular function (98%), biological process (85%) and cellular component (79%) using gene ontology. Some of the most interesting predicted roles were those associated with binding and transferase activity. We designed 94 primers pairs from the 74 immunity-related genes (IRGs) to amplify the corresponding genomic regions on six genotypes that included resistant and susceptible materials. From these, we selected 17 single band amplicons and sequenced them in 14 F. oxysporum resistant and susceptible genotypes. Sequence polymorphisms were analyzed through preliminary candidate gene association, which allowed the detection of one SNP at the PpIRG-63 marker revealing a nonsynonymous mutation in the predicted LRR domain suggesting functional roles for resistance.

Functional differentiation in the leucine-rich repeat domains of closely related plant virus-resistance proteins that recognize common avr proteins

The N' gene of Nicotiana sylvestris and L genes of Capsicum plants confer the resistance response accompanying the hypersensitive response (HR) elicited by tobamovirus coat proteins (CP) but with different viral specificities. Here, we report the identification of the N' gene. We amplified and cloned an N' candidate using polymerase chain reaction primers designed from L gene sequences. The N' candidate gene was a single 4143 base pairs fragment encoding a coiled-coil nucleotide-binding leucine-rich repeat (LRR)-type resistance protein of 1,380 amino acids. The candidate gene induced the HR in response to the coexpression of tobamovirus CP with the identical specificity as reported for N'. Analysis of N'-containing and tobamovirus-susceptible N. tabacum accessions supported the hypothesis that the candidate is the N' gene itself. Chimera analysis between N' and L(3) revealed that their LRR domains determine the spectrum of their tobamovirus CP recognition. Deletion and mutation analyses of N' and L(3) revealed that the conserved sequences in their C-terminal regions have important roles but contribute differentially to the recognition of common avirulence proteins. The results collectively suggest that Nicotiana N' and Capsicum L genes, which most likely evolved from a common ancestor, differentiated in their recognition specificity through changes in the structural requirements for LRR function.

Leucine-Rich Repeat (LRR) Domains Containing Intervening Motifs in Plants

LRRs (leucine rich repeats) are present in over 14,000 proteins. Non-LRR, island regions (IRs) interrupting LRRs are widely distributed. The present article reviews 19 families of LRR proteins having non-LRR IRs (LRR@IR proteins) from various plant species. The LRR@IR proteins are LRR-containing receptor-like kinases (LRR-RLKs), LRR-containing receptor-like proteins (LRR-RLPs), TONSOKU/BRUSHY1, and MJK13.7; the LRR-RLKs are homologs of TMK1/Rhg4, BRI1, PSKR, PSYR1, Arabidopsis At1g74360, and RPK2, while the LRR-RLPs are those of Cf-9/Cf-4, Cf-2/Cf-5, Ve, HcrVf, RPP27, EIX1, clavata 2, fascinated ear2, RLP2, rice Os10g0479700, and putative soybean disease resistance protein. The LRRs are intersected by single, non-LRR IRs; only the RPK2 homologs have two IRs. In most of the LRR-RLKs and LRR-RLPs, the number of repeat units in the preceding LRR block (N₁) is greater than the number of the following block (N₂); N₁ » N₂ in which N₁ is variable in the homologs of individual families, while N₂ is highly conserved. The five families of the LRR-RLKs except for the RPK2 family show N₁ = 8 − 18 and N₂ = 3 − 5. The nine families of the LRR-RLPs show N₁ = 12 − 33 and N₂ = 4; while N₁ = 6 and N₂ = 4 for the rice Os10g0479700 family and the N₁ = 4 − 28 and N₂ = 4 for the soybean protein family. The rule of N₁ » N₂ might play a common, significant role in ligand interaction, dimerization, and/or signal transduction of the LRR-RLKs and the LRR-RLPs. The structure and evolution of the LRR domains with non-LRR IRs and their proteins are also discussed.

Evolutionary meta-analysis of solanaceous resistance gene and solanum resistance gene analog sequences and a practical framework for cross-species comparisons

Cross-species comparative genomics approaches have been employed to map and clone many important disease resistance (R) genes from Solanum species-especially wild relatives of potato and tomato. These efforts will increase with the recent release of potato genome sequence and the impending release of tomato genome sequence. Most R genes belong to the prominent nucleotide binding site-leucine rich repeat (NBS-LRR) class and conserved NBS-LRR protein motifs enable survey of the R gene space of a plant genome by generation of resistance gene analogs (RGA), polymerase chain reaction fragments derived from R genes. We generated a collection of 97 RGA from the disease-resistant wild potato S. bulbocastanum, complementing smaller collections from other Solanum species. To further comparative genomics approaches, we combined all known Solanum RGA and cloned solanaceous NBS-LRR gene sequences, nearly 800 sequences in total, into a single meta-analysis. We defined R gene diversity bins that reflect both evolutionary relationships and DNA cross-hybridization results. The resulting framework is amendable and expandable, providing the research community with a common vocabulary for present and future study of R gene lineages. Through a series of sequence and hybridization experiments, we demonstrate that all tested R gene lineages are of ancient origin, are shared between Solanum species, and can be successfully accessed via comparative genomics approaches.

A novel protein RLS1 with NB-ARM domains is involved in chloroplast degradation during leaf senescence in rice

Leaf senescence, a type of programmed cell death (PCD) characterized by chlorophyll degradation, is important to plant growth and crop productivity. It emerges that autophagy is involved in chloroplast degradation during leaf senescence. However, the molecular mechanism(s) involved in the process is not well understood. In this study, the genetic and physiological characteristics of the rice rls1 (rapid leaf senescence 1) mutant were identified. The rls1 mutant developed small, yellow-brown lesions resembling disease scattered over the whole surfaces of leaves that displayed earlier senescence than those of wild-type plants. The rapid loss of chlorophyll content during senescence was the main cause of accelerated leaf senescence in rls1. Microscopic observation indicated that PCD was misregulated, probably resulting in the accelerated degradation of chloroplasts in rls1 leaves. Map-based cloning of the RLS1 gene revealed that it encodes a previously uncharacterized NB (nucleotide-binding site)-containing protein with an ARM (armadillo) domain at the carboxyl terminus. Consistent with its involvement in leaf senescence, RLS1 was up-regulated during dark-induced leaf senescence and down-regulated by cytokinin. Intriguingly, constitutive expression of RLS1 also slightly accelerated leaf senescence with decreased chlorophyll content in transgenic rice plants. Our study identified a previously uncharacterized NB-ARM protein involved in PCD during plant growth and development, providing a unique tool for dissecting possible autophagy-mediated PCD during senescence in plants.

Molecular evolution of a family of resistance gene analogs of nucleotide-binding site sequences in Solanum lycopersicum

Nucleotide-binding site-leucine-rich repeats (NBS-LRR) gene families are one of the major plant resistance genes. Genomic NBS evolution was studied in many plant species for diverse arrays of NBS gene families. In this study, we focused on one family of NBS sequences in an attempt to understand how closely related NBS sequences evolved in the light of selection in domesticated plant species. A phylogenetic analysis revealed five major clades (A-E) and five subclades (A1-A5) within clade A of cloned NBS sequences. Positive selection was only detected in newly evolved NBS lineages in subclades of clade A. Positively selected codon sites were found among NBS sequences of clade A. A sliding-window analysis revealed that regions with Ka/Ks ratios of >1 were in the inter-motifs when paired clades were compared, but regions with Ka/Ks ratios of >1 were found across NBS sequences when subclades of clade A were compared. Our results based on a family of closely related NBS sequences showed that positive selection was first exerted on specific lineages across all NBS sequences after selective constraints. Subsequently, sequences with mutations in commonly conserved motifs were scrutinized by purifying selection. In the long term, conserved high frequency alleles in commonly conserved motifs and changes in inter-motifs were maintained in the investigated family of NBS sequences. Moreover, codons identified to be under positive selection in the inter-motifs were mainly located in regions involved in functions of ATP binding or hydrolysis.

A nested leucine rich repeat (LRR) domain: the precursor of LRRs is a ten or eleven residue motif

Background

Leucine rich repeats (LRRs) are present in over 60,000 proteins that have been identified in viruses, bacteria, archae, and eukaryotes. All known structures of repeated LRRs adopt an arc shape. Most LRRs are 20-30 residues long. All LRRs contain LxxLxLxxNxL, in which "L" is Leu, Ile, Val, or Phe and "N" is Asn, Thr, Ser, or Cys and "x" is any amino acid. Seven classes of LRRs have been identified. However, other LRR classes remains to be characterized. The evolution of LRRs is not well understood.

Results

Here we describe a novel LRR domain, or nested repeat observed in 134 proteins from 54 bacterial species. This novel LRR domain has 21 residues with the consensus sequence of LxxLxLxxNxLxxLDLxx(N/L/Q/x)xx or LxxLxCxxNxLxxLDLxx(N/L/x)xx. This LRR domain is characterized by a nested periodicity; it consists of alternating 10- and 11- residues units of LxxLxLxxNx(x/-). We call it "IRREKO" LRR, since the Japanese word for "nested" is "IRREKO". The first unit of the "IRREKO" LRR domain is frequently occupied by an "SDS22-like" LRR with the consensus of LxxLxLxxNxLxxLxxLxxLxx or a "Bacterial" LRR with the consensus of LxxLxLxxNxLxxLPxLPxx. In some proteins an "SDS22-like" LRR intervenes between "IRREKO" LRRs.

Conclusion

Proteins having "IRREKO" LRR domain are almost exclusively found in bacteria. It is suggested that IRREKO@LRR evolved from a common ancestor with "SDS22-like" and "Bacterial" classes and that the ancestor of IRREKO@LRR is 10 or 11 residues of LxxLxLxxNx(x/-). The "IRREKO" LRR is predicted to adopt an arc shape with smaller curvature in which β-strands are formed on both concave and convex surfaces.

Analyses of non-leucine-rich repeat (non-LRR) regions intervening between LRRs in proteins

BACKGROUND

Many proteins have LRR (leucine-rich repeat) units interrupted by non-LRRs which we call IR (non-LRR island region).

METHODS

We identified proteins containing LRR@IRs (LRRs having IR) by using a new method and then analyzed their natures and distributions.

RESULTS

LRR@IR proteins were found in over two hundred proteins from prokaryotes and from eukaryotes. These are divided into twenty-one different protein families. The IRs occur one to four times in LRR regions and range in length from 5 to 11,265 residues. The IR lengths in Fungi adenylate cyclases (acys) range from 5 to 116 residues; there are 22 LRR repeats. The IRs in Leishmania proteophosphoglycans (ppgs) vary from 105 to 11,265 residues. These results indicate that the IRs evolved rapidly. A group of LRR@IR proteins-LRRC17, chondroadherin-like protein, ppgs, and four Pseudomonas proteins-have a super motif consisting of an LRR block and its adjacent LRR@IR region. This indicates that the entire super motif experienced duplication. The sequence analysis of IRs offers functional similarity in some LRR@IR protein families.

GENERAL SIGNIFICANCE

This study suggests that various IRs and super motifs provide a great variety of structures and functions for LRRs.

Silencing of the N family of resistance genes in Nicotiana edwardsonii compromises the hypersensitive response to tombusviruses

The nontarget effects associated with silencing of the N gene in Nicotiana edwardsonii, an amphidiploid species derived from N. glutinosa and N. clevelandii, have been characterized in this study. The N protein confers resistance to Tobacco mosaic virus (TMV), and is representative of a family of nucleotide-binding site leucine-rich repeat proteins present in N. glutinosa. Previous studies have shown that silencing of the N gene or of other plant genes associated with N-mediated defenses abolishes host resistance to TMV, and this effect can be measured through enhancements in movement or replication of TMV in the N-silenced plants. However, the nontarget effects of gene silencing have not been investigated thoroughly. Notably, are the functions of other resistance (R) genes also affected in experiments designed to silence the N gene? To investigate whether heterologous sequences could silence the N gene, we selected an R gene homolog from N. glutinosa that differed from the N gene by approximately 17%, created a hairpin transgene, and developed transgenic N. edwardsonii plants. Expression of this hairpin in the transgenic N. edwardsonii plants compromised the hypersensitive response to TMV, demonstrating that a single hairpin transgene could silence a block of R genes related by sequence similarity. We then investigated whether the response of N-silenced plants to other viruses would be altered, and found that the hypersensitive response triggered against the tombusviruses Tomato bushy stunt virus and Cymbidium ringspot virus also was compromised. This study indicates that a Tombusvirus R gene shares some homology with the N gene, which could facilitate the cloning of this gene.