|
|
||||||||
Systematics and Phytogeography |
2Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; 3Graduate School, Chinese Academy of Sciences, Beijing 100039, China
Received for publication December 27, 2004. Accepted for publication May 16, 2005.
ABSTRACT
The phylogeny and evolutionary history of the rice tribe (Oryzeae) were explored using sequences of five DNA fragments (matK, trnL, nad1, Adh2, and GPA1) from chloroplast, mitochondrial, and nuclear genomes. Results indicate that (1) Oryzeae is monophyletic and falls into two main clades corresponding to the traditionally recognized subtribes; (2) previous recognition of three monotypic genera (Hydrochloa, Porteresia, and Prosphytochloa) is not justified; and (3) close affinities of the monoecious genera are not supported, suggesting the possibility of multiple origins of unisexual florets. Based on the magnitude of matK and GPA1 sequence divergence, we suggest that Oryza and Leersia branched off from the remaining genera of Oryzeae
20 million years ago (mya), and separated from each other
14 mya. A divergence time of
9 mya is obtained for the most basal split within Oryza. These estimates suggest that Oryzeae diverged during the Miocene, and thus imply that long-distance dispersal appears to be one of the important factors in the diversification of the tribe.
Key Words: biogeography divergence time multiple sequences Oryzeae phylogeny
The tribe Oryzeae (Poaceae), as conventionally delimited, includes approximately 12 genera and more than 70 species distributed throughout the tropical and temperate regions of the world (Clayton and Renvoize, 1986
; Vaughan, 1994
). As the largest tribe in the subfamily Ehrhartoideae, Oryzeae contains more than half of both genera and species of the subfamily (Watson and Dallwitz, 1999
; GPWG, 2001
). In this tribe, Oryza L. with approximately 23 species and Leersia Sw. with approximately 18 species are the two largest genera and distributed worldwide (Watson and Dallwitz, 1999
; Terrell et al., 2001
). Except for Zizania L., which is disjunctly distributed in eastern Asia and eastern North America, the remaining nine genera are distributed regionally or confined to a specific continent (Watson and Dallwitz, 1999
). Chikusichloa Koidz and two monotypic genera (Hygroryza Nees and Porteresia Tateoka) are distributed in Asia. Zizaniopsis Doell ex Asch. and Luziola Juss., as well as the monotypic Rhynchoryza Baill., occur in North and South America, while monotypic Potamophila R. Br. is endemic to Australia and Prosphytochloa Schweick. to Africa.
In the genus Oryza, the Asian cultivated rice (O. sativa) is one of the world's most important crops and a primary food source for more than one-half of the world's population (Chandler and Wessler, 2001
; Ge et al., 2001
). Other members of Oryzeae have important economic value, such as the wild species of Oryza that contain agronomically useful traits for rice genetic improvement, Zizania species that are a well-known part of cuisine in China (Z. latifolia) and North America (Z. palustris L.), and Leersia species that are useful forage (Vaughan, 1994
; Kennard et al., 1999
). Importantly, rice is becoming an excellent research model due to its small genome size, efficient genetic transformation, and extensive colinearity with other grass species (Gale and Devos, 1998
; Schmidt, 2000; Shimamoto and Kyozuka, 2002
). In particular, the completion of the rice genome sequence provides tremendous value of rice as a model plant (Bennetzen, 2002
; Rensink and Buell, 2004
). Therefore, a clear picture of phylogenetic relationships of rice and its relatives will provide an important foundation for future breeding programs and for biological studies of rice and grasses in general (Kellogg, 1998
; Ge et al., 2001
; Rensink and Buell, 2004
).
Because of the economical and theoretical importance of rice, classification and phylogeny of Oryza and its closely related genera have been studied through various approaches including morphology (Weatherwax, 1929
; Launert, 1965
; Terrell and Robinson, 1974
; Kellogg and Watson, 1993
; Terrell et al., 2001
), cytology (Ramanujam, 1938
; Li et al., 1963
; Nayar, 1973
), restriction site analysis (Zhang and Second, 1989
; Duvall et al., 1993
), and DNA sequencing (Ge et al., 2002
; Hass et al., 2003
). Although evidence showed that this tribe was a distinct and monophyletic group (Duistermaat, 1987
; Zhang and Second, 1989
; Kellogg and Watson, 1993
; Ge et al., 2002
), its subdivision and phylogenetic relationships among genera have been inconsistent among studies (Hitchcock, 1920
; Hubbard, 1934
, 1959
; Hitchcock and Chase, 1951
; Stebbins and Crampton, 1961
; Pyrah, 1969
; Terrell and Robinson, 1974
; Duvall et al., 1993
). Previous morphological investigations treated the oryzoid grasses as a taxonomic bifurcation of monoecious vs. bisexual groups and recognized the two groups variously at tribal (Oryzeae vs. Zizanieae) (Hitchcock, 1920
; Hitchcock and Chase, 1951
; Stebbins and Crampton, 1961
) or subtribal (Oryzinae vs. Zizaniinae) (Hubbard, 1934
, 1959
; Pyrah, 1969
) levels. Through morphological and anatomical investigations, Terrell and Robinson (1974)
established a new subtribe Luziolinae and proposed a three-subtribe system, but this treatment was not supported by later studies of restriction sites and sequence data (Zhang and Second, 1989
; Duvall et al., 1993
; Ge et al., 2002
). In addition, the delimitation and taxonomic position of some genera in this tribe have long been problematic and are still controversial mainly due to limited samples or insufficient character information used in previous studies. For example, the mutually exclusive hypotheses that Zizania is more closely related to the monocious genera or to the bisexual genera have been in dispute for decades (Terrell and Robinson, 1974
; Duvall et al., 1993
; Ge et al., 2002
). Whether or not to place the monotypic genus Hydrochloa Beauv. in synonymy under Luziola and to treat the monotypic genus Porteresia as a member of Oryza are also not resolved with certainty (Terrell and Robinson, 1974
; Duvall et al., 1993
).
Based on sequences of the chloroplast matK gene, Ge et al. (2002)
studied the phylogeny of Oryzeae. They verified that the tribe was a monophyletic group consisting of two strongly supported monophyletic groups corresponding to the two traditionally recognized subtribes, and they also suggested treating Porteresia as a member of Oryza. Although the matK work was so far the most comprehensive phylogenetic study of this tribe using sequence data in terms of species and generic inclusion, many of the questions just mentioned remain unsolved or the answers uncertain because of the relatively low resolution of the chloroplast gene. Growing examples in the recent literature have demonstrated that low-copy nuclear genes have great potential to compensate cpDNA and nrDNA for improvement of resolution and robustness of plant phylogenetic reconstruction (Soltis and Soltis, 1998
; Doyle and Doyle, 1999
; Sang, 2002
; Grob et al., 2004
). The alcohol dehydrogenase (Adh) gene is the most widely used low-copy nuclear gene, whereas the chloroplast matK gene and trnL intron as well as the mitochondrial nad1 intron 2 are also used in phylogenetic studies (Soltis and Soltis, 1998
; Ge et al., 1999
, 2002
; Gugerli et al., 2001
; Sang, 2002
; Hass et al., 2003
). Because nuclear GPA1 that encodes the G protein
subunit is present as a single copy in higher plants and is well characterized in function and structure (Seo et al., 1995
; Fujisawa et al., 1999
), we successfully used it to investigate reticulate evolution of tetraploids in Oryza (Bao and Ge, 2004
). In the present study, we reconstruct the phylogeny of Oryzeae using sequences of five DNA fragments from three genomes, including two chloroplast (matK and trnL intron), one mitochondrial (nad1 intron 2), and two nuclear (Adh2 and GPA1) fragments. Our specific objectives are: (1) to evaluate the circumscription of the groups at tribal and subtribal levels and the subtribal classification systems of the Oryzeae, (2) to investigate generic relationships within the tribe with particular emphasis on the systematic positions of some genera, and (3) to gain insights into the origin and biogeographic history of the tribe from a molecular phylogenetic perspective.
MATERIALS AND METHODS
Plant materials and DNA isolation
Twenty-four species were used in the present study, representing all 12 genera of the tribe Oryzeae (Clayton and Renvoize, 1986
; Vaughan, 1994
), except for Maltebrunia that was placed in synonymy under the genus Potamophila by Duistermaat (1987)
. Nine Oryza species representing all six diploid genomes and four Leersia species including two diploids and two tetraploids were sampled in this study. Two species from each of Luziola and Zizania and one species from each of the remaining seven genera, including five monotypic genera (Hygroryza, Porteresia, Potamophila, Prosphytochloa, Rhynchoryza), were sampled. Based on a recent comprehensive study of the subfamilial classification of the grass family (GPWG, 2001
), the subfamily Ehrhartoideae includes three tribes: Ehrharteae, Oryzeae, and Phyllorachideae, and the tribe Ehrharteae is the most closely related tribe to Oryzeae. Therefore, we chose one species of the genus Ehrharta Thunb. as an outgroup. One additional species of the closely related subfamily Bambusoideae was also included in the phylogenetic analysis as an outgroup. The details of the sampled species with their scientific names, accession numbers or vouchers, origins, chromosome numbers, and GenBank accession numbers are listed in the Appendix. Total DNA was isolated from silica-gel dried or fresh leaves using the cetyltrimethylammonium bromide (CTAB) method as described previous by Ge et al. (1999)
.
Amplification, cloning, and sequencing
To amplify the chloroplast matK and trnL intron, mitochondrial nad1 intron 2, and nuclear Adh2 and GPA1, the polymerase chain reactions (PCR) were carried out in a volume of 25 µL, containing 550 ng of genomic DNA, 5.0 pmol of each primer, 0.2 mM of each dNTP, 2 mM MgCl2 and 0.65 units of Taq DNA polymerase. Amplifications were performed in a Peltier thermal cycler (PTC-200, PerkinElmer Corporation, Wellesley, USA) .
The primers for amplifying and sequencing matK followed those in Ge et al. (1999)
. Most sequences of matK were used in Ge et al. (2002)
, but additional sequences from Leersia oryzoides, Luziola fluitans, and Ehrharta erecta were generated in this study. The chloroplast trnL intron was amplified and sequenced with primers trn-L1 and trn-L2. The two primers, located on the 5' exon and 3 exon of the trnL gene, respectively, were designed based on sequences of Oryza sativa (GenBank accession X15901), Triticum aestivum L. (AF148757), and Zea mays L. (X86563). The entire intron of mitochondrial nad1 intron 2 and a small part of the 3' exon was amplified and sequenced using the primers of Demesure et al. (1995)
. The PCR cycling parameters for the three loci were similar, with an initial 4 min at 70°C, followed by two cycles of 1 min at 94°C, 20 s at 52°C, and 1.5 min at 72°C; after that, 35 cycles of 20 s at 94°C, 20 s at 52°C, and 1.5 min at 72°C were conducted, with a final extension time of 10 min at 72°C.
The nuclear Adh2 gene was amplified and sequenced with primers described in Ge et al. (1999)
. Two PCR primers (AdhF1 and AdhR1) are located in exon 2 and exon 8, respectively, and thus the sequenced Adh2 region includes six introns and five exons (Ge et al., 1999
). Sequences of all Oryza species and Zizaniopsis villanensis were used in Ge et al. (1999)
. The nuclear GPA1 was amplified and sequenced with primers GPA1-FF and GPA1-14R (Bao and Ge, 2004
). These two primers were designed based on GPA1 sequences of Oryza sativa (L35844) and Zea mays (AF055471). The region amplified by the two primers spanned from exon 9 to exon 14, thus including four exons and five introns (Bao and Ge, 2004
). The temperature profile for thermocycling to amplify the two nuclear genes was an initial 4 min at 70°C, followed by two cycles of 1 min at 94°C, 20 s at 54°C, and 1.5 min at 72°C; After that, 35 cycles of 20 s at 94°C, 20 s at 54°C, and 1.5 min at 72°C, with a final extension time of 10 min at 72°C. All the amplifying and internal sequencing primers for these fragments are listed in Table 1.
|
Data analysis
Sequences were aligned using ClustalX version 1.81 (Thompson et al., 1997
) and refined manually. The possibility of sequence saturation was examined using DAMBE version 4.2.13 (Xia and Xie, 2001
). A plot of the number of transitions and transversions vs. divergence offers a visual display of substitution saturation, with an asymptotic relationship indicating the presence of saturation (Xia and Xie, 2001
). Phylogenetic analyses of the sequence data were performed using the parsimony and Bayesian Markov chain Monte Carlo (MCMC) methods. Maximum parsimony (MP) analyses were conducted using PAUP version 4.0b10 (Swofford, 2001
). All characters were equally weighted, gaps were treated as missing, and character states were treated as unordered. Heuristic search was performed with MULPARS option, tree-bisection-reconnection (TBR) branch swapping, and RANDOM stepwise addition with 1000 replicates. Topological robustness was assessed by bootstrap analysis with 1000 replicates using simple taxon addition (Felsenstein, 1985
). An appropriate nucleotide substitution model was determined using Modeltest version 3.5 (Posada and Crandall, 1998
) for each data set. The models were chosen according to the hierarchical likelihood ratio test (LRT) and then used for subsequent Bayesian analysis. Bayesian inference (BI) was conducted using MrBayes version 3.0b4 (Huelsenbeck and Ronquist, 2001
). One cold and three incrementally heated MCMC chains were run for 106 generations, with trees sampled every 10 generation, using the random tree as its starting point and a temperature parameter value of 0.2 (the default setting of MrBayes). For each data set, MCMC runs were repeated twice as a safeguard against spurious results. The first 104 trees were discarded as burn-in, and the remaining trees were used to construct Bayesian trees. Examination of the log-likelihoods and the observed consistency between runs suggested that the burn-in periods were sufficiently long enough for chains to have become stationary.
The partition homogeneity test, as implemented in PAUP version 4.0b10 (Swofford, 2001
), was used to evaluate the topological congruence between gene trees (Farris et al., 1995
). Replicates were analyzed using the parsimony criterion, and branch swapping using TBR was performed and one tree held at each step during the stepwise addition. The partition homogeneity test was performed on either cytoplasmic (chloroplast plus mitochondrial) or nuclear loci, repectively. The results were used to determine whether each combined data set had significant incongruence. Following Cunningham (1997)
, we took a P value of 0.01 as a significance criterion for this test.
For a rough estimate of the divergence time of Oryzeae and its major lineages, a simple molecular clock assumption was employed, and the clock calibration was carried out based on the assumption that maize and rice diverged 50 million years ago [mya] (Stebbins, 1981
; Gaut, 2002
). Divergence times (T) of the Oryzeae and its major lineages were calculated by dividing the the estimated numbers of substitutions per site between homologous sequences (K) by twice the rate of nucleotide substitution (r), namely, using the formula T = K/(2) (r). Relative rate tests were implemented by Hyphy (Kosakovsky Pond et al., 2005
) and RRTree (Robinson-Rechavi and Huchon, 2000
) and used to assess heterogeneity among different lineages. Rates of nucleotide substitution for individual genes were estimated by the Nei-Gojobori model with the Jukes-Cantor and Kimura two-parameter corrections as implemented in MEGA version 2.1 (Kumar et al., 2001
).
RESULTS
Amplification, sequencing, and sequence characteristics
Amplifications for three cytoplasmic loci from all the species included in this study were successful except for the matK sequence of E. erecta for which no amplified product was obtained by using primers matKF1 and matKR1. Therefore, we used another pair of primers (matK-AF and matK-8R) (Kato et al., 1998
) to amplify matK from E. erecta (Table 1). The resulting fragment was about 1.2 kilobases (kb), shorter than those of other species (1.5 kb). The sequence length of matK ranged from 1502 to 1562 bp (1143 bp for E. erecta) with an aligned length of 1574 bp. The characteristics of the three cytoplasmic loci are shown in Table 2. As expected, the mitochondrial nad1 intron 2 is less variable and informative than the chloroplast genes. The mean GC content of mitochondrial nad1 intron 2 (54.2%) is higher than those of chloroplast fragments (34.4% for matK and 34.9% for trnL). Overall, 3.6 kb of cytoplasmic DNA sequences are determined for all the species involved.
|
1.1-kb Adh2 fragment from P. parviflora. For multiple clones of diploid species, we obtained almost identical sequences but occasionally found several base pair differences among clones. However, initial phylogenetic analysis indicated that Adh2 sequences of different clones from the same diploid species formed a monophyletic group, and thus we only retained one clone in our final analyses because exclusion of these minor types of clones have no impact on relationships among species. For tetraploid species, two distinct types of sequences were identified at the Adh2 locus for the majority of tetraploid species except for two Leersia species (L. oryzoides and L. hexandra), in which only one type of sequence was sampled. The detection of only one Adh2 sequence for the tetraploid L. oryzoides and L. hexandra is either indication of autoploidy or due to screening an insufficient number of clones. An exhaustive search for the other copy was not carried out, because it does not seem to be critical for understanding the phylogenetic relationships of Oryzeae in this study. In contrast, we got two distinct types of sequences from L. perrieri, which was previously identified as a diploid species (2n = 24) (Katayama, 1995
A single product of c. 1 kb in size for the GPA1 gene was amplified in all species except for Hygroryza aristata and E. erecta, in which a single product of c.1.5 kb and 2.8 kb, respectively, was amplified. Therefore, we designed four additional primers (two forward and two reverse primers as listed in Table 1) to obtain the complete sequences of the two species. After comparing their sequences with those of other species, we found an
630-bp fragment inserted into intron 13 of GPA1 in H. aristata and three large segments inserted into intron 9 (
700 bp), intron 10 (
550 bp), and intron 13 (
620 bp), respectively, in E. erecta. Consequently, the sequence lengths for H. aristata and E. erecta are 1531 bp and 2860 bp, respectively, whereas the sequences for the remaining Oryzeae species vary from 871 to 990 bp. For the other outgroup, Phyllostachys aurea, a 280-bp fragment inserted into intron 10 was found, and the final sequence length of P. aurea is 1202 bp. Amplification products of GPA1 for all diploid species were sequenced directly except for four species (Luziola leiocarpa, Zizania latifolia, Zizaniopsis villanensis, and E. erecta) in which unreadable sequences were found, suggestive of some extent of heterogeneity. Therefore, products from the four species, along with those from the tetraploid species, were subsequently cloned and sequenced using a similar strategy to that mentioned earlier. For each tetraploid, we identified two different types of sequences at the GPA1 locus. Similar to the Adh2 gene, we also got two different types of GPA1 sequence for the diploid L. perrieri. After excluding the ambiguous regions, the final alignment of GPA1 sequences is 708 bp long (including 369 bp of exons), of which 322 sites (45.5%) are variable and 168 sites (23.7%) are informative.
As shown in Table 2, cytoplasmic DNA is much less variable than the nuclear genes. Of the two nuclear genes, Adh2 has higher parsimony information than GPA1, but the percentage of variable sites of GPA1 is higher than that of Adh2. Because the Adh2 and GPA1 sequences had more than 30% variable sites (excluding the outgroups), we examined possible saturation of the two data sets. A good linear relationship was found for both Adh2 and GPA1 (not shown), indicating that the two data sets are not saturated in Oryzeae. Each of the sampled sequences and their most appropriate models are also presented in Table 2.
Phylogenetic analysis of the cytoplasmic genes
Phylogenetic analyses of three individual data sets from cytoplasmic genomes (matK, trnL intron, and nad1 intron 2) provide approximately the same topology, in which the genera of the tribe fall into two main clades, but resolution was not sufficient to resolve relationships within clades (not shown). In the partition homogeneity test, there was no significant incongruence among the three data sets (P = 0.921). Therefore, we combined the three fragments with a total length of 3687 bp. The combined data set generates one most parsimonious tree (654 steps, consistency index [CI] = 0.806, retention index [RI] = 0.852) (Fig. 1). The MP analysis shows that the Oryzeae species form a clade with 100% bootstrap support (BS) and can be divided into two highly supported clades (100% and 80% BS, respectively). The first clade includes three genera, Leersia, Oryza, and Porteresia. In this clade, four Leersia species form a highly supported group (100% BS) with the diploid species L. tisserantii at the basal position. The species of Oryza and the monotypic genus Porteresia form a strongly supported group (83% BS). The second clade consists of the remaining eight genera. Although the resolution in this clade is slightly lower, there are still several clades with strong support, including the Prosphytochloa/Potamophila clade (100% BS), two Zizania species (100% BS), and the Zizaniopsis/Luziola clade (99% BS).
|
In addition, we conducted different rooting strategies using either Phyllostachys aurea or E. erecta or both as outgroups. In all analyses, ingroup topologies were not affected by outgroup selection, and the alternative outgroup was not nested within the ingroup, providing robust support for the monophyly of the tribe Oryzeae. For all analyses discussed here, both outgroups were used to root trees.
Phylogenetic analysis of the nuclear genes
Phylogenetic trees generated by individual data sets are essentially similar (not shown). Two and 18 most parsimonious trees were obtained for the Adh2 (911 steps, CI = 0.630, RI = 0.727) and the GPA1 (639 steps, CI = 0.714, RI = 0.789) data sets, respectively. The strict consensus of both gene sequences support the monophyly of the Oryzeae species, and two main clades within the tribe were recovered. The first main clade includes three genera, Leersia, Oryza, Porteresia (98% BS in Adh2 and 88% BS in GPA1). Four genera (Zizania, Rhynchoryza, Zizaniopsis, and Luziola) form a second main clade with varying degrees of bootstrap support. A major exception is the placement of the genus Hygroryza. This monotypic genus is sister to the Leersia/Oryza/Porteresia clade in the Adh2 tree but grouped with the Zizaniopsis/Luziola+Zizania/Rhynchoryza clade in the GPA1 tree (not shown but see Fig. 3). The Bayesian method generates similar topologies under the models of TrNef + G (Equal-frequency Tamura-Nei; G correction; Tamura and Nei, 1993
) (Adh2) and TrN + G (Tamura-Nei; G correction; Tamura and Nei, 1993
) (GPA1), respectively.
|
|
Monophyly and subdivision of the tribe Oryzeae
Since the tribe Oryzeae was proposed by Dumortier (1823)
, many genera have been added into Oryzeae and many of them were mistakenly treated. Pharus L. and Leptaspis R. Br. (including Scrotochloa Judziewicz) had been placed in Oryzeae by some researchers but later were treated as the separate tribe (Phareae) of Bambusoideae (Duistermaat, 1987
). Similarly, genera that were unreasonably treated as members of Oryzeae include Ehrharta Thunb., Microlaena R. Br., Tetrarrhena R. Br., Snowdenia C. E. Hubb., Achlaena Griseb., and Reynaudia Kunth (Duistermaat, 1987
). In addition, some authors (Avdulov, 1931
; Gardner, 1952
) had combined the tribes Oryzeae and Ehrharteae into a large rice tribe, namely Oryzeae sensu lato. Based on the study of anatomy and histology of leaf blades, Tateoka (1963)
suggested that the two tribes should be treated separately. Since then, the recognized genera in Oryzeae have varied from 10 to 13 (Tateoka, 1963
; Tzvelev, 1989
), and the classification of 12 genera by Clayton and Renvoize (1986)
and Vaughan (1994)
has been accepted by most researchers (Zhang and Second, 1989
; Duvall et al., 1993
; Kellogg and Watson, 1993
; Ge et al., 2002
). From a morphological point of view, monophyly of the tribe Oryzeae is supported by the following synapomorphies: one-flowered spikelets, which are compressed or terete, with a lemma and palea, and two well-developed bracts (degenerated flowers or sterile lemmas) (Pyrah, 1969
; Clayton and Renvoize, 1986
); the two tiny lobes below, sometimes referred to as cupules or glumes, appear to be expanded apices of the pedicels (Terrell et al., 2001
). Additionally, morphological and molecular studies also support the distinctiveness and monophyly of Oryzeae (Duistermaat, 1987
; Zhang and Second, 1989
; Kellogg and Watson, 1993
; Ge et al., 2002
). Nevertheless, because most previous researchers (Pyrah, 1969
; Clayton and Renvoize, 1986
; Duistermaat, 1987
; Zhang and Second, 1989
; Kellogg and Watson, 1993
) sampled limited species without some important genera included in their studies, the monophyly of Oryzeae was not determined with certainty. In this study, the analyses of both cytoplasmic and nuclear (either individual or combined) sequence data (Figs. 13) strongly support the monophyly of the tribe Oryzeae.
The subdivision of Oryzeae has varied among different researchers (Hitchcock, 1920
; Hubbard, 1934
, 1959
; Hitchcock and Chase, 1951
; Stebbins and Crampton, 1961
; Pyrah, 1969
; Terrell and Robinson, 1974
), and three main subdivision systems of Oryzeae are summaried in Table 3. Whether they were recognized as two groups at tribal (Oryzeae vs. Zizanieae) (Hitchcock, 1920
; Hitchcock and Chase, 1951
; Stebbins and Crampton, 1961
) or subtribal (Oryzinae vs. Zizaniinae) (Hubbard, 1934
, 1959
; Pyrah, 1969
) levels or treated as three subtribes (Terrell and Robinson, 1974
), they were all based on the monoecious vs. bisexual spikelet bifurcation. Based on cytoplasmic and nuclear sequence data, the present phylogenetic analyses of all genera in this tribe support recognition of two clades although the systematic position of Hygroryza remains ambiguous. The third subtribe Luziolinae proposed by Terrell and Robinson (1974)
is not justified because in the cytoplasmic and nuclear phylogenetic trees the Zizaniopsis/Luziola clade forms a clade with the Rhynchoryza/Zizania clade, which is sister to the Potamophila/Prosphytochloa/Chikusichloa clade (Figs. 1 and 2).
|
Delimitation and phylogenetic relationships of genera in Oryzeae
Within Oryzeae, the monophyly of some genera is supported by the present multiple gene analyses. Four Leersia species, including L. tisserantii and L. perrieri that were previously classified in Oryza (Tateoka, 1963
), form a highly supported monophyletic group, indicating that the transfer of L. tisserantii and L. perrieri from Oryza to Leersia by Launert (1965)
is justified (Zhang and Second, 1989
; Ge et al., 2002
). Luziola fluitans has been treated as the monotypic genus Hydrochloa and later considered a member of Luziola based on growth habitat, chromosome number, and simple racemose staminate inflorescences (Pohl and Davidse, 1971
). The work of Terrell and Robinson (1974)
also indicated that Hydrochloa could not be differentiated satisfactorily from Luziola as they appeared to differ mainly by the reduced inflorescence in Hydrochloa. In the present analysis, two Luziola species including Luziola fluitans form a monophyletic group with high statistical (100% BS and 1.00 PP) support, suggesting that Hydrochloa be retained in Luziola. Similarly, previous researchers (Hubbard, 1967
; Clayton, 1970
) divided the genus Potamophila R. Br. into three genera, including two monotypic genera (Prosphytochloa and Potamophila) and Maltebrunia Kunth. This treatment was questioned by Duistermaat (1987)
because these genera did not have fundamental differences in the structure of spikelets. In the present study, the close relationship between Potamophila and Prosphytochloa is highly supported in all the data sets (100 BS and 1.00 PP) (Figs. 1 3), and the sequence divergence between the two genera is much less than those between species of many other genera in the tribe (data not shown). Therefore, the present molecular data is in general accordance with the opinion that Prosphytochloa is within the generic limits of Potamophila (Duistermaat, 1987
; Vaughan, 1994
).
Whether the genus Zizania is more closely related to monoecious or to bisexual genera and its relationship with Zizaniopsis has not been resolved (Terrell and Robinson, 1974
; Zhang and Second, 1989
; Duvall et al., 1993
; Terrell et al., 2001
; Ge et al., 2002
). Zizaniopsis and Zizania are similar to one another in many characteristics including the lodicules, stamens, and pistil as well as the monoecious spikelets, and therefore Zizaniopsis was included under Zizania prior to 1871 (Weatherwax, 1929
). But subsequent morphological studies showed that the fruits of Zizania and Zizaniopsis were conspicuously different and suggested that the two genera belonged to separate genera or even separate tribes of grasses (Martin, 1946
; Terrell and Robinson, 1974
). Recent molecular studies also indicated that the two genera separated from each other (Zhang and Second, 1989
; Duvall et al., 1993
; Ge et al., 2002
). Our study of cytoplasmic and nuclear sequence data demonstrates that Zizaniopsis is more closely related to Luziola than to Zizania (Figs. 1 and 3). In addition, that Porteresia is nested within the clade of Oryza species further supports our previous statement that P. coarctata had a high affinity with Oryza species and should be treated as a member of Oryza rather than an independent monotypic genus (Ge et al., 2002
).
It is noteworthy that the monotypic genus Hygroryza forms a clade with the Oryzinae clade on the Adh2 tree but clusters with the Zizaniinae clade on the GPA1 tree (Fig. 3). The combined cytoplasmic phylogeny (Fig. 1) presents a topology approximately similar to that of GPA1. The incongruence among phylogenies may often be due to many different causes such as gene choice and insufficient data, convergent evolution, rapid diversification, hybridization and introgression, and lineage sorting as well as rate heterogeneity among taxa or among sites of the sequences used (for reviews, see Johnson and Soltis, 1998
; Wendel and Doyle, 1998
). Because the sequences we used represent multiple loci from three genomes and include the nuclear genes that evolve rapidly with their substitution rates distributed evenly over sites (not shown), insufficient data and rate heterogeneity among sites are not likely to be reasons for the incongruence. As indicated by many authors, lineage sorting or phylogenetic sorting is one of the important factors causing phylogenetic incongruence if a polymorphism transcends one or more organismal divergence events (Wendel and Doyle, 1998
). However, lineage sorting can be ruled out as a factor for discrepancy in the position of Hygroryza because it usually occurs at lower taxonomic ranks except in some exceptional cases (Rivers et al., 1993
; Wendel and Doyle, 1998
). Hybridization and subsequent polyploidization, which may result in reticulate evolution, would also be excluded for the cause of the incongruence because Hygroryza is a diploid (2n = 24) endemic to southeast Asia, and only one type of sequence was found for both Adh2 and GPA1 genes. Consequently, the ambiguous position of Hygroryza is probably the result of introgression and rapid speciation. Accumulating evidence shows that interspecific hybridization and introgression, though sometimes "cryptic," are exceptionally widespread and are powerful and creative forces in evolution in plants (Wendel and Doyle, 1998
; Wendel and Cronn, 2003
). Whether this incongruence reflects "real" underlying biological phenomena or just "spurious" insufficient data or some other unknown artifact is currently unclear and needs to be further investigated.
Implications for the origin and divergence of Oryzeae
Because the rice tribe is a diverse and geographically widespread group with many economically important species, its time and place of origin and subsequent divergence have been subjects of considerable discussion (Clayton, 1975
; Chang, 1985
; Second, 1985
; Duistermaat, 1987
). In the Poaceae, macrofossils provide greater taxonomic resolution than pollen because there are more diagnostic characters that facilitate broad comparisons in phylogenetic studies (Jacobs et al., 1999
). The oldest Oryzeae macrofossils have been reported from the end of the Eocene by Litke (1968)
who discovered fossilized bits of leaf epidermis in coal deposits of eastern Germany. Other reported macrofossils identified as Oryzeae include silicified anthoecia (fertile lemmas and paleas) found in Nebraska (North America) in the Miocene age (Thomasson, 1980
) and spikelets found in an excavation of Miocene age in Germany (Heer, 1855
). The former was described as Archaeoleersia nebraskensis Thomasson, which has features in common with those of the living Leersia ligularis Trin., and the latter as O. exasperata (A. Braun) Heer., which closely resembles O. meyeriana (Heer, 1855
). Therefore, the fossil reports suggest that Oryzeae had already differentiated by the end of the Eocene (36 to
40 mya).
Molecular tools provide an alternative source of information with which to estimate divergence times. Based on isozyme electrophoresis data, Second (1985)
estimated the divergence time of Oryza to be of the Miocene epoch (
15 mya) and the African A-genome species separated from the Asian A-genome species c. 7 mya. Although the divergence time of Poaceae has been estimated by molecular approaches (Bremer, 2002
; Gaut, 2002
), no attempt has been made to date the various events in Oryzeae using sequence data. Despite a number of limitations to the use of a molecular clock on sequence data, it is still useful for estimating divergence time if calibration can be made with confidence (Gaut, 2002
; Wendel and Cronn, 2003
).
To test the molecular clock hypothesis in Oryzeae, we performed relative rate tests for all the fragments used in this study. Results show that no rate heterogeneity exists at synonymous sites of matK sequences and at nonsynonymous sites of GPA1 sequences whereas either saturation or rate heterogeneity among lineages are found for other data sets. Therefore, the data sets of matK synonymous sites and GPA1 nonsynonymous sites are used for molecular dating. Another important aspect of divergence time estimation is clock calibration. Because the sequence divergence rates of these two loci either varied greatly among plant groups (matK) (Koch et al., 2001
) or are unknown (GPA1), we estimated the absolute substitution rates of matK synonymous sites and GPA1 nonsynonymous sites based on the assumption that maize and rice diverged 50 mya (Stebbins, 1981
; Gaut, 2002
). The estimated substitution rates for matK synonymous sites and GPA1 nonsynonymous sites are 1.96 x 109 and 1.02 x 109 substitutions per site per year, respectively. Using the molecular clocks of matK synonymous sites and GPA1 nonsynonymous sites, we calculated approximate divergence times for various lineages within the tribe and within Oryza. Based on matK synonymous sites, the mean sequence divergence between the two subtribes (the OryzaLeersia clade vs. the clade including the remaining genera) is 8.1%, which translates into a divergence of 20.5 mya, close to the 22.1 mya estimate based on GPA1 nonsynonymous sites. Interestingly, both matK and GPA1 estimates indicate that Oryza and Leersia separated from each other c. 14.2 mya. Within Oryza, we obtain the divergence time of 10.2 mya (matK) and 8.8 mya (GPA1), respectively, for the most basal split in the genus (O. granulata vs. the remaining species). These estimations suggest that the tribe Oryzeae diversified during the Miocene epoch, in good agreement with fossil evidence (Heer, 1855
; Litke, 1968
; Thomasson, 1980
). Therefore, the previously estimated divergence time among Oryza species of about 15 mya based on isozyme electrophoresis data (Second, 1985
) is probably inflated, and the assumption that the rice genus originated before the fragmentation of Gondwanaland (Chang, 1985
) is not justified.
Although molecular dating suffers from various errors including substitution rate heterogeneity among lineages, uncertainties of clock calibration, and other potential sources for errors (Gaut, 2002
; Wendel and Cronn, 2003
), the estimates in the present study provide a rough timeframe for understanding the evolutionary tempos and biogeographic history of this important group. To put the divergence of the rice tribe and rice genus into a temporal framework raises an interesting question regarding the biogeographic patterns in Oryzeae because both the tribe and particular genera, such as Oryza and Leersia, have achieved pantropical distribution. These distribution patterns imply that transoceanic dispersal or long-distance dispersal would contribute to the evolution and divergence of the Oryzeae given the fact that the two genera (Oryza and Leersia) branched from the remaining genera c. 20 mya and split from each other c. 14 mya. Long-distance dispersal events were also suggested for other recently diverged lineages such as Gossypium (Wendel and Cronn, 2003
). It is worthwhile to note that the closely related genera Chikusichloa, Potamophila, and Prosphytochloa (Fig. 3) that are distributed in eastern Asia, Australia, and Africa, respectively, are currently isolated geographically from one another by thousands of kilometers of open ocean. This similarly implicates oceanic dispersal as a major factor in their evolution and divergence and suggests that this clade is likely to originate in Asia because Chikusichola is at the basal position (Fig. 3). Compared with the Oryza/Leersia clade, low sequence divergence was found among the genera in the clade corresponding to the subtribe Zizaniinae (not shown), which can be translated into phylogenetically short interior branches (Fig. 3) and may be one of the main reasons for the low resolution within this clade. This phenomenon may reflect rapid divergence or radiation during the early divergence of the tribe. It is also evident that the dispersal from America to Asia occurred in the Rhynchoryza/ Zizania/Zizaniopsis/Luziola clade because Zizania is disjunctively distributed in eastern Asia and eastern North America while all the other genera are exclusively distributed in North and South America. Taken together, the rice tribe, in particular the worldwide genera Oryza and Leersia, is a successful group, as indicated by its extensive geographic ranges. Further studies will require a multidisciplinary approach and more extensive sampling to clarify their evolutionary history and biogeographic patterns.
Taxon (genome type of Oryza species); voucher or accession number; origin; chromosome number 2n; GenBank accession numbers: MatK; trnL; nad1; Adh2; GPA1.
O. sativa L. (AA); 30416c; Philippines; 24; AF148650a; AY792515; AY507930a; AF148602a; AY792541. O. glaberrima Steud. (AA); 104042c; Chad; 24; AF148654a; AY792516; AY792569; AF148606a; AY792542. O. meridionalis Ng (AA); 101147c; Australia; 24; AF148657a; AY792517; AY792570; AF148609a; AY792543. O. punctata Kotschy ex Steud. (BB); 104017c; Cameroon; 24; AF148661a; AY792518; AY507931a; AF148611a; AY299684a. O. officinalis Wall ex Watt (CC); 105085c; Philippines; 24; AF148658a; AY792519; AY507932a; AF148613a; AY188586a. O. rhizomatis Vaughan (CC); 105448c; Sri Lanka; 24; AF148660a; AY792520; AY792571; AF148614a; AY188587a. O. australiensis Domin (EE); 101144c; Australia; 24; AF148667a; AY792521; AY507937a; AF148623a; AY188596a. O. brachyantha A. Chev. et Roehr. (FF); 105151c; Sierra Leone; 24; AF148670a; AY792522; AY507934a; AF148632a; AY792546. O. granulata Nees et Arn. ex Watt. (GG); 2422; China; 24; AF148674a; AY792523; AY507938a; AF148631a; AY188597a. Porteresia coarctata Tateokab (HHKK); 104502c; Bangladesh; 48; AF148669a; AY792524; AY507935a; AF148627a, AF148628a; AY792544, AY792545. Leersia oryzoides (L.) SW. (); GS0203; China; 48; AY792566; AY792525; AY792572; AY792577; AY792547, AY792548. Leersia perrieri (A. Camus.) Launert (); 105164c; Madagascar; 24; AF148677a; AY792526; AY792573; AY792578, AY792579; AY792549, AY792550. Leersia hexandra Sw. (); 105252c; Philippines; 48; AF489909a; AY792527; AY507940a; AY792580; AY792551, AY792552. Leersia tisserantii (A. Chev.) Launert (); 105610c; Cameroon; 24; AF489910a; AY792528; AY507939a; AY792581; AY792553. Prosphytochloa prehensis (Nees) Schweickb (); unknownc; S. Africa; 24; AF489916a; AY792529; AY507943a; AY792582; AY792554. Potamophila parviflora R. Br.b (); 85424c; Australia; 24; AF489914a; AY792530; AY507944a; AY792583; AY792555. Chikusichloa aquatica Koidz. (); 106186c; Japan; 24; AF489912a; AY792531; AY507948a; AY792584; AY792556. Rhynchoryza subulata (Nees) Baillonb (); 100913c; Argentina; 24; AF148675a; AY792532; AY507945a; AY792585; AY792557. Zizania aquatica L. (); J. Alexander 200301; MA, USA; 30; AF164393a; AY792533; AY792574; AY792586; AY792558. Zizania latifolia (Griseb.) Turcz. ex Stapf (); GS0202; China; 34; AY092064a; AY792534; AY507941a; AY792587; AY792559. Zizaniopsis villanensis Quarin; 85425c (); Argentina; 24; AF489917a; AY792535; AY507942a; AF148635a; AY792560. Luziola leiocarpa Lindm. (); 82043c; Argentina; 24; AF489911a; AY792536; AY507946a; AY792588; AY792561. Luziola fluitans (Michx.) E. E. Terrell & H. Robinson (); 8434; USA; 24; AY792567; AY792537; AY792575; AY792589; AY792562. Hygroryza aristata (Retz.) Neesb (); 105460c; Sri Lanka; 24; AF489913a; AY792538; AY507947a; AY792590; AY792563. Ehrharta erecta Lam.; B. Bartholomeal 9130; USA; 48; AY792568; AY792539; AY792576; AY792591; AY792564. Phyllostachys aurea Riviére & C. Riviére (); GS0204; China; 48; AF1643901a; AY792540; AY507949a; AY792592; AY792565.
a used in our previous studies or downloaded from GenBank.
b monotypic genera.
c provided by the Genetic Resources Center of the International Rice Research Institute (IRRI) at Los Banos, Philippines and other accessions are obtained by the authors.
FOOTNOTES
1 The authors thank Bao-Rong Lu, Jun Wen, Bruce G. Baldwin, Curtis A. Jones, and Jian-Hua Li for help in collecting materials; Zhen-Yu Li and Fu-Wu Xing for identifying the samples; and Shou-Zhou Zhang, Fu-Min Zhang, Ying Bao, Qi-Hui Zhu, Fu-Hong An, and Yang Yu for technical assistance. We also thank De-Yuan Hong, Tao Sang, and Bao-Rong Lu for their encouragement and support. Comments and suggestions from Paul M. Peterson and an anonymous reviewer are sincerely appreciated. We are grateful to the International Rice Research Institute (Manila, Philippines) for providing most seed and leaf samples used in this study. This study was supported by grants from the Program for Key International Science and Technique Cooperation Project of P. R. China (2001CB711103), the National Natural Science Foundation of China (30025005, 30121003), and the Chinese Academy of Sciences (kscxz-sw-101A). ![]()
4 Author for correspondence (e-mail: gesong{at}ibcas.ac.cn
) ![]()
LITERATURE CITED
Avdulov N. 1931 Karyosystematische Untersuchung der Familie Gramineen. Bulletin Applied Botanical Genetics and Plant Breed 44: 1-428
Bao Y. S. Ge 2004 Origin and phylogeny of Oryza species with the CD genome based on multiple-gene sequence data. Plant Systematics and Evolution 249: 55-66[CrossRef][Web of Science]
Bennetzen J. 2002 Opening the door to comparative plant biology. Science 296: 60-62
Bremer K. 2002 Gondwanan evolution of the grass alliance of families (Poales). Evolution 56: 1374-1387[CrossRef][Web of Science][Medline]
Chandler V. L. S. Wessler 2001 Grasses. A collective model genetic system. Plant Physiology 125: 1155-1156
Chang T. T. 1985 Crop history and genetic conservation: ricea case study. Iowa State Journal of Research 59: 425-455[Web of Science]
Clayton W. D. 1970 Flora of tropical East Africa. Grasses 1: 23-33
Clayton W. D. 1975 Chorology of the genera of Gramineae. Kew Bulletin 30: 111-132[CrossRef]
Clayton W. D. S. A. Renvoize 1986 Genera Graminum. Kew Bulletin Additional Series XIII 1-389
Cunningham C. 1997 Can three incongruence tests predict when data should be combined?. Molecular Biology and Evolution 14: 733-740[Abstract]
Demesure B. N. Sodzi R. J. Petit 1995 A set of universal primers for amplification of polymorphic non-coding regions of mitochondrial and chloroplast DNA in plants. Molecular Ecology 4: 129-131[Medline]
Doyle J. J. J. L. Doyle 1999 Nuclear protein-coding genes in phylogeny reconstruction and homology assessment: some examples from Leguminosae. In P. M. Hollingsworth, R. M. Bateman, and R. J. Gornall [eds.], Molecular systematics and plant evolution, 229254. Taylor and Francis, London, UK
Duistermaat H. 1987 A revision of Oryza (Gramineae) in Malesia and Australia. Blumea 32: 157-193[Web of Science]
Dumortier B. C. 1823 Observation sur les graminées de la flore Belgique. J. Casterman, Tournay, Belgium
Duvall M. R. P. M. Peterson E. E. Terrell A. H. Christensen 1993 Phylogeny of North American oryzoid grasses as construed from maps of plastid DNA restriction sites. American Journal of Botany 80: 83-88[CrossRef][Web of Science]
Farris J. S. M. Källersjö A. G. Kluge C. Bult 1995 Testing significance of incongruence. Cladistics 10: 315-319[CrossRef][Web of Science]
Felsenstein J. 1985 Confidence limits on phylogeny: an approach using the bootstrap. Evolution 39: 783-791[CrossRef][Web of Science]
Fujisawa Y. T. Kato S. Ohki A. Ishikawa H. Kitano T. Sasaki T. Asahi Y. Iwasaki 1999 Suppression of the heterotrimeric G protein causes abnormal morphology, including dwarfism, in rice. Proceedings of the National Academy of Sciences, USA 96: 7575-7580
Gale M. D. K. M. Devos 1998 Comparative genetics in the grasses. Proceedings of the National Academy of Sciences, USA 95: 1971-1974
Gardner C. A. 1952 Flora of Western Australia, vol. 1, part 1. Gramineae. Government Press, Perth, Australia
Gaut B. S. 2002 Evolutionary dynamics of grass genomes. New Phytologist 154: 15-28[CrossRef][Web of Science]
Ge S. A. Li B. R. Lu S. Z. Zhang D. Y. Hong 2002 A phylogeny of the rice tribe Oryzeae (Poaceae) based on matK sequence data. American Journal of Botany 89: 1967-1972
Ge S. T. Sang B. R. Lu D. Y. Hong 1999 Phylogeny of rice genomes with emphasis on origins of allotetraploid species. Proceedings of the National Academy of Sciences, USA 96: 14400-14405
Ge S. T. Sang B. R. Lu D. Y. Hong 2001 Phylogeny of the genus Oryza as revealed by molecular approaches. In G. S. Khush, D. S. Brar, and B. Hardy [eds.], Rice genetics IV. Proceedings of the Fourth International Rice Genetics Symposium, 89105. IRRI, Los Banos, Philippines
GPWG. 2001 Phylogeny and subfamilial classification of the grasses (Poaceae). Annals of the Missouri Botanical Garden 88: 373-457[CrossRef][Web of Science]
Grob G. B. J. B. Gravendeel M. C. M. Eurlings 2004 Potential phylogenetic utility of the nuclear FLORICAULA/LEAFY second intron: comparison with three chloroplast DNA regions in Amorphophallus (Araceae). Molecular Phylogenetics and Evolution 30: 12-23
Gugerli F. J. Senn M. Anzidei A. Madaghiele U. Büchler C. Sperisen G. G. Vendramin 2001 Chloroplast microsatellites and mitochondrial nad1 intron 2 sequences indicate congruent phylogenetic relationships among Swiss stone pine (Pinus cembra), Siberian stone pine (Pinus sibirica), and Siberian dwarf pine (Pinus pumila). Molecular Ecology 10: 1489-1497[CrossRef][Medline]
Hasegawa M. H. Kishino T. Yano 1985 Dating the human-ape split by a molecular clock of mitochondrial DNA. Journal of Molecular Evolution 22: 160-174[CrossRef][Web of Science][Medline]
Hass B. L. J. C. Pires R. Porter R. L. Phillips S. A. Jackson 2003 Comparative genetics at the gene and chromosome levels between rice (Oryza sativa) and wildrice (Zizania palustris). Theoretical and Applied Genetics 107: 773-782
Heer D. 1855 Flora tertiaria helvetiae. Winterthur, Paris, France
Hitchcock A. S. 1920 The genera of grasses of the United States. Bulletin 772, Miscellaneous Publication. USDA, Washington, D.C., USA
Hitchcock A. S. A. Chase 1951 Manual of the grasses of the United States, 2nd ed. Miscellaneous Publication. USDA, Washington, D.C., USA
Hubbard C. E. 1934 Gramineae. In J. Hutchinson [ed.], The families of flowering plants, vol. II, Monocotyledons, 199229. Macmillan, London, UK
Hubbard C. E. 1959 Gramineae. In J. Hutchinson [ed.], The families of flowering plants, vol. II, Monocotyledons, 710741. Clarendon Press, London, UK
Hubbard C. E. 1967 Tabula 3647: Snowdenia petitiana (A. Rich.) C. E. Hubbard. Hooker's Icones Plantarum 37, t3647
Huelsenbeck J. P. F. R. Ronquist 2001 MrBayes: Bayesian inference of phylogeny. Bioinformatics 17: 754-755
Jacobs B. F. J. D. Kingston L. L. Jacobs 1999 The origin of grass-dominated ecosystems. Annals of the Missouri Botanical Garden 86: 590-643[CrossRef][Web of Science]
Johnson L. A. D. E. Soltis 1998 Assessing congruence: empirical examples from molecular data. In D. E. Soltis, P. S. Soltis, and J. J. Doyle [eds.], Molecular systematics of plants II: DNA sequencing, 297 348. Kluwer, Dordrecht, Netherlands
Katayama T. 1995 Cytogenetical studies on the genus Oryza. XIV. Intergeneric hybridizations between tetraploid Oryza species and diploid Leersia species. Japanese Journal of Genetics 70: 47-55[CrossRef][Medline]
Kato H. K. Oginuma Z. Gu B. Hammel H. Tobe 1998 Phylogenetic relationships of Betulaceae based on matK sequences with particular reference to the positions of Ostryopsis. Acta Phytotaxonomica et Geobotanica 49: 89-97
Kellogg E. A. 1998 Relationships of cereal crops and other grasses. Proceedings of the National Academy of Sciences, USA 95: 2005-2010
Kellogg E. A. L. Watson 1993 Phylogenetic studies of a large data set. I. Bambusoideae, Andropogonodae, and Pooideae (Gramineae). Botanical Review 59: 273-343
Kennard W. R. L. Phillips R. A. Porter A. W. Grombacher 1999 A comparative map of wild rice (Zizania palustris L. 2n = 2x = 30). Theoretical and Applied Genetics 101: 677-684[Web of Science]
Koch M. B. Haubold T. Mitchell-Olds 2001 Molecular systematics of the Brassicaceae: evidence from coding plastidic matK and nuclear Chs sequences. American Journal of Botany 88: 534-544
Kosakovsky Pond S. L. S. D. W. Frost S. V. Muse 2005 HyPhy: hypothesis testing using phylogenies. Bioinformatics 21: 676-679
Kumar S. K. Tamura I. B. Jakobsen M. Nei 2001 MEGA2: molecular evolutionary genetics analysis software. Bioinformatics 17: 1244-1245
Launert E. 1965 A survey of the genus Leersia in Africa. Senckenbergiana Biology 46: 129-153
Li H. W. C. C. Chen H. K. Wu K. C. L. Lu 1963 Cytogenetical studies of Oryza sativa L. and its related species. 5. Differential condensation and chromosome pairing in the hybrid O. sativa x O. australiensis. Cytologia 28: 264-277[Web of Science]
Litke R. 1968 Über den Nachweis tertiärer Gramineen. Monatsberichte der Deutschen Akademie der Wissenschaften zu Berlin 10: 462-471
Martin A. C. 1946 The comparative internal morphology of seeds. American Midland Naturalist 36: 513-660[CrossRef][Web of Science]
Nayar N. M. 1973 Origin and cytogenetics of rice. Advances in Genetics 17: 153-292
Pohl R. W. G. Davidse 1971 Chromosome numbers of Costa Rican grasses. Brittonia 23: 293-324
Posada D. K. A. Crandall 1998 MODELTEST: testing the model of DNA substitution. Bioinformatics 14: 817-818
Pyrah G. L. 1969 Taxonomic and distributional studies in Leersia (Gramineae). Iowa State Journal of Science 44: 215-270
Ramanujam S. 1938 Cytogenetical studies in the Oryzeae. Annals of Botany (London) 2: 107-125
Rensink W. A. C. R. Buell 2004 Arabidopsis to rice. Applying knowledge from a weed to enhance our understanding of a crop species. Plant Physiology 135: 622-629
Rivers B. A. R. Bernatzky S. J. Robinson W. Jahnen-Dechent 1993 Molecular diversity at the self-incompatibility locus is a salient feature in natural populations of wild tomato (Lycopersicon peruvianum). Molecular Genetics and Genomics 238: 419-427[CrossRef]
Robinson-Rechavi M. D. Huchon 2000 RRTree: relative-rate tests between groups of sequences on a phylogenetic tree. Bioinformatics 16: 296-297
Sang T. 2002 Utility of low-copy nuclear gene sequences in plant phylogenetics. Critical Reviews in Biochemistry and Molecular Biology 37: 121-147[CrossRef][Web of Science][Medline]
Schmidt R. 2002 Synteny: recent advances and future prospects. Current Opinion in Plant Biology 3: 97-102
Second G. 1985 Geographic origins, genetic diversity and the molecular clock hypothesis in the Oryzeae. In P. Jacquard, G. Heim, and J. Antonovics [eds.], Genetic differentiation and dispersial in plants, 4156. Springer-Verlag, Berlin, Germany
Seo H. S. H. Y. Kim J. Y. Jeong S. Y. Lee M. J. Cho J. D. Bahk 1995 Molecular cloning and characterization of RGA1 encoding a G protein alpha subunit from rice (Oryza sativa L. IR-36). Plant Molecular Biology 27: 1119-1131[CrossRef][Web of Science][Medline]
Shimamoto K. J. Kyozuka 2002 Rice as a model for comparative genomics of plants. Annual Review of Plant Biology 53: 399-419[CrossRef][Medline]
Soltis D. E. P. S. Soltis 1998 Choosing an approach and an appropriate gene for phylogenetic analysis. In D. E. Soltis, P. S. Soltis, and J. Doyle [eds.], Molecular systematics of plants II: DNA sequencing. Kluwer, Dordrecht, Netherlands
Stebbins G. B. Crampton 1961 A suggested revision of grass genera of temperate North America. Recent Advances in Botany 1: 133-145
Stebbins G. L. 1981 Coevolution of grasses and herbivores. Annals of the Missouri Botanical Garden 68: 75-86[CrossRef][Web of Science]
Swofford D. L. 2001 PAUP*: phylogenetic analysis using parsimony (* and other methods), version 4.0b10. Sinauer, Sunderland, Massachusetts, USA
Tamura K. M. Nei 1993 Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Molecular Biological and Evolution 10: 512-526
Tateoka T. 1963 Notes on some grasses. XIII. Relationship between Oryzeae and Ehrharteae, with special reference to leaf anatomy and histology. Botanical Gazette 124: 264-270
Terrell E. E. P. M. Peterson W. P. Wergin 2001 Epidermal features and spikelet micromorphology in Oryza and related genera (Poaceae: Oryzeae). Smithsonian Contributions to Botany 91: 1-50
Terrell E. E. H. Robinson 1974 Luziolinae, a new subtribe of oryzoid grass. Bulletin of the Torrey Botanical Club 101: 235-245[CrossRef][Web of Science]
Thomasson J. R. 1980 Archaeolersia nebraskensis gen. et sp. nov. (GramineaeOryzeae), a new fossil grass from the late Tertiary of Nebraska. American Journal of Botany 67: 876-882
Thompson J. D. T. J. Gibson F. Plewniak F. Jeanmougin D. G. Higgins 1997 The Clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research 25: 4876-4882
Tzvelev N. N. 1989 The system of grasses (Poaceae) and their evolution. Botanical Review 55: 141-204
Vaughan D. A. 1994 The wild relative of rice: a genetic resources handbook. International Rice Research Institute, Manila, Philippines
Watson L. M. J. Dallwitz 1999 Grass genera of the world: descriptions, illustrations, identification, and information retrieval; including synonyms, morphology, anatomy, physiology, phytochemistry, cytology, classification, pathogens, world and local distribution, and references. Version: 18 August 1999, website http://biodiversity.uno.edu/ delta/
Weatherwax P. 1929 The morphology of the spikelets of six genera of Oryzeae. American Journal of Botany 16: 547-555[CrossRef][Web of Science]
Wendel J. F. R. C. Cronn 2003 Polyploidy and the evolutionary history of cotton. Advances in Agronomy 78: 139-186
Wendel J. F. J. J. Doyle 1998 Phylogenetic incongruence: window into genomes history and molecular evolution. In D. E. Soltis, P. S. Soltis, and J. Doyle [eds.], Molecular systematics of plants II: DNA sequencing, 265296. Kluwer Academic, Dordrecht, Netherlands
Xia X. Z. Xie 2001 DAMBE: data analysis in molecular biology and evolution. Journal of Heredity 92: 371-373
Zhang S. H. G. Second 1989 Phylogenetic analysis of the tribe Oryzeae: total chloroplast DNA restriction fragment analysis (a preliminary report). Rice Genetics Newsletter 6: 76-80
This article has been cited by other articles:
![]() |
Y.-p. Zhang and S. Ge Molecular evolution study in China: progress and future promise Phil Trans R Soc B, June 29, 2007; 362(1482): 973 - 986. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |