Assembling the genome of the African wild rice Oryza longistaminata by exploiting synteny in closely related Oryza species

First draft reference genome and annotation of the alternative oil species Physaria fendleri

In the wake of increasing demand for renewable energy sources, plant-based sources including alternative oilseeds have come to the forefront of interest. Hydroxy fatty acids (HFAs), produced in a few oilseed species, are important chemical feed stocks for industrial applications. An integrated approach was taken to assemble the first draft genome of the alternative HFA producer Physaria fendleri (n = 6), an outcrossing species with high heterozygosity. Both de novo transcriptome assemblies and genome assemblies were produced with public and generated sequencing reads. Resulting intermediate assemblies were then scaffolded and patched with multiple data sources, followed by super-scaffolding onto a masked genome of Camelina laxa (n = 6). Despite a current lack of available resources for the physical mapping of genomic scaffolds of P. fendleri, topography of the genome with respect to repeat and gene content was preserved at the scaffold level and not significantly lost via super-scaffolding. Read representation, gene and genome completion statistics, and annotation results illustrated the creation of a functional draft genome and a tool for future research on alternative oil species.

Spatiotemporal transcriptomic atlas of rhizome formation in Oryza longistaminata

Rhizomes are modified stems that grow underground and produce new individuals genetically identical to the mother plant. Recently, a breakthrough has been made in efforts to convert annual grains into perennial ones by utilizing wild rhizomatous species as donors, yet the developmental biology of this organ is rarely studied. Oryza longistaminata, a wild rice species featuring strong rhizomes, provides a valuable model for exploration of rhizome development. Here, we first assembled a double-haplotype genome of O. longistaminata, which displays a 48-fold improvement in contiguity compared to the previously published assembly. Furthermore, spatiotemporal transcriptomics was performed to obtain the expression profiles of different tissues in O. longistaminata rhizomes and tillers. Two spatially reciprocal cell clusters, the vascular bundle 2 cluster and the parenchyma 2 cluster, were determined to be the primary distinctions between the rhizomes and tillers. We also captured meristem initiation cells in the sunken area of parenchyma located at the base of internodes, which is the starting point for rhizome initiation. Trajectory analysis further indicated that the rhizome is regenerated through de novo generation. Collectively, these analyses revealed a spatiotemporal transcriptional transition underlying the rhizome initiation, providing a valuable resource for future perennial crop breeding.

Temperature Effect on Rhizome Development in Perennial rice

Traditional agriculture is becoming increasingly not adapted to global climate change. Compared with annual rice, perennial rice has strong environmental adaptation and needs fewer natural resources and labor inputs. Rhizome, a kind of underground stem for rice to achieve perenniallity, can grow underground horizontally and then bend upward, developing into aerial stems. The temperature has a great influence on plant development. To date, the effect of temperature on rhizome development is still unknown. Fine temperature treatment of Oryza longistaminata (OL) proved that compared with higher temperatures (28-30 ℃), lower temperature (17-19 ℃) could promote the sprouting of axillary buds and enhance negative gravitropism of branches, resulting in shorter rhizomes. The upward growth of branches was earlier at low temperature than that at high temperature, leading to a high frequency of shorter rhizomes and smaller branch angles. Comparative transcriptome showed that plant hormones played an essential role in the response of OL to temperature. The expressions of ARF17, ARF25 and FucT were up-regulated at low temperature, resulting in prospectively asymmetric auxin distribution, which subsequently induced asymmetric expression of IAA20 and WOX11 between the upper and lower side of the rhizome, further leading to upward growth of the rhizome. Cytokinin and auxin are phytohormones that can promote and inhibit bud outgrowth, respectively. The auxin biosynthesis gene YUCCA1 and cytokinin oxidase/dehydrogenase gene CKX4 and CKX9 were up-regulated, while cytokinin biosynthesis gene IPT4 was down-regulated at high temperature. Moreover, the D3 and D14 in strigolactones pathways, negatively regulating bud outgrowth, were up-regulated at high temperature. These results indicated that cytokinin, auxins, and strigolactones jointly control bud outgrowth at different temperatures. Our research revealed that the outgrowth of axillary bud and the upward growth of OL rhizome were earlier at lower temperature, providing clues for understanding the rhizome growth habit under different temperatures, which would be helpful for cultivating perennial rice.

A haplotype-resolved genome provides insight into allele-specific expression in wild walnut (Juglans regia L.)

Wild germplasm resources are crucial for gene mining and molecular breeding because of their special trait performance. Haplotype-resolved genome is an ideal solution for fully understanding the biology of subgenomes in highly heterozygous species. Here, we surveyed the genome of a wild walnut tree from Gongliu County, Xinjiang, China, and generated a haplotype-resolved reference genome of 562.99 Mb (contig N50 = 34.10 Mb) for one haplotype (hap1) and 561.07 Mb (contig N50 = 33.91 Mb) for another haplotype (hap2) using PacBio high-fidelity (HiFi) reads and Hi-C technology. Approximately 527.20 Mb (93.64%) of hap1 and 526.40 Mb (93.82%) of hap2 were assigned to 16 pseudochromosomes. A total of 41039 and 39744 protein-coding gene models were predicted for hap1 and hap2, respectively. Moreover, 123 structural variations (SVs) were identified between the two haplotype genomes. Allele-specific expression genes (ASEGs) that respond to cold stress were ultimately identified. These datasets can be used to study subgenome evolution, for functional elite gene mining and to discover the transcriptional basis of specific traits related to environmental adaptation in wild walnut.

MCPtaggR: R package for accurate genotype calling in reduced representation sequencing data by eliminating error-prone markers based on genome comparison

Reduced representation sequencing (RRS) offers cost-effective, high-throughput genotyping platforms such as genotyping-by-sequencing (GBS). RRS reads are typically mapped onto a reference genome. However, mapping reads harbouring mismatches against the reference can potentially result in mismapping and biased mapping, leading to the detection of error-prone markers that provide incorrect genotype information. We established a genotype-calling pipeline named mappable collinear polymorphic tag genotyping (MCPtagg) to achieve accurate genotyping by eliminating error-prone markers. MCPtagg was designed for the RRS-based genotyping of a population derived from a biparental cross. The MCPtagg pipeline filters out error-prone markers prior to genotype calling based on marker collinearity information obtained by comparing the genome sequences of the parents of a population to be genotyped. A performance evaluation on real GBS data from a rice F2 population confirmed its effectiveness. Furthermore, our performance test using a genome assembly that was obtained by genome sequence polishing on an available genome assembly suggests that our pipeline performs well with converted genomes, rather than necessitating de novo assembly. This demonstrates its flexibility and scalability. The R package, MCPtaggR, was developed to provide functions for the pipeline and is available at https://github.com/tomoyukif/MCPtaggR.

Insights into the Transcriptomics of Crop Wild Relatives to Unravel the Salinity Stress Adaptive Mechanisms

The narrow genomic diversity of modern cultivars is a major bottleneck for enhancing the crop's salinity stress tolerance. The close relatives of modern cultivated plants, crop wild relatives (CWRs), can be a promising and sustainable resource to broaden the diversity of crops. Advances in transcriptomic technologies have revealed the untapped genetic diversity of CWRs that represents a practical gene pool for improving the plant's adaptability to salt stress. Thus, the present study emphasizes the transcriptomics of CWRs for salinity stress tolerance. In this review, the impacts of salt stress on the plant's physiological processes and development are overviewed, and the transcription factors (TFs) regulation of salinity stress tolerance is investigated. In addition to the molecular regulation, a brief discussion on the phytomorphological adaptation of plants under saline environments is provided. The study further highlights the availability and use of transcriptomic resources of CWR and their contribution to pangenome construction. Moreover, the utilization of CWRs' genetic resources in the molecular breeding of crops for salinity stress tolerance is explored. Several studies have shown that cytoplasmic components such as calcium and kinases, and ion transporter genes such as Salt Overly Sensitive 1 (SOS1) and High-affinity Potassium Transporters (HKTs) are involved in the signaling of salt stress, and in mediating the distribution of excess Na⁺ ions within the plant cells. Recent comparative analyses of transcriptomic profiling through RNA sequencing (RNA-Seq) between the crops and their wild relatives have unraveled several TFs, stress-responsive genes, and regulatory proteins for generating salinity stress tolerance. This review specifies that the use of CWRs transcriptomics in combination with modern breeding experimental approaches such as genomic editing, de novo domestication, and speed breeding can accelerate the CWRs utilization in the breeding programs for enhancing the crop's adaptability to saline conditions. The transcriptomic approaches optimize the crop genomes with the accumulation of favorable alleles that will be indispensable for designing salt-resilient crops.

Genomic diversity of aquaporins across genus Oryza provides a rich genetic resource for development of climate resilient rice cultivars

BACKGROUND

Plant aquaporins are critical genetic players performing multiple biological functions, especially climate resilience and water-use efficiency. Their genomic diversity across genus Oryza is yet to be explored.

RESULTS

This study identified 369 aquaporin-encoding genes from 11 cultivated and wild rice species and further categorized these into four major subfamilies, among which small basic intrinsic proteins are speculated to be ancestral to all land plant aquaporins. Evolutionarily conserved motifs in peptides of aquaporins participate in transmembrane transport of materials and their relatively complex gene structures provide an evolutionary playground for regulation of genome structure and transcription. Duplication and evolution analyses revealed higher genetic conservation among Oryza aquaporins and strong purifying selections are assisting in conserving the climate resilience associated functions. Promoter analysis highlighted enrichment of gene upstream regions with cis-acting regulatory elements involved in diverse biological processes, whereas miRNA target site prediction analysis unveiled substantial involvement of osa-miR2102-3p, osa-miR2927 and osa-miR5075 in post-transcriptional regulation of gene expression patterns. Moreover, expression patterns of japonica aquaporins were significantly perturbed in response to different treatment levels of six phytohormones and four abiotic stresses, suggesting their multifarious roles in plants survival under stressed environments. Furthermore, superior haplotypes of seven conserved orthologous aquaporins for higher thousand-grain weight are reported from a gold mine of 3,010 sequenced rice pangenomes.

CONCLUSIONS

This study unveils the complete genomic atlas of aquaporins across genus Oryza and provides a comprehensive genetic resource for genomics-assisted development of climate-resilient rice cultivars.

Can the Wild Perennial, Rhizomatous Rice Species Oryza longistaminata be a Candidate for De Novo Domestication?

As climate change intensifies, the development of resilient rice that can tolerate abiotic stresses is urgently needed. In nature, many wild plants have evolved a variety of mechanisms to protect themselves from environmental stresses. Wild relatives of rice may have abundant and virtually untapped genetic diversity and are an essential source of germplasm for the improvement of abiotic stress tolerance in cultivated rice. Unfortunately, the barriers of traditional breeding approaches, such as backcrossing and transgenesis, make it challenging and complex to transfer the underlying resilience traits between plants. However, de novo domestication via genome editing is a quick approach to produce rice with high yields from orphans or wild relatives. African wild rice, Oryza longistaminata, which is part of the AA-genome Oryza species has two types of propagation strategies viz. vegetative propagation via rhizome and seed propagation. It also shows tolerance to multiple types of abiotic stress, and therefore O. longistaminata is considered a key candidate of wild rice for heat, drought, and salinity tolerance, and it is also resistant to lodging. Importantly, O. longistaminata is perennial and propagates also via rhizomes both of which are traits that are highly valuable for the sustainable production of rice. Therefore, O. longistaminata may be a good candidate for de novo domestication through genome editing to obtain rice that is more climate resilient than modern elite cultivars of O. sativa.

Genome-Wide Characterization of the PIFs Family in Sweet Potato and Functional Identification of IbPIF3.1 under Drought and Fusarium Wilt Stresses

Phytochrome-interacting factors (PIFs) are essential for plant growth, development, and defense responses. However, research on the PIFs in sweet potato has been insufficient to date. In this study, we identified PIF genes in the cultivated hexaploid sweet potato (Ipomoea batatas) and its two wild relatives, Ipomoea triloba, and Ipomoea trifida. Phylogenetic analysis revealed that IbPIFs could be divided into four groups, showing the closest relationship with tomato and potato. Subsequently, the PIFs protein properties, chromosome location, gene structure, and protein interaction network were systematically analyzed. RNA-Seq and qRT-PCR analyses showed that IbPIFs were mainly expressed in stem, as well as had different gene expression patterns in response to various stresses. Among them, the expression of IbPIF3.1 was strongly induced by salt, drought, H₂O₂, cold, heat, Fusarium oxysporum f. sp. batatas (Fob), and stem nematodes, indicating that IbPIF3.1 might play an important role in response to abiotic and biotic stresses in sweet potato. Further research revealed that overexpression of IbPIF3.1 significantly enhanced drought and Fusarium wilt tolerance in transgenic tobacco plants. This study provides new insights for understanding PIF-mediated stress responses and lays a foundation for future investigation of sweet potato PIFs.

A super pan-genomic landscape of rice

Pan-genomes from large natural populations can capture genetic diversity and reveal genomic complexity. Using de novo long-read assembly, we generated a graph-based super pan-genome of rice consisting of a 251-accession panel comprising both cultivated and wild species of Asian and African rice. Our pan-genome reveals extensive structural variations (SVs) and gene presence/absence variations. Additionally, our pan-genome enables the accurate identification of nucleotide-binding leucine-rich repeat genes and characterization of their inter- and intraspecific diversity. Moreover, we uncovered grain weight-associated SVs which specify traits by affecting the expression of their nearby genes. We characterized genetic variants associated with submergence tolerance, seed shattering and plant architecture and found independent selection for a common set of genes that drove adaptation and domestication in Asian and African rice. This super pan-genome facilitates pinpointing of lineage-specific haplotypes for trait-associated genes and provides insights into the evolutionary events that have shaped the genomic architecture of various rice species.

Regulation of ammonium acquisition and use in Oryza longistaminata ramets under nitrogen source heterogeneity

Oryza longistaminata, a wild rice, vegetatively reproduces and forms a networked clonal colony consisting of ramets connected by rhizomes. Although water, nutrients, and other molecules can be transferred between ramets via the rhizomes, inter-ramet communication in response to spatially heterogeneous nitrogen availability is not well understood. We studied the response of ramet pairs to heterogeneous nitrogen availability using a split hydroponic system that allowed each ramet root to be exposed to different conditions. Ammonium uptake was compensatively enhanced in the sufficient-side root when roots of the ramet pairs were exposed to ammonium-sufficient and ammonium-deficient conditions. Comparative transcriptome analysis revealed that a gene regulatory network for effective ammonium assimilation and amino acid biosynthesis was activated in the sufficient-side roots. Allocation of absorbed nitrogen from the nitrogen-sufficient to the nitrogen-deficient ramets was rather limited. Nitrogen was preferentially used for newly growing axillary buds on the sufficient-side ramets. Biosynthesis of trans-zeatin (tZ), a cytokinin, was upregulated in response to the nitrogen supply, but tZ appeared not to target the compensatory regulation. Our results also implied that the O. longistaminata putative ortholog of rice (Oryza sativa) C-terminally encoded peptide1 plays a role as a nitrogen-deficient signal in inter-ramet communication, providing compensatory upregulation of nitrogen assimilatory genes. These results provide insights into the molecular basis for efficient growth strategies of asexually proliferating plants growing in areas where the distribution of ammonium ions is spatially heterogeneous.

Genome assembly, annotation, and comparative analysis of the cattail Typha latifolia

Cattails (Typha species) comprise a genus of emergent wetland plants with a global distribution. Typha latifolia and Typha angustifolia are two of the most widespread species, and in areas of sympatry can interbreed to produce the hybrid Typha × glauca. In some regions, the relatively high fitness of Typha × glauca allows it to outcompete and displace both parent species, while simultaneously reducing plant and invertebrate biodiversity, and modifying nutrient and water cycling. We generated a high-quality whole-genome assembly of T. latifolia using PacBio long-read and high coverage Illumina sequences that will facilitate evolutionary and ecological studies in this hybrid zone. Genome size was 287 Mb and consisted of 1158 scaffolds, with an N50 of 8.71 Mb; 43.84% of the genome were identified as repetitive elements. The assembly has a BUSCO score of 96.03%, and 27,432 genes and 2700 RNA sequences were putatively identified. Comparative analysis detected over 9000 shared orthologs with related taxa and phylogenomic analysis supporting T. latifolia as a divergent lineage within Poales. This high-quality scaffold-level reference genome will provide a useful resource for future population genomic analyses and improve our understanding of Typha hybrid dynamics.

Genome Size Variation and Evolution Driven by Transposable Elements in the Genus Oryza

Genome size variation and evolutionary forces behind have been long pursued in flowering plants. The genus Oryza, consisting of approximately 25 wild species and two cultivated rice, harbors eleven extant genome types, six of which are diploid (AA, BB, CC, EE, FF, and GG) and five of which are tetraploid (BBCC, CCDD, HHJJ, HHKK, and KKLL). To obtain the most comprehensive knowledge of genome size variation in the genus Oryza, we performed flow cytometry experiments and estimated genome sizes of 166 accessions belonging to 16 non-AA genome Oryza species. k-mer analyses were followed to verify the experimental results of the two accessions for each species. Our results showed that genome sizes largely varied fourfold in the genus Oryza, ranging from 279 Mb in Oryza brachyantha (FF) to 1,203 Mb in Oryza ridleyi (HHJJ). There was a 2-fold variation (ranging from 570 to 1,203 Mb) in genome size among the tetraploid species, while the diploid species had 3-fold variation, ranging from 279 Mb in Oryza brachyantha (FF) to 905 Mb in Oryza australiensis (EE). The genome sizes of the tetraploid species were not always two times larger than those of the diploid species, and some diploid species even had larger genome sizes than those of tetraploids. Nevertheless, we found that genome sizes of newly formed allotetraploids (BBCC-) were almost equal to totaling genome sizes of their parental progenitors. Our results showed that the species belonging to the same genome types had similar genome sizes, while genome sizes exhibited a gradually decreased trend during the evolutionary process in the clade with AA, BB, CC, and EE genome types. Comparative genomic analyses further showed that the species with different rice genome types may had experienced dissimilar amplification histories of retrotransposons, resulting in remarkably different genome sizes. On the other hand, the closely related rice species may have experienced similar amplification history. We observed that the contents of transposable elements, long terminal repeats (LTR) retrotransposons, and particularly LTR/Gypsy retrotransposons varied largely but were significantly correlated with genome sizes. Therefore, this study demonstrated that LTR retrotransposons act as an active driver of genome size variation in the genus Oryza.

Towards understanding the biological foundations of perenniality

Perennial life cycles enable plants to have remarkably long lifespans, as exemplified by trees that can live for thousands of years. For this, they require sophisticated regulatory networks that sense environmental changes and initiate adaptive responses in their growth patterns. Recent research has gradually elucidated fundamental mechanisms underlying the perennial life cycle. Intriguingly, several conserved components of the floral transition pathway in annuals such as Arabidopsis thaliana also participate in these regulatory mechanisms underpinning perenniality. Here, we provide an overview of perennials' physiological features and summarise their recently discovered molecular foundations. We also highlight the importance of deepening our understanding of perenniality in the development of perennial grain crops, which are promising elements of future sustainable agriculture.

Reap the crop wild relatives for breeding future crops

Crop wild relatives (CWRs) have provided breeders with several 'game-changing' traits or genes that have boosted crop resilience and global agricultural production. Advances in breeding and genomics have accelerated the identification of valuable CWRs for use in crop improvement. The enhanced genetic diversity of breeding pools carrying optimum combinations of favorable alleles for targeted crop-growing regions is crucial to sustain genetic gain. In parallel, growing sequence information on wild genomes in combination with precise gene-editing tools provide a fast-track route to transform CWRs into ideal future crops. Data-informed germplasm collection and management strategies together with adequate policy support will be equally important to improve access to CWRs and their sustainable use to meet food and nutrition security targets.

A Chromosome-level Assembly of A Wild Castor Genome Provides New Insights into the Adaptive Evolution in A Tropical Desert

Wild castor grows in the high-altitude tropical desert of the African Plateau, a region known for high ultraviolet radiation, strong light, and extremely dry condition. To investigate the potential genetic basis of adaptation to both highland and tropical deserts, we generated a chromosome-level genome sequence assembly of the wild castor accession WT05, with a genome size of 316 Mb, a scaffold N50 of 31.93 Mb, and a contig N50 of 8.96 Mb, respectively. Compared with cultivated castor and other Euphorbiaceae species, the wild castor exhibits positive selection and gene family expansion for genes involved in DNA repair, photosynthesis, and abiotic stress responses. Genetic variations associated with positive selection were identified in several key genes, such as LIG1, DDB2, and RECG1, involved in nucleotide excision repair. Moreover, a study of genomic diversity among wild and cultivated accessions revealed genomic regions containing selection signatures associated with the adaptation to extreme environments. The identification of the genes and alleles with selection signatures provides insights into the genetic mechanisms underlying the adaptation of wild castor to the high-altitude tropical desert and would facilitate direct improvement of modern castor varieties.

Perennial Rye: Genetics of Perenniality and Limited Fertility

Perenniality, the ability of plants to regrow after seed set, could be introgressed into cultivated rye by crossing with the wild relative and perennial Secale strictum. However, studies in the past showed that Secale cereale × Secale strictum-derived cultivars were also characterized by reduced fertility what was related to so called chromosomal multivalents, bulks of chromosomes that paired together in metaphase I of pollen mother cells instead of only two chromosomes (bivalents). Those multivalents could be caused by ancient translocations that occurred between both species. Genetic studies on perennial rye are quite old and especially the advent of molecular markers and genome sequencing paved the way for new insights and more comprehensive studies. After a brief review of the past research, we used a basic QTL mapping approach to analyze the genetic status of perennial rye. We could show that for the trait perennation 0.74 of the genetic variance in our population was explained by additively inherited QTLs on chromosome 2R, 3R, 4R, 5R and 7R. Fertility on the other hand was with 0.64 of explained genetic variance mainly attributed to a locus on chromosome 5R, what was most probably the self-incompatibility locus S5. Additionally, we could trace the Z locus on chromosome 2R by high segregation distortion of markers. Indications for chromosomal co-segregation, like multivalents, could not be found. This study opens new possibilities to use perennial rye as genetic resource and for alternative breeding methods, as well as a valuable resource for comparative studies of perennation across different species.

Collection, preservation and distribution of Oryza genetic resources by the National Bioresource Project RICE (NBRP-RICE)

Biological resources are the basic infrastructure of bioscience research. Rice (Oryza sativa L.) is a good experimental model for research in cereal crops and monocots and includes important genetic materials used in breeding. The availability of genetic materials, including mutants, is important for rice research. In addition, Oryza species are attractive to researchers for both finding useful genes for breeding and for understanding the mechanism of genome evolution that enables wild plants to adapt to their own habitats. NBRP-RICE contributes to rice research by promoting the usage of genetic materials, especially wild Oryza accessions and mutant lines. Our activity includes collection, preservation and distribution of those materials and the provision of basic information on them, such as morphological and physiological traits and genomic information. In this review paper, we introduce the activities of NBRP-RICE and our database, Oryzabase, which facilitates the access to NBRP-RICE resources and their genomic sequences as well as the current situation of wild Oryza genome sequencing efforts by NBRP-RICE and other institutes.

Nitrogen Nutrition Promotes Rhizome Bud Outgrowth via Regulation of Cytokinin Biosynthesis Genes and an Oryza longistaminata Ortholog of FINE CULM 1

Oryza longistaminata, a wild rice, can propagate vegetatively via rhizome formation and, thereby, expand its territory through horizontal growth of branched rhizomes. The structural features of rhizomes are similar to those of aerial stems; however, the physiological roles of the two organs are different. Nitrogen nutrition is presumed to be linked to the vegetative propagation activity of rhizomes, but the regulation of rhizome growth in response to nitrogen nutrition and the underlying biological processes have not been well characterized. In this study, we analyzed rhizome axillary bud growth in response to nitrogen nutrition and examined the involvement of cytokinin-mediated regulation in the promotion of bud outgrowth in O. longistaminata. Our results showed that nitrogen nutrition sufficiency promoted rhizome bud outgrowth to form secondary rhizomes. In early stages of the response to nitrogen application, glutamine accumulated rapidly, two cytokinin biosynthesis genes, isopentenyltransferase, and CYP735A, were up-regulated with accompanying cytokinin accumulation, and expression of an ortholog of FINE CULM1, a negative regulator of axillary bud outgrowth, was severely repressed in rhizomes. These results suggest that, despite differences in physiological roles of these organs, the nitrogen-dependent outgrowth of rhizome axillary buds in O. longistaminata is regulated by a mechanism similar to that of shoot axillary buds in O. sativa. Our findings provide a clue for understanding how branched rhizome growth is regulated to enhance nutrient acquisition strategies.

Development of a genome-wide InDel marker set for allele discrimination between rice (Oryza sativa) and the other seven AA-genome Oryza species

Wild relatives of rice in the genus Oryza (composed of 24 species with 11 different genome types) have been significantly contributing to the varietal improvement of rice (Oryza sativa). More than 4000 accessions of wild rice species are available and they are regarded as a "genetic reservoir" for further rice improvement. DNA markers are essential tools in genetic analysis and breeding. To date, genome-wide marker sets for wild rice species have not been well established and this is one of the major difficulties for the efficient use of wild germplasm. Here, we developed 541 genome-wide InDel markers for the discrimination of alleles between the cultivated species O. sativa and the other seven AA-genome species by positional multiple sequence alignments among five AA-genome species with four rice varieties. The newly developed markers were tested by PCR-agarose gel analysis of 24 accessions from eight AA genome species (three accessions per species) along with two representative cultivars (O. sativa subsp. indica cv. IR24 and subsp. japonica cv. Nipponbare). Marker polymorphism was validated for 475 markers. The number of polymorphic markers between IR24 and each species (three accessions) ranged from 338 (versus O. rufipogon) to 416 (versus O. longistaminata) and the values in comparison with Nipponbare ranged from 179 (versus O. glaberrima) to 323 (versus O. glumaepatula). These marker sets will be useful for genetic studies and use of the AA-genome wild rice species.

Evolutionary Genomics of Structural Variation in Asian Rice (Oryza sativa) Domestication

Structural variants (SVs) are a largely unstudied feature of plant genome evolution, despite the fact that SVs contribute substantially to phenotypes. In this study, we discovered SVs across a population sample of 347 high-coverage, resequenced genomes of Asian rice (Oryza sativa) and its wild ancestor (O. rufipogon). In addition to this short-read data set, we also inferred SVs from whole-genome assemblies and long-read data. Comparisons among data sets revealed different features of genome variability. For example, genome alignment identified a large (∼4.3 Mb) inversion in indica rice varieties relative to japonica varieties, and long-read analyses suggest that ∼9% of genes from the outgroup (O. longistaminata) are hemizygous. We focused, however, on the resequencing sample to investigate the population genomics of SVs. Clustering analyses with SVs recapitulated the rice cultivar groups that were also inferred from SNPs. However, the site-frequency spectrum of each SV type-which included inversions, duplications, deletions, translocations, and mobile element insertions-was skewed toward lower frequency variants than synonymous SNPs, suggesting that SVs may be predominantly deleterious. Among transposable elements, SINE and mariner insertions were found at especially low frequency. We also used SVs to study domestication by contrasting between rice and O. rufipogon. Cultivated genomes contained ∼25% more derived SVs and mobile element insertions than O. rufipogon, indicating that SVs contribute to the cost of domestication in rice. Peaks of SV divergence were enriched for known domestication genes, but we also detected hundreds of genes gained and lost during domestication, some of which were enriched for traits of agronomic interest.

A chromosome-level genome assembly of the wild rice Oryza rufipogon facilitates tracing the origins of Asian cultivated rice

Oryza rufipogon Griff. is a wild progenitor of the Asian cultivated rice Oryza sativa. To better understand the genomic diversity of the wild rice, high-quality reference genomes of O. rufipogon populations are needed, which also facilitate utilization of the wild genetic resources in rice breeding. In this study, we generated a chromosome-level genome assembly of O. rufipogon using a combination of short-read sequencing, single-molecule sequencing, BioNano and Hi-C platforms. The genome sequence (399.8 Mb) was assembled into 46 scaffolds on the 12 chromosomes, with contig N50 and scaffold N50 of 13.2 Mb and 20.3 Mb, respectively. The genome contains 36,520 protein-coding genes, and 49.37% of the genome consists of repetitive elements. The genome has strong synteny with those of the O. sativa subspecies indica and japonica, but containing some large structural variations. Evolutionary analysis unveiled the polyphyletic origins of O. sativa, in which the japonica and indica genome formations involved different divergent O. rufipogon (including O. nivara) lineages, accompanied by introgression of genomic regions between japonica and indica. This high-quality reference genome provides insight on the genome evolution of the wild rice and the origins of the O. sativa subspecies, and valuable information for basic research and rice breeding.

Novel allelic variation in the Phospholipase D alpha1 gene (OsPLDα1) of wild Oryza species implies to its low expression in rice bran

Rice bran, a by-product after milling, is a rich source of phytonutrients like oryzanols, tocopherols, tocotrienols, phytosterols, and dietary fibers. Moreover, exceptional properties of the rice bran oil make it unparalleled to other vegetable oils. However, a lipolytic enzyme Phospholipase D alpha1 (OsPLDα1) causes rancidity and ‘stale flavor’ in the oil, and thus limits the rice bran usage for human consumption. To improve the rice bran quality, sequence based allele mining at OsPLDα1 locus (3.6 Kb) was performed across 48 accessions representing 11 wild Oryza species, 8 accessions of African cultivated rice, and 7 Oryza sativa cultivars. From comparative sequence analysis, 216 SNPs and 30 InDels were detected at the OsPLDα1 locus. Phylogenetic analysis revealed 20 OsPLDα1 cDNA variants which further translated into 12 protein variants. The O. officinalis protein variant, when compared to Nipponbare, showed maximum variability comprising 22 amino acid substitutions and absence of two peptides and two β-sheets. Further, expression profiling indicated significant differences in transcript abundance within as well as between the OsPLDα1 variants. Also, a new OsPLDα1 transcript variant having third exon missing in it, Os01t0172400-06, has been revealed. An O. officinalis accession (IRGC101152) had lowest gene expression which suggests the presence of novel allele, named as OsPLDα1-1a (GenBank accession no. MF966931). The identified novel allele could be further deployed in the breeding programs to overcome rice bran rancidity in elite cultivars.

SMRT sequencing of the Oryza rufipogon genome reveals the genomic basis of rice adaptation

Asian cultivated rice is believed to have been domesticated from a wild progenitor, Oryza rufipogon, offering promising sources of alleles for world rice improvement. Here we first present a high-quality chromosome-scale genome of the typical O. rufipogon. Comparative genomic analyses of O. sativa and its two wild progenitors, O. nivara and O. rufipogon, identified many dispensable genes functionally enriched in the reproductive process. We detected millions of genomic variants, of which large-effect mutations could affect agronomically relevant traits. We demonstrate how lineage-specific expansion of gene families may have contributed to the formation of reproduction isolation. We document thousands of genes with signatures of positive selection that are mainly involved in the reproduction and response to biotic- and abiotic stresses. We show that selection pressures may serve as forces to govern substantial genomic alterations that form the genetic basis of rapid evolution of mating and reproductive systems under diverse habitats.

Li et al. provide a chromosome-based genome of the ancestral rice relative Oryza rufipogon using SMRT, 10x and Hi-C. Using comparative genomics analysis, they identify genes important for agronomical traits and provide insights for the evolution of reproductive systems which can facilitate rice breeding programs.

The Mitochondrial Iron-Regulated (MIR) gene is Oryza genus specific and evolved before speciation within the Oryza sativa complex

The MIR gene is not an Oryza sativa orphan gene, but an Oryza genus-specific gene that evolved before AA lineage speciation by a complex origination process. Rice (Oryza sativa L.) is a model species and an economically relevant crop. The Oryza genus comprises 25 species, with genomic data available for several Oryza species, making it a model for genetics and evolution. The Mitochondrial Iron-Regulated (MIR) gene was previously implicated in the O. sativa Fe deficiency response, and was considered an orphan gene present only in rice. Here we show that MIR is also found in other Oryza species that belong to the Oryza sativa complex, which have AA genome type and constitute the primary gene pool for O. sativa breeding. Our data suggest that MIR originated in a stepwise process, in which sequences derived from an exon fragment of the raffinose synthase gene were pseudogenized into non-coding, which in turn originated the MIR gene de novo. All species with a putative functional MIR gene conserve their regulation by Fe deficiency, with the exception of Oryza barthii. In O. barthii, the MIR coding sequence was translocated to a different chromosomal position and separated from its regulatory region, leading to a lack of Fe deficiency responsiveness. Moreover, the MIR co-expression subnetwork cluster in O. sativa is responsive to Fe deficiency, evidencing the importance of the newly originated gene in Fe uptake. This work establishes that MIR is not an orphan gene as previously proposed, but a de novo originated gene within the genus Oryza. We also showed that MIR is undergoing genomic changes in one species (O. barthii), with an impact on Fe deficiency response.

Suppression of Leaf Blade Development by BLADE-ON-PETIOLE Orthologs Is a Common Strategy for Underground Rhizome Growth

Rhizomes are modified stems that grow horizontally underground in various perennial species, a growth habit that is advantageous for vigorous asexual proliferation. In Oryza longistaminata, a rhizomatous wild relative of cultivated rice (Oryza sativa), leaves in the aerial shoots consist of a distal leaf blade and a proximal leaf sheath [1]. Leaf blade formation is, however, suppressed in rhizome leaves. In O. sativa, BLADE-ON-PETIOLE (BOP) genes are the main regulators of proximal-distal leaf patterning [2]. During the juvenile phase of O. sativa, BOP expression is maintained at high levels by the small regulatory RNA microRNA156 (miR156), leading to formation of leaves consisting predominantly of the sheath. Here, we show that in O. longistaminata, high expression of BOPs caused by miR156 was responsible for suppression of the blade in rhizomes and that bop loss-of-function mutants produced leaves consisting of the leaf blade only. Rhizome growth in soil was also hampered in the mutants due to a severe reduction in rhizome tip stiffness. Leaf blade formation is also suppressed in the stolons of Zoysia matrella, a monocot species, and in the rhizomes of Houttuynia cordata, a dicot species, indicating that leaf blade suppression is widely conserved. We also show that strong expression of BOP homologs in both rhizome and stolon leaves rather than in aerial leaves is another conserved feature. We propose that suppression of the leaf blade by BOP is an evolutionary strategy that has been commonly recruited by both rhizomatous and stoloniferous species to establish their unique growth habit.

Identification of Anther Length QTL and Construction of Chromosome Segment Substitution Lines of Oryza longistaminata

Life histories and breeding systems strongly affect the genetic diversity of seed plants, but the genetic architectures that promote outcrossing in Oryza longistaminata, a perennial wild species in Africa, are not understood. We conducted a genetic analysis of the anther length of O. longistaminata accession W1508 using advanced backcross quantitative trait locus (QTL) analysis and chromosomal segment substitution lines (CSSLs) in the genetic background of O. sativa Taichung 65 (T65), with simple sequence repeat markers. QTL analysis of the BC₃F₁ population (n = 100) revealed that four main QTL regions on chromosomes 3, 5, and 6 were associated to anther length. We selected a minimum set of BC₃F₂ plants for the development of CSSLs to cover as much of the W1508 genome as possible. The additional minor QTLs were suggested in the regional QTL analysis, using 21 to 24 plants in each of the selected BC₃F₂ population. The main QTLs found on chromosomes 3, 5, and 6 were validated and designated qATL3, qATL5, qATL6.1, and qATL6.2, as novel QTLs identified in O. longistaminata in the mapping populations of 94, 88, 70, and 95 BC₃F₄ plants. qATL3, qATL5, and qATL6.1 likely contributed to anther length by cell elongation, whereas qATL6.2 likely contributed by cell multiplication. The QTLs were confirmed again in an evaluation of the W1508ILs. In several chromosome segment substitution lines without the four validated QTLs, the anthers were also longer than those of T65, suggesting that other QTLs also increase anther length in W1508. The cloning and diversity analyses of genes conferring anther length QTLs promotes utilization of the genetic resources of wild species, and the understanding of haplotype evolution on the differentiation of annuality and perenniality in the genus Oryza.

Development of introgression lines of AA genome Oryza species, O. glaberrima, O. rufipogon, and O. nivara, in the genetic background of O. sativa L. cv. Taichung 65

To evaluate and utilize potentially valuable quantitative trait loci or genes of wild relatives in the genetic background of domesticated crop species, chromosome segment substitution lines (CSSLs) are a valuable tool. CSSLs can be constructed through the exchange of chromosome segments of AA genome species of the genus Oryza with cultivated rice, Oryza sativa L. Here we report the development of three sets of CSSLs carrying segments of AA genome species closely related to Oryza sativa-O. glaberrima (IRGC 103777 from Mali), O. rufipogon (W1962 from China), and O. nivara (IRGC 105715 from Cambodia)-in the genetic background of ssp. japonica cultivar Taichung 65 through the use of 101 to 121 simple-sequence-repeat markers in whole-genome genotyping and marker-assisted selection. The materials are available via the National Bioresource Project (Rice) Oryzabase Web page.

Fine-Tuning the Expression of Duplicate Genes by Translational Regulation in Arabidopsis and Maize

Plant genomes are extensively shaped by various types of gene duplication. However, in this active area of investigation, the vast majority of studies focus on the sequence and transcription of duplicate genes, leaving open the question of how translational regulation impacts the expression and evolution of duplicate genes. We explored this issue by analyzing the ribo- and mRNA-seq data sets across six tissue types and stress conditions in Arabidopsis thaliana and maize (Zea mays). We dissected the relative contributions of transcriptional and translational regulation to the divergence in the abundance of ribosome footprint (RF) for different types of duplicate genes. We found that the divergence in RF abundance was largely programmed at the transcription level and that translational regulation plays more of a modulatory role. Intriguingly, translational regulation is characterized by its strong directionality, with the divergence in translational efficiency (TE) globally counteracting the divergence in mRNA abundance, indicating partial buffering of the transcriptional divergence between paralogs by translational regulation. Divergence in TE was associated with several sequence features. The faster-evolving copy in a duplicate pair was more likely to show lower RF abundance, which possibly results from relaxed purifying selection compared with its paralog. A considerable proportion of duplicates displayed differential TE across tissue types and stress conditions, most of which were enriched in photosynthesis, energy production, and translation-related processes. Additionally, we constructed a database TDPDG-DB (http://www.plantdupribo.tk), providing an online platform for data exploration. Overall, our study illustrates the roles of translational regulation in fine-tuning duplicate gene expression in plants.