Dynamic evolution of oryza genomes is revealed by comparative genomic analysis of a genus-wide vertical data set

Fumarate reductase drives methane emissions in the genus Oryza through differential regulation of the rhizospheric ecosystem

The emergence of waterlogged Oryza species ∼15Mya (million years ago) supplied an anoxic warm bed for methane-producing microorganisms, and methane emissions have hence accompanied the entire evolutionary history of the genus Oryza. However, to date no study has addressed how methane emission has been altered during Oryza evolution. In this paper we used a diverse collection of wild and cultivated Oryza species to study the relation between Oryza evolution and methane emissions. Phylogenetic analyses and methane detection identified a co-evolutionary pattern between Oryza and methane emissions, mediated by the diversity of the rhizospheric ecosystems arising from different oxygen levels. Fumarate was identified as an oxygen substitute used to retain the electron transport/energy production in the anoxic rice root, and the contribution of fumarate reductase to Oryza evolution and methane emissions has also been assessed. We confirmed the between-species patterns using genetic dissection of the traits in a cross between a low and high methane-emitting species. Our findings provide novel insights on the evolutionary processes of rice paddy methane emissions: the evolution of wild rice produces different Oryza species with divergent rhizospheric ecosystem attributing to the different oxygen levels and fumarate reductase activities, methane emissions are comprehensively assessed by the rhizospheric environment of diversity Oryza species and result in a co-evolution pattern.

Fine-mapping and evolutionary history of R-BPMV, a dominant resistance gene to Bean pod mottle virus in Phaseolus vulgaris L

KEY MESSAGE

R-BPMV is located within a recently expanded TNL cluster in the Phaseolus genus with suppressed recombination and known for resistance to multiple pathogens including potyviruses controlled by the I gene. Bean pod mottle virus (BPMV) is a comovirus that infects common bean and legumes in general. BPMV is distributed throughout the world and is a major threat on soybean, a closely related species of common bean. In common bean, BAT93 was reported to carry the R-BPMV resistance gene conferring resistance to BPMV and linked with the I resistance gene. To fine map R-BPMV, 182 recombinant inbred lines (RILs) derived from the cross BAT93 × JaloEEP558 were genotyped with polymerase chain reaction (PCR)-based markers developed using genome assemblies from G19833 and BAT93, as well as BAT93 BAC clone sequences. Analysis of RILs carrying key recombination events positioned R-BPMV to a target region containing at least 16 TIR-NB-LRR (TNL) sequences in BAT93. Because the I cluster presents a suppression of recombination and a large number of repeated sequences, none of the 16 TNLs could be excluded as R-BPMV candidate gene. The evolutionary history of the TNLs for the I cluster were reconstructed using microsynteny and phylogenetic analyses within the legume family. A single I TNL was present in Medicago truncatula and lost in soybean, mirroring the absence of complete BPMV resistance in soybean. Amplification of TNLs in the I cluster predates the divergence of the Phaseolus species, in agreement with the emergence of R-BPMV before the separation of the common bean wild centers of diversity. This analysis provides PCR-based markers useful in marker-assisted selection (MAS) and laid the foundation for cloning of R-BPMV resistance gene in order to transfer the resistance into soybean.

Phylogenomic Analysis of micro-RNA Involved in Juvenile to Flowering-Stage Transition in Photophilic Rice and Its Sister Species

Vegetative to reproductive phase transition in phototropic plants is an important developmental process and is sequentially mediated by the expression of micro-RNA MIR172. To obtain insight into the evolution, adaptation, and function of MIR172 in photophilic rice and its wild relatives, we analyzed the genescape of a 100 kb segment harboring MIR172 homologs from 11 genomes. The expression analysis of MIR172 revealed its incremental accumulation from the 2-leaf to 10-leaf stage, with maximum expression coinciding with the flag-leaf stage in rice. Nonetheless, the microsynteny analysis of MIR172s revealed collinearity within the genus Oryza, but a loss of synteny was observed in (i) MIR172A in O. barthii (AA) and O. glaberima (AA); (ii) MIR172B in O. brachyantha (FF); and (iii) MIR172C in O. punctata (BB). Phylogenetic analysis of precursor sequences/region of MIR172 revealed a distinct tri-modal clade of evolution. The genomic information generated in this investigation through comparative analysis of MIRNA, suggests mature MIR172s to have evolved in a disruptive and conservative mode amongst all Oryza species with a common origin of descent. Further, the phylogenomic delineation provided an insight into the adaptation and molecular evolution of MIR172 to changing environmental conditions (biotic and abiotic) of phototropic rice through natural selection and the opportunity to harness untapped genomic regions from rice wild relatives (RWR).

The International Oryza Map Alignment Project (IOMAP): the Americas-past achievements and future directions

The wild relatives of rice hold unexplored genetic diversity that can be employed to feed an estimated population of 10 billion by 2050. The Oryza Map Alignment Project (OMAP) initiated in 2003 has provided comprehensive genomic resources for comparative, evolutionary, and functional characterization of the wild relatives of rice, facilitating the cloning of >600 rice genes, including those for grain width (GW5) and submergence tolerance (SUB1A). Following in the footsteps of the original project, the goal of 'IOMAP: the Americas' is to investigate the present and historic genetic diversity of wild Oryza species endemic to the Americas through the sequencing of herbaria and in situ specimens. The generation of a large diversity panel describing past and current genetic status and potential erosion of genetic variation in the populations will provide useful knowledge for the conservation of the biodiversity in these species. The wild relatives of rice in the Americas present a wide range of resistance traits useful for crop improvement and neodomestication approaches. In the race against time for a sustainable food future, the neodomestication of the first cereal species recently accomplished in O. alta opens the door to the potential neodomestication of the other wild Oryza species in Americas.

Analysis of Homologous Regions of Small RNAs MIR397 and MIR408 Reveals the Conservation of Microsynteny among Rice Crop-Wild Relatives

MIRNAs are small non-coding RNAs that play important roles in a wide range of biological processes in plant growth and development. MIR397 (involved in drought, low temperature, and nitrogen and copper (Cu) starvation) and MIR408 (differentially expressed in response to environmental stresses such as copper, light, mechanical stress, dehydration, cold, reactive oxygen species, and drought) belong to conserved MIRNA families that either negatively or positively regulate their target genes. In the present study, we identified the homologs of MIR397 and MIR408 in Oryza sativa and its six wild progenitors, three non-Oryza species, and one dicot species. We analyzed the 100 kb segments harboring MIRNA homologs from 11 genomes to obtain a comprehensive view of their community evolution around these loci in the farthest (distant) relatives of rice. Our study showed that mature MIR397 and MIR408 were highly conserved among all Oryza species. Comparative genomics analyses also revealed that the microsynteny of the 100 kb region surrounding MIRNAs was only conserved in Oryza spp.; disrupted in Sorghum, maize, and wheat; and completely lost in Arabidopsis. There were deletions, rearrangements, and translocations within the 100 kb segments in Oryza spp., but the overall microsynteny of the region was maintained. The phylogenetic analyses of the precursor regions of all MIRNAs under study revealed a bimodal clade of common origin. This comparative analysis of miRNA involved in abiotic stress tolerance in plants provides a powerful tool for future Oryza research. Crop wild relatives (CWRs) offer multiple traits with potential to decrease the amount of yield loss owing to biotic and abiotic stresses. Using a comparative genomics approach, the exploration of CWRs as a source of tolerance to these stresses by understanding their evolution can be further used to leverage their yield potential.

The ties of brotherhood between japonica and indica rice for regional adaptation

Selection of beneficial genomic variants was crucial for regional adaptation of crops during domestication, but the underlying genomic basis remains largely unexplored. Here we report a genome-wide selective-sweep analysis of 655 japonica and 1,205 indica accessions selected from 2,673 landraces through principal component analysis to identify 5,636 non-synonymous single nucleotide polymorphisms (SNPs) fixed in at least one subspecies. We classified these SNPs into three groups, jiS (japonica- and indica-selected), jS (japonica-selected only), and iS (indica-selected only), and documented evidence for selection acting on these groups, their relation to yield-related traits, such as heading date, and their practical value in cropping area prediction. We also demonstrated the role of a jiS-SNP-containing gene in temperature adaptability. Our study informs genes underpinning adaptation that may shape Green Super Rice and proposes a time-saving, cost-reducing selection strategy of genomic breeding, sweep-SNP-guided selection, for developing regionally-adapted heterosis.

Evolutionary alterations in gene expression and enzymatic activities of gibberellin 3-oxidase 1 in Oryza

Proper anther and pollen development are important for plant reproduction. The plant hormone gibberellin is important for anther development in rice, but its gametophytic functions remain largely unknown. Here, we report the functional and evolutionary analyses of rice gibberellin 3-oxidase 1 (OsGA3ox1), a gibberellin synthetic enzyme specifically expressed in the late developmental stages of anthers. Enzymatic and X-ray crystallography analyses reveal that OsGA3ox1 has a higher GA₇ synthesis ratio than OsGA3ox2. In addition, we generate an osga3ox1 knockout mutant by genome editing and demonstrate the bioactive gibberellic acid synthesis by the OsGA3ox1 action during starch accumulation in pollen via invertase regulation. Furthermore, we analyze the evolution of Oryza GA3ox1s and reveal that their enzyme activity and gene expression have evolved in a way that is characteristic of the Oryza genus and contribute to their male reproduction ability.

The authors solve the crystal structure of OsGA3ox2 and predict that of OsGA3ox1. These enzymes catalyze the final step in the biosynthesis of gibberellin, one of the plant hormones. Evolutionary analysis combined with the new structure reveal important aspects of the OsGA3ox1’s function in plant male reproduction.

StAR-Related Lipid Transfer (START) Domains Across the Rice Pangenome Reveal How Ontogeny Recapitulated Selection Pressures During Rice Domestication

The StAR-related lipid transfer (START) domain containing proteins or START proteins, encoded by a plant amplified family of evolutionary conserved genes, play important roles in lipid binding, transport, signaling, and modulation of transcriptional activity in the plant kingdom, but there is limited information on their evolution, duplication, and associated sub- or neo-functionalization. Here we perform a comprehensive investigation of this family across the rice pangenome, using 10 wild and cultivated varieties. Conservation of START domains across all 10 rice genomes suggests low dispensability and critical functional roles for this family, further supported by chromosomal mapping, duplication and domain structure patterns. Analysis of synteny highlights a preponderance of segmental and dispersed duplication among STARTs, while transcriptomic investigation of the main cultivated variety Oryza sativa var. japonica reveals sub-functionalization amongst genes family members in terms of preferential expression across various developmental stages and anatomical parts, such as flowering. Ka/Ks ratios confirmed strong negative/purifying selection on START family evolution, implying that ontogeny recapitulated selection pressures during rice domestication. Our findings provide evidence for high conservation of START genes across rice varieties in numbers, as well as in their stringent regulation of Ka/Ks ratio, and showed strong functional dependency of plants on START proteins for their growth and reproductive development. We believe that our findings advance the limited knowledge about plant START domain diversity and evolution, and pave the way for more detailed assessment of individual structural classes of START proteins among plants and their domain specific substrate preferences, to complement existing studies in animals and yeast.

Dynamic differential evolution schemes of WRKY transcription factors in domesticated and wild rice

WRKY transcription factors play key roles in stress responses, growth, and development. We previously reported on the evolution of WRKYs from unicellular green algae to land plants. To address recent evolution events, we studied three domesticated and eight wild species in the genus Oryza, an ideal model due to its long history of domestication, economic importance, and central role as a model system. We have identified prevalence of Group III WRKYs despite differences in breeding of cultivated and wild species. Same groups of WRKY genes tend to cluster together, suggesting recent, multiple duplication events. Duplications followed by divergence may result in neofunctionalizations of co-expressed WRKY genes that finely tune the regulation of target genes in a same metabolic or response pathway. WRKY genes have undergone recent rearrangements to form novel genes. Group Ib WRKYs, unique to AA genome type Oryza species, are derived from Group III genes dated back to 6.76 million years ago. Gene tree reconciliation analysis with the species tree revealed details of duplication and loss events in the 11 genomes. Selection analysis on single copy orthologs reveals the highly conserved nature of the WRKY domain and clusters of fast evolving sites under strong positive selection pressure. Also, the numbers of single copy orthologs under positive or negative selection almost evenly split. Our results provide valuable insights into the preservation and diversification of an important gene family under strong selective pressure for biotechnological improvements of the world's most valued food crop.

Xa7, a Small Orphan Gene Harboring Promoter Trap for AvrXa7, Leads to the Durable Resistance to Xanthomonas oryzae Pv. oryzae

BACKGROUND

The rice (Oryza sativa) gene Xa7 has been hypothesized to be a typical executor resistance gene against Xanthomonas oryzae pv. oryzae (Xoo), and has conferred durable resistance in the field for decades. Its identity and the molecular mechanisms underlying this resistance remain elusive.

RESULTS

Here, we filled in gaps of genome in Xa7 mapping locus via BAC library construction, revealing the presence of a 100-kb non-collinear sequence in the line IRBB7 compared with Nipponbare reference genomes. Complementary transformation with sequentially overlapping subclones of the BACs demonstrated that Xa7 is an orphan gene, encoding a small novel protein distinct from any other resistance proteins reported. A 27-bp effector binding element (EBE) in the Xa7 promoter is essential for AvrXa7-inducing expression model. XA7 is anchored in the endoplasmic reticulum membrane and triggers programmed cell death in rice and tobacco (Nicotiana benthamiana). The Xa7 gene is absent in most cultivars, landraces, and wild rice accessions, but highly homologs of XA7 were identified in Leersia perrieri, the nearest outgroup of the genus Oryza.

CONCLUSIONS

Xa7 acts as a trap to perceive AvrXa7 via EBEAvrXa7 in its promoter, leading to the initiation of resistant reaction. Since EBEAvrXa7 is ubiquitous in promoter of rice susceptible gene SWEET14, the elevated expression of which is conducive to the proliferation of Xoo, that lends a great benefit for the Xoo strains retaining AvrXa7. As a result, varieties harboring Xa7 would show more durable resistance in the field. Xa7 alleles analysis suggests that the discovery of new resistance genes could be extended beyond wild rice, to include wild grasses such as Leersia species.

Starch Synthesis-Related Genes (SSRG) Evolution in the Genus Oryza

Cooking quality is an important attribute in Common/Asian rice (Oryzasativa L.) varieties, being highly dependent on grain starch composition. This composition is known to be highly dependent on a cultivar’s genetics, but the way in which their genes express different phenotypes is not well understood. Further analysis of variation of grain quality genes using new information obtained from the wild relatives of rice should provide important insights into the evolution and potential use of these genetic resources. All analyses were conducted using bioinformatics approaches. The analysis of the protein sequences of grain quality genes across the Oryza suggest that the deletion/mutation of amino acids in active sites result in variations that can negatively affect specific steps of starch biosynthesis in the endosperm. On the other hand, the complete deletion of some genes in the wild species may not affect the amylose content. Here we present new insights for Starch Synthesis-Related Genes (SSRGs) evolution from starch-specific rice phenotypes.

Genome evolution in Oryza allopolyploids of various ages: Insights into the process of diploidization

The prevalence and recurrence of whole-genome duplication in plants and its major role in evolution have been well recognized. Despite great efforts, many aspects of genome evolution, particularly the temporal progression of genomic responses to allopolyploidy and the underlying mechanisms, remain poorly understood. The rice genus Oryza consists of both recently formed and older allopolyploid species, representing an attractive system for studying the genome evolution after allopolyploidy. In this study, through screening BAC libraries and sequencing and annotating the targeted BAC clones, we generated orthologous genomic sequences surrounding the DEP1 locus, a major grain yield QTL in cultivated rice, from four Oryza polyploids of various ages and their likely diploid genome donors or close relatives. Based on sequenced DEP1 region and published data from three other genomic regions, we investigated the temporal evolutionary dynamics of four polyploid genomes at both genetic and expression levels. In the recently formed BBCC polyploid, Oryza minuta, genome dominance was not observed and its short-term responses to allopolyploidy are mainly manifested as a high proportion of homoeologous gene pairs showing unequal expression. This could partly be explained by parental legacy, rewiring of divergent regulatory networks and epigenetic modulation. Moreover, we detected an ongoing diploidization process in this genus, and suggest that the expression divergence driven by changes of selective constraint probably plays a big role in the long-term diploidization. These findings add novel insights into our understanding of genome evolution after allopolyploidy, and could facilitate crop improvements through hybridization and polyploidization.

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.

Wild and Cultivated Species of Rice Have Distinctive Proteomic Responses to Drought

Drought often compromises yield in non-irrigated crops such as rainfed rice, imperiling the communities that depend upon it as a primary food source. In this study, two cultivated species (Oryza sativa cv. Nipponbare and Oryza glaberrima cv. CG14) and an endemic, perennial Australian wild species (Oryza australiensis) were grown in soil at 40% field capacity for 7 d (drought). The hypothesis was that the natural tolerance of O. australiensis to erratic water supply would be reflected in a unique proteomic profile. Leaves from droughted plants and well-watered controls were harvested for label-free quantitative shotgun proteomics. Physiological and gene ontology analysis confirmed that O. australiensis responded uniquely to drought, with superior leaf water status and enhanced levels of photosynthetic proteins. Distinctive patterns of protein accumulation in drought were observed across the O. australiensis proteome. Photosynthetic and stress-response proteins were more abundant in drought-affected O. glaberrima than O. sativa, and were further enriched in O. australiensis. In contrast, the level of accumulation of photosynthetic proteins decreased when O. sativa underwent drought, while a narrower range of stress-responsive proteins showed increased levels of accumulation. Distinctive proteomic profiles and the accumulated levels of individual proteins with specific functions in response to drought in O. australiensis indicate the importance of this species as a source of stress tolerance genes.

Back to the Wild: On a Quest for Donors Toward Salinity Tolerant Rice

Salinity stress affects global food producing areas by limiting both crop growth and yield. Attempts to develop salinity-tolerant rice varieties have had limited success due to the complexity of the salinity tolerance trait, high variation in the stress response and a lack of available donors for candidate genes for cultivated rice. As a result, finding suitable donors of genes and traits for salinity tolerance has become a major bottleneck in breeding for salinity tolerant crops. Twenty-two wild Oryza relatives have been recognized as important genetic resources for quantitatively inherited traits such as resistance and/or tolerance to abiotic and biotic stresses. In this review, we discuss the challenges and opportunities of such an approach by critically analyzing evolutionary, ecological, genetic, and physiological aspects of Oryza species. We argue that the strategy of rice breeding for better Na+ exclusion employed for the last few decades has reached a plateau and cannot deliver any further improvement in salinity tolerance in this species. This calls for a paradigm shift in rice breeding and more efforts toward targeting mechanisms of the tissue tolerance and a better utilization of the potential of wild rice where such traits are already present. We summarize the differences in salinity stress adaptation amongst cultivated and wild Oryza relatives and identify several key traits that should be targeted in future breeding programs. This includes: (1) efficient sequestration of Na+ in mesophyll cell vacuoles, with a strong emphasis on control of tonoplast leak channels; (2) more efficient control of xylem ion loading; (3) efficient cytosolic K+ retention in both root and leaf mesophyll cells; and (4) incorporating Na+ sequestration in trichrome. We conclude that while amongst all wild relatives, O. rufipogon is arguably a best source of germplasm at the moment, genes and traits from the wild relatives, O. coarctata, O. latifolia, and O. alta, should be targeted in future genetic programs to develop salt tolerant cultivated rice.

The unique epigenetic features of Pack-MULEs and their impact on chromosomal base composition and expression spectrum

Acquisition and rearrangement of host genes by transposable elements (TEs) is an important mechanism to increase gene diversity as exemplified by the ∼3000 Pack-Mutator-like TEs in the rice genome which have acquired gene sequences (Pack-MULEs), yet remain enigmatic. To identify signatures of functioning Pack-MULEs and Pack-MULE evolution, we generated transcriptome, translatome, and epigenome datasets and compared Pack-MULEs to genes and other TE families. Approximately 40% of Pack-MULEs were transcribed with 9% having translation evidence, clearly distinguishing them from other TEs. Pack-MULEs exhibited a unique expression profile associated with specificity in reproductive tissues that may be associated with seed traits. Expressed Pack-MULEs resemble regular protein-coding genes as exhibited by a low level of DNA methylation, association with active histone marks and DNase I hypersensitive sites, and an absence of repressive histone marks, suggesting that a substantial fraction of Pack-MULEs are potentially functional in vivo. Interestingly, the expression capacity of Pack-MULEs is independent of the local genomic environment, and the insertion and expression of Pack-MULEs may have altered the local chromosomal expression pattern as well as counteracted the impact of recombination on chromosomal base composition, which has profound consequences on the evolution of chromosome structure.

Introgression of the chromosomal region with the Pi-cd locus from Oryza meridionalis into O. sativa L. during rice domestication

The genotype of the Pi-cd locus found in blast-resistant rice variety Kitakurin, which is a cultivated rice from Japan belonging to Oryza sativa japonica, is identical to that of its wild relative O. meridionalis. Crop domestication from wild relatives to cultivated species has encompassed significant phenotypic changes. However, little is known about the genetic changes involved in domestication. Here, we surveyed the origin of the Pi-cd locus across Oryza species with AA genomes by comparison with the genome sequences of Hoshinoyume (HS), which does not carry the Pi-cd blast resistance gene, and Kitakurin (KK), which carries the Pi-cd blast resistance gene. We found that variety-specific transposons were enriched at the Pi-cd locus. The genotype of the Pi-cd locus characterized by transposons in HS and KK was specific to each Oryza species with the AA genome. The Kitaake (KT) genotype at the Pi-cd locus found in KK was identical only to that of O. meridionalis and distributed only in subgroups of japonica in the World Rice Collection and tropical japonica in the Japanese Rice Collection, whereas it was not present in O. rufipogon accessions. The distinct distributions of genotypes of the Pi-cd locus clearly demonstrated that the Pi-cd locus was introgressed from O. meridionalis into O. sativa, specific to tropical japonica.

The Genomics of Oryza Species Provides Insights into Rice Domestication and Heterosis

Here, we review recent progress in genetic and genomic studies of the diversity of Oryza species. In recent years, unlocking the genetic diversity of Oryza species has provided insights into the genomics of rice domestication, heterosis, and complex traits. Genome sequencing and analysis of numerous wild rice (Oryza rufipogon) and Asian cultivated rice (Oryza sativa) accessions have enabled the identification of genome-wide signatures of rice domestication and the unlocking of the origin of Asian cultivated rice. Moreover, similar studies on genome variations of African rice (Oryza glaberrima) cultivars and their closely related wild progenitor Oryza barthii accessions have provided strong evidence to support a theory of independent domestication in African rice. Integrated genomic approaches have efficiently investigated many heterotic loci in hybrid rice underlying yield heterosis advantages and revealed the genomic architecture of rice heterosis. We conclude that in-depth unlocking of genetic variations among Oryza species will further enhance rice breeding.

Lateral transfers of large DNA fragments spread functional genes among grasses

A fundamental tenet of multicellular eukaryotic evolution is that vertical inheritance is paramount, with natural selection acting on genetic variants transferred from parents to offspring. This lineal process means that an organism's adaptive potential can be restricted by its evolutionary history, the amount of standing genetic variation, and its mutation rate. Lateral gene transfer (LGT) theoretically provides a mechanism to bypass many of these limitations, but the evolutionary importance and frequency of this process in multicellular eukaryotes, such as plants, remains debated. We address this issue by assembling a chromosome-level genome for the grass Alloteropsis semialata, a species surmised to exhibit two LGTs, and screen it for other grass-to-grass LGTs using genomic data from 146 other grass species. Through stringent phylogenomic analyses, we discovered 57 additional LGTs in the A. semialata nuclear genome, involving at least nine different donor species. The LGTs are clustered in 23 laterally acquired genomic fragments that are up to 170 kb long and have accumulated during the diversification of Alloteropsis. The majority of the 59 LGTs in A. semialata are expressed, and we show that they have added functions to the recipient genome. Functional LGTs were further detected in the genomes of five other grass species, demonstrating that this process is likely widespread in this globally important group of plants. LGT therefore appears to represent a potent evolutionary force capable of spreading functional genes among distantly related grass species.

Potential of Oryza officinalis to augment the cold tolerance genetic mechanisms of Oryza sativa by network complementation

Oryza officinalis is an accessible alien donor for genetic improvement of rice. Comparison across a representative panel of Oryza species showed that the wild O. officinalis and cultivated O. sativa ssp. japonica have similar cold tolerance potentials. The possibility that either distinct or similar genetic mechanisms are involved in the low temperature responses of each species was addressed by comparing their transcriptional networks. General similarities were supported by shared transcriptomic signatures indicative of equivalent metabolic, hormonal, and defense status. However, O. officinalis has maintained an elaborate cold-responsive brassinosteroid-regulated BES1-network that appeared to have been fragmented in O. sativa. BES1-network is potentially important for integrating growth-related responses with physiological adjustments and defenses through the protection of photosynthetic machinery and maintenance of stomatal aperture, oxidative defenses, and osmotic adjustment. Equivalent physiological processes are functional in O. sativa but their genetic mechanisms are under the direct control of ABA-dependent, DREB-dependent and/or oxidative-mediated networks uncoupled to BES1. While O. officinalis and O. sativa represent long periods of speciation and domestication, their comparable cold tolerance potentials involve equivalent physiological processes but distinct genetic networks. BES1-network represents a novel attribute of O. officinalis with potential applications in diversifying or complementing other mechanisms in the cultivated germplasm.

Transcriptome Analysis and Functional Identification of Xa13 and Pi-ta Orthologs in Oryza granulata

Nees & Arn. ex Watt, a perennial wild rice species with a GG genome, preserves many important genes for cultivated rice ( L.) improvement. At present, however, no genetic resource is available for studying . Here, we report 91,562 high-quality transcripts of assembled de novo. Moreover, comparative transcriptome analysis revealed that 1311 single-copy orthologous pairs shared by and (Zoll. & Moritzi) Baill. that may have undergone adaptive evolution. We performed an analysis of the genes potentially involved in plant-pathogen interactions to explore the molecular basis of disease resistance, and isolated the full-length cDNAs of () and () orthologs from . The overexpression of in Nipponbare and functional characterization showed enhanced the resistance of transgenic Nipponbare to rice blast resulting from the presence of the gene. , an alternatively spliced transcript of the blast resistance gene in encodes a 1024-amino acid polypeptide with a C-terminal thioredoxin domain. This study provides an important resource for functional and evolutionary studies of the genus .

Monosomic alien addition lines (MAALs) of Oryza rhizomatis in Oryza sativa: production, cytology, alien trait introgression, molecular analysis and breeding application

Key message Development of MAALs and disomic introgression lines derived from the cross between O. sativa and O. rhizomatis to exploit and utilize the valuable traits for rice improvement. The CC genome wild species, Oryza rhizomatis, possesses valuable traits for rice improvement. Unlike other CC genome wild rice, O. rhizomatis is less studied and none of the research has focused on the utilization of this resource in rice breeding. The transfer of novel genes governing the valuable traits from O. rhizomatis is difficult due to high genome incompatibility with O. sativa. Here we report the development of backcross progenies and complete sets of monosomic alien addition lines (MAALs) for the first time from O. rhizomatis in O. sativa line IR31917-45-3-2. Autotetraploid IR31917-45-3-2 (4x = AAAA) was used to generate allotriploid F₁, and the F₁ plant was backcrossed to IR31917-45-3-2 (2x). Forty-seven BC₁F₁ and 73 BC₂F₁ plants were produced with chromosome numbers ranging from 24 to 33 (2x + 9) and 24 to 27 (2x + 3), respectively. A complete set of MAALs were identified by morphological, cytological and marker-based analysis. A total of 116 CC genome-specific InDel markers across the 12 chromosome of rice were used to detect O. rhizomatis chromosome segments in F₁, BC₁F₁, BC₂F₂, MAALs and disomic introgression lines (DILs). Expressions of major phenotypic traits inherited from O. rhizomatis were observed in MAAL-derived DILs. Small chromosomal segments of O. rhizomatis for chromosomes 1, 2, 4, 5, 6, 7, 10 and 12 were detected in DILs, and some of the introgression lines showed insect resistance against brown planthopper and green leafhopper. These newly developed MAALs and DILs will be useful for gene mining and more precise faster transfer of favorable genes to improve rice cultivars.

Antifungal and Biocontrol Evaluation of Four Lysobacter Strains Against Clubroot Disease

The effect of crude extract (Ce), seed coating agent (SCA) and whole bacterial broth culture (WBC) of Lysobacter strains was evaluated against the causal agent of clubroot formation in Cruciferous vegetables. The ability of four Lysobacter strains (L. antibioticus 6-B-1, L. antibioticus 6-T-4, L. antibioticus 13-B-1 and L. capsici ZST1-2) inhibited Plasmodiophora brassicae of resting spores and disease. Application of WBC of four Lysobacter strains inhibited clubroot disease, indicating that the disease suppression was due to antifungal compounds produced by the biocontrol bacterium in the culture. Development of clubroot on Chinese cabbage was inhibited when the WBC and SCA were applied before P. brassicae inoculation. Crude extract (Ce) of culture filtrate was effective in arresting the germination of resting spores of P. brassicae on slides. However, Lysobacter strains differed in their biocontrol effects, the strain L. capsci ZST1-2 recorded a high level of disease limiting effect.

Common Bean Subtelomeres Are Hot Spots of Recombination and Favor Resistance Gene Evolution

Subtelomeres of most eukaryotes contain fast-evolving genes usually involved in adaptive processes. In common bean (Phaseolus vulgaris), the Co-2 anthracnose resistance (R) locus corresponds to a cluster of nucleotide-binding-site leucine-rich-repeat (NL) encoding sequences, the prevalent class of plant R genes. To study the recent evolution of this R gene cluster, we used a combination of sequence, genetic and cytogenetic comparative analyses between common bean genotypes from two distinct gene pools (Andean and Mesoamerican) that diverged 0.165 million years ago. Co-2 is a large subtelomeric cluster on chromosome 11 comprising from 32 (Mesoamerican) to 52 (Andean) NL sequences embedded within khipu satellite repeats. Since the recent split between Andean and Mesoamerican gene pools, the Co-2 cluster has experienced numerous gene-pool specific NL losses, leading to distinct NL repertoires. The high proportion of solo-LTR retrotransposons indicates that the Co-2 cluster is located in a hot spot of unequal intra-strand homologous recombination. Furthermore, we observe large segmental duplications involving both Non-Homologous End Joining and Homologous Recombination double-strand break repair pathways. Finally, the identification of a Mesoamerican-specific subtelomeric sequence reveals frequent interchromosomal recombinations between common bean subtelomeres. Altogether, our results highlight that common bean subtelomeres are hot spots of recombination and favor the rapid evolution of R genes. We propose that chromosome ends could act as R gene incubators in many plant genomes.

Whole genome sequencing-based association study to unravel genetic architecture of cooked grain width and length traits in rice

In this study, we used 2.9 million single nucleotide polymorphisms (SNP) and 393,429 indels derived from whole genome sequences of 591 rice landraces to determine the genetic basis of cooked and raw grain length, width and shape using genome-wide association study (GWAS). We identified a unique fine-mapped genetic region GWi7.1 significantly associated with cooked and raw grain width. Additionally, GWi7.2 that harbors GL7/GW7 a cloned gene for grain dimension was found. Novel regions in chromosomes 10 and 11 were also found to be associated with cooked grain shape and raw grain width, respectively. The indel-based GWAS identified fine-mapped genetic regions GL3.1 and GWi5.1 that matched synteny breakpoints between indica and japonica. GL3.1 was positioned a few kilobases away from GS3, a cloned gene for cooked and raw grain lengths in indica. GWi5.1 found to be significantly associated with cooked and raw grain width. It anchors upstream of cloned gene GW5, which varied between indica and japonica accessions. GWi11.1 is present inside the 3'-UTR of a functional gene in indica that corresponds to a syntenic break in chromosome 11 of japonica. Our results identified novel allelic structural variants and haplotypes confirmed using single locus and multilocus SNP and indel-based GWAS.

Whole genome sequencing-based association study to unravel genetic architecture of cooked grain width and length traits in rice

In this study, we used 2.9 million single nucleotide polymorphisms (SNP) and 393,429 indels derived from whole genome sequences of 591 rice landraces to determine the genetic basis of cooked and raw grain length, width and shape using genome-wide association study (GWAS). We identified a unique fine-mapped genetic region GWi7.1 significantly associated with cooked and raw grain width. Additionally, GWi7.2 that harbors GL7/GW7 a cloned gene for grain dimension was found. Novel regions in chromosomes 10 and 11 were also found to be associated with cooked grain shape and raw grain width, respectively. The indel-based GWAS identified fine-mapped genetic regions GL3.1 and GWi5.1 that matched synteny breakpoints between indica and japonica. GL3.1 was positioned a few kilobases away from GS3, a cloned gene for cooked and raw grain lengths in indica. GWi5.1 found to be significantly associated with cooked and raw grain width. It anchors upstream of cloned gene GW5, which varied between indica and japonica accessions. GWi11.1 is present inside the 3'-UTR of a functional gene in indica that corresponds to a syntenic break in chromosome 11 of japonica. Our results identified novel allelic structural variants and haplotypes confirmed using single locus and multilocus SNP and indel-based GWAS.

Comprehensive survey and evolutionary analysis of genome-wide miRNA genes from ten diploid Oryza species

BACKGROUND

MicroRNAs (miRNAs) are non-coding RNAs that play versatile roles in post-transcriptional gene regulation. Although much is known about their biogenesis, and gene regulation very little is known about their evolutionary relation among the closely related species.

RESULT

All the orthologous miRNA genes of Oryza sativa (japonica) from 10 different Oryza species were identified, and the evolutionary changes among these genes were analysed. Significant differences in the expansion of miRNA gene families were observed across the Oryza species. Analysis of the nucleotide substitution rates indicated that the mature sequences show the least substitution rates among the different regions of miRNA genes, and also show a very much less substitution rates as compared to that of all protein-coding genes across the Oryza species. Evolution of miRNA genes was also found to be contributed by transposons. A non-neutral selection was observed at 80 different miRNA loci across Oryza species which were estimated to have lost ~87% of the sequence diversity during the domestication. The phylogenetic analysis revealed that O. longistaminata diverged first among the AA-genomes, whereas O. brachyantha and O. punctata appeared as the eminent out-groups. The miR1861 family organised into nine distinct compact clusters in the studied Oryza species except O. brachyantha. Further, the expression analysis showed that 11 salt-responsive miRNAs were differentially regulated between O. coarctata and O. glaberrima.

CONCLUSION

Our study provides the evolutionary dynamics in the miRNA genes of 10 different Oryza species which will support more investigations about the structural and functional organization of miRNA genes of Oryza species.

Rice: The First Crop Genome

Rice was the first sequenced crop genome, paving the way for the sequencing of additional and more complicated crop genomes. The impact that the genome sequence made on rice genetics and breeding research was immediate, as evidence by citations and DNA marker use. The impact on other crop genomes was evident too, particularly for those within the grass family. As we celebrate 10 years since the completion of the rice genome sequence, we look forward to new empowering tool sets that will further revolutionize research in rice genetics and breeding and result in varieties that will continue to feed a growing population.

Genome biogeography reveals the intraspecific spread of adaptive mutations for a complex trait

Physiological novelties are often studied at macro-evolutionary scales such that their micro-evolutionary origins remain poorly understood. Here, we test the hypothesis that key components of a complex trait can evolve in isolation and later be combined by gene flow. We use C4 photosynthesis as a study system, a derived physiology that increases plant productivity in warm, dry conditions. The grass Alloteropsis semialata includes C4 and non-C4 genotypes, with some populations using laterally acquired C4 -adaptive loci, providing an outstanding system to track the spread of novel adaptive mutations. Using genome data from C4 and non-C4 A. semialata individuals spanning the species' range, we infer and date past migrations of different parts of the genome. Our results show that photosynthetic types initially diverged in isolated populations, where key C4 components were acquired. However, rare but recurrent subsequent gene flow allowed the spread of adaptive loci across genetic pools. Indeed, laterally acquired genes for key C4 functions were rapidly passed between populations with otherwise distinct genomic backgrounds. Thus, our intraspecific study of C4 -related genomic variation indicates that components of adaptive traits can evolve separately and later be combined through secondary gene flow, leading to the assembly and optimization of evolutionary innovations.

The Evolutionary Basis of Naturally Diverse Rice Leaves Anatomy

Rice contains genetically and ecologically diverse wild and cultivated species that show a wide variation in plant and leaf architecture. A systematic characterization of leaf anatomy is essential in understanding the dynamics behind such diversity. Therefore, leaf anatomies of 24 Oryza species spanning 11 genetically diverse rice genomes were studied in both lateral and longitudinal directions and possible evolutionary trends were examined. A significant inter-species variation in mesophyll cells, bundle sheath cells, and vein structure was observed, suggesting precise genetic control over these major rice leaf anatomical traits. Cellular dimensions, measured along three growth axes, were further combined proportionately to construct three-dimensional (3D) leaf anatomy models to compare the relative size and orientation of the major cell types present in a fully expanded leaf. A reconstruction of the ancestral leaf state revealed that the following are the major characteristics of recently evolved rice species: fewer veins, larger and laterally elongated mesophyll cells, with an increase in total mesophyll area and in bundle sheath cell number. A huge diversity in leaf anatomy within wild and domesticated rice species has been portrayed in this study, on an evolutionary context, predicting a two-pronged evolutionary pathway leading to the 'sativa leaf type' that we see today in domesticated species.

Development of Oryza sativa L. by Oryza punctata Kotschy ex Steud. monosomic addition lines with high value traits by interspecific hybridization

This paper describes the development of monosomic alien addition and disomic introgression lines through a cross between autotetraploid indica rice and Oryza punctata toward tapping valuable traits for rice improvement. Oryza punctata is a distantly related wild Oryza species having BB genome with untapped genetic resources for rice improvement. Low crossability between the cultivated O. sativa and O. punctata restricts the success of transferring many desirable traits into cultivated rice. Artificially induced autotetraploids of an elite breeding line, IR31917-45-3-2, were produced and crossed with O. punctata. Allotriploid F1 plants were backcrossed to IR31917-45-3-2 and generated progenies with extra chromosomes from O. punctata. Twenty BC1F1 and 59 BC2F1 plants were produced with chromosome numbers ranging from 24 (2n) to 29 (2n + 5) and 2n (24) to 26 (2n + 2), respectively. Eleven monosomic alien addition lines (MAALs) were characterized morphologically and cytologically and designated as MAAL 1-12. MAALs were genotyped using O. punctata genome-specific molecular markers and detected chromosome segments inherited from O. punctata. O. punctata introgressions across all the chromosomes of O. sativa were identified except for chromosome 8. The most frequent introgressions were observed in chromosomes 4, 6, 10, and 11, which could be the recombination hotspots between A and B genomes. Some of the qualitative traits such as black hull, purple coleoptile base, purple stigma, long awn, and short grain size from O. punctata were inherited in some disomic introgression lines (DILs). Several DILs inherited genes from O. punctata conferring resistance to brown planthopper, green leafhopper, and diseases such as bacterial blight and blast. This is the first report on successful gene transfer from O. punctata into O. sativa.

Ricebase: a breeding and genetics platform for rice, integrating individual molecular markers, pedigrees and whole-genome-based data

Ricebase (http://ricebase.org) is an integrative genomic database for rice (Oryza sativa) with an emphasis on combining datasets in a way that maintains the key links between past and current genetic studies. Ricebase includes DNA sequence data, gene annotations, nucleotide variation data and molecular marker fragment size data. Rice research has benefited from early adoption and extensive use of simple sequence repeat (SSR) markers; however, the majority of rice SSR markers were developed prior to the latest rice pseudomolecule assembly. Interpretation of new research using SNPs in the context of literature citing SSRs requires a common coordinate system. A new pipeline, using a stepwise relaxation of stringency, was used to map SSR primers onto the latest rice pseudomolecule assembly. The SSR markers and experimentally assayed amplicon sizes are presented in a relational database with a web-based front end, and are available as a track loaded in a genome browser with links connecting the browser and database. The combined capabilities of Ricebase link genetic markers, genome context, allele states across rice germplasm and potentially user curated phenotypic interpretations as a community resource for genetic discovery and breeding in rice.

Genome-wide survey and expression analysis of the OSCA gene family in rice

BACKGROUND

Reception of and response to exogenous and endogenous osmotic changes is important to sustain plant growth and development, as well as reproductive formation. Hyperosmolality-gated calcium-permeable channels (OSCA) were first characterised as an osmosensor in Arabidopsis and are involved in the perception of extracellular changes to trigger hyperosmolality-induced [Ca(2+)]i increases (OICI). To explore the potential biological functions of OSCAs in rice, we performed a bioinformatics and expression analysis of the OsOSCA gene family.

RESULTS

A total of 11 OsOSCA genes were identified from the genome database of Oryza sativa L. Japonica. Based on their sequence composition and phylogenetic relationship, the OsOSCA family was classified into four clades. Gene and protein structure analysis indicated that the 11 OsOSCAs shared similar structures with their homologs in Oryza sativa L. ssp. Indica, Oryza glaberrima, and Oryza brachyantha. Multiple sequence alignment analysis revealed a conserved DUF221 domain in these members, in which the first three TMs were conserved, while the others were not. The expression profiles of OsOSCA genes were analysed at different stages of vegetative growth, reproductive development, and under osmotic-associated abiotic stresses. We found that four and six OsOSCA genes showed a clear correlation between the expression profile and osmotic changes during caryopsis development and seed imbibition, respectively. Orchestrated transcription of three OsOSCAs was strongly associated with the circadian clock. Moreover, osmotic-related abiotic stress differentially induced the expression of 10 genes.

CONCLUSION

The entire OSCA family is characterised by the presence of a conserved DUF221 domain, which functions as an osmotic-sensing calcium channel. The phylogenetic tree of OSCA genes showed that two subspecies of cultivated rice, Oryza sativa L. ssp. Japonica and Oryza sativa L. ssp. Indica, are more closely related than wild rice Oryza glaberrima, while Oryza brachyantha was less closely related. OsOSCA expression is organ- and tissue-specific and regulated by different osmotic-related abiotic stresses in rice. These findings will facilitate further research in this gene family and provide potential target genes for generation of genetically modified osmotic-stress-resistant plants.

Range expansion and habitat shift triggered elevated diversification of the rice genus (Oryza, Poaceae) during the Pleistocene

Background

The rice genus (Oryza) contains many wild genetic resources that are vital to the well-being of humans. However, little is known about the process by which the genus diversified or the factors that drove its speciation. Here, we integrated the phylogenetic, molecular dating and biogeographic methods to investigate the spatial-temporal patterns of Oryza diversification, and used a series of model tests to examine whether intercontinental migrations and/or key innovations were associated with significant changes in diversification rates in the genus.

Results

Oryza became differentiated in tropical Asia in the Miocene. There were two migrations from the ancestral area into Africa and Australia during the Miocene. We inferred at least 10 migration events out of tropical Asia since the Pleistocene, mainly involving the species adapting open habitat. A rapid increase in diversification rates of the whole Oryza occurred during the Pleistocene. Intercontinental migrations from tropical Asia to other tropical regions were positively correlated with shift in habitat, but not with changes in life history. A habitat preference shift from shade tolerant to open habitat predated the burst in diversification rates.

Conclusions

Rice species may have been pre-adapted to invade open habitat. Significant increase in diversification rates occurred during the Pleistocene and is associated with range expansion and habitat shift, but not with life history. The rice genus provides an excellent case supporting the idea that range expansion and invasion of novel habitats can drive the diversification of a group.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0459-1) contains supplementary material, which is available to authorized users.

Genetic shift in local rice populations during rice breeding programs in the northern limit of rice cultivation in the world

The rapid accumulation of pre-existing mutations may play major roles in the establishment and shaping of adaptability for local regions in current rice breeding programs. The cultivated rice, Oryza sativa L., which originated from tropical regions, is now grown worldwide due to the concerted efforts of breeding programs. However, the process of establishing local populations and their origins remain unclear. In the present study, we characterized DNA polymorphisms in the rice variety KITAAKE from Hokkaido, one of the northern limits of rice cultivation in the world. Indel polymorphisms were attributed to transposable element-like insertions, tandem duplications, and non-TE deletions as the original mutation events in the NIPPONBARE and KITAAKE genomes. The allele frequencies of the KITAAKE alleles markedly shifted to the current variety types among the local population from Hokkaido in the last two decades. The KITAAKE alleles widely distributed throughout wild rice and cultivated rice over the world. These have accumulated in the local population from Hokkaido via Japanese landraces as the ancestral population of Hokkaido. These results strongly suggested that combinations of pre-existing mutations played a role in the establishment of adaptability. This approach using the re-sequencing of local varieties in unique environmental conditions will be useful as a genetic resource in plant breeding programs in local regions.

Rapid diversification of five Oryza AA genomes associated with rice adaptation

Comparative genomic analyses among closely related species can greatly enhance our understanding of plant gene and genome evolution. We report de novo-assembled AA-genome sequences for Oryza nivara, Oryza glaberrima, Oryza barthii, Oryza glumaepatula, and Oryza meridionalis. Our analyses reveal massive levels of genomic structural variation, including segmental duplication and rapid gene family turnover, with particularly high instability in defense-related genes. We show, on a genomic scale, how lineage-specific expansion or contraction of gene families has led to their morphological and reproductive diversification, thus enlightening the evolutionary process of speciation and adaptation. Despite strong purifying selective pressures on most Oryza genes, we documented a large number of positively selected genes, especially those genes involved in flower development, reproduction, and resistance-related processes. These diversifying genes are expected to have played key roles in adaptations to their ecological niches in Asia, South America, Africa and Australia. Extensive variation in noncoding RNA gene numbers, function enrichment, and rates of sequence divergence might also help account for the different genetic adaptations of these rice species. Collectively, these resources provide new opportunities for evolutionary genomics, numerous insights into recent speciation, a valuable database of functional variation for crop improvement, and tools for efficient conservation of wild rice germplasm.

Diversification of the C-TERMINALLY ENCODED PEPTIDE (CEP) gene family in angiosperms, and evolution of plant-family specific CEP genes

BACKGROUND

Small, secreted signaling peptides work in parallel with phytohormones to control important aspects of plant growth and development. Genes from the C-TERMINALLY ENCODED PEPTIDE (CEP) family produce such peptides which negatively regulate plant growth, especially under stress, and affect other important developmental processes. To illuminate how the CEP gene family has evolved within the plant kingdom, including its emergence, diversification and variation between lineages, a comprehensive survey was undertaken to identify and characterize CEP genes in 106 plant genomes.

RESULTS

Using a motif-based system developed for this study to identify canonical CEP peptide domains, a total of 916 CEP genes and 1,223 CEP domains were found in angiosperms and for the first time in gymnosperms. This defines a narrow band for the emergence of CEP genes in plants, from the divergence of lycophytes to the angiosperm/gymnosperm split. Both CEP genes and domains were found to have diversified in angiosperms, particularly in the Poaceae and Solanaceae plant families. Multispecies orthologous relationships were determined for 22% of identified CEP genes, and further analysis of those groups found selective constraints upon residues within the CEP peptide and within the previously little-characterized variable region. An examination of public Oryza sativa RNA-Seq datasets revealed an expression pattern that links OsCEP5 and OsCEP6 to panicle development and flowering, and CEP gene trees reveal these emerged from a duplication event associated with the Poaceae plant family.

CONCLUSIONS

The characterization of the plant-family specific CEP genes OsCEP5 and OsCEP6, the association of CEP genes with angiosperm-specific development processes like panicle development, and the diversification of CEP genes in angiosperms provides further support for the hypothesis that CEP genes have been integral to the evolution of novel traits within the angiosperm lineage. Beyond these findings, the comprehensive set of CEP genes and their properties reported here will be a resource for future research on CEP genes and peptides.

Identification and diversity of functional centromere satellites in the wild rice species Oryza brachyantha

The centromere is a key chromosomal component for sister chromatid cohesion and is the site for kinetochore assembly and spindle fiber attachment, allowing each sister chromatid to faithfully segregate to each daughter cell during cell division. It is not clear what types of sequences act as functional centromeres and how centromere sequences are organized in Oryza brachyantha, an FF genome species. In this study, we found that the three classes of centromere-specific CentO-F satellites (CentO-F1, CentO-F2, and CentOF3) in O. brachyantha share no homology with the CentO satellites in Oryza sativa. The three classes of CentO-F satellites are all located within the chromosomal regions to which the spindle fibers attach and are characterized by megabase tandem arrays that are flanked by centromere-specific retrotransposons, CRR-F, in the O. brachyantha centromeres. Although these CentO-F satellites are quantitatively variable among 12 O. brachyantha centromeres, immunostaining with an antibody specific to CENH3 indicates that they are colocated with CENH3 in functional centromere regions. Our results demonstrate that the three classes of CentO-F satellites may be the major components of functional centromeres in O. brachyantha.

The genome sequence of African rice (Oryza glaberrima) and evidence for independent domestication

Mingsheng Chen, Klaus Mayer, Steve Rounsley, Rod Wing and colleagues report the genome sequence of African rice (Oryza glaberrima), a different species than Asian rice. The authors resequenced 20 O. glaberrima accessions and 94 Oryza barthii accessions (the putative progenitor species of O. glaberrima), and their analyses support the hypothesis that O. glaberrima was domesticated in a single region along the upper Niger river.

The cultivation of rice in Africa dates back more than 3,000 years. Interestingly, African rice is not of the same origin as Asian rice (Oryza sativa L.) but rather is an entirely different species (i.e., Oryza glaberrima Steud.). Here we present a high-quality assembly and annotation of the O. glaberrima genome and detailed analyses of its evolutionary history of domestication and selection. Population genomics analyses of 20 O. glaberrima and 94 Oryza barthii accessions support the hypothesis that O. glaberrima was domesticated in a single region along the Niger river as opposed to noncentric domestication events across Africa. We detected evidence for artificial selection at a genome-wide scale, as well as with a set of O. glaberrima genes orthologous to O. sativa genes that are known to be associated with domestication, thus indicating convergent yet independent selection of a common set of genes during two geographically and culturally distinct domestication processes.

GC content evolution in coding regions of angiosperm genomes: a unifying hypothesis

In angiosperms (as in other species), GC content varies along and between genes, within a genome, and between genomes of different species, but the reason for this distribution is still an open question. Grass genomes are particularly intriguing because they exhibit a strong bimodal distribution of genic GC content and a sharp 5'-3' decreasing GC content gradient along most genes. Here, we propose a unifying model to explain the main patterns of GC content variation at the gene and genome scale. We argue that GC content patterns could be mainly determined by the interactions between gene structure, recombination patterns, and GC-biased gene conversion. Recent studies on fine-scale recombination maps in angiosperms support this hypothesis and previous results also fit this model. We propose that our model could be used as a null hypothesis to search for additional forces that affect GC content in angiosperms.

Rice MtN3/saliva/SWEET gene family: Evolution, expression profiling, and sugar transport

The rice MtN3/saliva/SWEET gene family consists of 21 paralogs. However, their functions in physiological processes are largely unknown, although at least three of the 21 paralogs are used by pathogenic bacteria to infect rice. Here, we report the evolutionary features, transcriptional characteristics, and putative functions in sugar transport of this gene family. The wild rice accessions in this study included those with AA, BB, CC, BBCC, CCDD, EE, and GG genomes, which appeared approximately 0.58-14.6 million years ago. The structures, chromosomal locations, phylogenetic relationships, and homologous distribution among the accessions suggest that the number of rice MtN3/saliva/SWEET paralogs gradually increased as the Oryza genus evolved, and one third of the paralogs may have originated recently. These paralogs are differentially expressed in vegetative and reproductive tissues, in the leaf senescence process, and in signaling dependent on gibberellic acid, cytokinin, or 1-naphthalene acetic acid (an analog of auxin), suggesting that they may be associated with multiple physiological processes. Four paralogs could transport galactose in yeast, which suggests that they may have a similar function in rice. These results will help to elucidate their roles and biochemical functions in rice development, adaptation to environment, host-pathogen interaction, and so forth.

Evolution of centromeric retrotransposons in grasses

Centromeric retrotransposons (CRs) constitute a family of plant retroelements, some of which have the ability to target their insertion almost exclusively to the functional centromeres. Our exhaustive analysis of CR family members in four grass genomes revealed not only horizontal transfer (HT) of CR elements between the oryzoid and panicoid grass lineages but also their subsequent recombination with endogenous elements that in some cases created prolific recombinants in foxtail millet and sorghum. HT events are easily identifiable only in cases where host genome divergence significantly predates HT, thus documented HT events likely represent only a fraction of the total. If the more difficult to detect ancient HT events occurred at frequencies similar to those observable in present day grasses, the extant long terminal repeat retrotransposons represent the mosaic products of HT and recombination that are optimized for retrotransposition in their host genomes. This complicates not only phylogenetic analysis but also the establishment of a meaningful retrotransposon nomenclature, which we have nevertheless attempted to implement here. In contrast to the plant-centric naming convention used currently for CR elements, we classify elements primarily based on their phylogenetic relationships regardless of host plant, using the exhaustively studied maize elements assigned to six different subfamilies as a standard. The CR2 subfamily is the most widely distributed of the six CR subfamilies discovered in grass genomes to date and thus the most likely to play a functional role at grass centromeres.

Fifteen million years of evolution in the Oryza genus shows extensive gene family expansion

In analyzing gene families in the whole-genome sequences available for O. sativa (AA), O. glaberrima (AA), and O. brachyantha (FF), we observed large size expansions in the AA genomes compared to FF genomes for the super-families F-box and NB-ARC, and five additional families: the Aspartic proteases, BTB/POZ proteins (BTB), Glutaredoxins, Trypsin α-amylase inhibitor proteins, and Zf-Dof proteins. Their evolutionary dynamic was investigated to understand how and why such important size variations are observed between these closely related species. We show that expansions resulted from both amplification, largely by tandem duplications, and contraction by gene losses. For the F-box and NB-ARC gene families, the genes conserved in all species were under strong purifying selection while expanded orthologous genes were under more relaxed purifying selection. In F-box, NB-ARC, and BTB, the expanded groups were enriched in genes with little evidence of expression, in comparison with conserved groups. We also detected 87 loci under positive selection in the expanded groups. These results show that most of the duplicated copies in the expanded groups evolve neutrally after duplication because of functional redundancy but a fraction of these genes were preserved following neofunctionalization. Hence, the lineage-specific expansions observed between Oryza species were partly driven by directional selection.

Boom-Bust Turnovers of Megabase-Sized Centromeric DNA in Solanum Species: Rapid Evolution of DNA Sequences Associated with Centromeres

Centromeres are composed of long arrays of satellite repeats in most multicellular eukaryotes investigated to date. The satellite repeat-based centromeres are believed to have evolved from "neocentromeres" that originally contained only single- or low-copy sequences. However, the emergence and evolution of the satellite repeats in centromeres has been elusive. Potato (Solanum tuberosum) provides a model system for studying centromere evolution because each of its 12 centromeres contains distinct DNA sequences, allowing comparative analysis of homoeologous centromeres from related species. We conducted genome-wide analysis of the centromeric sequences in Solanum verrucosum, a wild species closely related to potato. Unambiguous homoeologous centromeric sequences were detected in only a single centromere (Cen9) between the two species. Four centromeres (Cen2, Cen4, Cen7, and Cen10) in S. verrucosum contained distinct satellite repeats that were amplified from retrotransposon-related sequences. Strikingly, the same four centromeres in potato contain either different satellite repeats (Cen2 and Cen7) or exclusively single- and low-copy sequences (Cen4 and Cen10). Our sequence comparison of five homoeologous centromeres in two Solanum species reveals rapid divergence of centromeric sequences among closely related species. We propose that centromeric satellite repeats undergo boom-bust cycles before a favorable repeat is fixed in the population.

Major repeat components covering one-third of the ginseng (Panax ginseng C.A. Meyer) genome and evidence for allotetraploidy

Ginseng (Panax ginseng) is a famous medicinal herb, but the composition and structure of its genome are largely unknown. Here we characterized the major repeat components and inspected their distribution in the ginseng genome. By analyzing three repeat-rich bacterial artificial chromosome (BAC) sequences from ginseng, we identified complex insertion patterns of 34 long terminal repeat retrotransposons (LTR-RTs) and 11 LTR-RT derivatives accounting for more than 80% of the BAC sequences. The LTR-RTs were classified into three Ty3/gypsy (PgDel, PgTat and PgAthila) and two Ty1/Copia (PgTork and PgOryco) families. Mapping of 30-Gbp Illumina whole-genome shotgun reads to the BAC sequences revealed that these five LTR-RT families occupy at least 34% of the ginseng genome. The Ty3/Gypsy families were predominant, comprising 74 and 33% of the BAC sequences and the genome, respectively. In particular, the PgDel family accounted for 29% of the genome and presumably played major roles in enlargement of the size of the ginseng genome. Fluorescence in situ hybridization (FISH) revealed that the PgDel1 elements are distributed throughout the chromosomes along dispersed heterochromatic regions except for ribosomal DNA blocks. The intensity of the PgDel2 FISH signals was biased toward 24 out of 48 chromosomes. Unique gene probes showed two pairs of signals with different locations, one pair in subtelomeric regions on PgDel2-rich chromosomes and the other in interstitial regions on PgDel2-poor chromosomes, demonstrating allotetraploidy in ginseng. Our findings promote understanding of the evolution of the ginseng genome and of that of related species in the Araliaceae.

Global genomic diversity of Oryza sativa varieties revealed by comparative physical mapping

Bacterial artificial chromosome (BAC) physical maps embedding a large number of BAC end sequences (BESs) were generated for Oryza sativa ssp. indica varieties Minghui 63 (MH63) and Zhenshan 97 (ZS97) and were compared with the genome sequences of O. sativa spp. japonica cv. Nipponbare and O. sativa ssp. indica cv. 93-11. The comparisons exhibited substantial diversities in terms of large structural variations and small substitutions and indels. Genome-wide BAC-sized and contig-sized structural variations were detected, and the shared variations were analyzed. In the expansion regions of the Nipponbare reference sequence, in comparison to the MH63 and ZS97 physical maps, as well as to the previously constructed 93-11 physical map, the amounts and types of the repeat contents, and the outputs of gene ontology analysis, were significantly different from those of the whole genome. Using the physical maps of four wild Oryza species from OMAP (http://www.omap.org) as a control, we detected many conserved and divergent regions related to the evolution process of O. sativa. Between the BESs of MH63 and ZS97 and the two reference sequences, a total of 1532 polymorphic simple sequence repeats (SSRs), 71,383 SNPs, 1767 multiple nucleotide polymorphisms, 6340 insertions, and 9137 deletions were identified. This study provides independent whole-genome resources for intra- and intersubspecies comparisons and functional genomics studies in O. sativa. Both the comparative physical maps and the GBrowse, which integrated the QTL and molecular markers from GRAMENE (http://www.gramene.org) with our physical maps and analysis results, are open to the public through our Web site (http://gresource.hzau.edu.cn/resource/resource.html).

High occurrence of functional new chimeric genes in survey of rice chromosome 3 short arm genome sequences

In an effort to identify newly evolved genes in rice, we searched the genomes of Asian-cultivated rice Oryza sativa ssp. japonica and its wild progenitors, looking for lineage-specific genes. Using genome pairwise comparison of approximately 20-Mb DNA sequences from the chromosome 3 short arm (Chr3s) in six rice species, O. sativa, O. nivara, O. rufipogon, O. glaberrima, O. barthii, and O. punctata, combined with synonymous substitution rate tests and other evidence, we were able to identify potential recently duplicated genes, which evolved within the last 1 Myr. We identified 28 functional O. sativa genes, which likely originated after O. sativa diverged from O. glaberrima. These genes account for around 1% (28/3,176) of all annotated genes on O. sativa's Chr3s. Among the 28 new genes, two recently duplicated segments contained eight genes. Fourteen of the 28 new genes consist of chimeric gene structure derived from one or multiple parental genes and flanking targeting sequences. Although the majority of these 28 new genes were formed by single or segmental DNA-based gene duplication and recombination, we found two genes that were likely originated partially through exon shuffling. Sequence divergence tests between new genes and their putative progenitors indicated that new genes were most likely evolving under natural selection. We showed all 28 new genes appeared to be functional, as suggested by Ka/Ks analysis and the presence of RNA-seq, cDNA, expressed sequence tag, massively parallel signature sequencing, and/or small RNA data. The high rate of new gene origination and of chimeric gene formation in rice may demonstrate rice's broad diversification, domestication, its environmental adaptation, and the role of new genes in rice speciation.