Dacheng Tian

Fitness costs of R-gene-mediated resistance in Arabidopsis thaliana.

Resistance genes (R-genes) act as an immune system in plants by recognizing pathogens and inducing defensive pathways. Many R-gene loci are present in plant genomes, presumably reflecting the need to maintain a large repertoire of resistance alleles. These loci also often segregate for resistance and susceptibility alleles that natural selection has maintained as polymorphisms within a species for millions of years. Given the obvious advantage to an individual of being disease resistant, what prevents these resistance alleles from being driven to fixation by natural selection? A cost of resistance is one potential explanation; most models require a lower fitness of resistant individuals in the absence of pathogens for long-term persistence of susceptibility alleles. Here we test for the presence of a cost of resistance at the RPM1 locus of Arabidopsis thaliana. Results of a field experiment comparing the fitness of isogenic strains that differ in the presence or absence of RPM1 and its natural promoter reveal a large cost of RPM1, providing the first evidence that costs contribute to the maintenance of an ancient R-gene polymorphism.

Genome-wide identification of NBS genes in japonica rice reveals significant expansion of divergent non-TIR NBS-LRR genes

A complete set of candidate disease resistance ( R) genes encoding nucleotide-binding sites (NBSs) was identified in the genome sequence of japonica rice ( Oryza sativa L. var. Nipponbare). These putative R genes were characterized with respect to structural diversity, phylogenetic relationships and chromosomal distribution, and compared with those in Arabidopsis thaliana. We found 535 NBS-coding sequences, including 480 non-TIR (Toll/IL-1 receptor) NBS-LRR (Leucine Rich Repeat) genes. TIR NBS-LRR genes, which are common in A. thaliana, have not been identified in the rice genome. The number of non-TIR NBS-LRR genes in rice is 8.7 times higher than that in A. thaliana, and they account for about 1% of all of predicted ORFs in the rice genome. Some 76% of the NBS genes were located in 44 gene clusters or in 57 tandem arrays, and 16 apparent gene duplications were detected in these regions. Phylogenetic analyses based both NBS and N-terminal regions classified the genes into about 200 groups, but no deep clades were detected, in contrast to the two distinct clusters found in A. thaliana. The structural and genetic diversity that exists among NBS-LRR proteins in rice is remarkable, and suggests that diversifying selection has played an important role in the evolution of R genes in this agronomically important species. (Supplemental material is available online at http://gattaca.nju.edu.cn.)

Signature of balancing selection in Arabidopsis

Natural selection and genetic linkage cause DNA segments to have genealogical histories resembling those of the selected sites. When a polymorphism maintained by selection is old, it will have an island of enhanced sequence variability surrounding it, which represents a detectable “signature of selection.” We investigate the structure of single-nucleotide polymorphisms (SNPs) in a 20-kb interval containing the Arabidopsis thaliana disease resistance gene RPS5, a locus containing common alleles for the presence/absence of the entire locus. The alleles are considerably diverged at surrounding sites, indicative of an old polymorphism maintained by selection. The island of “enhanced” variability extends several kilobases to either side of the RPS5 deletion junction, and these SNPs are in nearly complete linkage disequilibrium with the RPS5 insertion/deletion. At a distance of 10 kb to either side of the locus, however, we find low levels of polymorphism and the absence of linkage disequilibrium between individual SNPs and RPS5 alleles. Our results show that the interval of enhanced variability surrounding this balanced polymorphism in Arabidopsis is large enough to be readily detected, but small enough to span the focal gene and few others. For this species it should be possible to identify the complete set of genes with long-lived polymorphisms, a potentially important subset of genes segregating for functional variants.

Recent duplications dominate NBS-encoding gene expansion in two woody species

Most disease resistance genes in plants encode NBS-LRR proteins. However, in woody species, little is known about the evolutionary history of these genes. Here, we identified 459 and 330 respective NBS-LRRs in grapevine and poplar genomes. We subsequently investigated protein motif composition, phylogenetic relationships and physical locations. We found significant excesses of recent duplications in perennial species, compared with those of annuals, represented by rice and Arabidopsis. Consequently, we observed higher nucleotide identity among paralogs and a higher percentage of NBS-encoding genes positioned in numerous clusters in the grapevine and poplar. These results suggested that recent tandem duplication played a major role in NBS-encoding gene expansion in perennial species. These duplication events, together with a higher probability of recombination revealed in this study, could compensate for the longer generation time in woody perennial species e.g. duplication and recombination could serve to generate novel resistance specificities. In addition, we observed extensive species-specific expansion in TIR-NBS-encoding genes. Non-TIR-NBS-encoding genes were poly- or paraphyletic, i.e. genes from three or more plant species were nested in different clades, suggesting different evolutionary patterns between these two gene types.

Single-nucleotide mutation rate increases close to insertions/deletions in eukaryotes

Mutation hotspots are commonly observed in genomic sequences and certain human disease loci, but general mechanisms for their formation remain elusive. Here we investigate the distribution of single-nucleotide changes around insertions/deletions (indels) in six independent genome comparisons, including primates, rodents, fruitfly, rice and yeast. In each of these genomic comparisons, nucleotide divergence (D) is substantially elevated surrounding indels and decreases monotonically to near-background levels over several hundred bases. D is significantly correlated with both size and abundance of nearby indels. In comparisons of closely related species, derived nucleotide substitutions surrounding indels occur in significantly greater numbers in the lineage containing the indel than in the one containing the ancestral (non-indel) allele; the same holds within species for single-nucleotide mutations surrounding polymorphic indels. We propose that heterozygosity for an indel is mutagenic to surrounding sequences, and use yeast genome-wide polymorphism data to estimate the increase in mutation rate. The consistency of these patterns within and between species suggests that indel-associated substitution is a general mutational mechanism.

Tracing the origin and evolutionary history of plant nucleotide‐binding site–leucine‐rich repeat (NBS‐LRR) genes

Plant disease resistance genes (R genes) encode proteins that function to monitor signals indicating pathogenic infection, thus playing a critical role in the plants defense system. Although many studies have been performed to explore the functional details of these important genes, their origin and evolutionary history remain unclear. In this study, focusing on the largest group of R genes, the nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes, we conducted an extensive genome-wide survey of 38 representative model organisms and obtained insights into the evolutionary stage and timing of NBS-LRR genes. Our data show that the two major domains, NBS and LRR, existed before the split of prokaryotes and eukaryotes but their fusion was observed only in land plant lineages. The Toll/interleukin-1 receptor (TIR) class of NBS-LRR genes probably had an earlier origin than its nonTIR counterpart. The similarities of the innate immune systems of plants and animals are likely to have been shaped by convergent evolution after their independent origins. Our findings start to unravel the evolutionary history of these important genes from the perspective of comparative genomics and also highlight the important role of reorganizing pre-existing building blocks in generating evolutionary novelties.

Unique Evolutionary Mechanism in R-Genes Under the Presence/Absence Polymorphism in Arabidopsis thaliana

While the presence/absence polymorphism is commonly observed in disease resistance (R-) genes in Arabidopsis, only a few R-genes under the presence/absence polymorphism (R-P/A) have been investigated. To understand the mechanism of the molecular evolution of R-P/A, we investigated genetic variation of nine R-P/A in A. thaliana from worldwide populations. The number of possessed R-genes varied widely among accessions (two to nine, on average 4.3 ± 1.6/accession). No pair of accessions shared the same haplotype, and no clear geographic differentiation was observed with respect to the pattern of presence/absence of the R-genes investigated. Presence allele frequencies also varied among loci (25–70%), and no linkage disequilibrium was detected among them. Although the LRR region in regular R-genes is known to be highly polymorphic and has a high Ka/Ks ratio in A. thaliana, nucleotide sequences of this region in the R-P/A showed a relatively low level of genetic variation (π = 0.0002–0.016) and low Ka/Ks (0.03–0.70, <1). In contrast, the nucleotide diversities around the deletion junction of R-P/A were constantly high between presence and absence accessions for the R-genes (Dxy = 0.031–0.103). Our results suggest that R-P/A loci evolved differently from other R-gene loci and that balancing selection plays an important role in molecular evolution of R-P/A.

Patterns of exon-intron architecture variation of genes in eukaryotic genomes

BackgroundThe origin and importance of exon-intron architecture comprises one of the remaining mysteries of gene evolution. Several studies have investigated the variations of intron length, GC content, ordinal position in a gene and divergence. However, there is little study about the structural variation of exons and introns.ResultsWe investigated the length, GC content, ordinal position and divergence in both exons and introns of 13 eukaryotic genomes, representing plant and animal. Our analyses revealed that three basic patterns of exon-intron variation were present in nearly all analyzed genomes (P < 0.001 in most cases): an ordinal reduction of length and divergence in both exon and intron, a co-variation between exon and its flanking introns in their length, GC content and divergence, and a decrease of average exon (or intron) length, GC content and divergence as the total exon numbers of a gene increased. In addition, we observed that the shorter introns had either low or high GC content, and the GC content of long introns was intermediate.ConclusionAlthough the factors contributing to these patterns have not been identified, our results provide three important clues: common factor(s) exist and may shape both exons and introns; the ordinal reduction patterns may reflect a time-orderly evolution; and the larger first and last exons may be splicing-required. These clues provide a framework for elucidating mechanisms involved in the organization of eukaryotic genomes and particularly in building exon-intron structures.

Genome-wide investigation on the genetic variations of rice disease resistance genes

Exploitation of plant disease resistance (R) gene in breeding programs has been proven to be the most efficient strategy for coping with the threat of pathogens. An understanding of R-gene variation is the basis for this strategy. Here we report a genome-wide investigation on the variation of NBS-LRR-encoding genes, the common type of R genes, between two sequenced rice genomes, Oryza sativa L. var. Nipponbare and 93–11. We show that the allelic nucleotide diversity in 65.0% of 397 least-divergent pairs is not high (0.344% on average), while the remaining 35% display a greater diversity (5.4% on average). The majority of conserved R genes is single-copy and/or located as a singleton. The clustered, particularly the complex-clustered, R-genes contribute greatly to the rich genetic variation. Surprisingly only 11.2% of R-genes have remarkably high ratios of non-synonymous to synonymous rates, which is much less than the 17.4% observed between Arabidopsis genomes. Noticeable “artificially selective sweeping” could be detected in a large proportion of the conserved R-genes, a scenario described in the “arms race” co-evolutionary model. Based on our study, a variation pattern of R-genes is proposed and confirmed by the analysis of R-genes from other rice lines, indicating that the observed variation pattern may be common in all rice lines.

Variation in the Ratio of Nucleotide Substitution and Indel Rates across Genomes in Mammals and Bacteria

Rates of nucleotide substitution and insertion/deletion (indel) are known to vary across the functional components of a genome. Little attention has been paid, however, to the quantitative relationship between the two. Here we investigate the ratio of nucleotide substitutions to indels (S/I) in different regions of 4 primates, 70 bacteria, and 8 other genomes. We find that the ratio differs at 5.4-times between coding and noncoding, 3.3-times between conserved and less conserved coding sequences, and 1.46-times between nonrepeat and repeat regions. The S/I ratio is also positively correlated with the level of divergence between the genomes compared. Our results suggest that the S/I ratio may reflect differences in the efficacy of selection against indels. Due to the sensitivity of indel density in different regions, this ratio varies over a much larger range. With the recent discovery suggesting that indels act as local enhancers of mutation in surrounding sequences, nucleotide substitution rates are expected to be accelerated in regions of low constraint, where indels tend to accumulate, but will otherwise be modulated in proportion to the level of a sequences functional constraint. Indels, therefore, may play a nontrivial role in controlling differences in genetic variation and divergence across functional regions of a genome.