Jian-Qun Chen

Exogenous plant MIR168a specifically targets mammalian LDLRAP1: evidence of cross-kingdom regulation by microRNA

Our previous studies have demonstrated that stable microRNAs (miRNAs) in mammalian serum and plasma are actively secreted from tissues and cells and can serve as a novel class of biomarkers for diseases, and act as signaling molecules in intercellular communication. Here, we report the surprising finding that exogenous plant miRNAs are present in the sera and tissues of various animals and that these exogenous plant miRNAs are primarily acquired orally, through food intake. MIR168a is abundant in rice and is one of the most highly enriched exogenous plant miRNAs in the sera of Chinese subjects. Functional studies in vitro and in vivo demonstrated that MIR168a could bind to the human/mouse low-density lipoprotein receptor adapter protein 1 (LDLRAP1) mRNA, inhibit LDLRAP1 expression in liver, and consequently decrease LDL removal from mouse plasma. These findings demonstrate that exogenous plant miRNAs in food can regulate the expression of target genes in mammals.

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.)

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.

Parent-progeny sequencing indicates higher mutation rates in heterozygotes

Mutation rates vary within genomes, but the causes of this remain unclear. As many prior inferences rely on methods that assume an absence of selection, potentially leading to artefactual results, we call mutation events directly using a parent–offspring sequencing strategy focusing on Arabidopsis and using rice and honey bee for replication. Here we show that mutation rates are higher in heterozygotes and in proximity to crossover events. A correlation between recombination rate and intraspecific diversity is in part owing to a higher mutation rate in domains of high recombination/diversity. Implicating diversity per se as a cause, we find an ∼3.5-fold higher mutation rate in heterozygotes than in homozygotes, with mutations occurring in closer proximity to heterozygous sites than expected by chance. In a genome that is a patchwork of heterozygous and homozygous domains, mutations occur disproportionately more often in the heterozygous domains. If segregating mutations predispose to a higher local mutation rate, clusters of genes dominantly under purifying selection (more commonly homozygous) and under balancing selection (more commonly heterozygous), might have low and high mutation rates, respectively. Our results are consistent with this, there being a ten times higher mutation rate in pathogen resistance genes, expected to be under positive or balancing selection. Consequently, we do not necessarily need to evoke extremely weak selection on the mutation rate to explain why mutational hot and cold spots might correspond to regions under positive/balancing and purifying selection, respectively.