Genome Wide Association Study Pinpoints Key Agronomic QTLs in African Rice Oryza glaberrima

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.

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.

Chlorophyll fluorescence-based high-throughput phenotyping facilitates the genetic dissection of photosynthetic heat tolerance in African (Oryza glaberrima) and Asian (Oryza sativa) rice

Rising temperatures and extreme heat events threaten rice production. Half of the global population relies on rice for basic nutrition, and therefore developing heat-tolerant rice is essential. During vegetative development, reduced photosynthetic rates can limit growth and the capacity to store soluble carbohydrates. The photosystem II (PSII) complex is a particularly heat-labile component of photosynthesis. We have developed a high-throughput chlorophyll fluorescence-based screen for photosynthetic heat tolerance capable of screening hundreds of plants daily. Through measuring the response of maximum PSII efficiency to increasing temperature, this platform generates data for modelling the PSII-temperature relationship in large populations in a small amount of time. Coefficients from these models (photosynthetic heat tolerance traits) demonstrated high heritabilities across African (Oryza glaberrima) and Asian (Oryza sativa, Bengal Assam Aus Panel) rice diversity sets, highlighting valuable genetic variation accessible for breeding. Genome-wide association studies were performed across both species for these traits, representing the first documented attempt to characterize the genetic basis of photosynthetic heat tolerance in any species to date. A total of 133 candidate genes were highlighted. These were significantly enriched with genes whose predicted roles suggested influence on PSII activity and the response to stress. We discuss the most promising candidates for improving photosynthetic heat tolerance in rice.

GWAS analysis to elucidate genetic composition underlying a photoperiod-insensitive rice population, North Korea

Heading date (Hd) is one of the main factors determining rice production and regional adaptation. To identify the genetic factors involved in the wide regional adaptability of rice, we conducted a genome-wide association study (GWAS) with 190 North Korean rice accessions selected for non-precocious flowering in the Philippines, a low-latitude region. Using both linear mixed models (LMM) and fixed and random model circulating probability unification (FarmCPU), we identified five significant loci for Hd in trials in 2018 and 2019. Among the five lead single nucleotide polymorphisms (SNPs), three were located adjacent to the known Hd genes, Heading date 3a (Hd3a), Heading date 5 (Hd5), and GF14-c. In contrast, three SNPs were located in novel loci with minor effects on heading. Further GWAS analysis for photoperiod insensitivity (PS) revealed no significant genes associated with PS, supporting that this North Korean (NK) population is largely photoperiod-insensitive. Haplotyping analysis showed that more than 80% of the NK varieties harbored nonfunctional alleles of major Hd genes investigated, of which a nonfunctional allele of Heading date 1 (Hd1) was observed in 66% of the varieties. Geographical distribution analysis of Hd allele combination types showed that nonfunctional alleles of floral repressor Hd genes enabled rice cultivation in high-latitude regions. In contrast, Hd1 alleles largely contributed to the wide regional adaptation of rice varieties. In conclusion, an allelic combination of Hd genes is critical for rice cultivation across wide areas.

Physiological and genetic control of transpiration efficiency in African rice, Oryza glaberrima Steud

Improving crop water use efficiency, the amount of carbon assimilated as biomass per unit of water used by a plant, is of major importance as water for agriculture becomes scarcer. In rice, the genetic bases of transpiration efficiency, the derivation of water use efficiency at the whole-plant scale, and its putative component trait transpiration restriction under high evaporative demand remain unknown. These traits were measured in 2019 in a panel of 147 African rice (Oryza glaberrima) genotypes known to be potential sources of tolerance genes to biotic and abiotic stresses. Our results reveal that higher transpiration efficiency is associated with transpiration restriction in African rice. Detailed measurements in a subset of highly contrasted genotypes in terms of biomass accumulation and transpiration confirmed these associations and suggested that root to shoot ratio played an important role in transpiration restriction. Genome wide association studies identified marker-trait associations for transpiration response to evaporative demand, transpiration efficiency, and its residuals, with links to genes involved in water transport and cell wall patterning. Our data suggest that root-shoot partitioning is an important component of transpiration restriction that has a positive effect on transpiration efficiency in African rice. Both traits are heritable and define targets for breeding rice with improved water use strategies.

Out of Africa: characterizing the natural variation in dynamic photosynthetic traits in a diverse population of African rice (Oryza glaberrima)

African rice (Oryza glaberrima) has adapted to challenging environments and is a promising source of genetic variation. We analysed dynamics of photosynthesis and morphology in a reference set of 155 O. glaberrima accessions. Plants were grown in an agronomy glasshouse to late tillering stage. Photosynthesis induction from darkness and the decrease in low light was measured by gas exchange and chlorophyll fluorescence along with root and shoot biomass, stomatal density, and leaf area. Steady-state and kinetic responses were modelled. We describe extensive natural variation in O. glaberrima for steady-state, induction, and reduction responses of photosynthesis that has value for gene discovery and crop improvement. Principal component analyses indicated key clusters of plant biomass, kinetics of photosynthesis (CO2 assimilation, A), and photoprotection induction and reduction (measured by non-photochemical quenching, NPQ), consistent with diverse adaptation. Accessions also clustered according to countries with differing water availability, stomatal conductance (gs), A, and NPQ, indicating that dynamic photosynthesis has adaptive value in O. glaberrima. Kinetics of NPQ, A, and gs showed high correlation with biomass and leaf area. We conclude that dynamic photosynthetic traits and NPQ are important within O. glaberrima, and we highlight NPQ kinetics and NPQ under low light.

Out of Africa: characterizing the natural variation in dynamic photosynthetic traits in a diverse population of African rice (Oryza glaberrima)

African rice (Oryza glaberrima) has adapted to challenging environments and is a promising source of genetic variation. We analysed dynamics of photosynthesis and morphology in a reference set of 155 O. glaberrima accessions. Plants were grown in an agronomy glasshouse to late tillering stage. Photosynthesis induction from darkness and the decrease in low light was measured by gas exchange and chlorophyll fluorescence along with root and shoot biomass, stomatal density, and leaf area. Steady-state and kinetic responses were modelled. We describe extensive natural variation in O. glaberrima for steady-state, induction, and reduction responses of photosynthesis that has value for gene discovery and crop improvement. Principal component analyses indicated key clusters of plant biomass, kinetics of photosynthesis (CO2 assimilation, A), and photoprotection induction and reduction (measured by non-photochemical quenching, NPQ), consistent with diverse adaptation. Accessions also clustered according to countries with differing water availability, stomatal conductance (gs), A, and NPQ, indicating that dynamic photosynthesis has adaptive value in O. glaberrima. Kinetics of NPQ, A, and gs showed high correlation with biomass and leaf area. We conclude that dynamic photosynthetic traits and NPQ are important within O. glaberrima, and we highlight NPQ kinetics and NPQ under low light.

Deciphering the Genetic Architecture of Plant Virus Resistance by GWAS, State of the Art and Potential Advances

Growing virus resistant varieties is a highly effective means to avoid yield loss due to infection by many types of virus. The challenge is to be able to detect resistance donors within plant species diversity and then quickly introduce alleles conferring resistance into elite genetic backgrounds. Until now, mainly monogenic forms of resistance with major effects have been introduced in crops. Polygenic resistance is harder to map and introduce in susceptible genetic backgrounds, but it is likely more durable. Genome wide association studies (GWAS) offer an opportunity to accelerate mapping of both monogenic and polygenic resistance, but have seldom been implemented and described in the plant-virus interaction context. Yet, all of the 48 plant-virus GWAS published so far have successfully mapped QTLs involved in plant virus resistance. In this review, we analyzed general and specific GWAS issues regarding plant virus resistance. We have identified and described several key steps throughout the GWAS pipeline, from diversity panel assembly to GWAS result analyses. Based on the 48 published articles, we analyzed the impact of each key step on the GWAS power and showcase several GWAS methods tailored to all types of viruses.

Dissection of the Genetic Basis of Rice Panicle Architecture Using a Genome-wide Association Study

Panicle architecture is one of the major factors influencing productivity of rice crops. The regulatory mechanisms underlying this complex trait are still unclear and genetic resources for rice breeders to improve panicle architecture are limited. Here, we have performed a genome-wide association study (GWAS) to analyze and identify genetic determinants underlying three panicle architecture traits. A population of 340 rice accessions from the 3000 Rice Genomes Project was phenotyped for panicle length, primary panicle number and secondary branch number over two years; GWAS was performed across the whole panel, and also across the japonica and indica sub-panels. A total of 153 quantitative trait loci (QTLs) were detected, of which 5 were associated with multiple traits, 8 were unique to either indica or japonica sub-panels, while 37 QTLs were stable across both years. Using haplotype and expression analysis, we reveal that genetic variations in the OsSPL18 promoter significantly affect gene expression and correlate with panicle length phenotypes. Three new candidate genes with putative roles in determining panicle length were also identified. Haplotype analysis of OsGRRP and LOC_Os03g03480 revealed high association with panicle length variation. Gene expression of DSM2, involved in abscisic acid biosynthesis, was up-regulated in long panicle accessions. Our results provide valuable information and resources for further unravelling the genetic basis determining rice panicle architecture. Identified candidate genes and molecular markers can be used in marker-assisted selection to improve rice panicle architecture through molecular breeding.

Introgression among North American wild grapes (Vitis) fuels biotic and abiotic adaptation

BACKGROUND

Introgressive hybridization can reassort genetic variants into beneficial combinations, permitting adaptation to new ecological niches. To evaluate evolutionary patterns and dynamics that contribute to introgression, we investigate six wild Vitis species that are native to the Southwestern United States and useful for breeding grapevine (V. vinifera) rootstocks.

RESULTS

By creating a reference genome assembly from one wild species, V. arizonica, and by resequencing 130 accessions, we focus on identifying putatively introgressed regions (pIRs) between species. We find six species pairs with signals of introgression between them, comprising up to ~ 8% of the extant genome for some pairs. The pIRs tend to be gene poor, located in regions of high recombination and enriched for genes implicated in disease resistance functions. To assess potential pIR function, we explore SNP associations to bioclimatic variables and to bacterial levels after infection with the causative agent of Pierce's disease (Xylella fastidiosa). pIRs are enriched for SNPs associated with both climate and bacterial levels, suggesting that introgression is driven by adaptation to biotic and abiotic stressors.

CONCLUSIONS

Altogether, this study yields insights into the genomic extent of introgression, potential pressures that shape adaptive introgression, and the evolutionary history of economically important wild relatives of a critical crop.

Genome-wide association study and its applications in the non-model crop Sesamum indicum

BACKGROUND

Sesame is a rare example of non-model and minor crop for which numerous genetic loci and candidate genes underlying features of interest have been disclosed at relatively high resolution. These progresses have been achieved thanks to the applications of the genome-wide association study (GWAS) approach. GWAS has benefited from the availability of high-quality genomes, re-sequencing data from thousands of genotypes, extensive transcriptome sequencing, development of haplotype map and web-based functional databases in sesame.

RESULTS

In this paper, we reviewed the GWAS methods, the underlying statistical models and the applications for genetic discovery of important traits in sesame. A novel online database SiGeDiD ( http://sigedid.ucad.sn/ ) has been developed to provide access to all genetic and genomic discoveries through GWAS in sesame. We also tested for the first time, applications of various new GWAS multi-locus models in sesame.

CONCLUSIONS

Collectively, this work portrays steps and provides guidelines for efficient GWAS implementation in sesame, a non-model crop.