Phylogenetic Analysis of AA-genome Oryza Species (Poaceae) Based on Chloroplast, Mitochondrial, and Nuclear DNA Sequences

African Cultivated, Wild and Weedy Rice (Oryza spp.): Anticipating Further Genomic Studies

Rice is a staple crop in sub-Saharan Africa, and it is mostly produced by Asian cultivars of Oryza sativa that were introduced to the continent around the fifteenth or sixteenth century. O. glaberrima, the native African rice, has also been planted due to its valuable traits of insect and drought tolerance. Due to competition and resistance evolution, weedy rice has evolved from O. sativa and O. glaberrima, posing an increasing threat to rice production. This paper provides an overview of current knowledge on the introduction and domestication history of cultivated rice in Africa, as well as the genetic properties of African weedy rice that invades paddy fields. Recent developments in genome sequencing have made it possible to uncover findings about O. glaberrima's population structure, stress resilience genes, and domestication bottleneck. Future rice genomic research in Africa should prioritize producing more high-quality reference genomes, quantifying the impact of crop-wild hybridization, elucidating weed adaptation mechanisms through resequencing, and establishing a connection between genomic variation and stress tolerance phenotypes to accelerate breeding efforts.

Solving the mystery of Obake rice in Africa: population structure analyses of Oryza longistaminata reveal three genetic groups and evidence of both recent and ancient introgression with O. sativa

The undomesticated rice relative Oryza longistaminata is a valuable genetic resource for the improvement of the domesticated Asian rice, Oryza sativa. To facilitate the conservation, management, and use of O. longistaminata germplasm, we sought to quantify the population structure and diversity of this species across its geographic range, which includes most of sub-Saharan Africa, and to determine phylogenetic relationships to other AA-genome species of rice present in Africa, including the prevalence of interspecific hybridization between O. longistaminata and O. sativa. Though past plant breeding efforts to introgress genes from O. longistaminata have improved biotic stress resistance, ratooning ability, and yield in O. sativa, progress has been limited by substantial breeding barriers. Nevertheless, despite the strong breeding barriers observed by plant breeders who have attempted this interspecific cross, there have been multiple reports of spontaneous hybrids of O. sativa and O. longistaminata (aka "Obake") obtained from natural populations in Africa. However, the frequency and extent of such natural introgressions and their effect on the evolution of O. longistaminata had not been previously investigated. We studied 190 O. longistaminata accessions, primarily from the International Rice Research Institute genebank collection, along with 309 O. sativa, 25 Oryza barthii, and 83 Oryza glaberrima control outgroups, and 17 control interspecific O. sativa/O. longistaminata hybrids. We analyzed the materials using 178,651 single-nucleotide polymorphisms (SNPs) and seven plastid microsatellite markers. This study identified three genetic subpopulations of O. longistaminata, which correspond geographically to Northwestern Africa, Pan-Africa, and Southern Africa. We confirmed that O. longistaminata is, perhaps counterintuitively, more closely related to the Asian species, O. sativa, than the African species O. barthii and O. glaberrima. We identified 19 recent spontaneous interspecific hybrid individuals between O. sativa and O. longistaminata in the germplasm sampled. Notably, the recent introgression between O. sativa and O. longistaminata has been bidirectional. Moreover, low levels of O. sativa alleles admixed in many predominantly O. longistaminata accessions suggest that introgression also occurred in the distant past, but only in Southern Africa.

A large deletion within intron 20 sequence of single-copy PolA1 gene as a useful marker for the speciation in Oryza AA-genome species

Oryza AA-genome complex comprises five wild species, O. rufipogon, O. barthii, O. longistaminata, O. glumaepatula, and O. meridionalis. Evolutionary relationships among these five wild species have remained contentious and inconclusive. We found that intron 20 of PolA1, a single-copy nuclear gene, was short (S-type: 141-142 bp) in O. rufipogon, O. barthii, and O. glumaepatula, while long (L-type: ca. 1.5 kb) introns were apparent in O. longistaminata and O. meridionalis. Because Oryza species containing BB, CC, EE, FF, and GG genome showed L-type introns, the S-type intron was probably derived from the L-type intron by the deletion of a 1.4 kb fragment through intramolecular homologous recombination between two tandem TTTTGC repeats. Excluding the large deletion sequence, intron 20 sequence of O. barthii was identical to that of O. longistaminata. As more than 3,470 accessions of O. rufipogon and O. sativa also contained the same intron 20 sequence with O. longistaminata except for single T-nucleotide deletion, which was shared with O. glumaepatuala, the deletion of the T-nucleotide probably occurred in the L-type intron 20 of O. logistaminata. Deletions of a large 1.4 kb fragment and single T-nucleotide within the intron 20 of PolA1 gene were considered as useful DNA markers to study the evolutionary relationships among Oryza AA-genome species.

Comparative genomic study on the complete plastomes of four officinal Ardisia species in China

Ardisia Sw. (Primulaceae) is naturally distributed in tropical and subtropical areas. Most of them possess edible and medicinal values and are popular in clinical and daily use in China. However, ambiguous species delineation and genetic information limit the development and utilization of this genus. In this study, the chloroplast genomes of four Ardisia species, namely A. gigantifolia Stapf, A. crenata Sims, A. villosa Roxb. and A. mamillata Hance, were sequenced, annotated, and analyzed comparatively. All the four chloroplast genomes possess a typical quadripartite structure, and each of the genomes is about 156 Kb in size. The structure and gene content of the Ardisia plastomes were conservative and showed low sequence divergence. Furthermore, we identified five mutation hotspots as candidate DNA barcodes for Ardisia, namely, trnT-psbD, ndhF-rpl32, rpl32-ccsA, ccsA-ndhD and ycf1. Phylogenetic analysis based on the whole-chloroplast genomes data showed that Ardisia was sister to Tapeinosperma Hook. f. In addition, the results revealed a great topological profile of Ardisia's with strong support values, which matches their geographical distribution patterns. Summarily, our results provide useful information for investigations on taxonomic differences, molecular identification, and phylogenetic relationships of Ardisia plants.

Genetics and Genomics of African Rice (Oryza glaberrima Steud) Domestication

African rice (Oryza glaberrima Steud) is one of the two independently domesticated rice species, the other one being Asian rice (Oryza sativa L.). Despite major progress being made in understanding the evolutionary and domestication history of African rice, key outstanding issues remain controversial. There appears to be an underlying difficulty in identifying the domestication centre and number of times the crop has been domesticated. Advances in genomics have provided unprecedented opportunities for understanding the genetic architecture of domestication related traits. For most of the domestication traits, the underlying genes and mutations have been identified. Comparative analysis of domestication genes between Asian and African rice has revealed that the two species went through an independent but convergent evolution process. The genetic and developmental basis of some of the domestic traits are conserved not only between Asian and African rice but also with other domesticated crop species. Analysis of genome data and its interpretation is emerging as a major challenge facing studies of domestication in African rice as key studies continue giving contradictory findings and conclusions. Insights obtained on the domestication of this species are vital for guiding crop improvement efforts.

Comparisons of molecular diversity indices, selective sweeps and population structure of African rice with its wild progenitor and Asian rice

The extent of molecular diversity parameters across three rice species was compared using large germplasm collection genotyped with genomewide SNPs and SNPs that fell within selective sweep regions. Previous studies conducted on limited number of accessions have reported very low genetic variation in African rice (Oryza glaberrima Steud.) as compared to its wild progenitor (O. barthii A. Chev.) and to Asian rice (O. sativa L.). Here, we characterized a large collection of African rice and compared its molecular diversity indices and population structure with the two other species using genomewide single nucleotide polymorphisms (SNPs) and SNPs that mapped within selective sweeps. A total of 3245 samples representing African rice (2358), Asian rice (772) and O. barthii (115) were genotyped with 26,073 physically mapped SNPs. Using all SNPs, the level of marker polymorphism, average genetic distance and nucleotide diversity in African rice accounted for 59.1%, 63.2% and 37.1% of that of O. barthii, respectively. SNP polymorphism and overall nucleotide diversity of the African rice accounted for 20.1-32.1 and 16.3-37.3% of that of the Asian rice, respectively. We identified 780 SNPs that fell within 37 candidate selective sweeps in African rice, which were distributed across all 12 rice chromosomes. Nucleotide diversity of the African rice estimated from the 780 SNPs was 8.3 × 10^-4, which is not only 20-fold smaller than the value estimated from all genomewide SNPs (π = 1.6 × 10^-2), but also accounted for just 4.1%, 0.9% and 2.1% of that of O. barthii, lowland Asian rice and upland Asian rice, respectively. The genotype data generated for a large collection of rice accessions conserved at the AfricaRice genebank will be highly useful for the global rice community and promote germplasm use.

Morphological and molecular dissection of wild rices from eastern India suggests distinct speciation between O. rufipogon and O. nivara populations

The inter relationships between the two progenitors is interesting as both wild relatives are known to be the great untapped gene reservoirs. The debate continues on granting a separate species status to Oryza nivara. The present study was conducted on populations of Oryza rufipogon and Oryza nivara from Eastern India employing morphological and molecular characteristics. The cluster analysis of the data on morphological traits could clearly classify the two wild forms into two separate discrete groups without any overlaps i.e. lack of intermediate forms, suggesting the non-sympatric existence of the wild forms. Amplification of hyper variable regions of the genome could reveal 144 alleles suggesting high genetic diversity values (average He = 0.566). Moreover, with 42.37% of uncommon alleles between the two wild relatives, the molecular variance analysis (AMOVA) could detect only 21% of total variation (p < 0.001) among them and rest 59% was within them. The population structure analysis clearly classified these two wild populations into two distinct sub-populations (K = 2) without any overlaps i.e. lack of intermediate forms, suggesting the non-sympatric existence of the wild forms. Clear differentiation into two distinct groups indicates that O. rufipogon and O. nivara could be treated as two different species.

Complete chloroplast and ribosomal sequences for 30 accessions elucidate evolution of Oryza AA genome species

Cytoplasmic chloroplast (cp) genomes and nuclear ribosomal DNA (nR) are the primary sequences used to understand plant diversity and evolution. We introduce a high-throughput method to simultaneously obtain complete cp and nR sequences using Illumina platform whole-genome sequence. We applied the method to 30 rice specimens belonging to nine Oryza species. Concurrent phylogenomic analysis using cp and nR of several of specimens of the same Oryza AA genome species provides insight into the evolution and domestication of cultivated rice, clarifying three ambiguous but important issues in the evolution of wild Oryza species. First, cp-based trees clearly classify each lineage but can be biased by inter-subspecies cross-hybridization events during speciation. Second, O. glumaepatula, a South American wild rice, includes two cytoplasm types, one of which is derived from a recent interspecies hybridization with O. longistminata. Third, the Australian O. rufipogan-type rice is a perennial form of O. meridionalis.

Exploring natural selection to guide breeding for agriculture

Climate change threatens reduced crop production and poses major challenges to food security. The breeding of climate-resilient crop varieties is increasingly urgent. Wild plant populations evolve to cope with changes in their environment due to the forces of natural selection. This adaptation may be followed over time in populations at the same site or explored by examining differences between populations growing in different environments or across an environmental gradient. Survival in the wild has important differences to the objective of agriculture to maximize crop yields. However, understanding the nature of adaptation in wild populations at the whole genome level may suggest strategies for crop breeding to deliver agricultural production with more resilience to climate variability.

Chloroplast Genome of Novel Rice Germplasm Identified in Northern Australia

Rice (Oryza sativa L.) was probably domesticated from O. rufipogon in Asia in the last 10,000 years. Relatives of cultivated rice (A genome species of Oryza) are found in South America, Africa, Australia and Asia. These A genome species are the close relatives of cultivated rice and represent the effective gene pool for rice improvement. Members of this group in Northern Australia include, an annual species, O. meridionalis, and two recently distinguished perennial taxa, to one of which the name O. rufipogon has been applied and the other a perennial form of O. meridionalis. Comparison of whole chloroplast genome sequences of these taxa has now been used to determine the relationships between the wild taxa and cultivated rice. The chloroplast genomes of the perennials were both found to be distinguished from O. rufipogon from Asia by 124 or 125 variations and were distinguished from each other by 53 variations. These populations have remained isolated from the overwhelming genetic impact of the large domesticated rice populations in Asia and may be unique descendants of the gene pool from which domesticated rice arose. The conservation of this wild genetic resource may be critical for global food security.

Conservation and utilization of African Oryza genetic resources

Africa contains a huge diversity of both cultivated and wild rice species. The region has eight species representing six of the ten known genome types. Genetic resources of these species are conserved in various global germplasm repositories but they remain under collected and hence underrepresented in germplasm collections. Moreover, they are under characterized and therefore grossly underutilized. The lack of in situ conservation programs further exposes them to possible genetic erosion or extinction. In order to obtain maximum benefits from these resources, it is imperative that they are collected, efficiently conserved and optimally utilized. High throughput molecular approaches such as genome sequencing could be employed to more precisely study their genetic diversity and value and thereby enhance their use in rice improvement. Oryza sativa was the first crop plant to have its reference genome sequence released marking a major milestone that opened numerous opportunities for functional characterization of the entire rice genome. Studies have however demonstrated that one reference genome sequence is not enough to fully explore the genetic variation in the Oryza genus, hence the need to have reference sequences for other species in the genus. An overview of the state of conservation and utilization of African Oryza is hereby presented. Progress in the release of reference genome sequences for these species is also highlighted.

Phylogenetic relationships and genome divergence among the AA- genome species of the genus Oryza as revealed by 53 nuclear genes and 16 intergenic regions

Rapid radiations have long been regarded as the most challenging issue for elucidating poorly resolved phylogenies in evolutionary biology. The eight diploid AA- genome species in the genus Oryza represent a typical example of a closely spaced series of recent speciation events in plants. However, questions regarding when and how they diversified have long been an issue of extensive interest but remain a mystery. Here, a data set comprising >60 kb of 53 singleton fragments and 16 intergenic regions is used to perform phylogenomic analyses of all eight AA- genome species plus four diploid Oryza species with BB-, CC-, EE- and GG- genomes. We fully reconstruct phylogenetic relationships of AA- genome species with confidence. Oryza meridionalis, native to Australia, is found to be the earliest divergent lineage around 2.93 mya, whereas O. punctata, a BB- genome species, serves as the best outgroup to distinguish their phylogenetic relationships. They separated from O. punctata approximately 9.11 mya during the Miocene epoch, and subsequently radiated to generate the entire AA- genome lineage diversity. The success in resolving the phylogeny of AA- genome species highlights the potential of phylogenomics to determine their divergence and evolutionary histories.

Local differentiation amidst extensive allele sharing in Oryza nivara and O. rufipogon

Genetic variation patterns within and between species may change along geographic gradients and at different spatial scales. This was revealed by microsatellite data at 29 loci obtained from 119 accessions of three Oryza series Sativae species in Asia Pacific: Oryza nivara Sharma and Shastry, O. rufipogon Griff., and O. meridionalis Ng. Genetic similarities between O. nivara and O. rufipogon across their distribution are evident in the clustering and ordination results and in the large proportion of shared alleles between these taxa. However, local-level species separation is recognized by Bayesian clustering and neighbor-joining analyses. At the regional scale, the two species seem more differentiated in South Asia than in Southeast Asia as revealed by F ST analysis. The presence of strong gene flow barriers in smaller spatial units is also suggested in the analysis of molecular variance (AMOVA) results where 64% of the genetic variation is contained among populations (as compared to 26% within populations and 10% among species). Oryza nivara (H E = 0.67) exhibits slightly lower diversity and greater population differentiation than O. rufipogon (H E = 0.70). Bayesian inference identified four, and at a finer structural level eight, genetically distinct population groups that correspond to geographic populations within the three taxa. Oryza meridionalis and the Nepalese O. nivara seemed diverged from all the population groups of the series, whereas the Australasian O. rufipogon appeared distinct from the rest of the species.

Rubisco activity is associated with photosynthetic thermotolerance in a wild rice (Oryza meridionalis)

Oryza meridionalis is a wild species of rice, endemic to tropical Australia. It shares a significant genome homology with the common domesticated rice Oryza sativa. Exploiting the fact that the two species are highly related but O. meridionalis has superior heat tolerance, experiments were undertaken to identify the impact of temperature on key events in photosynthesis. At an ambient CO(2) partial pressure of 38 Pa and irradiance of 1500 µmol quanta m(-2) s(-1), the temperature optimum of photosynthesis was 33.7 ± 0.8°C for O. meridionalis, significantly higher than the 30.6 ± 0.7°C temperature optimum of O. sativa. To understand the basis for this difference, we measured gas exchange and rubisco activation state between 20 and 42°C and modeled the response to determine the rate-limiting steps of photosynthesis. The temperature response of light respiration (R(light)) and the CO(2) compensation point in the absence of respiration (Γ(*)) were determined and found to be similar for the two species. C3 photosynthesis modeling showed that despite the difference in susceptibility to high temperature, both species had a similar temperature-dependent limitation to photosynthesis. Both rice species were limited by ribulose-1,5-bisphosphate (RuBP) regeneration at temperatures of 25 and 30°C but became RuBP carboxylation limited at 35 and 40°C. The activation state of rubisco in O. meridionalis was more stable at higher temperatures, explaining its greater heat tolerance compared with O. sativa.

Temperature response of mesophyll conductance in cultivated and wild Oryza species with contrasting mesophyll cell wall thickness

A critical component of photosynthetic capacity is the conductance of CO(2) from intercellular airspaces to the sites of CO(2) fixation in the stroma of chloroplasts, termed mesophyll conductance (g(m)). Leaf anatomy has been identified as an important determinant of g(m). There are few studies of the temperature response of g(m) and none has examined the implications of leaf anatomy. Hence, we compared a cultivar of Oryza sativa with two wild Oryza relatives endemic to the hot northern savannah of Australia, namely Oryza meridionalis and Oryza australiensis. All three species had similar leaf anatomical properties, except that the wild relatives had significantly thicker mesophyll cell walls than O. sativa. Thicker mesophyll cell walls in the wild rice species are likely to have contributed to the reduction in g(m) , which was associated with a greater drawdown of CO(2) into chloroplasts (C(i) -C(c) ) compared with O. sativa. Mesophyll conductance increased at higher temperatures, whereas the rate of CO(2) assimilation was relatively stable between 20 and 40 °C. Consequently, C(i) -C(c) decreased for all three species as temperature increased.

A next-generation sequencing method for overcoming the multiple gene copy problem in polyploid phylogenetics, applied to Poa grasses

BACKGROUND

Polyploidy is important from a phylogenetic perspective because of its immense past impact on evolution and its potential future impact on diversification, survival and adaptation, especially in plants. Molecular population genetics studies of polyploid organisms have been difficult because of problems in sequencing multiple-copy nuclear genes using Sanger sequencing. This paper describes a method for sequencing a barcoded mixture of targeted gene regions using next-generation sequencing methods to overcome these problems.

RESULTS

Using 64 3-bp barcodes, we successfully sequenced three chloroplast and two nuclear gene regions (each of which contained two gene copies with up to two alleles per individual) in a total of 60 individuals across 11 species of Australian Poa grasses. This method had high replicability, a low sequencing error rate (after appropriate quality control) and a low rate of missing data. Eighty-eight percent of the 320 gene/individual combinations produced sequence reads, and >80% of individuals produced sufficient reads to detect all four possible nuclear alleles of the homeologous nuclear loci with 95% probability.We applied this method to a group of sympatric Australian alpine Poa species, which we discovered to share an allopolyploid ancestor with a group of American Poa species. All markers revealed extensive allele sharing among the Australian species and so we recommend that the current taxonomy be re-examined. We also detected hypermutation in the trnH-psbA marker, suggesting it should not be used as a land plant barcode region. Some markers indicated differentiation between Tasmanian and mainland samples. Significant positive spatial genetic structure was detected at <100 km with chloroplast but not nuclear markers, which may be a result of restricted seed flow and long-distance pollen flow in this wind-pollinated group.

CONCLUSIONS

Our results demonstrate that 454 sequencing of barcoded amplicon mixtures can be used to reliably sample all alleles of homeologous loci in polyploid species and successfully investigate phylogenetic relationships among species, as well as to investigate phylogeographic hypotheses. This next-generation sequencing method is more affordable than and at least as reliable as bacterial cloning. It could be applied to any experiment involving sequencing of amplicon mixtures.

Orthologous comparisons of the Hd1 region across genera reveal Hd1 gene lability within diploid Oryza species and disruptions to microsynteny in Sorghum

Heading date is one of the most important quantitative traits responsible for the domestication of rice. We compared a 155-kb reference segment of the Oryza sativa ssp. japonica cv. Nipponbare genome surrounding Hd1, a major heading date gene in rice, with orthologous regions from nine diploid Oryza species that diverged over a relatively short time frame (∼16 My) to study sequence evolution around a domestication locus. The orthologous Hd1 region from Sorghum bicolor was included to compare and contrast the evolution in a more distant relative of rice. Consistent with other observations at the adh1/adh2, monoculm1, and sh2/a1 loci in grass species, we found high gene colinearity in the Hd1 region amidst size differences that were lineage specific and long terminal repeat retrotransposon driven. Unexpectedly, the Hd1 gene was deleted in O. glaberrima, whereas the O. rufipogon and O. punctata copies had degenerative mutations, suggesting that other heading date loci might compensate for the loss or nonfunctionality of Hd1 in these species. Compared with the japonica Hd1 region, the orthologous region in sorghum exhibited micro-rearrangements including gene translocations, seven additional genes, and a gene triplication and truncation event predating the divergence from Oryza.

Paleogenomic analysis of the short arm of chromosome 3 reveals the history of the African and Asian progenitors of cultivated rices

Rice is one of the most important crops, feeding more than half of the world population. There are two cultivated species, the African rice Oryza glaberrima and the Asian rice O. sativa. Although the African species is gradually replaced by O. sativa in most of African rice agrosystems, this species represents an important reservoir of genes of agronomical interest. Their exploitation for the development of modern African rice varieties requires a good understanding of the genetic relationships between the two cultivated species. We took advantage of the recent availability of the sequence of the chromosome 3 short arm of O. glaberrima to estimate the date of radiation between O. glaberrima and O. sativa lineages, using all the long terminal repeat (LTR)-retrotransposons as paleogenomic markers. We first demonstrated that in two distinct lineages, LTR-retrotransposons mutate at the same rate. Based on LTR-retrotransposons shared by both species in orthologous position, we then estimated that O. glaberrima and O. sativa progenitors diverged 1.2 Ma. This constitutes one of the first studies using such a large sample of transposable elements to reconstruct the phylogeny of species. Given the number of genome sequencing projects, there is no doubt that such approach will allow to resolve phylogenetic incongruities. The application of this method to other plant genomes will also facilitate further understanding of evolution of LTR-retrotransposons and eventually of the whole genome in divergent plant lineages.

Relationships among genetic makeup, active ingredient content, and place of origin of the medicinal plant Gastrodia tuber

Gastrodia tuber and its component gastrodin have many pharmacological effects. The chemical fingerprints and gastrodin contents of eight Gastrodia populations were determined, and the genomic DNA polymorphism of the populations was investigated. Genetic distance coefficients among the populations were calculated using the DNA polymorphism data. A dendrogram of the genetic similarities between the populations was constructed using the genetic distance coefficients. The results indicated that the genomic DNA of Gastrodia tubers was highly polymorphic; the eight populations clustered into three major groups, and the gastrodin content varied greatly among these groups. There were obvious correlations among genetic makeup, gastrodin content, and place of origin. The ecological environments in Guizhou and Shanxi may be conducive to evolution and to gastrodin biosynthesis, and more suitable for cultivation of Gastrodia tubers. These findings may provide a scientific basis for overall genetic resource management and for the selection of locations for cultivating Gastrodia tubers.

Retrotransposon insertion polymorphisms in six rice genes and their evolutionary history

Retrotransposons are abundant in higher plant genomes. Although retrotransposons associated with plant genes have been identified, little is known about their evolutionary conservation at the level of species and subspecies. In the present study, we investigated the phylogenetic distribution of long terminal repeat (LTR) retrotransposon, long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) insertions in six genes in 95 cultivated and wild rice genotypes. These six genes are likely to be functional based on nonsynonymous (Ka) to synonymous (Ks) substitution ratios which were found to be significantly <1. Different conservation patterns of these retrotransposons in genes were observed in cultivated and wild rice species. Four out of seven retrotransposon insertions appear to predate the ancestral Oryza AA genome. Two of these insertions in genes 4 and 5 occurred early in the evolutionary history of Oryza. Two retrotransposon insertions in gene 1 arose after the divergence of Asian cultivated rice from its wild ancestor. Furthermore, the retrotransposon insertion in gene 3 appears to have occurred in the ancestral lineage leading to temperate japonicas. Conservation of retrotransposon insertions in genes in specific groups, species, and lineages might be related to their specific function.