Genomic introgressions from African rice (Oryza glaberrima) in Asian rice (O. sativa) lead to the identification of key QTLs for panicle architecture.

BACKGROUND: Developing high yielding varieties is a major challenge for breeders tackling the challenges of climate change in agriculture. The panicle (inflorescence) architecture of rice is one of the key components of yield potential and displays high inter- and intra-specific variability. The genus Oryza features two different crop species: Asian rice (Oryza sativa L.) and the African rice (O. glaberrima Steud.). One of the main morphological differences between the two independently domesticated species is the structure (or complexity) of the panicle, with O. sativa displaying a highly branched panicle, which in turn produces a larger number of grains than that of O. glaberrima. The gene regulatory network that governs intra- and interspecific panicle diversity is still under-studied. RESULTS: To identify genetic factors linked to panicle architecture diversity in the two species, we used a set of 60 Chromosome Segment Substitution Lines (CSSLs) issued from third generation backcross (BC3DH) and carrying genomic segments from O. glaberrima cv. MG12 in the genetic background of O. sativa Tropical Japonica cv. Caiapó. Phenotypic data were collected for rachis and primary branch length, primary, secondary and tertiary branch number and spikelet number. A total of 15 QTLs were localized on chromosomes 1, 2, 3, 7, 11 and 12, QTLs associated with enhanced secondary and tertiary branch numbers were detected in two CSSLs. Furthermore, BC4F3:5 lines carrying different combinations of substituted segments were produced to decipher the effects of the identified QTL regions on variations in panicle architecture. A detailed analysis of phenotypes versus genotypes was carried out between the two parental genomes within these regions in order to understand how O. glaberrima introgression events may lead to alterations in panicle traits. CONCLUSION: Our analysis led to the detection of genomic variations between O. sativa cv. Caiapó and O. glaberrima cv. MG12 in regions associated with enhanced panicle traits in specific CSSLs. These regions contain a number of key genes that regulate panicle development in O. sativa and their interspecific genomic variations may explain the phenotypic effects observed.

Genetic bases of resistance to the rice hoja blanca disease deciphered by a quantitative trait locus approach.

Rice hoja blanca (RHB) is one of the most serious diseases in rice-growing areas in tropical Americas. Its causal agent is RHB virus (RHBV), transmitted by the planthopper Tagosodes orizicolus Müir. Genetic resistance is the most effective and environment-friendly way of controlling the disease. So far, only 1 major quantitative trait locus (QTL) of Oryza sativa ssp. japonica origin, qHBV4.1, that alters the incidence of the virus symptoms in 2 Colombian cultivars has been reported. This resistance has already started to be broken, stressing the urgent need for diversifying the resistance sources. In the present study, we performed a search for new QTLs of O. sativa indica origin associated with RHB resistance. We used 4 F2:3-segregating populations derived from indica-resistant varieties crossed with a highly susceptible japonica pivot parent. Besides the standard method for measuring disease incidence, we developed a new method based on computer-assisted image processing to determine the affected leaf area (ALA) as a measure of symptom severity. Based on the disease severity and incidence scores in the F3 families under greenhouse conditions and SNP genotyping of the F2 individuals, we identified 4 new indica QTLs for RHB resistance on rice chromosomes 4, 6, and 11, namely, qHBV4.2WAS208, qHBV6.1PTB25, qHBV11.1, and qHBV11.2, respectively. We also confirmed the wide-range action of qHBV4.1. Among the 5 QTLs, qHBV4.1 and qHBV11.1 had the largest effects on incidence and severity, respectively. These results provide a more complete understanding of the genetic bases of RHBV resistance in the cultivated rice gene pool and can be used to develop marker-aided breeding strategies to improve RHB resistance. The power of joint- and meta-analyses allowed precise mapping and candidate gene identification, providing the basis for positional cloning of the 2 major QTLs qHBV4.1 and qHBV11.1.

Genome-wide association analysis identifies natural allelic variants associated with panicle architecture variation in African rice, Oryza glaberrima Steud.

African rice (Oryza glaberrima Steud), a short-day cereal crop closely related to Asian rice (Oryza sativa L.), has been cultivated in Sub-Saharan Africa for ∼ 3,000 years. Although less cultivated globally, it is a valuable genetic resource in creating high-yielding cultivars that are better adapted to diverse biotic and abiotic stresses. While inflorescence architecture, a key trait for rice grain yield improvement, has been extensively studied in Asian rice, the morphological and genetic determinants of this complex trait are less understood in African rice. In this study, using a previously developed association panel of 162 O. glaberrima accessions and new SNP variants characterized through mapping to a new version of the O. glaberrima reference genome, we conducted a genome-wide association study of four major morphological panicle traits. We have found a total of 41 stable genomic regions that are significantly associated with these traits, of which 13 co-localized with previously identified QTLs in O. sativa populations and 28 were unique for this association panel. Additionally, we found a genomic region of interest on chromosome 3 that was associated with the number of spikelets and primary and secondary branches. Within this region was localized the O. sativa ortholog of the PHYTOCHROME B gene (Oglab_006903/OgPHYB). Haplotype analysis revealed the occurrence of natural sequence variants at the OgPHYB locus associated with panicle architecture variation through modulation of the flowering time phenotype, whereas no equivalent alleles were found in O. sativa. The identification in this study of genomic regions specific to O. glaberrima indicates panicle-related intra-specific genetic variation in this species, increasing our understanding of the underlying molecular processes governing panicle architecture. Identified candidate genes and major haplotypes may facilitate the breeding of new African rice cultivars with preferred panicle traits.

Pan-genome inversion index reveals evolutionary insights into the subpopulation structure of Asian rice.

Understanding and exploiting genetic diversity is a key factor for the productive and stable production of rice. Here, we utilize 73 high-quality genomes that encompass the subpopulation structure of Asian rice (Oryza sativa), plus the genomes of two wild relatives (O. rufipogon and O. punctata), to build a pan-genome inversion index of 1769 non-redundant inversions that span an average of ~29% of the O. sativa cv. Nipponbare reference genome sequence. Using this index, we estimate an inversion rate of ~700 inversions per million years in Asian rice, which is 16 to 50 times higher than previously estimated for plants. Detailed analyses of these inversions show evidence of their effects on gene expression, recombination rate, and linkage disequilibrium. Our study uncovers the prevalence and scale of large inversions (≥100 bp) across the pan-genome of Asian rice and hints at their largely unexplored role in functional biology and crop performance.

Prediction of crossover recombination using parental genomes.

Meiotic recombination is a crucial cellular process, being one of the major drivers of evolution and adaptation of species. In plant breeding, crossing is used to introduce genetic variation among individuals and populations. While different approaches to predict recombination rates for different species have been developed, they fail to estimate the outcome of crossings between two specific accessions. This paper builds on the hypothesis that chromosomal recombination correlates positively to a measure of sequence identity. It presents a model that uses sequence identity, combined with other features derived from a genome alignment (including the number of variants, inversions, absent bases, and CentO sequences) to predict local chromosomal recombination in rice. Model performance is validated in an inter-subspecific indica x japonica cross, using 212 recombinant inbred lines. Across chromosomes, an average correlation of about 0.8 between experimental and prediction rates is achieved. The proposed model, a characterization of the variation of the recombination rates along the chromosomes, can enable breeding programs to increase the chances of creating novel allele combinations and, more generally, to introduce new varieties with a collection of desirable traits. It can be part of a modern panel of tools that breeders can use to reduce costs and execution times of crossing experiments.

Methylation in the CHH Context Allows to Predict Recombination in Rice.

DNA methylation is the most studied epigenetic trait. It is considered a key factor in regulating plant development and physiology, and has been associated with the regulation of several genomic features, including transposon silencing, regulation of gene expression, and recombination rates. Nonetheless, understanding the relation between DNA methylation and recombination rates remains a challenge. This work explores the association between recombination rates and DNA methylation for two commercial rice varieties. The results show negative correlations between recombination rates and methylated cytosine counts for all contexts tested at the same time, and for CG and CHG contexts independently. In contrast, a positive correlation between recombination rates and methylated cytosine count is reported in CHH contexts. Similar behavior is observed when considering only methylated cytosines within genes, transposons, and retrotransposons. Moreover, it is shown that the centromere region strongly affects the relationship between recombination rates and methylation. Finally, machine learning regression models are applied to predict recombination using the count of methylated cytosines in the CHH context as the entrance feature. These findings shed light on the understanding of the recombination landscape of rice and represent a reference framework for future studies in rice breeding, genetics, and epigenetics.

The ÓMICAS alliance, an international research program on multi-omics for crop breeding optimization.

The OMICAS alliance is part of the Colombian government's Scientific Ecosystem, established between 2017-2018 to promote world-class research, technological advancement and improved competency of higher education across the nation. Since the program's kick-off, OMICAS has focused on consolidating and validating a multi-scale, multi-institutional, multi-disciplinary strategy and infrastructure to advance discoveries in plant science and the development of new technological solutions for improving agricultural productivity and sustainability. The strategy and methods described in this article, involve the characterization of different crop models, using high-throughput, real-time phenotyping technologies as well as experimental tissue characterization at different levels of the omics hierarchy and under contrasting conditions, to elucidate epigenome-, genome-, proteome- and metabolome-phenome relationships. The massive data sets are used to derive in-silico models, methods and tools to discover complex underlying structure-function associations, which are then carried over to the production of new germplasm with improved agricultural traits. Here, we describe OMICAS' R&D trans-disciplinary multi-project architecture, explain the overall strategy and methods for crop-breeding, recent progress and results, and the overarching challenges that lay ahead in the field.

Insights Into the Genetics of the Zhonghua 11 Resistance to Meloidogyne graminicola and Its Molecular Determinism in Rice.

Meloidogyne graminicola is a widely spread nematode pest of rice that reduces crop yield up to 20% on average in Asia, with devastating consequences for local and global rice production. Due to the ban on many chemical nematicides and the recent changes in water management practices in rice agriculture, an even greater impact of M. graminicola can be expected in the future, stressing the demand for the development of new sustainable nematode management solutions. Recently, a source of resistance to M. graminicola was identified in the Oryza sativa japonica rice variety Zhonghua 11 (Zh11). In the present study, we examine the genetics of the Zh11 resistance to M. graminicola and provide new insights into its cellular and molecular mechanisms. The segregation of the resistance in F2 hybrid populations indicated that two dominant genes may be contributing to the resistance. The incompatible interaction of M. graminicola in Zh11 was distinguished by a lack of swelling of the root tips normally observed in compatible interactions. At the cellular level, the incompatible interaction was characterised by a rapid accumulation of reactive oxygen species in the vicinity of the nematodes, accompanied by extensive necrosis of neighbouring cells. The expression profiles of several genes involved in plant immunity were analysed at the early stages of infection during compatible (susceptible plant) and incompatible (resistant plant) interactions. Notably, the expression of OsAtg4 and OsAtg7, significantly increased in roots of resistant plants in parallel with the cell death response, suggesting that autophagy is activated and may contribute to the resistance-mediated hypersensitive response. Similarly, transcriptional regulation of genes involved in hormonal pathways in Zh11 indicated that salicylate signalling may be important in the resistance response towards M. graminicola. Finally, the nature of the resistance to M. graminicola and the potential exploitation of the Zh11 resistance for breeding are discussed.

Corrigendum: The ÓMICAS alliance, an international research program on multi-omics for crop breeding optimization.

[This corrects the article DOI: 10.3389/fpls.2022.992663.].

A platinum standard pan-genome resource that represents the population structure of Asian rice.

As the human population grows from 7.8 billion to 10 billion over the next 30 years, breeders must do everything possible to create crops that are highly productive and nutritious, while simultaneously having less of an environmental footprint. Rice will play a critical role in meeting this demand and thus, knowledge of the full repertoire of genetic diversity that exists in germplasm banks across the globe is required. To meet this demand, we describe the generation, validation and preliminary analyses of transposable element and long-range structural variation content of 12 near-gap-free reference genome sequences (RefSeqs) from representatives of 12 of 15 subpopulations of cultivated Asian rice. When combined with 4 existing RefSeqs, that represent the 3 remaining rice subpopulations and the largest admixed population, this collection of 16 Platinum Standard RefSeqs (PSRefSeq) can be used as a template to map resequencing data to detect virtually all standing natural variation that exists in the pan-genome of cultivated Asian rice.

A set of AP2-like genes is associated with inflorescence branching and architecture in domesticated rice.

Rice yield is influenced by inflorescence size and architecture, and inflorescences from domesticated rice accessions produce more branches and grains. Neither the molecular control of branching nor the developmental differences between wild and domesticated rice accessions are fully understood. We surveyed phenotypes related to branching, size, and grain yield across 91 wild and domesticated African and Asian accessions. Characteristics related to axillary meristem identity were the main phenotypic differences between inflorescences from wild and domesticated accessions. We used whole transcriptome sequencing in developing inflorescences to measure gene expression before and after the transition from branching axillary meristems to determinate spikelet meristems. We identified a core set of genes associated with axillary meristem identity in Asian and African rice, and another set associated with phenotypic variability between wild and domesticated accessions. AP2/EREBP-like genes were enriched in both sets, suggesting that they are key factors in inflorescence branching and rice domestication. Our work has identified new candidates in the molecular control of inflorescence development and grain yield, and provides a detailed description of the effects of domestication on phenotype and gene expression.

A Hypersensitivity-Like Response to Meloidogyne graminicola in Rice (Oryza sativa).

Meloidogyne graminicola is a major plant-parasitic nematode affecting rice cultivation in Asia. Resistance to this nematode was found in the African rice genotypes Oryza glaberrima and O. longistaminata; however, due to interspecific hybrid sterility, the introgression of resistance genes in the widely consumed O. sativa varieties remains challenging. Recently, resistance was found in O. sativa and, here, we report for the first time the histological and genetic characterization of the resistance to M. graminicola in Zhonghua 11, an O. sativa variety. Bright-light microscopy and fluorescence observations of the root tissue of this variety revealed that the root cells surrounding the nematode displayed a hypersensitivity-like reaction with necrotic cells at early stages of infection when nematodes are migrating in the root's mesoderm. An accumulation of presumably phenolic compounds in the nematodes' neighboring root cells was also observed. In addition, at a later stage of infection, not only were few feeding sites observed but also the giant cells were underdeveloped, underlining an incompatible interaction. Furthermore, we generated a hybrid O. sativa population by crossing Zhonghua 11 with the susceptible O. sativa variety IR64 in order to describe the genetic background of this resistance. Our data suggested that the resistance to M. graminicola infection was qualitative rather than quantitative and, therefore, major resistance genes must be involved in this infection process. The full characterization of the defense mechanism and the preliminary study of the genetic inheritance of novel sources of resistance to Meloidogyne spp. in rice constitute a major step toward their use in crop breeding.

Constructing linkage maps in the genomics era with MapDisto 2.0.

MOTIVATION: Genotyping by sequencing (GBS) generates datasets that are challenging to handle by current genetic mapping software with graphical interface. Geneticists need new user-friendly computer programs that can analyze GBS data on desktop computers. This requires improvements in computation efficiency, both in terms of speed and use of random-access memory (RAM). RESULTS: MapDisto v.2.0 is a user-friendly computer program for construction of genetic linkage maps. It includes several new major features: (i) handling of very large genotyping datasets like the ones generated by GBS; (ii) direct importation and conversion of Variant Call Format (VCF) files; (iii) detection of linkage, i.e. construction of linkage groups in case of segregation distortion; (iv) data imputation on VCF files using a new approach, called LB-Impute. Features i to iv operate through inclusion of new Java modules that are used transparently by MapDisto; (v) QTL detection via a new R/qtl graphical interface. AVAILABILITY AND IMPLEMENTATION: The program is available free of charge at mapdisto.free.fr. CONTACT: mapdisto@gmail.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genetic Architecture of a Rice Nested Association Mapping Population.

Describing the genetic diversity in the gene pool of crops will provide breeders with novel resources for varietal improvement. Nested Association Mapping (NAM) populations are uniquely suited for characterizing parental diversity through the shuffling and fixation of parental haplotypes. Here, we describe a set of 1879 rice NAM lines created through the selfing and single-seed descent of F1 hybrids derived from elite IR64 indica crossed with 10 diverse tropical japonica lines. Genotyping data indicated tropical japonica alleles were captured at every queried locus despite the presence of segregation distortion factors. Several distortion loci were mapped, both shared and unique, among the 10 populations. Using two-point and multi-point genetic map calculations, our datasets achieved the ∼1500 cM expected map size in rice. Finally, we highlighted the utility of the NAM lines for QTL mapping, including joint analysis across the 10 populations, by confirming known QTL locations for the trait days to heading.

Identification of novel major and minor QTLs associated with Xanthomonas oryzae pv. oryzae (African strains) resistance in rice (Oryza sativa L.).

BACKGROUND: Xanthomonas oryzae pv. oryzae (Xoo) is the causal agent of Bacterial Leaf Blight (BB), an emerging disease in rice in West-Africa which can induce up to 50 % of yield losses. So far, no specific resistance gene or QTL to African Xoo were mapped. The objectives of this study were to identify and map novels and specific resistance QTLs to African Xoo strains. RESULTS: The reference recombinant inbred lines (RIL) mapping population derived from the cross between IR64 and Azucena was used to investigate Xoo resistance. Resistance to African and Philippine Xoo strains representing different races was assessed on the RIL population under greenhouse conditions. Five major quantitative trait loci (QTL) for resistance against African Xoo were located on different chromosomes. Loci on chromosomes 1, 7, 9, 10 and 11 explained as much as 13 %, 37 %, 13 %, 11 % and 15 % of resistance variation, respectively. A major novel QTL located on chromosome 7 explained 37 % of the phenotypic variance to the African Xoo corresponding to race A3 whereas that on chromosome 11 is effective to all African races tested. Together with genes and QTLs for resistance to bacterial blight previously described, the QTLs described here were mapped onto the reference O. sativa subs japonica (var. Nipponbare) physical map. CONCLUSION: We characterized new resistance QTLs. While some co-localize with known resistance genes/QTLs to Asian strains, others are specific to African strains. We result with new information on genes and QTLs for resistance to bacterial blight that will be useful for controlling the disease.

Whole genome sequencing of elite rice cultivars as a comprehensive information resource for marker assisted selection.

Current advances in sequencing technologies and bioinformatics revealed the genomic background of rice, a staple food for the poor people, and provided the basis to develop large genomic variation databases for thousands of cultivars. Proper analysis of this massive resource is expected to give novel insights into the structure, function, and evolution of the rice genome, and to aid the development of rice varieties through marker assisted selection or genomic selection. In this work we present sequencing and bioinformatics analyses of 104 rice varieties belonging to the major subspecies of Oryza sativa. We identified repetitive elements and recurrent copy number variation covering about 200 Mbp of the rice genome. Genotyping of over 18 million polymorphic locations within O. sativa allowed us to reconstruct the individual haplotype patterns shaping the genomic background of elite varieties used by farmers throughout the Americas. Based on a reconstruction of the alleles for the gene GBSSI, we could identify novel genetic markers for selection of varieties with high amylose content. We expect that both the analysis methods and the genomic information described here would be of great use for the rice research community and for other groups carrying on similar sequencing efforts in other crops.

Development and GBS-genotyping of introgression lines (ILs) using two wild species of rice, O. meridionalis and O. rufipogon, in a common recurrent parent, O. sativa cv. Curinga.

Two populations of interspecific introgression lines (ILs) in a common recurrent parent were developed for use in pre-breeding and QTL mapping. The ILs were derived from crosses between cv Curinga, a tropical japonica upland cultivar, and two different wild donors, Oryza meridionalis Ng. accession (W2112) and Oryza rufipogon Griff. accession (IRGC 105491). The lines were genotyped using genotyping-by-sequencing (GBS) and SSRs. The 32 Curinga/O. meridionalis ILs contain 76.73 % of the donor genome in individual introgressed segments, and each line has an average of 94.9 % recurrent parent genome. The 48 Curinga/O. rufipogon ILs collectively contain 97.6 % of the donor genome with an average of 89.9 % recurrent parent genome per line. To confirm that these populations were segregating for traits of interest, they were phenotyped for pericarp color in the greenhouse and for four agronomic traits-days to flowering, plant height, number of tillers, and number of panicles-in an upland field environment. Seeds from these IL libraries and the accompanying GBS datasets are publicly available and represent valuable genetic resources for exploring the genetics and breeding potential of rice wild relatives.

An extensive analysis of the African rice genetic diversity through a global genotyping.

KEY MESSAGE: We present here the first curated collection of wild and cultivated African rice species. For that, we designed specific SNPs and were able to structure these very low diverse species. Oryza glaberrima, the cultivated African rice, is endemic from Africa. This species and its direct ancestor, O. barthii, are valuable tool for improvement of Asian rice O. sativa in terms of abiotic and biotic stress resistance. However, only a few limited studies about the genetic diversity of these species were performed. In the present paper, and for the first time at such extend, we genotyped 279 O. glaberrima, selected both for their impact in current breeding and for their geographical distribution, and 101 O. barthii, chosen based on their geographic origin, using a set of 235 SNPs specifically designed for African rice diversity. Using those data, we were able to structure the individuals from our sample in three populations for O. barthii, related to geography, and two populations in O. glaberrima; these two last populations cannot be linked however to any currently phenotyped trait. Moreover, we were also able to identify misclassification in O. glaberrima as well as in O. barthii and identified new form of O. sativa from the set of African varieties.

Major QTLs control resistance to rice hoja blanca virus and its vector Tagosodes orizicolus.

Rice hoja blanca (white leaf) disease can cause severe yield losses in rice in the Americas. The disease is caused by the rice hoja blanca virus (RHBV), which is transmitted by the planthopper vector Tagosodes orizicolus. Because classical breeding schemes for this disease rely on expensive, time-consuming screenings, there is a need for alternatives such as marker-aided selection. The varieties Fedearroz 2000 and Fedearroz 50, which are resistant to RHBV and to the feeding damage caused by T. orizicolus, were crossed with the susceptible line WC366 to produce segregating F2:3 populations. The F3 families were scored for their resistance level to RHBV and T. orizicolus. The F2:3 lines of both crosses were genotyped using microsatellite markers. One major QTL on the short arm of chromosome 4 was identified for resistance to RHBV in the two populations. Two major QTL on chromosomes 5 and 7 were identified for resistance to T. orizicolus in the Fd2000 × WC366 and Fd50 × WC366 crosses, respectively. This comparative study using two distinct rice populations allowed for a better understanding of how the resistance to RHBV and its vector are controlled genetically. Simple marker-aided breeding schemes based on QTL information can be designed to improve rice germplasm to reduce losses caused by this important disease.