Sequence variation and selection of small RNAs in domesticated rice
- Yu Wang†1,
- Dan Shen†1,
- Shiping Bo1,
- Huan Chen1,
- Jian Zheng1,
- Qian-Hao Zhu2,
- Daguang Cai3,
- Chris Helliwell2 and
- Longjiang Fan1Email author
© Wang et al; licensee BioMed Central Ltd. 2010
Received: 8 July 2009
Accepted: 30 April 2010
Published: 30 April 2010
Endogenous non-coding small RNAs (21-24 nt) play an important role in post-transcriptional gene regulation in plants. Domestication selection is the most important evolutionary force in shaping crop genomes. The extent of polymorphism at small RNA loci in domesticated rice and whether small RNA loci are targets of domestication selection have not yet been determined.
A polymorphism survey of 94 small RNA loci (88 MIRNAs, four TAS3 loci and two miRNA-like long hairpins) was conducted in domesticated rice, generating 2 Mb of sequence data. Many mutations (substitution or insertion/deletion) were observed at small RNA loci in domesticated rice, e.g. 12 mutation sites were observed in the mature miRNA sequences of 11 MIRNAs (12.5% of the investigated MIRNAs). Several small RNA loci showed significant signals for positive selection and/or potential domestication selection.
Sequence variation at miRNAs and other small RNAs is higher than expected in domesticated rice. Like protein-coding genes, non-coding small RNA loci could be targets of domestication selection and play an important role in rice domestication and improvement.
Endogenous non-coding small RNAs (21-24 nt) play an important role in post-transcriptional gene regulation in plants. In general, small RNAs are grouped into two major classes, microRNAs (miRNAs) and short-interfering RNAs (siRNAs) based on the mechanisms by which they are synthesized and function . miRNAs are generated from a stem-loop primary transcript by the endonuclease Dicer-like 1 (DCL1) and its partners. To date, 414 miRNAs from 151 miRNA families (miRBase, http://microrna.sanger.ac.uk/sequences/, Release 14.0), including several natural antisense miRNAs (nat-miRNAs) and a mitron, have been identified in rice [2, 3]. siRNAs derive from double-stranded RNA precursors that are processed mainly by the endonucleases DCL2, DCL3 and/or DCL4 . Several types of siRNAs, such as trans-acting siRNAs (ta-siRNAs), natural antisense transcript-derived siRNAs (nat-siRNAs), repeat-associated siRNAs (ra-siRNAs) and miRNA-like long hairpin siRNAs, have been identified in rice and other plants [3, 5–13].
mRNAs are regulated by several classes of small RNAs through base complementarity leading to cleavage or repression of translation of the mRNA. It has been observed that single nucleotide polymorphisms (SNP) within a miRNA binding site can cause significant phenotypic changes . Genome-wide investigation of nucleotide variation at small RNA loci can determine the prevalence of functional variations in a species; this was first assessed in humans . In plants, a large-scale survey of miRNA polymorphism has been performed in Arabidopsis . This survey showed that mutations in mature miRNA sequences are rare and that the level of polymorphism in miRNAs is lower than in their flanking regions, suggesting purifying selection in the miRNA sequences. Evidence of positive selection on some miRNA loci has been found in Arabidopsis [16, 17]. Our previous studies on rice miRNAs and miRNA binding sites suggested that positive selection and nucleotide mutations play an important role in co-evolution of miRNAs and their targets [18, 19].
Despite the importance of small RNAs in regulation of development in crops, investigation of polymorphism of miRNAs and other small RNAs in rice and other crops is still in its infancy. The extent of sequence variation at functional regions of small RNAs and their potential contribution to phenotypic evolution during crop domestication are unclear. Unlike Arabidopsis, rice has experienced domestication and subsequent artificial genetic improvement and therefore has been shaped by an additional bottleneck effect. Rice is believed to have been domesticated approximately 10,000 years ago [20, 21]. Several protein-coding genes controlling domestication traits have been identified and a notable feature of these domestication genes is that they generally encode transcription factors that regulate other protein-coding genes by direct binding to their DNA . Most conserved miRNAs target transcription factors with an important role in plant development , but whether miRNAs and other types of non-coding small RNAs are targets of domestication selection in crops is yet to be determined.
Nucleotide divergence and neutrality of small RNAs
A total of 94 small RNA loci, including 88 MIRNAs from 40 miRNA families, four members of the TAS3 family and two miRNA-like long hairpins, were selected for a polymorphism survey of rice small RNAs (Additional file 1). With the exception of the two miRNA-like long hairpins, genomic fragments containing the intact precursor sequences of the small RNAs were amplified and sequenced from a population of 33 domesticated rice accessions (17 indica and 16 japonica cultivars), which were collected from diverse geographic locations (Additional file 2). For the two miRNA-like long hairpins which have large loop regions, the two stem regions were amplified and sequenced. Fragments from each rice accession were aligned to each other for each small RNA locus. The average fragment length of all alignments was 579.0 bp, which covered the precursor sequences of the investigated small RNAs (Additional file 3).
In addition to the mature miRNA sequences, many SNPs were also observed in the flanking regions of miRNAs. Compared to the mature miRNA sequences, significantly higher SNP densities or diversities were observed in the precursor miRNA (pre-miRNA) sequences, but the highest polymorphism was found in the upstream or downstream region of pre-miRNAs (Figure 2A and 2B). These results suggest that purifying selection was the predominant evolutionary force acting on miRNA sequences in rice as reported for Arabidopsis . Some miRNAs are conserved in both dicot and monocot plants while others are species-specific. Of the 88 MIRNA genes investigated in this study, 64 are conserved and 17 are rice-specific (Additional file 1). A dramatic difference in nucleotide diversities in the mature miRNAs and their flanking regions was observed between the conserved and rice-specific MIRNA loci (Figure 2C and 2D). The nucleotide diversity of the rice-specific MIRNAs was nearly twice that of the conserved MIRNAs. However, no significant difference in nucleotide diversity could be found in the miRNAs from the two subspecies of cultivated rice although the average number of SNPs per site of the indica subgroup was slightly lower than that of the japonica subgroup in both the mature miRNA and its flanking regions (Figure 2E and 2F).
Polymorphisms were also found in the functional regions of siRNA-generating loci in rice, for example, at phases P5'_5 and P5'_9 of a miRNA-like long hairpin (AK120922), and in ta-siARF and the 3' miR390 binding site of TAS3a2 (Additional file 5).
To examine neutrality of small RNAs in rice, Tajima's D , a widely used neutrality test, was employed. Tajima's D test measures the frequency distribution of polymorphisms and selection is expected to skew the population frequency of genetic variants relative to the neutral equilibrium model (NE). Under NE, the mean Tajima's D is expected to be zero and a negative value indicates an excess of rare sequence variants relative to NE expectation and a recent positive selection is thus inferred. In contrast to positive selection, balanced selection retains genetic differences and elevates the Tajima's D statistic towards a positive value .
Potential domestication selection
Indica and japonica rice were domesticated independently and their domestication genes usually show different evolutionary histories of selection between cultivated and wild rice populations [26, 27]. To find candidate small RNA loci that experienced domestication selection, 20 loci were selected from the set of 94 small RNA loci for further analysis (Figure 1) using a cultivated rice population including 29 indica and 25 japonica cultivars, a wild population including 15 accessions of O. rufipogon and an outgroup (Africa wild rice, O. barthii) (Additional file 2). Of the 20 loci selected, three (miR395a/b, TAS3a2 and MIR399d) had significant negative Tajima's D values and one (miR390) had extremely low divergence according to the above neutral test results (Additional file 3); the remaining loci were randomly chosen from the small RNA list shown in Additional file 1.
Sequence diversity, Tajima's D and HKA test results
To further investigate positive selection in these small RNA loci, in addition to the Tajima's D test, another independent test, the Hudson-Kreitman-Aguade (HKA) test , was performed. The HKA test is able to measure the degree of nucleotide diversity reduced by positive selection, a feature of the sequence data that cannot be measured by the Tajima's D test. Rejection of the multi-locus HKA test that measures the difference in diversity within species relative to divergence between species is often considered as evidence of a positive selection. The multi-locus HKA test requires an outgroup and reference loci that are believed to have not been affected by selection. For this purpose, African wild rice (O. barthii) that is the nearest phylogenetic node to O. sativa and O. rufipogon  was used as the outgroup and the aforementioned four neutral genes were used as the reference loci. Of the 20 small RNA loci, four miRNA genes (MIR164e, MIR390, MIR395a/b and MIR399d) and two siRNA genes (TAS3a2 and AK120922) had signatures of positive selection in indica and/or japonica subgroup according to Tajima's D or HKA test (Table 1). No positive selection was detected by the two tests for other 14 MIRNA loci. In the wild rice population, MIR399d, MIR1862d and the miRNA-like long hairpin (AK120922) presented evidence of positive selection (Tajima's D value of -2.04, -1.88 and -2.05, respectively), suggesting the non-neutral evolution of these three genes in the wild rice population. A gene can be considered as a candidate domestication gene if a significantly positive selection signal is detected in the cultivated population but not in the wild population . Based on this criterion, our results thus suggest that MIR164e, MIR390, MIR395a/b and TAS3a2 are potential candidates of small RNA loci that have experienced direct selection during rice domestication.
To date, large-scale studies of nucleotide variation at miRNAs have only been carried out in human  and Arabidopsis . Overall, our results in rice support the findings in human and Arabidopsis that purifying selection is one of the main evolutionary forces acting on rice MIRNA genes, maintaining lower levels of sequence divergence in mature miRNAs than in their flanking regions. In Arabidopsis, only a few substitution events were observed in the mature miRNA sequences (2 of the 66 investigated miRNAs) in a population including 24 diverse wild accessions. However, nucleotide mutations in the mature miRNA sequences were found in 11 of the 88 investigated rice miRNAs in this study. We also found polymorphisms in at least seven miRNAs in the wild rice population (data not shown). The higher percentage of SNPs or indels found in rice miRNAs implies that the polymorphism levels of miRNAs and other miRNA-like small RNAs may vary among species.
Our results suggest selection, and probably domestication selection, on some small RNA loci (Table 1 and Additional file 3). Small RNAs regulate genes that control a wide range of traits; the functions of these small RNA regulated genes might provide clues for the traits under selection via small RNA. One interesting result in this study is that both MIR390 and TAS3 are potential candidates of selection. miR390 acts as a guide for processing of TAS3 and biogenesis of ta-siRNAs, including the functional ta-siARF [5, 10, 11]. ARFs (auxin response factors) are encoded by a large gene family (25 members in rice) and are involved in regulation of a wide range of biological functions in rice . ta-siARF downregulates ARF2, ARF3 and ARF4 [5, 12] that affect gynoecium patterning , proper timing of vegetative shoot development and establishment of leaf polarity . At least two ARFs have been detected as selected genes in maize [27, 35]. It is unclear whether rice ARF genes have experienced selection during evolution. Our results suggest that genes involved in the same regulatory pathway could all be targeted by domestication selection for enhanced growth response and productivity. miR164e targets NAC domain transcription factors that play an important role in development, stress-tolerance and disease-resistance in rice [36, 37]. miR395 targets ATP-sulfurylases that are involved in sulfate assimilation. Expression of miR395 is greatly increased under sulfate starvation conditions . Sulfur is an essential macronutrient required by plants. miR395 regulated sulfate metabolism might have been selected for high sulfate usage efficiency. Our results suggest that miRNAs and other small RNA loci should not be ignored in the endeavor to identify the molecular basis of domestication and improvement of crops.
Note that there are three limitations in our effort to identify candidate small RNA loci experienced domestication selection. First, a bottleneck effect or demographic effect can also cause significant reduction of genetic diversity and it could not be excluded in our analysis. Second, the positive selection detected in this study could have arisen by selective sweep rather than direct selection on the small RNA loci themselves. Third, although we identified several mutations in the functional regions (mature miRNA, miRNA precursor or promoter) of the candidate small RNA genes in domesticated rice, such as in the mature miR395a/b sequences and in tasi-ARF and the miR390 binding sites of TAS3 (Additional file 4 and 5), and found distinct expression levels of miR164e in the cultivated and wild rice (data not shown), it does not necessarily indicate a direct consequence of the domestication event. Further experiments are needed to reveal the consequence of positive selection, such as phenotypic changes.
Our sequencing results revealed sequence polymorphism in miR166e and its target gene Os03 g16320 (Figure 3). By checking the recently released rice SNP dataset (http://www.oryzasnp.org), we found a similar situation with miR443 and the target gene Os01 g49940 (a C → T mutation at the fifth position from the 3' end of miR443 and a G → A mutation at the eighth position from the 3' end of the miR443 binding site in Os01 g49940). We also observed SNPs in the miRNA binding sites of targets of miR164d/e and miR530, but no SNP was observed in the mature sequences of the two miRNAs. These results suggest that evolution of miRNAs and their targets could be independent, although a co-evolution relationship between miRNAs and their target sequences does happen in cultivated rice as demonstrated previously .
At least one member each of four miRNA families (miR166, miR167, miR171 and miR395) were found to have experienced positive selection based on Tajima's D test in the natural Arabidopsis populations [16, 17]. These four miRNA families are conserved in dicot and monocot plants. In our study, significant signals of positive selection were also detected by Tajima's D test in at least one member of the miR166, miR167 and miR395 family in indica or japonica subgroup (Additional file 3), implying a potential conservation of adaptive evolution between rice and Arabidopsis for these miRNA families. A common feature of these three families is that they have relatively large number of family members (e.g. 23 and 6 miR395 members in rice and Arabidopsis, respectively). These homologous genes might have experienced rapid functional divergence after gene expansion. It has been suggested that positive selection acts on gene copy number variations . Positive selection could promote adaptive evolution but whether plants enhance their fitness in the changing environment by adding copy number of miRNAs requires further investigation.
Our survey of polymorphisms in 94 small RNA loci in rice, the first such effort in crops, revealed a higher level of sequence variations in miRNAs in domesticated rice than in Arabidopsis. In agreement with previous studies, our investigation suggests that purifying selection dominated the evolution of small RNAs in rice. A neutrality test revealed non-neutral evolution of small RNA loci in which both positive selection and balanced selection are involved. Our results suggest that domestication selection on small RNA loci could play an important role in rice domestication and/or improvement.
Fifty-four cultivated rice accessions (Oryza sativa, 29 indica and 25 japonica) from a wide range of geographical locations and 15 accessions of the wild ancestor, O. rufipogon, were selected for investigation of nucleotide diversity in small RNA loci. Two accessions of the Africa wild rice relative, O. barthii, were used as an outgroup species in a neutrality test. Details of the 71 rice accessions are in Additional file 2.
Small RNAs selected for analyses of nucleotide diversity and domestication selection
The majority of the MIRNAs selected are conserved in other species to ensure they are bona fide MIRNAs. For all the rice-specific MIRNAs, authenticity is supported by experimental evidence and BLAST analysis has been used to remove those potentially derived from repeat regions of the rice genome. The siRNAs used in this study have been experimentally identified and confirmed by our previous study 
PCR and DNA Sequencing
Primers were designed based on the genomic sequence of japonica cultivar Nipponbare of the Rice Genome Annotation Project (http://rice.plantbiology.msu.edu/; [40, 41] using the Primer3 program http://fokker.wi.mit.edu/primer3/input.htm and were compared to the rice genome sequence to ensure their specificity. Primer pairs were designed to amplify products with a length of 600 - 1300 bp, including the pre-miRNA and its flanking regions (for MIRNA loci). Details of the primers and their products were provided in Additional file 6. Polymerase chain reaction (PCR) amplification was carried out in a DNAEngine Peltier Thermal Cycler (BIO-RAD) in a total volume of 50 μl. The reaction mixture contained 10-100 ng of template DNA, 2 μl of each primer (10-20 pmol), 0.6 μl of Taq DNA Polymerase (5 U/μl), 1 μl of dNTPs (10 mM), 3 μl of MgCl2 (25 mM), 5 μl of 10× PCR buffer (Sangon, China). PCR was performed using the following conditions: 5 min at 95°C followed by 35 cycles of 30 s at 95°C, 30 s at 54°C, and 90 s at 72°C, and a final extension at 72°C for 10 minutes. When this program failed to amplify the expected products, specific conditions were experimentally determined. For O. sativa cultivars, purified PCR products were directly sequenced using the forward or reverse primer. For wild rice, in which either homozygous or heterozygous individuals might exist, PCR fragments were cloned into pGEM T-Easy vector (PMD19-T, Takara) and sequenced using the forward or reverse primer. The sequencing quality files were checked manually to avoid the sequencing errors. Sequences with poor quality were re-sequenced. All sequences have been deposited into GenBank under the accession numbers GQ418390-GQ420345 and HM138917-HM140183. In total, approximately 2.0 Mb of sequences were generated in this study.
CLUSTALW version 1.82-UNIX was used for multiple sequence alignments  and segregating sites were identified manually using GeneDoc . Two summary statistics parameters (θ and π) for nucleotide diversity and Tajima's D statistic  were calculated with a perl script. θ and π estimate the population mutation rate per locus based on the number of segregating sites and the mean value of pairwise divergence per locus, respectively [45, 46]. Multi-locus HKA tests were performed using the HKA program available at http://lifesci.rutgers.edu/hey/home. The entire segments that contain the small RNAs were amplified by the primers shown in Additional file 6 and used for a neutrality test. For MIRNA loci, nucleotide diversities were separately estimated for the mature miRNAs, the pre-miRNAs, and the upstream and downstream regions of the mature miRNA (the miR395 family that is located at four clusters in which most members are less than 100 bp apart from each other was not included in this analysis). Sequences of the mature miRNAs and the pre-miRNAs were retrieved from miRBase http://microrna.sanger.ac.uk/sequences/. Conservation of miRNAs in dicots and monocots is based on Sunkar and Jagadeeswaran . Targets of miRNA were predicted using Patscan  by searching the full length cDNA sequence dataset (release 5) of the Rice Genome Annotation http://rice.plantbiology.msu.edu with a penalty score up to four.
short inferring RNA
- TAS3 :
trans-acting siRNA locus 3
the Hudson-Kreitman-Aguade test
insertion and deletion
single nucleotide polymorphism.
This work was supported by the National Basic Research Program of China (2006CB101700/2008CB117002), National High Technology Research and Development Program of China (2006AA10A102) and the Program of Scientific Methodology by Minister of Science and Technology (2008IM020800).
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.View ArticlePubMedGoogle Scholar
- Lu C, Jeong DH, Kulkarni K, Pillay M, Nobuta K, German R, Thatcher SR, Maher C, Zhang L, Ware D, et al: Genome-wide analysis for discovery of rice microRNAs reveals natural antisense microRNAs (nat-miRNAs). Proc Natl Acad Sci USA. 2008, 105 (12): 4951-4956. 10.1073/pnas.0708743105.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhu QH, Spriggs A, Matthew L, Fan L, Kennedy G, Gubler F, Helliwell C: A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res. 2008, 18 (9): 1456-1465. 10.1101/gr.075572.107.PubMed CentralView ArticlePubMedGoogle Scholar
- Vazquez F: Arabidopsis endogenous small RNAs: highways and byways. Trends Plant Sci. 2006, 11 (9): 460-468. 10.1016/j.tplants.2006.07.006.View ArticlePubMedGoogle Scholar
- Allen E, Xie Z, Gustafson AM, Carrington JC: microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell. 2005, 121 (2): 207-221. 10.1016/j.cell.2005.04.004.View ArticlePubMedGoogle Scholar
- Axtell MJ, Jan C, Rajagopalan R, Bartel DP: A two-hit trigger for siRNA biogenesis in plants. Cell. 2006, 127 (3): 565-577. 10.1016/j.cell.2006.09.032.View ArticlePubMedGoogle Scholar
- Borsani O, Zhu J, Verslues PE, Sunkar R, Zhu JK: Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell. 2005, 123 (7): 1279-1291. 10.1016/j.cell.2005.11.035.PubMed CentralView ArticlePubMedGoogle Scholar
- Heisel SE, Zhang Y, Allen E, Guo L, Reynolds TL, Yang X, Kovalic D, Roberts JK: Characterization of unique small RNA populations from rice grain. PLoS One. 2008, 3 (8): e2871-10.1371/journal.pone.0002871.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu B, Chen Z, Song X, Liu C, Cui X, Zhao X, Fang J, Xu W, Zhang H, Wang X, et al: Oryza sativa dicer-like4 reveals a key role for small interfering RNA silencing in plant development. Plant Cell. 2007, 19 (9): 2705-2718. 10.1105/tpc.107.052209.PubMed CentralView ArticlePubMedGoogle Scholar
- Peragine A, Yoshikawa M, Wu G, Albrecht HL, Poethig RS: SGS3 and SGS2/SDE1/RDR6 are required for juvenile development and the production of trans-acting siRNAs in Arabidopsis. Genes Dev. 2004, 18 (9): 2368-2379. 10.1101/gad.1231804.PubMed CentralView ArticlePubMedGoogle Scholar
- Vazquez F, Vaucheret H, Rajagopalan R, Lepers C, Gasciolli V, Mallory AC, Hilbert JL, Bartel DP, Crete P: Endogenous trans-acting siRNAs regulate the accumulation of Arabidopsis mRNAs. Mol Cell. 2004, 16 (1): 69-79. 10.1016/j.molcel.2004.09.028.View ArticlePubMedGoogle Scholar
- Williams L, Carles CC, Osmont KS, Fletcher JC: A database analysis method identifies an endogenous trans-acting short-interfering RNA that targets the Arabidopsis ARF2, ARF3, and ARF4 genes. Proc Natl Acad Sci USA. 2005, 102 (27): 9703-9708. 10.1073/pnas.0504029102.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou X, Sunkar R, Jin H, Zhu JK, Zhang W: Genome-wide identification and analysis of small RNAs originated from natural antisense transcripts in Oryza sativa. Genome Res. 2009, 19 (1): 70-78. 10.1101/gr.084806.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Clop A, Marcq F, Takeda H, Pirottin D, Tordoir X, Bibe B, Bouix J, Caiment F, Elsen JM, Eychenne F, et al: A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nat Genet. 2006, 38 (7): 813-818. 10.1038/ng1810.View ArticlePubMedGoogle Scholar
- Saunders MA, Liang H, Li WH: Human polymorphism at microRNAs and microRNA target sites. Proc Natl Acad Sci USA. 2007, 104 (9): 3300-3305. 10.1073/pnas.0611347104.PubMed CentralView ArticlePubMedGoogle Scholar
- Ehrenreich IM, Purugganan MD: Sequence variation of MicroRNAs and their binding sites in Arabidopsis. Plant Physiol. 2008, 146 (4): 1974-1982. 10.1104/pp.108.116582.PubMed CentralView ArticlePubMedGoogle Scholar
- de Meaux J, Hu JY, Tartler U, Goebel U: Structurally different alleles of the ath-MIR824 microRNA precursor are maintained at high frequency in Arabidopsis thaliana. Proc Natl Acad Sci USA. 2008, 105 (26): 8994-8999. 10.1073/pnas.0803218105.PubMed CentralView ArticlePubMedGoogle Scholar
- Guo X, Gui Y, Wang Y, Zhu QH, Helliwell C, Fan L: Selection and mutation on microRNA target sequences during rice evolution. BMC Genomics. 2008, 9: 454-10.1186/1471-2164-9-454.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang S, Zhu QH, Guo X, Gui Y, Bao J, Helliwell C, Fan L: Molecular evolution and selection of a gene encoding two tandem microRNAs in rice. FEBS Lett. 2007, 581 (24): 4789-4793. 10.1016/j.febslet.2007.09.002.View ArticlePubMedGoogle Scholar
- Kovach MJ, Sweeney MT, McCouch SR: New insights into the history of rice domestication. Trends Genet. 2007, 23 (1): 578-587. 10.1016/j.tig.2007.08.012.View ArticlePubMedGoogle Scholar
- Vaughan DA, Lu B, Tomooka N: The evolving story of rice evolution. Plant science. 2008, 174: 394-408.View ArticleGoogle Scholar
- Doebley JF, Gaut BS, Smith BD: The molecular genetics of crop domestication. Cell. 2006, 127 (7): 1309-1321. 10.1016/j.cell.2006.12.006.View ArticlePubMedGoogle Scholar
- McNally KL, Childs KL, Bohnert R, Davidson RM, Zhao K, Ulat VJ, Zeller G, Clark RM, Hoen DR, Bureau TE, et al: Genomewide SNP variation reveals relationships among landraces and modern varieties of rice. P Natl Acad Sci USA. 2009, 106 (30): 12273-12278. 10.1073/pnas.0900992106.View ArticleGoogle Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.PubMed CentralPubMedGoogle Scholar
- Wright SI, Gaut BS: Molecular population genetics and the search for adaptive evolution in plants. Mol Biol Evol. 2005, 22 (3): 506-519. 10.1093/molbev/msi035.View ArticlePubMedGoogle Scholar
- Londo JP, Chiang YC, Hung KH, Chiang TY, Schaal BA: Phylogeography of Asian wild rice, Oryza rufipogon, reveals multiple independent domestications of cultivated rice, Oryza sativa. Proc Natl Acad Sci USA. 2006, 103 (25): 9578-9583. 10.1073/pnas.0603152103.PubMed CentralView ArticlePubMedGoogle Scholar
- Yamasaki M, Tenaillon MI, Bi IV, Schroeder SG, Sanchez-Villeda H, Doebley JF, Gaut BS, McMullen MD: A large-scale screen for artificial selection in maize identifies candidate agronomic loci for domestication and crop improvement. Plant Cell. 2005, 17 (1): 2859-2872. 10.1105/tpc.105.037242.PubMed CentralView ArticlePubMedGoogle Scholar
- Olsen KM, Purugganan MD: Molecular evidence on the origin and evolution of glutinous rice. Genetics. 2002, 162 (2): 941-950.PubMed CentralPubMedGoogle Scholar
- Zhu Q, Zheng X, Luo J, Gaut BS, Ge S: Multilocus analysis of nucleotide variation of Oryza sativa and its wild relatives: severe bottleneck during domestication of rice. Mol Biol Evol. 2007, 24 (3): 875-888. 10.1093/molbev/msm005.View ArticlePubMedGoogle Scholar
- Hudson RR, Kreitman M, Aguade M: A test of neutral molecular evolution based on nucleotide data. Genetics. 1987, 116 (1): 153-159.PubMed CentralPubMedGoogle Scholar
- Zou XH, Zhang FM, Zhang JG, Zang LL, Tang L, Wang J, Sang T, Ge S: Analysis of 142 genes resolves the rapid diversification of the rice genus. Genome Biol. 2008, 9 (3): R49-10.1186/gb-2008-9-3-r49.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang D, Pei K, Fu Y, Sun Z, Li S, Liu H, Tang K, Han B, Tao Y: Genome-wide analysis of the auxin response factors (ARF) gene family in rice (Oryza sativa). Gene. 2007, 394 (1-2): 13-24. 10.1016/j.gene.2007.01.006.View ArticlePubMedGoogle Scholar
- Nemhauser JL, Feldman LJ, Zambryski PC: Auxin and ETTIN in Arabidopsis gynoecium morphogenesis. Development. 2000, 127 (8): 3877-3888.PubMedGoogle Scholar
- Hunter C, Willmann MR, Wu G, Yoshikawa M, de la Luz Gutierrez-Nava M, Poethig SR: Trans-acting siRNA-mediated repression of ETTIN and ARF4 regulates heteroblasty in Arabidopsis. Development. 2006, 133 (5): 2973-2981. 10.1242/dev.02491.PubMed CentralView ArticlePubMedGoogle Scholar
- Wright SI, Bi IV, Schroeder SG, Yamasaki M, Doebley JF, McMullen MD, Gaut BS: The effects of artificial selection on the maize genome. Science. 2005, 308 (5726): 1310-1314. 10.1126/science.1107891.View ArticlePubMedGoogle Scholar
- Lin RM, Zhao WS, Meng XB, Wang M, Peng YL: Rice gene OsNAC19 encodes a novel NAC-domain transcription factor and responds to infection by Magnaporthe grisea. Plant science. 2007, 172 (1): 120-130. 10.1016/j.plantsci.2006.07.019.View ArticleGoogle Scholar
- Zheng X, Chen B, Lu G, Han B: Overexpression of a NAC transcription factor enhances rice drought and salt tolerance. Biochem Biophys Res Commun. 2009, 379 (4): 985-989. 10.1016/j.bbrc.2008.12.163.View ArticlePubMedGoogle Scholar
- Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14 (6): 787-799. 10.1016/j.molcel.2004.05.027.View ArticlePubMedGoogle Scholar
- Kelley JL, Swanson WJ: Positive selection in the human genome: from genome scans to biological significance. Annu Rev Genomics Hum Genet. 2008, 9: 143-160. 10.1146/annurev.genom.9.081307.164411.View ArticlePubMedGoogle Scholar
- Ouyang S, Zhu W, Hamilton J, Lin H, Campbell M, Childs K, Thibaud-Nissen F, Malek RL, Lee Y, Zheng L, et al: The TIGR Rice Genome Annotation Resource: improvements and new features. Nucleic Acids Res. 2007, D883-887. 10.1093/nar/gkl976. 35 Database
- Yuan Q, Ouyang S, Wang A, Zhu W, Maiti R, Lin H, Hamilton J, Haas B, Sultana R, Cheung F, et al: The institute for genomic research Osa1 rice genome annotation database. Plant Physiol. 2005, 138 (1): 18-26. 10.1104/pp.104.059063.PubMed CentralView ArticlePubMedGoogle Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.PubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Nicholas KB, Nicholas JHB: GeneDoc: a tool for editing and annotating multiple sequence alignments. Distributed by the authors. 1997, [http://www.psc.edu/biomed/genedoc]Google Scholar
- Tajima F: Evolutionary relationship of DNA sequences in finite populations. Genetics. 1983, 105 (2): 437-460.PubMed CentralPubMedGoogle Scholar
- Watterson GA: On the number of segregating sites in genetical models without recombination. Theor Popul Biol. 1975, 7 (2): 256-276. 10.1016/0040-5809(75)90020-9.View ArticlePubMedGoogle Scholar
- Sunkar R, Jagadeeswaran G: In silico identification of conserved microRNAs in large number of diverse plant species. BMC Plant Biol. 2008, 8: 37-10.1186/1471-2229-8-37.PubMed CentralView ArticlePubMedGoogle Scholar
- Dsouza M, Larsen N, Overbeek R: Searching for patterns in genomic data. Trends Genet. 1997, 13 (2): 497-498. 10.1016/S0168-9525(97)01347-4.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.