Point mutations with positive selection were a major force during the evolution of a receptor-kinase resistance gene family of rice

The Role of Liriodendron Dof Gene Family in Abiotic Stress Response

The DOF (DNA-binding with one finger) transcription factors are exclusive to plants and play crucial roles in plant growth, development, and environmental adaptation. Although extensive research has been conducted on the Dof gene family in Arabidopsis, maize, and Solanum, investigations concerning the role of this gene family in Liriodendron remain unreported, leaving its biological function largely unknown. In this study, we performed a comprehensive genome-wide identification of the Dof gene family based on the Liriodendron genome, resulting in the discovery of a total of 17 LcDof gene members. Based on the results of phylogenetic analysis, the 17 LcDof proteins were classified into eight subfamilies. The motif analysis revealed the diverse nature of motifs within the D1 subfamily, which includes a distinct type of Dof transcription factor known as CDF (Cycling Dof Factor). We further characterized the chromosomal distribution, gene structure, conserved protein motifs, and cis-elements in the promoter regions. Additionally, utilizing transcriptome data from Liriodendron hybrids and conducting RT-qPCR experiments, we investigated the expression patterns of LhDofs under various abiotic stresses such as drought, cold, and heat stress. Notably, we found that several LhDofs, particularly LhDof4 and LhDof6, were significantly upregulated in response to abiotic stress. Furthermore, we cloned LhDof4 and LhDof6 genes and found that its encoding protein was mainly located in the nucleus by transient transformation in Liriodendron hybrids protoplast. Subsequently, we used LhDof6-overexpressing Liriodendron hybrid seedlings. We found that overexpression of LhDof6 enhanced the cold tolerance of the plants, increasing their survival rate at -20 °C. This result was further validated by changes in physiological indicators.

The First Genome-Wide Mildew Locus O Genes Characterization in the Lamiaceae Plant Family

Powdery mildew (PM) is a widespread plant disease that causes significant economic losses in thousands crops of temperate climates, including Lamiaceae species. Multiple scientific studies describe a peculiar form of PM-resistance associated at the inactivation of specific members of the Mildew Locus O (MLO) gene family, referred to as mlo-resistance. The characterization of Lamiaceae MLO genes, at the genomic level, would be a first step toward their potential use in breeding programs. We carried out a genome-wide characterization of the MLO gene family in 11 Lamiaceae species, providing a manual curated catalog of 324 MLO proteins. Evolutionary history and phylogenetic relationships were studied through maximum likelihood analysis and motif patter reconstruction. Our approach highlighted seven different clades diversified starting from an ancestral MLO domain pattern organized in 18 highly conserved motifs. In addition, 74 Lamiaceae putative PM susceptibility genes, clustering in clade V, were identified. Finally, we performed a codon-based evolutionary analysis, revealing a general high level of purifying selection in the eleven Lamiaceae MLO gene families, and the occurrence of few regions under diversifying selection in candidate susceptibility factors. The results of this work may help to address further biological questions concerning MLOs involved in PM susceptibility. In follow-up studies, it could be investigated whether the silencing or loss-of-function mutations in one or more of these candidate genes may lead to PM resistance.

Integrated Molecular and Bioinformatics Approaches for Disease-Related Genes in Plants

Modern plant pathology relies on bioinformatics approaches to create novel plant disease diagnostic tools. In recent years, a significant amount of biological data has been generated due to rapid developments in genomics and molecular biology techniques. The progress in the sequencing of agriculturally important crops has made it possible to develop a better understanding of plant-pathogen interactions and plant resistance. The availability of host-pathogen genome data offers effective assistance in retrieving, annotating, analyzing, and identifying the functional aspects for characterization at the gene and genome levels. Physical mapping facilitates the identification and isolation of several candidate resistance (R) genes from diverse plant species. A large number of genetic variations, such as disease-causing mutations in the genome, have been identified and characterized using bioinformatics tools, and these desirable mutations were exploited to develop disease resistance. Moreover, crop genome editing tools, namely the CRISPR (clustered regulatory interspaced short palindromic repeats)/Cas9 (CRISPR-associated) system, offer novel and efficient strategies for developing durable resistance. This review paper describes some aspects concerning the databases, tools, and techniques used to characterize resistance (R) genes for plant disease management.

Genome-wide characterization and expression analysis of the MLO gene family sheds light on powdery mildew resistance in Lagenaria siceraria

MLO (mildew locus O) genes play a vital role in plant disease defense system, especially powdery mildew (PM). Lagenaria siceraria is a distinct Cucurbitaceae crop, and PM is one of the most serious diseases threatening crop production and quality. Although MLOs have been exploited in many Cucurbitaceae species, genome-wide mining of MLO gene family in bottle gourd has not been surveyed yet. Here we identified 16 MLO genes in our recently assembled L. siceraria genome. A total of 343 unique MLO protein sequences from 20 species were characterized and compared to deduce a generally high level of purifying selection and the occurrence of regions related to candidate susceptibility factors in the evolutional divergence. LsMLOs were clustered in six clades containing seven conserved transmembrane domains and 10 clade-specific motifs along with deletion and variation. Three genes (LsMLO3, LsMLO6, and LsMLO13) in clade V showed high sequence identity with orthologues involved in PM susceptibility. The expression pattern of LsMLOs was tissue-specific but not cultivar-specific. Furthermore, it was indicated by qRT-PCR and RNA-seq that LsMLO3 and LsMLO13 were highly upregulated in response to PM stress. Subsequent sequence analysis revealed the structural deletion of LsMLO13 and a single nonsynonymous substitution of LsMLO3 in the PM-resistant genotype. Taken all together, it is speculated that LsMLO13 is likely a major PM susceptibility factor. The results of this study provide new insights into MLO family genes in bottle gourd and find a potential candidate S gene for PM tolerance breeding.

Molecular characterization of bacterial leaf streak resistance in hard winter wheat

Bacterial leaf streak (BLS) caused by Xanthomonas campestris pv. translucens is one of the major bacterial diseases threatening wheat production in the United States Northern Great Plains (NGP) region. It is a sporadic but widespread wheat disease that can cause significant loss in grain yield and quality. Identification and characterization of genomic regions in wheat that confer resistance to BLS will help track resistance genes/QTLs in future wheat breeding. In this study, we evaluated a hard winter wheat association mapping panel (HWWAMP) containing 299 hard winter wheat lines from the US hard winter wheat growing region for their reactions to BLS. We observed a range of BLS responses among the lines, importantly, we identified ten genotypes that showed a resistant reaction both in greenhouse and field evaluation. -Genome-wide association analysis with 15,990 SNPs was conducted using an exponentially compressed mixed linear model. Five genomic regions (p < 0.001) that regulate the resistance to BLS were identified on chromosomes 1AL, 1BS, 3AL, 4AL, and 7AS. The QTLs Q.bls.sdsu-1AL, Q.bls.sdsu-1BS, Q.bls.sdsu-3AL, Q.bls.sdsu-4AL, and Q.bls.sdsu-7AS explain a total of 42% of the variation. In silico analysis of sequences in the candidate regions on chromosomes 1AL, 1BS, 3AL, 4AL, and 7AS identified 10, 25, 22, eight, and nine genes, respectively with known plant defense-related functions. Comparative analysis with rice showed two syntenic regions in rice that harbor genes for bacterial leaf streak resistance. The ten BLS resistant genotypes and SNP markers linked to the QTLs identified in our study could facilitate breeding for BLS resistance in winter wheat.

Improvement of multiple agronomic traits by a disease resistance gene via cell wall reinforcement

The major disease resistance gene Xa4 confers race-specific durable resistance against Xanthomonas oryzae pv. oryzae, which causes the most damaging bacterial disease in rice worldwide. Although Xa4 has been one of the most widely exploited resistance genes in rice production worldwide, its molecular nature remains unknown. Here we show that Xa4, encoding a cell wall-associated kinase, improves multiple traits of agronomic importance without compromising grain yield by strengthening the cell wall via promoting cellulose synthesis and suppressing cell wall loosening. Strengthening of the cell wall by Xa4 enhances resistance to bacterial infection, and also increases mechanical strength of the culm with slightly reduced plant height, which may improve lodging resistance of the rice plant. The simultaneous improvement of multiple agronomic traits conferred by Xa4 may account for its widespread and lasting utilization in rice breeding programmes globally.

Mapping and Genetic Structure Analysis of the Anthracnose Resistance Locus Co-1HY in the Common Bean (Phaseolus vulgaris L.)

Anthracnose is a destructive disease of the common bean (Phaseolus vulgaris L.). The Andean cultivar Hongyundou has been demonstrated to possess strong resistance to anthracnose race 81. To study the genetics of this resistance, the Hongyundou cultivar was crossed with a susceptible genotype Jingdou. Segregation of resistance for race 81 was assessed in the F2 population and F2:3 lines under controlled conditions. Results indicate that Hongyundou carries a single dominant gene for anthracnose resistance. An allele test by crossing Hongyundou with another resistant cultivar revealed that the resistance gene is in the Co-1 locus (therefore named Co-1HY). The physical distance between this locus and the two flanking markers was 46 kb, and this region included four candidate genes, namely, Phvul.001G243500, Phvul.001G243600, Phvul.001G243700 and Phvul.001G243800. These candidate genes encoded serine/threonine-protein kinases. Expression analysis of the four candidate genes in the resistant and susceptible cultivars under control condition and inoculated treatment revealed that all the four candidate genes are expressed at significantly higher levels in the resistant genotype than in susceptible genotype. Phvul.001G243600 and Phvul.001G243700 are expressed nearly 15-fold and 90-fold higher in the resistant genotype than in the susceptible parent before inoculation, respectively. Four candidate genes will provide useful information for further research into the resistance mechanism of anthracnose in common bean. The closely linked flanking markers identified here may be useful for transferring the resistance allele Co-1HY from Hongyundou to elite anthracnose susceptible common bean lines.

Structure, evolution and functional inference on the Mildew Locus O (MLO) gene family in three cultivated Cucurbitaceae spp

BACKGROUND

The powdery mildew disease affects thousands of plant species and arguably represents the major fungal threat for many Cucurbitaceae crops, including melon (Cucumis melo L.), watermelon (Citrullus lanatus L.) and zucchini (Cucurbita pepo L.). Several studies revealed that specific members of the Mildew Locus O (MLO) gene family act as powdery mildew susceptibility factors. Indeed, their inactivation, as the result of gene knock-out or knock-down, is associated with a peculiar form of resistance, referred to as mlo resistance.

RESULTS

We exploited recently available genomic information to provide a comprehensive overview of the MLO gene family in Cucurbitaceae. We report the identification of 16 MLO homologs in C. melo, 14 in C. lanatus and 18 in C. pepo genomes. Bioinformatic treatment of data allowed phylogenetic inference and the prediction of several ortholog pairs and groups. Comparison with functionally characterized MLO genes and, in C. lanatus, gene expression analysis, resulted in the detection of candidate powdery mildew susceptibility factors. We identified a series of conserved amino acid residues and motifs that are likely to play a major role for the function of MLO proteins. Finally, we performed a codon-based evolutionary analysis indicating a general high level of purifying selection in the three Cucurbitaceae MLO gene families, and the occurrence of regions under diversifying selection in candidate susceptibility factors.

CONCLUSIONS

Results of this study may help to address further biological questions concerning the evolution and function of MLO genes. Moreover, data reported here could be conveniently used by breeding research, aiming to select powdery mildew resistant cultivars in Cucurbitaceae.

Prolonged Adaptive Evolution of a Defensive Gene in the Solanaceae

Although plants and their natural enemies may coevolve for prolonged periods, little is known about how long individual plant defensive genes are involved in the coevolutionary process. We address this issue by examining patterns of selection on the defensive gene threonine deaminase (TD). Tomato (Solanum lycopersicum) has two copies of this gene. One performs the canonical housekeeping function in amino acid metabolism of catalyzing the first reaction in the conversion of threonine to isoleucine. The second copy functions as an antinutritive defense against lepidopteran herbivores by depleting threonine in the insect gut. Wild tobacco (Nicotiana attenuata) also contains a defensive copy. We show that a single copy of TD underwent two or three duplications near the base of the Solanaceae. One copy retains the housekeeping function, whereas a second copy evolved defensive functions. Positive selection occurred on the branch of the TD2 gene tree subtending the common ancestor of the Nicotianoideae and Solanoideae. It also occurred within the Solanoideae clade but not within the Nicotianoideae clade. Finally, it occurred on most branches leading from the common ancestor to S. lycopersicum. Based on recent calibrations of the Solanaceae phylogeny, TD2 experienced adaptive substitutions for a period of 30-50 My. We suggest that the most likely explanation for this result is fluctuating herbivore abundances: When herbivores are rare, relaxed selection increases the likelihood that slightly disadvantageous mutations will be fixed by drift; when herbivores are common, increased selection causes the evolution of compensatory adaptive mutations. Alternative explanations are also discussed.

Disease Resistance Gene Analogs (RGAs) in Plants

Plants have developed effective mechanisms to recognize and respond to infections caused by pathogens. Plant resistance gene analogs (RGAs), as resistance (R) gene candidates, have conserved domains and motifs that play specific roles in pathogens' resistance. Well-known RGAs are nucleotide binding site leucine rich repeats, receptor like kinases, and receptor like proteins. Others include pentatricopeptide repeats and apoplastic peroxidases. RGAs can be detected using bioinformatics tools based on their conserved structural features. Thousands of RGAs have been identified from sequenced plant genomes. High-density genome-wide RGA genetic maps are useful for designing diagnostic markers and identifying quantitative trait loci (QTL) or markers associated with plant disease resistance. This review focuses on recent advances in structures and mechanisms of RGAs, and their identification from sequenced genomes using bioinformatics tools. Applications in enhancing fine mapping and cloning of plant disease resistance genes are also discussed.

Kinase domain-targeted isolation of defense-related receptor-like kinases (RLK/Pelle) in Platanus×acerifolia: phylogenetic and structural analysis

BACKGROUND

Plant receptor-like kinase (RLK/Pelle) family regulates growth and developmental processes and interaction with pathogens and symbionts.Platanaceae is one of the earliest branches of Eudicots temporally located before the split which gave rise to Rosids and Asterids. Thus investigations into the RLK family in Platanus can provide information on the evolution of this gene family in the land plants.Moreover RLKs are good candidates for finding genes that are able to confer resistance to Platanus pathogens.

RESULTS

Degenerate oligonucleotide primers targeting the kinase domain of stress-related RLKs were used to isolate for the first time 111 RLK gene fragments in Platanus×acerifolia. Sequences were classified as candidates of the following subfamilies: CrRLK1L, LRR XII, WAK-like, and LRR X-BRI1 group. All the structural features typical of the RLK kinase domain were identified, including the non-RD motif which marks potential pathogen recognition receptors (PRRs). The LRR XII candidates, whose counterpart in Arabidopsis and rice comprises non-RD PRRs, were mostly non-RD kinases, suggesting a group of PRRs. Region-specific signatures of a relaxed purifying selection in the LRR XII candidates were also found, which is novel for plant RLK kinase domain and further supports the role of LRR XII candidates as PRRs. As we obtained CrRLK1L candidates using primers designed on Pto of tomato, we analysed the phylogenetic relationship between CrRLK1L and Pto-like of plant species. We thus classified all non-solanaceous Pto-like genes as CrRLK1L and highlighted for the first time the close phylogenetic vicinity between CrRLK1L and Pto group. The origins of Pto from CrRLK1L is proposed as an evolutionary mechanism.

CONCLUSIONS

The signatures of relaxed purifying selection highlight that a group of RLKs might have been involved in the expression of phenotypic plasticity and is thus a good candidate for investigations into pathogen resistance.Search of Pto-like genes in Platanus highlighted the close relationship between CrRLK1L and Pto group. It will be exciting to verify if sensu strictu Pto are present in taxonomic groups other than Solanaceae, in order to further clarify the evolutionary link with CrRLK1L.We obtained a first valuable resource useful for an in-depth study on stress perception systems.

Diversity in receptor-like kinase genes is a major determinant of quantitative resistance to Fusarium oxysporum f.sp. matthioli

Resistance to wilt fungus Fusarium oxysporum f.sp. matthioli (FOM) is a polygenic trait in Arabidopsis thaliana. RFO3 is one of six quantitative trait loci accounting for the complete resistance of accession Columbia-0 (Col-0) and susceptibility of accession Taynuilt-0 (Ty-0). We find that Col-0 and Ty-0 alleles of RFO3 are representative of two common variants in wild Arabidopsis accessions, that resistance and susceptibility to FOM are ancestral features of the two variants and that resistance from RFO3 is unrivalled by other genes in a genome-wide survey of diversity in accessions. A single receptor-like kinase (RLK) gene in Col-0 is responsible for the resistance of RFO3, although the susceptible Ty-0 allele codes for two RLK homologs. Expression of RFO3 is highest in vascular tissue, which F. oxysporum infects, and root-expressed RFO3 restricts FOM infection of the vascular system. RFO3 confers specific resistance to FOM and provides no resistance to two other crucifer-infecting F. oxysporum pathogens. RFO3's identity, expression and specificity suggest that RFO3 represents diversity in pattern-recognition receptor (PRR) genes. The characteristics of RFO3 and the previously published RFO1 suggest that diversity in RLK PRRs is a major determinant of quantitative resistance in wild plant populations.

Involvement of OsJAZ8 in jasmonate-induced resistance to bacterial blight in rice

The plant hormone jasmonic acid (JA) has a crucial role in both host immunity and development in plants. Here, we report the importance of JA signaling in the defense system of rice. Exogenous application of JA conferred resistance to bacterial blight caused by Xanthomonas oryzae pv. oryzae (Xoo) in rice. Expression of OsJAZ8, a rice jasmonate ZIM-domain protein, was highly up-regulated by JA. OsJAZ8 interacted with a putative OsCOI1, which is a component of the SCF(COI1) E3 ubiquitin ligase complex, in a coronatine-dependent manner. OsJAZ8 also formed heterodimers with other OsJAZ proteins but did not form homodimer. JA treatment caused OsJAZ8 degradation and this degradation was dependent on the 26S proteasome pathway. Furthermore, the JA-dependent OsJAZ8 degradation was mediated by the Jas domain. Transgenic rice plants overexpressing OsJAZ8ΔC, which lacks the Jas domain, exhibited a JA-insensitive phenotype. A large-scale analysis using a rice DNA microarray revealed that overexpression of OsJAZ8ΔC altered the expression of JA-responsive genes, including defense-related genes, in rice. Furthermore, OsJAZ8ΔC negatively regulated the JA-induced resistance to Xoo in rice. On the basis of these data, we conclude that JA plays an important role in resistance to Xoo, and OsJAZ8 acts as a repressor of JA signaling in rice.

Contrasted patterns of selective pressure in three recent paralogous gene pairs in the Medicago genus (L.)

Background

Gene duplications are a molecular mechanism potentially mediating generation of functional novelty. However, the probabilities of maintenance and functional divergence of duplicated genes are shaped by selective pressures acting on gene copies immediately after the duplication event. The ratio of non-synonymous to synonymous substitution rates in protein-coding sequences provides a means to investigate selective pressures based on genic sequences. Three molecular signatures can reveal early stages of functional divergence between gene copies: change in the level of purifying selection between paralogous genes, occurrence of positive selection, and transient relaxed purifying selection following gene duplication. We studied three pairs of genes that are known to be involved in an interaction with symbiotic bacteria and were recently duplicated in the history of the Medicago genus (Fabaceae). We sequenced two pairs of polygalacturonase genes (Pg11-Pg3 and Pg11a-Pg11c) and one pair of auxine transporter-like genes (Lax2-Lax4) in 17 species belonging to the Medicago genus, and sought for molecular signatures of differentiation between copies.

Results

Selective histories revealed by these three signatures of molecular differentiation were found to be markedly different between each pair of paralogs. We found sites under positive selection in the Pg11 paralogs while Pg3 has mainly evolved under purifying selection. The most recent paralogs examined Pg11a and Pg11c, are both undergoing positive selection and might be acquiring new functions. Lax2 and Lax4 paralogs are both under strong purifying selection, but still underwent a temporary relaxation of purifying selection immediately after duplication.

Conclusions

This study illustrates the variety of selective pressures undergone by duplicated genes and the effect of age of the duplication. We found that relaxation of selective constraints immediately after duplication might promote adaptive divergence.

Genetical genomics for evolutionary studies

Genetical genomics combines acquired high-throughput genomic data with genetic analysis. In this chapter, we discuss the application of genetical genomics for evolutionary studies, where new high-throughput molecular technologies are combined with mapping quantitative trait loci (QTL) on the genome in segregating populations.The recent explosion of high-throughput data--measuring thousands of proteins and metabolites, deep sequencing, chromatin, and methyl-DNA immunoprecipitation--allows the study of the genetic variation underlying quantitative phenotypes, together termed xQTL. At the same time, mining information is not getting easier. To deal with the sheer amount of information, powerful statistical tools are needed to analyze multidimensional relationships. In the context of evolutionary computational biology, a well-designed experiment may help dissect a complex evolutionary trait using proven statistical methods for associating phenotypical variation with genomic locations.Evolutionary expression QTL (eQTL) studies of the last years focus on gene expression adaptations, mapping the gene expression landscape, and, tentatively, eQTL networks. Here, we discuss the possibility of introducing an evolutionary prior, in the form of gene families displaying evidence of positive selection, and using that in the context of an eQTL experiment for elucidating host-pathogen protein-protein interactions. Through the example of an experimental design, we discuss the choice of xQTL platform, analysis methods, and scope of results. The resulting eQTL can be matched, resulting in putative interacting genes and their regulators. In addition, a prior may help distinguish QTL causality from reactivity, or independence of traits, by creating QTL networks.

Ortholog alleles at Xa3/Xa26 locus confer conserved race-specific resistance against Xanthomonas oryzae in rice

The rice disease resistance (R) gene Xa3/Xa26 (having also been named Xa3 and Xa26) against Xanthomonas oryzae pv. oryzae (Xoo), which causes bacterial blight disease, belongs to a multiple gene family clustered in chromosome 11 and is from an AA genome rice cultivar (Oryza sativa L.). This family encodes leucine-rich repeat (LRR) receptor kinase-type proteins. Here, we show that the orthologs (alleles) of Xa3/Xa26, Xa3/Xa26-2, and Xa3/Xa26-3, from wild Oryza species O. officinalis (CC genome) and O. minuta (BBCC genome), respectively, were also R genes against Xoo. Xa3/Xa26-2 and Xa3/Xa26-3 conferred resistance to 16 of the 18 Xoo strains examined. Comparative sequence analysis of the Xa3/Xa26 families in the two wild Oryza species showed that Xa3/Xa26-3 appeared to have originated from the CC genome of O. minuta. The predicted proteins encoded by Xa3/Xa26, Xa3/Xa26-2, and Xa3/Xa26-3 share 91-99% sequence identity and 94-99% sequence similarity. Transgenic plants carrying a single copy of Xa3/Xa26, Xa3/Xa26-2, or Xa3/Xa26-3, in the same genetic background, showed a similar resistance spectrum to a set of Xoo strains, although plants carrying Xa3/Xa26-2 or Xa3/Xa26-3 showed lower resistance levels than the plants carrying Xa3/Xa26. These results suggest that the Xa3/Xa26 locus predates the speciation of A and C genome, which is approximately 7.5 million years ago. Thus, the resistance specificity of this locus has been conserved for a long time.

A pair of orthologs of a leucine-rich repeat receptor kinase-like disease resistance gene family regulates rice response to raised temperature

BACKGROUND

Rice Xa3/Xa26 disease-resistance gene encodes a leucine-rich repeat (LRR) receptor kinase-type protein against Xanthomonas oryzae pv. oryzae (Xoo) and belongs to a multigene family. However, the functions of most genes in this family are unknown.

RESULTS

Here we report that two orthologs of this family, the NRKe from rice variety Nipponbare and 9RKe from variety 93-11 at the RKe locus, have similar functions although they encode different proteins. This pair of orthologs could not mediate resistance to Xoo, but they were transcriptionally induced by raised temperature. Transcriptional activation of NRKe or 9RKe resulted in the formation of temperature-sensitive lesion mimics, which were spots of dead cells associated with accumulation of superoxides, in different organs of the transgenic plants. These plants were more sensitive to high temperature shock than wild-type controls. Transgenic plants carrying a chimeric protein consisting of the LRR domain of NRKe and the kinase domain of Xa3/Xa26 developed the same lesion mimics as the NRKe-transgenic plants, whereas transgenic plants carrying another chimeric protein consisting of the LRR domain of Xa3/Xa26 and the kinase domain of NRKe were free of lesion mimic. All the transgenic plants carrying a chimeric protein were susceptible to Xoo.

CONCLUSION

These results suggest that the RKe locus is involved in rice response to raised temperature. The LRR domain of RKe protein appears to be important to sense increased temperature. The RKe-involved temperature-related pathway and Xa3/Xa26-mediated disease-resistance pathway may partially overlap.

Selection-driven divergence after gene duplication in Arabidopsis thaliana

Gene duplications are one of the most important mechanisms for the origin of evolutionary novelties. Even though various models of the fate of duplicated genes have been established, current knowledge about the role of divergent selection after gene duplication is rather limited. In this study, we analyzed sequence divergence in response to neo- and subfunctionalization of segmentally duplicated genes in the genome of Arabidopsis thaliana. We compared the genomes of A. thaliana and the poplar Populus trichocarpa to identify orthologous pairs of genes and their corresponding inparalogs. Maximum-likelihood analyses of the nonsynonymous and synonymous substitution rate ratio [Formula: see text] of pairs of A. thaliana inparalogs were used to detect differences in the evolutionary rates of protein coding sequences. We analyzed 1,924 A. thaliana paralogous pairs and our results indicate that around 6.9% show divergent ω values between the lineages for a fraction of sites. We observe an enrichment of regulatory sequences, a reduced level of co-expression and an increased number of substitutions that can be attributed to positive selection based on an McDonald-Kreitman type of analysis. Taken together, these results show that selection after duplication contributes substantially to gene novelties and hence functional divergence in plants.

Genome-wide identification, characterization and phylogenetic analysis of the rice LRR-kinases

LRR-kinases constitute the largest subfamily of receptor-like kinases in plants and regulate a wide variety of processes related to development and defense. Through a reiterative process of sequence analysis and re-annotation, we identified 309 LRR-kinase genes in the rice genome (Nipponbare). Among them, 127 genes in the Rice Annotation Project Database and 85 in Refseq of NCBI were amended (in addition, 62 LRR-kinase genes were not annotated in Refseq). The complete set of LRR-kinases was characterized. These LRR-kinases were classified into five groups according to phylogenetic analysis, and the genes in groups 1, 2, 3 and 4 usually have fewer introns than those in group 5. The introns in the LRR domain, which are highly conserved in regards to their positions and configurations, split the first Leu or other amino residues at this position of the ‘xxLxLxx’ motif with phase 2 and usually separate one or more LRR repeats exactly. Tandemly repeated LRR motifs have evolved from exon duplication, mutation and exon shuffling. The extensive distribution and diversity of the LRR-kinase genes have been mainly generated by tandem duplication and mutation after whole genome duplication. Positive selection has made a limited contribution to the sequence diversity after duplication, but positively selected sites located in the LRR domain are thought to involve in the protein-protein interaction.

A pair of allelic WRKY genes play opposite roles in rice-bacteria interactions

Although allelic diversity of genes has been reported to play important roles in different physiological processes, information on allelic diversity of defense-responsive genes in host-pathogen interactions is limited. Here, we report that a pair of allelic genes, OsWRKY45-1 and OsWRKY45-2, which encode proteins with a 10-amino acid difference, play opposite roles in rice (Oryza sativa) resistance against bacterial pathogens. Bacterial blight caused by Xanthomonas oryzae pv oryzae (Xoo), bacterial streak caused by Xanthomonas oryzae pv oryzicola (Xoc), and fungal blast caused by Magnaporthe grisea are devastating diseases of rice worldwide. OsWRKY45-1-overexpressing plants showed increased susceptibility and OsWRKY45-1-knockout plants showed enhanced resistance to Xoo and Xoc. In contrast, OsWRKY45-2-overexpressing plants showed enhanced resistance and OsWRKY45-2-suppressing plants showed increased susceptibility to Xoo and Xoc. Interestingly, both OsWRKY45-1- and OsWRKY45-2-overexpressing plants showed enhanced resistance to M. grisea. OsWRKY45-1-regulated Xoo resistance was accompanied by increased accumulation of salicylic acid and jasmonic acid and induced expression of a subset of defense-responsive genes, while OsWRKY45-2-regulated Xoo resistance was accompanied by increased accumulation of jasmonic acid but not salicylic acid and induced expression of another subset of defense-responsive genes. These results suggest that both OsWRKY45-1 and OsWRKY45-2 are positive regulators in rice resistance against M. grisea, but the former is a negative regulator and the latter is a positive regulator in rice resistance against Xoo and Xoc. The opposite roles of the two allelic genes in rice-Xoo interaction appear to be due to their mediation of different defense signaling pathways.

Dissection of the factors affecting development-controlled and race-specific disease resistance conferred by leucine-rich repeat receptor kinase-type R genes in rice

Development-controlled resistance and resistance specificity frequently restrict the application of a disease resistance (R) gene in crop breeding programs. Xa3/Xa26 and Xa21, encoding leucine-rich repeat (LRR)-kinase type plasma membrane proteins, mediate race-specific resistance to Xanthomonas oryzae pv. oryzae (Xoo), which causes bacterial blight, one of the most devastating rice diseases. Plants carrying Xa3/Xa26 and plants carrying Xa21 have different resistance spectra and the functions of the two R genes are regulated by developmental stage. Four chimeric genes encoding proteins consisting of different parts of XA3/XA26 and XA21 were constructed by domain swapping and transformed into a susceptible rice variety. The resistance spectra and development-regulated resistance of the transgenic plants carrying Xa3/Xa26, Xa21, or chimeric gene to different Xoo strains were analyzed in the same genetic background. The results suggest that the gradually increased expression of Xa3/Xa26 and Xa21 plays an important role in the progressively enhanced Xoo resistance during rice development. In addition, the LRR domains of XA3/XA26 and XA21 are important determinants of race-specific recognition during rice-Xoo interaction, but juxtamembrane regions of the two R proteins also appear to contribute to resistance specificity.

Molecular evolution of the CPP-like gene family in plants: insights from comparative genomics of Arabidopsis and rice

CPP-like genes are members of a small family which features the existence of two similar Cys-rich domains termed CXC domains in their protein products and are distributed widely in plants and animals but do not exist in yeast. The members of this family in plants play an important role in development of reproductive tissue and control of cell division. To gain insights into how CPP-like genes evolved in plants, we conducted a comparative phylogenetic and molecular evolutionary analysis of the CPP-like gene family in Arabidopsis and rice. The results of phylogeny revealed that both gene loss and species-specific expansion contributed to the evolution of this family in Arabidopsis and rice. Both intron gain and intron loss were observed through intron/exon structure analysis for duplicated genes. Our results also suggested that positive selection was a major force during the evolution of CPP-like genes in plants, and most amino acid residues under positive selection were disproportionately located in the region outside the CXC domains. Further analysis revealed that two CXC domains and sequences connecting them might have coevolved during the long evolutionary period.

Molecular evolution of the CPP-like gene family in plants: insights from comparative genomics of Arabidopsis and rice

CPP-like genes are members of a small family which features the existence of two similar Cys-rich domains termed CXC domains in their protein products and are distributed widely in plants and animals but do not exist in yeast. The members of this family in plants play an important role in development of reproductive tissue and control of cell division. To gain insights into how CPP-like genes evolved in plants, we conducted a comparative phylogenetic and molecular evolutionary analysis of the CPP-like gene family in Arabidopsis and rice. The results of phylogeny revealed that both gene loss and species-specific expansion contributed to the evolution of this family in Arabidopsis and rice. Both intron gain and intron loss were observed through intron/exon structure analysis for duplicated genes. Our results also suggested that positive selection was a major force during the evolution of CPP-like genes in plants, and most amino acid residues under positive selection were disproportionately located in the region outside the CXC domains. Further analysis revealed that two CXC domains and sequences connecting them might have coevolved during the long evolutionary period.

Evidence for positive selection in the TLR9 gene of teleosts

Toll-like receptors (TLRs) have been identified as key sensors of invading microbes by identifying pathogen-associated molecular patterns and activating innate immune responses. Whereas purifying selection has been suggested in mammalian TLR9, evolutionary features of TLR9 in teleosts have not been investigated in detail. We therefore analysed TLR9 DNA sequences of eight teleost species, including zebrafish (Danio rerio), Japanese flounder (Paralichthys olivaceus), pufferfish (Takifugu rubripes), and five seabreams. Eleven sites subjected to positive selection were identified using the codon-substitution models of PAML 3.15. Ten of these 11 sites were found to be associated with leucine-rich repeats (LRRs). Seven of these 10 positively selected sites were associated with the convex surface of the LRR solenoids, leading to variations of the structures of the LRRs possibly by the introduction of flexibility into the LRR solenoids. The positive selection of LRRs in TLR9 may indicate the adaptation of teleosts to different oligodeoxynucleotides present in different bacterial species.

Molecular adaptation and expression evolution following duplication of genes for organellar ribosomal protein S13 in rosids

BACKGROUND

Gene duplication has been a fundamental process in the evolution of eukaryotic genomes. After duplication one copy (or both) can undergo divergence in sequence, expression pattern, and function. Two divergent copies of the ribosomal protein S13 gene (rps13) of chloroplast origin are found in the nucleus of the rosids Arabidopsis, Gossypium, and Glycine. One encodes chloroplast-imported RPS13 (nucp rps13), while the other encodes mitochondria-imported RPS13 (numit rps13). The rps13 gene has been lost from mitochondrial DNA (mt rps13) of many rosids.

RESULTS

We studied sequence evolution of numit rps13 in comparison with nucp rps13 in seven rosid genera. Ka/Ks analysis and likelihood ratio tests showed considerably higher Ka values and Ka/Ks ratios in numit rps13 than in nucp rps13, indicating increased amino acid sequence divergence in numit rps13. Two positively selected codons were detected in numit RPS13 in regions that are inferred to interact with the 16S rRNA. Several amino acids in numit RPS13 have changed from the one present in nucp RPS13 to the one present in mt RPS13, showing that numit rps13 is becoming more like mt rps13. Comparison of expression patterns and levels of numit rps13 and nucp rps13 in Arabidopsis using microarray data indicated divergence in gene expression. We discovered that in addition to numit rps13, Malus (apple) contains a transcribed mt rps13 gene. To determine if partitioning of expression takes place between numit rps13 and mt rps13, expression of both copies and RNA editing of mt rps13 were examined by RT-PCR, qRT-PCR, and sequencing from 14 different organ types plus seedlings subjected to five different abiotic stresses. Co-expression of numit rps13 and mt rps13 was observed in all the organs and various stress treatments. We determined that purifying selection is acting on both numit rps13 and mt rps13 in Malus.

CONCLUSION

Our data provide evidence that numit rps13 genes in rosids have experienced adaptive sequence evolution and convergent evolution with mt rps13. Co-expression of numit rps13 and mt rps13 and purifying selection on both genes in Malus suggest that both are functional. The three organellar rps13 genes in rosids provide a distinctive case of gene duplication involving the co-evolution of the nuclear and cytoplasmic genomes.

Adaptive evolution of the symbiotic gene NORK is not correlated with shifts of rhizobial specificity in the genus Medicago

Background

The NODULATION RECEPTOR KINASE (NORK) gene encodes a Leucine-Rich Repeat (LRR)-containing receptor-like protein and controls the infection by symbiotic rhizobia and endomycorrhizal fungi in Legumes. The occurrence of numerous amino acid changes driven by directional selection has been reported in this gene, using a limited number of messenger RNA sequences, but the functional reason of these changes remains obscure. The Medicago genus, where changes in rhizobial associations have been previously examined, is a good model to test whether the evolution of NORK is influenced by rhizobial interactions.

Results

We sequenced a region of 3610 nucleotides (encoding a 392 amino acid-long region of the NORK protein) in 32 Medicago species. We confirm that positive selection in NORK has occurred within the Medicago genus and find that the amino acid positions targeted by selection occur in sites outside of solvent-exposed regions in LRRs, and other sites in the N-terminal region of the protein. We tested if branches of the Medicago phylogeny where changes of rhizobial symbionts occurred displayed accelerated rates of amino acid substitutions. Only one branch out of five tested, leading to M. noeana, displays such a pattern. Among other branches, the most likely for having undergone positive selection is not associated with documented shift of rhizobial specificity.

Conclusion

Adaptive changes in the sequence of the NORK receptor have involved the LRRs, but targeted different sites than in most previous studies of LRR proteins evolution. The fact that positive selection in NORK tends not to be associated to changes in rhizobial specificity indicates that this gene was probably not involved in evolving rhizobial preferences. Other explanations (e.g. coevolutionary arms race) must be tested to explain the adaptive evolution of NORK.

Functional analysis of Xa3/Xa26 family members in rice resistance to Xanthomonas oryzae pv. oryzae

Plant disease resistant (R) genes are frequently clustered in the genome. The diversity of members in a complex R-gene family may provide variation in resistance specificity. Rice Xa3/Xa26, conferring resistance to Xanthomonas oryzae pv. oryzae (Xoo) encodes a leucine-rich repeat (LRR) receptor kinase-type protein and belongs to a multigene family, consisting of Xa3/Xa26, MRKa, MRKc and MRKd in rice cultivar Minghui 63. MRKa and MRKc are intact genes, while MRKd is a pseudogene. Complementary analyses showed that MRKa and MRKc could not mediate resistance to Xoo when regulated by their native promoters, but MRKa not MRKc conferred partial resistance to Xoo when regulated by a strong constitutive promoter. Plants carrying truncated XA3/XA26, which lacked the kinase domain, were compromised in their resistance to Xoo. However, the kinase domain of MRKa could partially restore the function of the truncated XA3/XA26 in resistance. MRKa and MRKc showed similar expression pattern as Xa3/Xa26, which expressed only in the vascular systems of different tissues. The expressional characteristic of MRKa and MRKc perfectly fits the function of genes conferring resistance to Xoo, a vascular pathogen. These results suggest that although MRKa and MRKc cannot mediate bacterial blight resistance nowadays, they may be once effective genes for Xoo resistance. Their expressional characteristic and sequence similarity to Xa3/Xa26 will provide templates for generating novel recognition specificity to face the evolution of Xoo. In addition, both LRR and kinase domains encoded by Xa3/Xa26 and MRKa are the functional determinants and MRKa-mediated resistance is dosage-dependent.

Introgressed and endogenous Mi-1 gene clusters in tomato differ by complex rearrangements in flanking sequences and show sequence exchange and diversifying selection among homologues

Many plant disease resistance genes (R-genes) encode proteins characterized by the presence of a nucleotide-binding site (NBS) and a leucine-rich repeat (LRR) region and occur in clusters of related genes in plant genomes. One such gene, Mi-1, confers isolate-specific resistance against root-knot nematodes, aphids and whiteflies in cultivated tomato, Solanum lycopersicon. The DNA region carrying Mi-1 and six closely related sequences was introgressed into tomato from Solanum peruvianum in the 1940s. For both susceptible and resistant tomato, Mi-1 homologues are present in two clusters with 3 and 4 copies each on the short arm of chromosome 6. Two homologues from each source are pseudogenes, and one homologue from each source encodes a truncated product. DNA sequence identity among the homologues including the truncated genes, but excluding the pseudogenes, ranges from 92.9 to 96.7%. All the non-pseudogene homologues are transcribed. Comparison of homologues suggests that extensive sequence exchange has occurred. Regions of diversifying selection are present in the ARC2 domain of the NBS region and dispersed throughout the LRR region, suggesting that these regions are possible locations of specificity determinants. Other sequences in the introgressed region have similarity to the Arabidopsis auxin-receptor protein TIR1, a jumonji-like transcription factor and a Na(+)/H(+) antiporter. Analysis of sequences flanking the Mi-1-homologues reveals blocks of homology, but complex differences in arrangement of these blocks when susceptible and resistant genotypes are compared indicating that the region has undergone considerable rearrangement during evolution, perhaps contributing to evolution of specificity.

Xanthomonas oryzae pathovars: model pathogens of a model crop

SUMMARY Xanthomonas oryzae pv. oryzae and Xanthomonas oryzae pv. oryzicola cause bacterial blight and bacterial leaf streak of rice (Oryza sativa), which constrain production of this staple crop in much of Asia and parts of Africa. Tremendous progress has been made in characterizing the diseases and breeding for resistance. X. oryzae pv. oryzae causes bacterial blight by invading the vascular tissue, while X. oryzae pv. oryzicola causes bacterial leaf streak by colonizing the parenchyma. In rice there are 29 major genes for resistance to bacterial blight, but so far only a few quantitative resistance loci for bacterial leaf streak. Over 30 races of X. oryzae pv. oryzae have been reported. Both pathogens exhibit genetic variation among isolates. Mechanisms of pathogenesis and resistance have begun to be elucidated. Members of the AvrBs3/PthA family of transcription activator-like effectors play a major role in the virulence of X. oryzae pv. oryzae and possibly X. oryzae pv. oryzicola. Cloning of six rice resistance genes for bacterial blight and one from maize effective against bacterial leaf streak has uncovered a diversity of structure and function, some shared by genes involved in defence in animals. This article reviews research that spans a century. It also presents a perspective on challenges for sustainable control, and opportunities that interactions of X. oryzae pathovars with rice present as models for understanding fundamental aspects of bacterial pathogenesis of plants and plant disease resistance, as well as other aspects of plant and microbial biology, with implications also for animal innate immunity.