Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa

Advances and Challenges for QTL Analysis and GWAS in the Plant-Breeding of High-Yielding: A Focus on Rapeseed

Yield is one of the most important agronomic traits for the breeding of rapeseed (Brassica napus L), but its genetic dissection for the formation of high yield remains enigmatic, given the rapid population growth. In the present review, we review the discovery of major loci underlying important agronomic traits and the recent advancement in the selection of complex traits. Further, we discuss the benchmark summary of high-throughput techniques for the high-resolution genetic breeding of rapeseed. Biparental linkage analysis and association mapping have become powerful strategies to comprehend the genetic architecture of complex agronomic traits in crops. The generation of improved crop varieties, especially rapeseed, is greatly urged to enhance yield productivity. In this sense, the whole-genome sequencing of rapeseed has become achievable to clone and identify quantitative trait loci (QTLs). Moreover, the generation of high-throughput sequencing and genotyping techniques has significantly enhanced the precision of QTL mapping and genome-wide association study (GWAS) methodologies. Furthermore, this study demonstrates the first attempt to identify novel QTLs of yield-related traits, specifically focusing on ovule number per pod (ON). We also highlight the recent breakthrough concerning single-locus-GWAS (SL-GWAS) and multi-locus GWAS (ML-GWAS), which aim to enhance the potential and robust control of GWAS for improved complex traits.

Genomic insight into balancing high yield, good quality, and blast resistance of japonica rice

BACKGROUND

Balancing the yield, quality and resistance to disease is a daunting challenge in crop breeding due to the negative relationship among these traits. Large-scale genomic landscape analysis of germplasm resources is considered to be an efficient approach to dissect the genetic basis of the complex traits. Central China is one of the main regions where the japonica rice is produced. However, dozens of high-yield rice varieties in this region still exist with low quality or susceptibility to blast disease, severely limiting their application in rice production.

RESULTS

Here, we re-sequence 200 japonica rice varieties grown in central China over the past 30 years and analyze the genetic structure of these cultivars using 2.4 million polymorphic SNP markers. Genome-wide association mapping and selection scans indicate that strong selection for high-yield and taste quality associated with low-amylose content may have led to the loss of resistance to the rice blast fungus Magnaporthe oryzae. By extensive bioinformatic analyses of yield components, resistance to rice blast, and taste quality, we identify several superior alleles for these traits in the population. Based on this information, we successfully introduce excellent taste quality and blast-resistant alleles into the background of two high-yield cultivars and develop two elite lines, XY99 and JXY1, with excellent taste, high yield, and broad-spectrum of blast resistance.

CONCLUSIONS

This is the first large-scale genomic landscape analysis of japonica rice varieties grown in central China and we demonstrate a balancing of multiple agronomic traits by genomic-based strategy.

A genome-wide association study reveals that the 2-oxoglutarate/malate translocator mediates seed vigor in rice

Seed vigor is an important trait for the direct seeding of rice (Oryza sativa L.). In this study, we examined the genetic architecture of variation in the germination rate using a diverse panel of rice accessions. Four quantitative trait loci for germination rate were identified using a genome-wide association study during early germination. One candidate gene, encoding the 2-oxoglutarate/malate translocator (OsOMT), was validated for qGR11. Disruption of this gene (Osomt mutants) reduced seed vigor, including seed germination and seedling growth, in rice. Functional analysis revealed that OsOMT influences seed vigor mainly by modulating amino acid levels and glycolysis and tricarboxylic acid cycle processes. The levels of most amino acids, including the Glu family (Glu, Pro, Arg, and GABA), Asp family (Asp, Thr, Lys, Ile, and Met), Ser family (Ser, Gly, and Cys), and others (His, Ala, Leu, and Val), were significantly reduced in the mature grains and the early germinating seeds of Osomt mutants compared to wild type (WT). The glucose and soluble sugar contents, as well as adenosine triphosphate levels, were significantly decreased in germinating seeds of Osomt mutants compared to WT. These results provide important insights into the role of OsOMT in seed vigor in rice.

Genetic Dissection of Grain Yield Component Traits Under High Nighttime Temperature Stress in a Rice Diversity Panel

To dissect the genetic complexity of rice grain yield (GY) and quality in response to heat stress at the reproductive stage, a diverse panel of 190 rice accessions in the United States Department of Agriculture (USDA) rice mini-core collection (URMC) diversity panel were treated with high nighttime temperature (HNT) stress at the reproductive stage of panicle initiation. The quantifiable yield component response traits were then measured. The traits, panicle length (PL), and number of spikelets per panicle (NSP) were evaluated in subsets of the panel comprising the rice subspecies Oryza sativa ssp. Indica and ssp. Japonica. Under HNT stress, the Japonica ssp. exhibited lower reductions in PL and NSP and a higher level of genetic variation compared with the other subpopulations. Whole genome sequencing identified 6.5 million single nucleotide polymorphisms (SNPs) that were used for the genome-wide association studies (GWASs) of the PL and NSP traits. The GWAS analysis in the Combined, Indica, and Japonica populations under HNT stress identified 83, 60, and 803 highly significant SNPs associated with PL, compared to the 30, 30, and 11 highly significant SNPs associated with NSP. Among these trait-associated SNPs, 140 were coincident with genomic regions previously reported for major GY component quantitative trait loci (QTLs) under heat stress. Using extents of linkage disequilibrium in the rice populations, Venn diagram analysis showed that the highest number of putative candidate genes were identified in the Japonica population, with 20 putative candidate genes being common in the Combined, Indica and Japonica populations. Network analysis of the genes linked to significant SNPs associated with PL and NSP identified modules that were involved in primary and secondary metabolisms. The findings in this study could be useful to understand the pathways/mechanisms involved in rice GY and its components under HNT stress for the acceleration of rice-breeding programs and further functional analysis by molecular geneticists.

Overexpression of OsPHR3 improves growth traits and facilitates nitrogen use efficiency under low phosphate condition

Phosphorus (P) and nitrogen (N) are both essential macronutrients for maintaining plant growth and development. In rice (Oryza sativa L.), OsPHR3 is one of the four paralogs of PHR1, which acts as a central regulator of phosphate (Pi) homeostasis, as well being involved in N homeostasis. However, the functions of OsPHR3 in N utilization under different Pi conditions have yet to be fully studied. In this study, we aimed to dissect the effect of OsPHR3-overexpression on N utilization under Pi deficient regimes. Biochemical, molecular and physiological assays were performed to determine the N-influx, translocation, and accumulation in OsPHR3-overexpressing rice lines, grown under Pi-sufficient and -deficient conditions, in both hydroponic and soil systems. Furthermore, important agronomic traits of these plants were also evaluated. The overexpression of OsPHR3 increased N uptake under Pi stress regimes. Increased N uptake also elevated total N concentrations in these plants by inducing N transporter genes expression. Furthermore, overexpression of OsPHR3 increased N use efficiency, 1000-grain weight and grain yield under different Pi conditions. We established new findings that OsPHR3-overexpression facilitates N utilization under Pi deficient conditions. This will help achieving higher yields by coordinating the utilization of N and P.

Reciprocal adaptation of rice and Xanthomonas oryzae pv. oryzae: cross-species 2D GWAS reveals the underlying genetics

A 1D/2D genome-wide association study strategy was adopted to investigate the genetic systems underlying the reciprocal adaptation of rice (Oryza sativa) and its bacterial pathogen, Xanthomonas oryzae pv. oryzae (Xoo) using the whole-genome sequencing and large-scale phenotyping data of 701 rice accessions and 23 diverse Xoo strains. Forty-seven Xoo virulence-related genes and 318 rice quantitative resistance genes (QR-genes) mainly located in 41 genomic regions, and genome-wide interactions between the detected virulence-related genes and QR genes were identified, including well-known resistance genes/virulence genes plus many previously uncharacterized ones. The relationship between rice and Xoo was characterized by strong differentiation among Xoo races corresponding to the subspecific differentiation of rice, by strong shifts toward increased resistance/virulence of rice/Xoo populations and by rich genetic diversity at the detected rice QR-genes and Xoo virulence genes, and by genome-wide interactions between many rice QR-genes and Xoo virulence genes in a multiple-to-multiple manner, presumably resulting either from direct protein-protein interactions or from genetic epistasis. The observed complex genetic interaction system between rice and Xoo likely exists in other crop-pathogen systems that would maintain high levels of diversity at their QR-loci/virulence-loci, resulting in dynamic coevolutionary consequences during their reciprocal adaptation.

Genome-Wide Association Mapping of Mixed Linkage (1,3;1,4)-β-Glucan and Starch Contents in Rice Whole Grain

The glucan content of rice is a key factor defining its nutritional and economic value. Starch and its derivatives have many industrial applications such as in fuel and material production. Non-starch glucans such as (1,3;1,4)-β-D-glucan (mixed-linkage β-glucan, MLG) have many benefits in human health, including lowering cholesterol, boosting the immune system, and modulating the gut microbiome. In this study, the genetic variability of MLG and starch contents were analyzed in rice (Oryza sativa L.) whole grain, by performing a new quantitative analysis of the polysaccharide content of rice grains. The 197 rice accessions investigated had an average MLG content of 252 μg/mg, which was negatively correlated with the grain starch content. A new genome-wide association study revealed seven significant quantitative trait loci (QTLs) associated with the MLG content and two QTLs associated with the starch content in rice whole grain. Novel genes associated with the MLG content were a hexose transporter and anthocyanidin 5,3-O-glucosyltransferase. Also, the novel gene associated with the starch content was a nodulin-like domain. The data pave the way for a better understanding of the genes involved in determining both MLG and starch contents in rice grains and should facilitate future plant breeding programs.

Genome-Wide Association Studies Provide Insights into the Genetic Determination of Flower and Leaf Traits of Actinidia eriantha

Kiwifruit (Actinidia eriantha) is a dioecious vine, and the pollen of its male cultivars has a direct effect on the quality of its fruits. In this study, to facilitate molecular breeding and gene identification, we performed genome-wide association studies (GWAS) on 11 traits of flower and leaf. A total of 946,337 highly consistent SNP markers were obtained in the whole genome. Phylogenetic tree analysis and population structure analysis showed that the 143 germplasms can be divided into two groups. The linkage disequilibrium analysis showed that A. eriantha have a relatively fast attenuation rate, and that the average attenuation distance of LD was 0.1-0.3 Kb. The MLM (QK) model was determined as best for correlation analysis, and eight and three SNPs associated with flower- and leaf-related traits were identified, respectively, at 0.01 significance level. However, SNP markers associated with stamen number per flower, pollen viability, total chlorophyll content, and total flavonoid content were not identified at the 0.01 significant level, although it is worth noting that one, one, five, and two SNPs were identified to be associated with these traits at the 0.05 significant level. This study provides insights into the complex flower- and leaf-related biology, and identifies genes controlling important traits in A. eriantha through GWAS, which extends the genetic resources and basis for facilitating molecular breeding in kiwifruits.

Novel Sources of Pre-Harvest Sprouting Resistance for Japonica Rice Improvement

Pre-harvest sprouting (PHS), induced by unexpected weather events, such as typhoons, at the late seed maturity stage, is becoming a serious threat to rice production, especially in the state of California, USA, Japan, and the Republic of Korea, where japonica varieties (mostly susceptible to PHS) are mainly cultivated. A projected economic loss by severe PHS in these three countries could range between 8-10 billion USD per year during the next 10 years. Here, we present promising rice germplasm with strong resistance to PHS that were selected from a diverse rice panel of accessions held in the International Rice Genebank (IRG) at the International Rice Research Institute (IRRI). To induce PHS, three panicle samples per accession were harvested at 20 and 30 days after flowering (DAF), respectively, and incubated at 100% relative humidity (RH), 30 °C in a growth chamber for 15 days. A genome-wide association (GWA) analysis using a 4.8 million single nucleotide polymorphisms (SNP) marker set was performed to identify loci and candidate genes conferring PHS resistance. Interestingly, two tropical japonica and four temperate japonica accessions showed outstanding PHS resistance as compared to tolerant indica accessions. Two major loci on chromosomes 1 and 4 were associated with PHS resistance. A priori candidate genes interactions with rice gene networks, which are based on the gene ontology (GO), co-expression, and other evidence, suggested that a key resistance mechanism is related to abscisic acid (ABA), gibberellic acid (GA), and auxin mediated signaling pathways.

Loci and natural alleles for cadmium-mediated growth responses revealed by a genome wide association study and transcriptome analysis in rice

BACKGROUND

Cadmium (Cd) is a toxic heavy metal that is harmful to the environment and human health. Cd pollution threatens the cultivation of rice (Oryza sativa L.) in many countries. Improving rice performance under Cd stress could potentially improve rice productivity.

RESULTS

In this study, 9 growth traits of 188 different cultivated rice accessions under normal and Cd stress conditions were found to be highly variable during the seedling stage. Based on ~3.3 million single nucleotide polymorphisms (SNPs), 119 Cd-mediated growth response (CGR) quantitative trait loci (QTL) were identified by a genome-wide association study (GWAS), 55 of which have been validated by previously reported QTL and 64 were new CGR loci. Combined with the data from the GWAS, transcriptome analysis, gene annotations from the gene ontology (GO) Slim database, and annotations and functions of homologous genes, 148 CGR candidate genes were obtained. Additionally, several reported genes have been found to play certain roles in CGRs. Seven Cd-related cloned genes were found among the CGR genes. Natural elite haplotypes/alleles in these genes that increased Cd tolerance were identified by a haplotype analysis of a diverse mini core collection. More importantly, this study was the first to uncover the natural variations of 5 GST genes that play important roles in CGRs.

CONCLUSION

The exploration of Cd-resistant rice germplasm resources and the identification of elite natural variations related to Cd-resistance will help improve the tolerance of current major rice varieties to Cd, as well as provide raw materials and new genes for breeding Cd-resistant varieties.

Loci and natural alleles for cadmium-mediated growth responses revealed by a genome wide association study and transcriptome analysis in rice

BACKGROUND

Cadmium (Cd) is a toxic heavy metal that is harmful to the environment and human health. Cd pollution threatens the cultivation of rice (Oryza sativa L.) in many countries. Improving rice performance under Cd stress could potentially improve rice productivity.

RESULTS

In this study, 9 growth traits of 188 different cultivated rice accessions under normal and Cd stress conditions were found to be highly variable during the seedling stage. Based on ~3.3 million single nucleotide polymorphisms (SNPs), 119 Cd-mediated growth response (CGR) quantitative trait loci (QTL) were identified by a genome-wide association study (GWAS), 55 of which have been validated by previously reported QTL and 64 were new CGR loci. Combined with the data from the GWAS, transcriptome analysis, gene annotations from the gene ontology (GO) Slim database, and annotations and functions of homologous genes, 148 CGR candidate genes were obtained. Additionally, several reported genes have been found to play certain roles in CGRs. Seven Cd-related cloned genes were found among the CGR genes. Natural elite haplotypes/alleles in these genes that increased Cd tolerance were identified by a haplotype analysis of a diverse mini core collection. More importantly, this study was the first to uncover the natural variations of 5 GST genes that play important roles in CGRs.

CONCLUSION

The exploration of Cd-resistant rice germplasm resources and the identification of elite natural variations related to Cd-resistance will help improve the tolerance of current major rice varieties to Cd, as well as provide raw materials and new genes for breeding Cd-resistant varieties.

Identification and allele mining of new candidate genes underlying rice grain weight and grain shape by genome-wide association study

BACKGROUND

Grain weight and grain shape are important agronomic traits that affect the grain yield potential and grain quality of rice. Both grain weight and grain shape are controlled by multiple genes. The 3,000 Rice Genomes Project (3 K RGP) greatly facilitates the discovery of agriculturally important genetic variants and germplasm resources for grain weight and grain shape.

RESULTS

Abundant natural variations and distinct phenotic differentiation among the subgroups in grain weight and grain shape were observed in a large population of 2,453 accessions from the 3 K RGP. A total of 21 stable quantitative trait nucleotides (QTNs) for the four traits were consistently identified in at least two of 3-year trials by genome-wide association study (GWAS), including six new QTNs (qTGW3.1, qTGW9, qTGW11, qGL4/qRLW4, qGL10, and qRLW1) for grain weight and grain shape. We further predicted seven candidate genes (Os03g0186600, Os09g0544400, Os11g0163600, Os04g0580700, Os10g0399700, Os10g0400100 and Os01g0171000) for the six new QTNs by high-density association and gene-based haplotype analyses. The favorable haplotypes of the seven candidate genes and five previously cloned genes in elite accessions with high TGW and RLW are also provided.

CONCLUSIONS

Our results deepen the understanding of the genetic basis of grain weight and grain shape in rice and provide valuable information for improving rice grain yield and grain quality through molecular breeding.

Potential of rice landraces with strong culms as genetic resources for improving lodging resistance against super typhoons

It is generally believed that rice landraces with long culms are susceptible to lodging, and have not been utilized for breeding to improve lodging resistance. However, little is known about the structural culm strength of landraces and their beneficial genetic loci. Therefore, in this study, genome-wide association studies (GWAS) were performed using a rice population panel including Japanese rice landraces to identify beneficial loci associated with strong culms. As a result, the landraces were found to have higher structural culm strength and greater diversity than the breeding varieties. Genetic loci associated with strong culms were identified, and it was demonstrated that haplotypes with positive effects of those loci were present in a high proportion of these landraces. These results indicated that the utilization of the strong culm-associated loci present in Japanese rice landraces may further improve the lodging resistance of modern breeding varieties that have relied on semi-dwarfism.

Genetic architecture of root and shoot ionomes in rice (Oryza sativa L.)

KEY MESSAGE

Association analysis for ionomic concentrations of 20 elements identified independent genetic factors underlying the root and shoot ionomes of rice, providing a platform for selecting and dissecting causal genetic variants. Understanding the genetic basis of mineral nutrient acquisition is key to fully describing how terrestrial organisms interact with the non-living environment. Rice (Oryza sativa L.) serves both as a model organism for genetic studies and as an important component of the global food system. Studies in rice ionomics have primarily focused on above ground tissues evaluated from field-grown plants. Here, we describe a comprehensive study of the genetic basis of the rice ionome in both roots and shoots of 6-week-old rice plants for 20 elements using a controlled hydroponics growth system. Building on the wealth of publicly available rice genomic resources, including a panel of 373 diverse rice lines, 4.8 M genome-wide single-nucleotide polymorphisms, single- and multi-marker analysis pipelines, an extensive tome of 321 candidate genes and legacy QTLs from across 15 years of rice genetics literature, we used genome-wide association analysis and biparental QTL analysis to identify 114 genomic regions associated with ionomic variation. The genetic basis for root and shoot ionomes was highly distinct; 78 loci were associated with roots and 36 loci with shoots, with no overlapping genomic regions for the same element across tissues. We further describe the distribution of phenotypic variation across haplotypes and identify candidate genes within highly significant regions associated with sulfur, manganese, cadmium, and molybdenum. Our analysis provides critical insight into the genetic basis of natural phenotypic variation for both root and shoot ionomes in rice and provides a comprehensive resource for dissecting and testing causal genetic variants.

Genomic signatures of natural selection at phenology-related genes in a widely distributed tree species Fagus sylvatica L

BACKGROUND

Diversity among phenology-related genes is predicted to be a contributing factor in local adaptations seen in widely distributed plant species that grow in climatically variable geographic areas, such as forest trees. European beech (Fagus sylvatica L.) is widespread, and is one of the most important broadleaved tree species in Europe; however, its potential for adaptation to climate change is a matter of uncertainty, and little is known about the molecular basis of climate change-relevant traits like bud burst.

RESULTS

We explored single nucleotide polymorphisms (SNP) at candidate genes related to bud burst in beech individuals sampled across 47 populations from Europe. SNP diversity was monitored for 380 candidate genes using a sequence capture approach, providing 2909 unlinked SNP loci. We used two complementary analytical methods to find loci significantly associated with geographic variables, climatic variables (expressed as principal components), or phenotypic variables (spring and autumn phenology, height, survival). Redundancy analysis (RDA) was used to detect candidate markers across two spatial scales (entire study area and within subregions). We revealed 201 candidate SNPs at the broadest scale, 53.2% of which were associated with phenotypic variables. Additive polygenic scores, which provide a measure of the cumulative signal across significant candidate SNPs, were correlated with a climate variable (first principal component, PC1) related to temperature and precipitation availability, and spring phenology. However, different genotype-environment associations were identified within Southeastern Europe as compared to the entire geographic range of European beech.

CONCLUSIONS

Environmental conditions play important roles as drivers of genetic diversity of phenology-related genes that could influence local adaptation in European beech. Selection in beech favors genotypes with earlier bud burst under warmer and wetter habitats within its range; however, selection pressures may differ across spatial scales.

Validation of a novel associative transcriptomics pipeline in Brassica oleracea: identifying candidates for vernalisation response

BACKGROUND

Associative transcriptomics has been used extensively in Brassica napus to enable the rapid identification of markers correlated with traits of interest. However, within the important vegetable crop species, Brassica oleracea, the use of associative transcriptomics has been limited due to a lack of fixed genetic resources and the difficulties in generating material due to self-incompatibility. Within Brassica vegetables, the harvestable product can be vegetative or floral tissues and therefore synchronisation of the floral transition is an important goal for growers and breeders. Vernalisation is known to be a key determinant of the floral transition, yet how different vernalisation treatments influence flowering in B. oleracea is not well understood.

RESULTS

Here, we present results from phenotyping a diverse set of 69 B. oleracea accessions for heading and flowering traits under different environmental conditions. We developed a new associative transcriptomics pipeline, and inferred and validated a population structure, for the phenotyped accessions. A genome-wide association study identified miR172D as a candidate for the vernalisation response. Gene expression marker association identified variation in expression of BoFLC.C2 as a further candidate for vernalisation response.

CONCLUSIONS

This study describes a new pipeline for performing associative transcriptomics studies in B. oleracea. Using flowering time as an example trait, it provides insights into the genetic basis of vernalisation response in B. oleracea through associative transcriptomics and confirms its characterisation as a complex G x E trait. Candidate leads were identified in miR172D and BoFLC.C2. These results could facilitate marker-based breeding efforts to produce B. oleracea lines with more synchronous heading dates, potentially leading to improved yields.

A statistical framework for QTL hotspot detection

Quantitative trait loci (QTL) hotspots (genomic locations enriched in QTL) are a common and notable feature when collecting many QTL for various traits in many areas of biological studies. The QTL hotspots are important and attractive since they are highly informative and may harbor genes for the quantitative traits. So far, the current statistical methods for QTL hotspot detection use either the individual-level data from the genetical genomics experiments or the summarized data from public QTL databases to proceed with the detection analysis. These methods may suffer from the problems of ignoring the correlation structure among traits, neglecting the magnitude of LOD scores for the QTL, or paying a very high computational cost, which often lead to the detection of excessive spurious hotspots, failure to discover biologically interesting hotspots composed of a small-to-moderate number of QTL with strong LOD scores, and computational intractability, respectively, during the detection process. In this article, we describe a statistical framework that can handle both types of data as well as address all the problems at a time for QTL hotspot detection. Our statistical framework directly operates on the QTL matrix and hence has a very cheap computational cost and is deployed to take advantage of the QTL mapping results for assisting the detection analysis. Two special devices, trait grouping and top γn,α profile, are introduced into the framework. The trait grouping attempts to group the traits controlled by closely linked or pleiotropic QTL together into the same trait groups and randomly allocates these QTL together across the genomic positions separately by trait group to account for the correlation structure among traits, so as to have the ability to obtain much stricter thresholds and dismiss spurious hotspots. The top γn,α profile is designed to outline the LOD-score pattern of QTL in a hotspot across the different hotspot architectures, so that it can serve to identify and characterize the types of QTL hotspots with varying sizes and LOD-score distributions. Real examples, numerical analysis, and simulation study are performed to validate our statistical framework, investigate the detection properties, and also compare with the current methods in QTL hotspot detection. The results demonstrate that the proposed statistical framework can effectively accommodate the correlation structure among traits, identify the types of hotspots, and still keep the notable features of easy implementation and fast computation for practical QTL hotspot detection.

Genetic relationships and identification of core germplasm among rice photoperiod- and thermo-sensitive genic male sterile lines

BACKGROUND

Harnessing heterosis is one of the major approaches to increase rice yield and has made a great contribution to food security. The identification and selection of outstanding parental genotypes especially among male sterile lines is a key step for exploiting heterosis. Two-line hybrid system is based on the discovery and application of photoperiod- and thermo-sensitive genic sensitive male sterile (PTGMS) materials. The development of wide-range of male sterile lines from a common gene pool leads to a narrower genetic diversity, which is vulnerable to biotic and abiotic stress. Hence, it is valuable to ascertain the genetic background of PTGMS lines and to understand their relationships in order to select and design a future breeding strategy.

RESULTS

A collection of 118 male sterile rice lines and 13 conventional breeding lines from the major rice growing regions of China was evaluated and screened against the photosensitive (pms3) and temperature sensitive male sterility (tms5) genes. The total gene pool was divided into four major populations as P1 possessing the pms3, P2 possessing tms5, P3 possessing both pms3 and tms5 genes, and P4 containing conventional breeding lines without any male sterility allele. The high genetic purity was revealed by homozygous alleles in all populations. The population admixture, principle components and the phylogenetic analysis revealed the close relations of P2 and P3 with P4. The population differentiation analysis showed that P1 has the highest differentiation coefficient. The lines from P1 were observed as the ancestors of other three populations in a phylogenetic tree, while the lines in P2 and P3 showed a close genetic relation with conventional lines. A core collection of top 10% lines with maximum within and among populations genetic diversity was constructed for future research and breeding efforts.

CONCLUSION

The low genetic diversity and close genetic relationship among PTGMS lines in P2, P3 and P4 populations suggest a selection sweep and they might result from a backcrossing with common ancestors including the pure lines of P1. The core collection from PTGMS panel updated with new diverse germplasm will serve best for further two-line hybrid breeding.

Genome-Wide Association Mapping and Gene Expression Analysis Reveal the Negative Role of OsMYB21 in Regulating Bacterial Blight Resistance in Rice

BACKGROUND

Bacterial blight (BB), caused by Xanthomonas oryzae pv. oryzae (Xoo), is one of the most devastating diseases in rice all over the world. Due to the diversity and rapid evolution of Xoo, identification and use of the non-race specific quantitative resistance QTLs has been considered the preferred strategy for effective control of this disease. Although numerous QTLs for BB resistance have been identified, they haven't been effectively used for improvement of BB resistance in rice due to their small effects and lack of knowledge on the function of genes underlying the QTLs.

RESULTS

In the present study, a genome-wide association study of BB resistance was performed in a rice core collection from South China. A total of 17 QTLs were identified to be associated with BB resistance. Among them, 13 QTLs were newly identified in the present study and the other 4 QTLs were co-localized with the previously reported QTLs or Xa genes that confer qualitative resistance to Xoo strains. Particularly, the qBBR11-4 on chromosome 11 explained the largest phenotypic variation in this study and was co-localized with the previously identified QTLs for BB and bacterial leaf streak (BLS) resistance against diverse strains in three studies, suggesting its broad-spectrum resistance and potential value in rice breeding. Through combined analysis of differential expression and annotations of the predicted genes within qBBR11-4 between two sets of rice accessions selected based on haplotypes and disease phenotypes, we identified the transcription factor OsMYB21 as the candidate gene for qBBR11-4. The OsMYB21 overexpressing plants exhibited decreased resistance to bacterial blight, accompanied with down-regulation of several defense-related genes compared with the wild-type plants.

CONCLUSION

The results suggest that OsMYB21 negatively regulates bacterial blight resistance in rice, and this gene can be a promising target in rice breeding by using the gene editing method. In addition, the potential candidate genes for the 13 novel QTLs for BB resistance were also analyzed in this study, providing a new source for cloning of genes associated with BB resistance and molecular breeding in rice.

Robustification of GWAS to explore effective SNPs addressing the challenges of hidden population stratification and polygenic effects

Genome-wide association studies (GWAS) play a vital role in identifying important genes those is associated with the phenotypic variations of living organisms. There are several statistical methods for GWAS including the linear mixed model (LMM) which is popular for addressing the challenges of hidden population stratification and polygenic effects. However, most of these methods including LMM are sensitive to phenotypic outliers that may lead the misleading results. To overcome this problem, in this paper, we proposed a way to robustify the LMM approach for reducing the influence of outlying observations using the β-divergence method. The performance of the proposed method was investigated using both synthetic and real data analysis. Simulation results showed that the proposed method performs better than both linear regression model (LRM) and LMM approaches in terms of powers and false discovery rates in presence of phenotypic outliers. On the other hand, the proposed method performed almost similar to LMM approach but much better than LRM approach in absence of outliers. In the case of real data analysis, our proposed method identified 11 SNPs that are significantly associated with the rice flowering time. Among the identified candidate SNPs, some were involved in seed development and flowering time pathways, and some were connected with flower and other developmental processes. These identified candidate SNPs could assist rice breeding programs effectively. Thus, our findings highlighted the importance of robust GWAS in identifying candidate genes.

Importance of correcting genomic relationships in single-locus QTL mapping model with an advanced backcross population

Advanced backcross (AB) populations have been widely used to identify and utilize beneficial alleles in various crops such as rice, tomato, wheat, and barley. For the development of an AB population, a controlled crossing scheme is used and this controlled crossing along with the selection (both natural and artificial) of agronomically adapted alleles during the development of AB population may lead to unbalanced allele frequencies in the population. However, it is commonly believed that interval mapping of traits in experimental crosses such as AB populations is immune to the deviations from the expected frequencies under Mendelian segregation. Using two AB populations and simulated data sets as examples, we describe the severity of the problem caused by unbalanced allele frequencies in quantitative trait loci mapping and demonstrate how it can be corrected using the linear mixed model having a polygenic effect with the covariance structure (genomic relationship matrix) calculated from molecular markers.

Genotyping-by-Sequencing-Based Genome-Wide Association Studies of Fusarium Wilt Resistance in Radishes (Raphanus sativus L.)

Fusarium wilt (FW) is a fungal disease that causes severe yield losses in radish production. The most effective method to control the FW is the development and use of resistant varieties in cultivation. The identification of marker loci linked to FW resistance are expected to facilitate the breeding of disease-resistant radishes. In the present study, we applied an integrated framework of genome-wide association studies (GWAS) using genotyping-by-sequencing (GBS) to identify FW resistance loci among a panel of 225 radish accessions, including 58 elite breeding lines. Phenotyping was conducted by manual inoculation of seedlings with the FW pathogen, and scoring for the disease index was conducted three weeks after inoculation during two constitutive years. The GWAS analysis identified 44 single nucleotide polymorphisms (SNPs) and twenty putative candidate genes that were significantly associated with FW resistance. In addition, a total of four QTLs were identified from F2 population derived from a FW resistant line and a susceptible line, one of which was co-located with the SNPs on chromosome 7, detected in GWAS study. These markers will be valuable for molecular breeding programs and marker-assisted selection to develop FW resistant varieties of R. sativus.

Coping with inclement weather conditions due to high temperature and water deficit in rice: An insight from genetic and biochemical perspectives

Climatic fluctuations, temperature extremes, and water scarcity are becoming increasingly unpredictable with the passage of time. Such environmental atrocities have been the scourge of agriculture over the ages, bringing with them poor harvests and threat of famine. Rice production, owing to its high-water requirement for cultivation, is highly vulnerable to the threat of changing climate, particularly prolonged drought and high temperature, individually or in combination. Amidst all the abiotic stresses, heat and drought are considered as the most important concurrent stressors, largely affecting rice yield and productivity under the current scenario. Such threats heighten the need for new breeding and cultivation strategies in generating abiotic stress-resilient rice varieties with better yield potential. Responses of rice to these stresses can be categorized at the morphological, physiological and biochemical levels. This review examines the physiological and molecular mechanism, in the form of up regulation of several defense machineries of rice varieties to cope with drought stress (DS), high temperature stress (HTS), and their combination (DS-HTS). Genotypic differences among rice varieties in their tolerance ability have also been addressed. The review also appraises research studies conducted in rice regarding various phenotypic traits, genetic loci and response mechanisms to stress conditions to help craft new breeding strategies for improved tolerance to DS and HTS, singly or in combination. The review also encompasses the gene regulatory networks and transcription factors, and their cross-talks in mediating tolerance to such stresses. Understanding the epigenetic regulation, involving DNA methylation and histone modification during such hostile situations, will also play a crucial role in our comprehensive understanding of combinatorial stress responses. Taken together, this review consolidates current research and available information on promising rice cultivars with desirable traits as well as advocates synergistic and complementary approaches in molecular and systems biology to develop new rice breeds that favorably respond to DS-HTS-induced abiotic stress.

Natural variation in plant telomere length is associated with flowering time

Telomeres are highly repetitive DNA sequences found at the ends of chromosomes that protect the chromosomes from deterioration duringcell division. Here, using whole-genome re-sequencing and terminal restriction fragment assays, we found substantial natural intraspecific variation in telomere length in Arabidopsis thaliana, rice (Oryza sativa), and maize (Zea mays). Genome-wide association study (GWAS) mapping in A. thaliana identified 13 regions with GWAS-significant associations underlying telomere length variation, including a region that harbors the telomerase reverse transcriptase (TERT) gene. Population genomic analysis provided evidence for a selective sweep at the TERT region associated with longer telomeres. We found that telomere length is negatively correlated with flowering time variation not only in A. thaliana, but also in maize and rice, indicating a link between life-history traits and chromosome integrity. Our results point to several possible reasons for this correlation, including the possibility that longer telomeres may be more adaptive in plants that have faster developmental rates (and therefore flower earlier). Our work suggests that chromosomal structure itself might be an adaptive trait associated with plant life-history strategies.

Chromosomal characteristics of salt stress heritable gene expression in the rice genome

BACKGROUND

Gene expression is potentially an important heritable quantitative trait that mediates between genetic variation and higher-level complex phenotypes through time and condition-dependent regulatory interactions. Therefore, we sought to explore both the genomic and condition-specific characteristics of gene expression heritability within the context of chromosomal structure.

RESULTS

Heritability was estimated for biological gene expression using a diverse, 84-line, Oryza sativa (rice) population under optimal and salt-stressed conditions. Overall, 5936 genes were found to have heritable expression regardless of condition and 1377 genes were found to have heritable expression only during salt stress. These genes with salt-specific heritable expression are enriched for functional terms associated with response to stimulus and transcription factor activity. Additionally, we discovered that highly and lowly expressed genes, and genes with heritable expression are distributed differently along the chromosomes in patterns that follow previously identified high-throughput chromosomal conformation capture (Hi-C) A/B chromatin compartments. Furthermore, multiple genomic hot-spots enriched for genes with salt-specific heritability were identified on chromosomes 1, 4, 6, and 8. These hotspots were found to contain genes functionally enriched for transcriptional regulation and overlaps with a previously identified major QTL for salt-tolerance in rice.

CONCLUSIONS

Investigating the heritability of traits, and in-particular gene expression traits, is important towards developing a basic understanding of how regulatory networks behave across a population. This work provides insights into spatial patterns of heritable gene expression at the chromosomal level.

Genome-wide comparative analyses of GATA transcription factors among 19 Arabidopsis ecotype genomes: Intraspecific characteristics of GATA transcription factors

GATA transcription factors (TFs) are widespread eukaryotic regulators whose DNA-binding domain is a class IV zinc finger motif (CX2CX17-20CX2C) followed by a basic region. Due to the low cost of genome sequencing, multiple strains of specific species have been sequenced: e.g., number of plant genomes in the Plant Genome Database (http://www.plantgenome.info/) is 2,174 originated from 713 plant species. Thus, we investigated GATA TFs of 19 Arabidopsis thaliana genome-widely to understand intraspecific features of Arabidopsis GATA TFs with the pipeline of GATA database (http://gata.genefamily.info/). Numbers of GATA genes and GATA TFs of each A. thaliana genome range from 29 to 30 and from 39 to 42, respectively. Four cases of different pattern of alternative splicing forms of GATA genes among 19 A. thaliana genomes are identified. 22 of 2,195 amino acids (1.002%) from the alignment of GATA domain amino acid sequences display variations across 19 ecotype genomes. In addition, maximally four different amino acid sequences per each GATA domain identified in this study indicate that these position-specific amino acid variations may invoke intraspecific functional variations. Among 15 functionally characterized GATA genes, only five GATA genes display variations of amino acids across ecotypes of A. thaliana, implying variations of their biological roles across natural isolates of A. thaliana. PCA results from 28 characteristics of GATA genes display the four groups, same to those defined by the number of GATA genes. Topologies of bootstrapped phylogenetic trees of Arabidopsis chloroplasts and common GATA genes are mostly incongruent. Moreover, no relationship between geographical distribution and their phylogenetic relationships was found. Our results present that intraspecific variations of GATA TFs in A. thaliana are conserved and evolutionarily neutral along with 19 ecotypes, which is congruent to the fact that GATA TFs are one of the main regulators for controlling essential mechanisms, such as seed germination and hypocotyl elongation.

Tests of association based on genomic windows can lead to spurious associations when using genotype panels with heterogeneous SNP densities

Dense single nucleotide polymorphism (SNP) panels are widely used for genome-wide association studies (GWAS). In these panels, SNPs within a genomic segment tend to be highly correlated. Thus, association studies based on testing the significance of single SNPs are not very effective, and genomic-window based tests have been proposed to address this problem. However, when the SNP density on the genotype panel is not homogeneous, genomic-window based tests can lead to the detection of spurious associations by declaring effects of genomic windows that explain a large proportion of genetic variance as significant. We propose two methods to solve this problem.

Genomic Approaches to Identify Molecular Bases of Crop Resistance to Diseases and to Develop Future Breeding Strategies

Plant diseases are responsible for substantial crop losses each year and affect food security and agricultural sustainability. The improvement of crop resistance to pathogens through breeding represents an environmentally sound method for managing disease and minimizing these losses. The challenge is to breed varieties with a stable and broad-spectrum resistance. Different approaches, from markers to recent genomic and 'post-genomic era' technologies, will be reviewed in order to contribute to a better understanding of the complexity of host-pathogen interactions and genes, including those with small phenotypic effects and mechanisms that underlie resistance. An efficient combination of these approaches is herein proposed as the basis to develop a successful breeding strategy to obtain resistant crop varieties that yield higher in increasing disease scenarios.

Evolutionary Processes Involved in the Emergence and Expansion of an Atypical O. sativa Group in Madagascar

Understanding crops genetic diversity and the evolutionary processes that accompanied their worldwide spread is useful for designing effective breeding strategies. Madagascar Island was one of the last major Old World areas where human settlement brought the introduction of Oryza sativa. Early studies in the island had reported the presence of a rice group specific to Madagascar. Using 24 K SNP, we compared diversity patterns at the whole genome and at haplotype (30 SNP-long segments along the genome) levels, between 620 Malagasy and 1929 Asian rice accessions. The haplotype level analysis aimed at identifying local genotypic variations, relative to the whole genome level, using a group assignment method that relies on kernel density estimation in a Principal Component Analysis feature space. Migration bottleneck had resulted in 10-25% reduction of diversity among the Malagasy representatives of indica and japonica populations. Compared to their Asian counterpart, they showed slightly lower indica and japonica introgressions, suggesting the two populations had undergone less recombination when migration to the island occurred. The origins of the Malagasy indica and japonica groups were delineated to indica subpopulation from the Indian subcontinent and to tropical japonica from the Malay Archipelago, respectively. The Malagasy-specific group (Gm) had a rather high gene diversity and an original haplotype pattern: much lower share of indica haplotypes, and much higher share of Aus and japonica haplotypes than indica. Its emergence and expansion are most probably due to inter-group recombination facilitated by sympatry between indica-Aus admixes and "Bulu" type landraces of japonica in the central high plateaux of Madagascar, and to human selection for adaptation to the lowland rice cultivation. Pattern of rice genetic diversity was also tightly associated with the history of human settlement in the island. Emergence of the Gm group is associated with the latest arrivals of Austronesians, who founded the Merina kingdom in the high plateaux and developed lowland rice cultivation. As an intermediary form between Aus, indica and japonica, the three pillars of O. sativa domestication, Gm represents a very valuable genetic resource in breeding for adaptation to cold tolerance in tropical highlands. We proposed the name Rojo for this new rice group.

Superior haplotypes towards development of low glycemic index rice with preferred grain and cooking quality

Increasing trends in the occurrence of diabetes underline the need to develop low glycemic index (GI) rice with preferred grain quality. In the current study, a diverse set of 3 K sub-panel of rice consisting of 150 accessions was evaluated for resistant starch and predicted glycemic index, including nine other quality traits under transplanted situation. Significant variations were noticed among the accessions for the traits evaluated. Trait associations had shown that amylose content possess significant positive and negative association with resistant starch and predicted glycemic index. Genome-wide association studies with 500 K SNPs based on MLM model resulted in a total of 41 marker-trait associations (MTAs), which were further confirmed and validated with mrMLM multi-locus model. We have also determined the allelic effect of identified MTAs for 11 targeted traits and found favorable SNPs for 8 traits. A total of 11 genes were selected for haplo-pheno analysis to identify the superior haplotypes for the target traits where haplotypes ranges from 2 (Os10g0469000-GC) to 15 (Os06g18720-AC). Superior haplotypes for RS and PGI, the candidate gene Os06g11100 (H4-3.28% for high RS) and Os08g12590 (H13-62.52 as intermediate PGI). The identified superior donors possessing superior haplotype combinations may be utilized in Haplotype-based breeding to developing next-generation tailor-made high quality healthier rice varieties suiting consumer preference and market demand.

FLOURY ENDOSPERM19 encoding a class I glutamine amidotransferase affects grain quality in rice

As a staple food for more than half of the world's population, the importance of rice is self-evident. Compared with ordinary rice, rice cultivars with superior eating quality and appearance quality are more popular with consumers due to their unique taste and ornamental value, even if their price is much higher. Appearance quality and CEQ (cooking and eating quality) are two very important aspects in the evaluation of rice quality. Here, we performed a genome-wide association study on floury endosperm in a diverse panel of 533 cultivated rice accessions. We identified a batch of potential floury genes and prioritize one (LOC_Os03g48060) for functional analyses. Two floury outer endosperm mutants (flo19-1 and flo19-2) were generated through editing LOC_Os03g48060 (named as FLO19 in this study), which encodes a class I glutamine amidotransferase. The different performances of the two mutants in various storage substances directly led to completely different changes in CEQ. The mutation of FLO19 gene caused the damage of carbon and nitrogen metabolism in rice, which affected the normal growth and development of rice, including decreased plant height and yield loss by decreased grain filling rate. Through haplotype analysis, we identified a haplotype of FLO19 that can improve both CEQ and appearance quality of rice, Hap2, which provides a selection target for rice quality improvement, especially for high-yield indica rice varieties.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01226-z.

Genome-wide discovery of natural variation in pre-mRNA splicing and prioritising causal alternative splicing to salt stress response in rice

Pre-mRNA splicing is an essential step for the regulation of gene expression. In order to specifically capture splicing variants in plants for genome-wide association studies (GWAS), we developed a software tool to quantify and visualise Variations of Splicing in Population (VaSP). VaSP can quantify splicing variants from short-read RNA-seq datasets and discover genotype-specific splicing (GSS) events, which can be used to prioritise causal pre-mRNA splicing events in GWAS. We applied our method to an RNA-seq dataset with 328 samples from 82 genotypes from a rice diversity panel exposed to optimal and saline growing conditions. In total, 764 significant GSS events were identified in salt stress conditions. GSS events were used as markers for a GWAS with the shoot Na⁺ accumulation, which identified six GSS events in five genes significantly associated with the shoot Na⁺ content. Two of these genes, OsNUC1 and OsRAD23 emerged as top candidate genes with splice variants that exhibited significant divergence between the variants for shoot growth under salt stress conditions. VaSP is a versatile tool for alternative splicing analysis in plants and a powerful tool for prioritising candidate causal pre-mRNA splicing and corresponding genomic variations in GWAS.

On the use of GBLUP and its extension for GWAS with additive and epistatic effects

Genomic best linear unbiased prediction (GBLUP) is the most widely used model for genome-wide predictions. Interestingly, it is also possible to perform genome-wide association studies (GWAS) based on GBLUP. Although the estimated marker effects in GBLUP are shrunken and the conventional test based on such effects has low power, it was observed that a modified test statistic can be produced and the result of test was identical to a standard GWAS model. Later, a mathematical proof was given for the special case that there is no fixed covariate in GBLUP. Since then, the new approach has been called “GWAS by GBLUP”. Nevertheless, covariates such as environmental and subpopulation effects are very common in GBLUP. Thus, it is necessary to confirm the equivalence in the general case. Recently, the concept was generalized to GWAS for epistatic effects and the new approach was termed rapid epistatic mixed-model association analysis (REMMA) because it greatly improved the computational efficiency. However, the relationship between REMMA and the standard GWAS model has not been investigated. In this study, we first provided a general mathematical proof of the equivalence between “GWAS by GBLUP” and the standard GWAS model for additive effects. Then, we compared REMMA with the standard GWAS model for epistatic effects by a theoretical investigation and by empirical data analyses. We hypothesized that the similarity of the two models is influenced by the relative contribution of additive and epistatic effects to the phenotypic variance, which was verified by empirical and simulation studies.

Genome-Wide Association Study Reveals Marker-Trait Associations for Early Vegetative Stage Salinity Tolerance in Rice

Rice germplasm is a rich resource for discovering genes associated with salt tolerance. In the current study, a set of 96 accessions were evaluated for seedling stage salinity tolerance and its component traits. Significant phenotypic variation was observed among the genotypes for all the measured traits and eleven accessions with high level of salt tolerance at seedling stage were identified. The germplasm set comprised of three sub-populations and genome-wide association study (GWAS) identified a total of 23 marker-trait associations (MTAs) for traits studied. These MTAs were located on rice chromosomes 1, 2, 5, 6, 7, 9, and 12 and explained the trait phenotypic variances ranging from 13.98 to 29.88 %. Twenty-one MTAs identified in this study were located either in or near the previously reported quantitative trait loci (QTLs), while two MTAs namely, qSDW2.1 and qSNC5 were novel. A total of 18 and 13 putative annotated candidate genes were identified in a genomic region spanning ~200 kb around the MTAs qSDW2.1 and qSNC5, respectively. Some of the important genes underlying the novel MTAs were OsFBA1,OsFBL7, and mTERF which are known to be associated with salinity tolerance in crops. These MTAs pave way for combining salinity tolerance with high yield in rice genotypes through molecular breeding.

Genome-Wide Association Study on Resistance to Rice Black-Streaked Dwarf Disease Caused by Rice black-streaked dwarf virus

Rice black-streaked dwarf disease caused by Rice black-streaked dwarf virus (RBSDV) is one of the most destructive viral diseases of rice. Thus, it is imperative that resistant rice germplasms are screened for novel RBSDV-resistant genes. RBSDV resistance of a diverse global collection comprising 1,953 rice accessions was evaluated under natural conditions across 3 years. The average disease incidences of the Xian/indica (XI) subgroup were significantly lower than those of the Geng/japonica (GJ) subgroup. Interestingly, most XI-1A accessions in the Xian subgroup were significantly more susceptible than XI-1B accessions, even though XI-1A and XI-1B have a close phylogenetic relationship. Four Xian accessions stably and highly resistant to RBSDV were consistently identified in 2 years. Ten genomic regions (GRs) with 147 single nucleotide polymorphisms associated with RBSDV resistance were detected by a single-locus genome-wide association study (GWAS), of which five were repeatedly identified in a multilocus GWAS. Two previously reported GRs, grRBSDV-6.1 and grRBSDV-6.3, which were repeatedly detected as stably and highly associated with RBSDV resistance, contained 17 and seven genes, respectively, with significant differences of resistance among haplotypes. Haplotype analyses of the candidate genes LOC_Os06g03150 in grRBSDV-6.1 and LOC_Os06g31190 in grRBSDV-6.3 suggested that the former gene is mainly associated with the differentiation of resistance within the Xian subgroup and the latter gene mainly explains the difference in the resistance between Xian and Geng. Another three novel resistance GRs (grRBSDV-1.1, grRBSDV-7.1, and grRBSDV-9.1) were identified. Our findings may enhance the application of disease-resistant rice germplasms for breeding RBSDV-resistant varieties.

Population genetic structure and association mapping for iron toxicity tolerance in rice

Iron (Fe) toxicity is a major abiotic stress which severely reduces rice yield in many countries of the world. Genetic variation for this stress tolerance exists in rice germplasms. Mapping of gene(s)/QTL controlling the stress tolerance and transfer of the traits into high yielding rice varieties are essential for improvement against the stress. A panel population of 119 genotypes from 352 germplasm lines was constituted for detecting the candidate gene(s)/QTL through association mapping. STRUCTURE, GenAlEx and Darwin softwares were used to classify the population. The marker-trait association was detected by considering both the Generalized Linear Model (GLM) and Mixed Linear Model (MLM) analyses. Wide genetic variation was observed among the genotypes present in the panel population for the stress tolerance. Linkage disequilibrium was detected in the population for iron toxicity tolerance. The population was categorized into three genetic structure groups. Marker-trait association study considering both the Generalized Linear Model (GLM) and Mixed Linear Model (MLM) showed significant association of leaf browning index (LBI) with markers RM471, RM3, RM590 and RM243. Three novel QTL controlling Fe-toxicity tolerance were detected and designated as qFeTox4.3, qFeTox6.1 and qFeTox10.1. A QTL reported earlier in the marker interval of C955-C885 on chromosome 1 is validated using this panel population. The present study showed that QTL controlling Fe-toxicity tolerance to be co-localized with the QTL for Fe-biofortification of rice grain indicating involvement of common pathway for Fe toxicity tolerance and Fe content in rice grain. Fe-toxicity tolerance QTL qFeTox6.1 was co-localized with grain Fe-biofortification QTLs qFe6.1 and qFe6.2 on chromosome 6, whereas qFeTox10.1 was co-localized with qFe10.1 on chromosome 10. The Fe-toxicity tolerance QTL detected from this mapping study will be useful in marker-assisted breeding programs.

Genome-Wide Association Study Dissects Resistance Loci against Bacterial Blight in a Diverse Rice Panel from the 3000 Rice Genomes Project

BACKGROUND

Bacterial blight (BB), caused by Xanthomonas oryzae pv. oryzae (Xoo) is one of the most devastating bacterial diseases of rice in temperate and tropical regions. Breeding and deployment of resistant cultivars carrying major resistance (R) genes has been the most effective approach for BB management. However, because of specific interaction of each R gene with the product of the corresponding pathogen avirulence or effector gene, new pathogen strains that can overcome the deployed resistance often emerge rapidly. To deal with ever-evolving Xoo, it is necessary to identify novel R genes and resistance quantitative trait loci (QTL).

RESULTS

BB resistance of a diverse panel of 340 accessions from the 3000 Rice Genomes Project (3 K RGP) was evaluated by artificial inoculation with four representative Xoo strains, namely Z173 (C4), GD1358 (C5), V from China and PXO339 (P9a) from Philippines. Using the 3 K RG 4.8mio filtered SNP Dataset, a total of 11 QTL associated with BB resistance on chromosomes 4, 5, 11 and 12 were identified through a genome-wide association study (GWAS). Among them, eight resistance loci, which were narrowed down to relatively small genomic intervals, coincided with previously reported QTL or R genes, e.g. xa5, xa25, xa44(t). The other three QTL were putative novel loci associated with BB resistance. Linear regression analysis showed a dependence of BB lesion length on the number of favorable alleles, suggesting that pyramiding QTL using marker-assisted selection would be an effective approach for improving resistance. In addition, the Hap2 allele of LOC_Os11g46250 underlying qC5-11.1 was validated as positively regulating resistance against strain C5.

CONCLUSIONS

Our findings provide valuable information for the genetic improvement of BB resistance and application of germplasm resources in rice breeding programs.

A genome-wide association study reveals that the glucosyltransferase OsIAGLU regulates root growth in rice

Good root growth in the early post-germination stages is an important trait for direct seeding in rice, but its genetic control is poorly understood. In this study, we examined the genetic architecture of variation in primary root length using a diverse panel of 178 accessions. Four QTLs for root length (qRL3, qRL6, qRL7, and qRL11) were identified using genome-wide association studies. One candidate gene was validated for the major QTL qRL11, namely the glucosyltransferase OsIAGLU. Disruption of this gene in Osiaglu mutants reduced the primary root length and the numbers of lateral and crown roots. The natural allelic variations of OsIAGLU contributing to root growth were identified. Functional analysis revealed that OsIAGLU regulates root growth mainly via modulating multiple hormones in the roots, including levels of auxin, jasmonic acid, abscisic acid, and cytokinin. OsIAGLU also influences the expression of multiple hormone-related genes associated with root growth. The regulation of root growth through multiple hormone pathways by OsIAGLU makes it a potential target for future rice breeding for crop improvement.

Genome-wide association of rice response to blast fungus identifies loci for robust resistance under high nitrogen

BACKGROUND

Nitrogen fertilization is known to increase disease susceptibility, a phenomenon called Nitrogen-Induced Susceptibility (NIS). In rice, this phenomenon has been observed in infections with the blast fungus Magnaporthe oryzae. A previous classical genetic study revealed a locus (NIS1) that enhances susceptibility to rice blast under high nitrogen fertilization. In order to further address the underlying genetics of plasticity in susceptibility to rice blast after fertilization, we analyzed NIS under greenhouse-controlled conditions in a panel of 139 temperate japonica rice strains. A genome-wide association analysis was conducted to identify loci potentially involved in NIS by comparing susceptibility loci identified under high and low nitrogen conditions, an approach allowing for the identification of loci validated across different nitrogen environments. We also used a novel NIS Index to identify loci potentially contributing to plasticity in susceptibility under different nitrogen fertilization regimes.

RESULTS

A global NIS effect was observed in the population, with the density of lesions increasing by 8%, on average, under high nitrogen fertilization. Three new QTL, other than NIS1, were identified. A rare allele of the RRobN1 locus on chromosome 6 provides robust resistance in high and low nitrogen environments. A frequent allele of the NIS2 locus, on chromosome 5, exacerbates blast susceptibility under the high nitrogen condition. Finally, an allele of NIS3, on chromosome 10, buffers the increase of susceptibility arising from nitrogen fertilization but increases global levels of susceptibility. This allele is almost fixed in temperate japonicas, as a probable consequence of genetic hitchhiking with a locus involved in cold stress adaptation.

CONCLUSIONS

Our results extend to an entire rice subspecies the initial finding that nitrogen increases rice blast susceptibility. We demonstrate the usefulness of estimating plasticity for the identification of novel loci involved in the response of rice to the blast fungus under different nitrogen regimes.

Genetic loci regulating cadmium content in rice grains

It has been estimated that up to 90% of human exposure to cadmium is through food, and that cadmium within rice grains can be a major contributor to that dietary source. In this study genome wide association mapping was conducted on the Bengal and Assam Aus Panel (BAAP) of rice to identify quantitative trait loci and candidate genes for lowering grain cadmium. Field experiments were conducted over two years under two different irrigation systems: continually flooded and alternate wetting and drying (AWD). There was significant effects of water treatment, genotype, and genotype by water treatment interaction. Importantly, AWD increased grain cadmium, on average, by 49.6% and 108.8% in year 1 and 2 respectively. There was between 4.6 and 28 fold variation in cadmium concentration. A total of 58 QTLs were detected but no loci are clearly specific to one water regime despite approximately 20% of variation attributable to genotype by water regime interaction. A number of QTLs were consistent across most water treatments and years. These included QTLs on chromosome 7 (7.23-7.61, 8.93-9.04, and 29.12-29.14 Mbp), chromosome 5 (8.66-8.72 Mbp), and chromosome 9 (11.46-11.64 Mbp). Further analysis of the loci on chromosome 7 (8.93-9.04 Mbp), identified the candidate gene OsNRAMP1, where cultivars with a deletion upstream of the gene had higher concentrations of cadmium compared to the cultivars that did not have the deletion. The distribution of alleles within the BAAP suggest this QTL is easily detected in this population because it is composed of aus cultivars. Local genome cluster analysis suggest high Cd alleles are uncommon, but should be avoided in breeding.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at (10.1007/s10681-020-02752-1).

A quantitative genomics map of rice provides genetic insights and guides breeding

Extensive allelic variation in agronomically important genes serves as the basis of rice breeding. Here, we present a comprehensive map of rice quantitative trait nucleotides (QTNs) and inferred QTN effects based on eight genome-wide association study cohorts. Population genetic analyses revealed that domestication, local adaptation and heterosis are all associated with QTN allele frequency changes. A genome navigation system, RiceNavi, was developed for QTN pyramiding and breeding route optimization, and implemented in the improvement of a widely cultivated indica variety. This work presents an efficient platform that bridges ever-increasing genomic knowledge and diverse improvement needs in rice.

Genomic resources in plant breeding for sustainable agriculture

Climate change during the last 40 years has had a serious impact on agriculture and threatens global food and nutritional security. From over half a million plant species, cereals and legumes are the most important for food and nutritional security. Although systematic plant breeding has a relatively short history, conventional breeding coupled with advances in technology and crop management strategies has increased crop yields by 56 % globally between 1965-85, referred to as the Green Revolution. Nevertheless, increased demand for food, feed, fiber, and fuel necessitates the need to break existing yield barriers in many crop plants. In the first decade of the 21st century we witnessed rapid discovery, transformative technological development and declining costs of genomics technologies. In the second decade, the field turned towards making sense of the vast amount of genomic information and subsequently moved towards accurately predicting gene-to-phenotype associations and tailoring plants for climate resilience and global food security. In this review we focus on genomic resources, genome and germplasm sequencing, sequencing-based trait mapping, and genomics-assisted breeding approaches aimed at developing biotic stress resistant, abiotic stress tolerant and high nutrition varieties in six major cereals (rice, maize, wheat, barley, sorghum and pearl millet), and six major legumes (soybean, groundnut, cowpea, common bean, chickpea and pigeonpea). We further provide a perspective and way forward to use genomic breeding approaches including marker-assisted selection, marker-assisted backcrossing, haplotype based breeding and genomic prediction approaches coupled with machine learning and artificial intelligence, to speed breeding approaches. The overall goal is to accelerate genetic gains and deliver climate resilient and high nutrition crop varieties for sustainable agriculture.

SeedExtractor: An Open-Source GUI for Seed Image Analysis

Accurate measurement of seed size parameters is essential for both breeding efforts aimed at enhancing yields and basic research focused on discovering genetic components that regulate seed size. To address this need, we have developed an open-source graphical user interface (GUI) software, SeedExtractor that determines seed size and shape (including area, perimeter, length, width, circularity, and centroid), and seed color with capability to process a large number of images in a time-efficient manner. In this context, our application takes ∼2 s for analyzing an image, i.e., significantly less compared to the other tools. As this software is open-source, it can be modified by users to serve more specific needs. The adaptability of SeedExtractor was demonstrated by analyzing scanned seeds from multiple crops. We further validated the utility of this application by analyzing mature-rice seeds from 231 accessions in Rice Diversity Panel 1. The derived seed-size traits, such as seed length, width, were used for genome-wide association analysis. We identified known loci for regulating seed length (GS3) and width (qSW5/GW5) in rice, which demonstrates the accuracy of this application to extract seed phenotypes and accelerate trait discovery. In summary, we present a publicly available application that can be used to determine key yield-related traits in crops.

Plant-ImputeDB: an integrated multiple plant reference panel database for genotype imputation

Genotype imputation is a process that estimates missing genotypes in terms of the haplotypes and genotypes in a reference panel. It can effectively increase the density of single nucleotide polymorphisms (SNPs), boost the power to identify genetic association and promote the combination of genetic studies. However, there has been a lack of high-quality reference panels for most plants, which greatly hinders the application of genotype imputation. Here, we developed Plant-ImputeDB (http://gong_lab.hzau.edu.cn/Plant_imputeDB/), a comprehensive database with reference panels of 12 plant species for online genotype imputation, SNP and block search and free download. By integrating genotype data and whole-genome resequencing data of plants from various studies and databases, the current Plant-ImputeDB provides high-quality reference panels of 12 plant species, including ∼69.9 million SNPs from 34 244 samples. It also provides an easy-to-use online tool with the option of two popular tools specifically designed for genotype imputation. In addition, Plant-ImputeDB accepts submissions of different types of genomic variations, and provides free and open access to all publicly available data in support of related research worldwide. In general, Plant-ImputeDB may serve as an important resource for plant genotype imputation and greatly facilitate the research on plant genetic research.

Explore the genetics of weedy traits using rice 3K database

BACKGROUND

Weedy rice, a conspecific weedy counterpart of the cultivated rice (Oryza sativa L.), has been problematic in rice-production area worldwide. Although we started to know about the origin of some weedy traits for some rice-growing regions, an overall assessment of weedy trait-related loci was not yet available. On the other hand, the advances in sequencing technologies, together with community efforts, have made publicly available a large amount of genomic data. Given the availability of public data and the need of "weedy" allele mining for a better management of weedy rice, the objective of the present study was to explore the genetic architecture of weedy traits based on publicly available data, mainly from the 3000 Rice Genome Project (3K-RGP).

RESULTS

Based on the results of population structure analysis, we have selected 1378 individuals from four sub-populations (aus, indica, temperate japonica, tropical japonica) without admixed genomic composition for genome-wide association analysis (GWAS). Five traits were investigated: awn color, seed shattering, seed threshability, seed coat color, and seedling height. GWAS was conducted for each sub-population × trait combination and we have identified 66 population-specific trait-associated SNPs. Eleven significant SNPs fell into an annotated gene and four other SNPs were close to a putative candidate gene (± 25 kb). SNPs located in or close to Rc were particularly predictive of the occurrence of seed coat color and our results showed that different sub-populations required different SNPs for a better seed coat color prediction. We compared the data of 3K-RGP to a publicly available weedy rice dataset. The profile of allele frequency, phenotype-genotype segregation of target SNP, as well as GWAS results for the presence and absence of awns diverged between the two sets of data.

CONCLUSIONS

The genotype of trait-associated SNPs identified in this study, especially those located in or close to Rc, can be developed to diagnostic SNPs to trace the origin of weedy trait occurred in the field. The difference of results from the two publicly available datasets used in this study emphasized the importance of laboratory experiments to confirm the allele mining results based on publicly available data.

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

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

Dynamic genome-wide association analysis and identification of candidate genes involved in anaerobic germination tolerance in rice

BACKGROUND

In Asian rice production, an increasing number of countries now choose the direct seeding mode because of rising costs, labour shortages and water shortages. The ability of rice seeds to undergo anaerobic germination (AG) plays an important role in the success of direct seeding.

RESULTS

In this study, we used 2,123,725 single nucleotide polymorphism (SNP) markers based on resequencing to conduct a dynamic genome-wide association study (GWAS) of coleoptile length (CL) and coleoptile diameter (CD) in 209 natural rice populations. A total of 26 SNP loci were detected in these two phenotypes, of which 5 overlapped with previously reported loci (S1_ 39674301, S6_ 20797781, S7_ 18722403, S8_ 9946213, S11_ 19165397), and two sites were detected repeatedly at different time points (S3_ 24689629 and S5_ 27918754). We suggest that these 7 loci (-log₁₀ (P) value > 7.3271) are the key sites that affect AG tolerance. To screen the candidate genes more effectively, we sequenced the transcriptome of the flooding-tolerant variety R151 in six key stages, including anaerobic (AN) and the oxygen conversion point (AN-A), and obtained high-quality differential expression profiles. Four reliable candidate genes were identified: Os01g0911700 (OsVP1), Os05g0560900 (OsGA2ox8), Os05g0562200 (OsDi19-1) and Os06g0548200. Then qRT-PCR and LC-MS/ MS targeting metabolite detection technology were used to further verify that the up-regulated expression of these four candidate genes was closely related to AG.

CONCLUSION

The four novel candidate genes were associated with gibberellin (GA) and abscisic acid (ABA) regulation and cell wall metabolism under oxygen-deficiency conditions and promoted coleoptile elongation while avoiding adverse effects, allowing the coleoptile to obtain oxygen, escape the low-oxygen environment and germinate rapidly. The results of this study improve our understanding of the genetic basis of AG in rice seeds, which is conducive to the selection of flooding-tolerant varieties suitable for direct seeding.

Genome-Wide Association Study of Root and Shoot Related Traits in Spring Soybean (Glycine max L.) at Seedling Stages Using SLAF-Seq

Root systems can display variable genetic architectures leading to nutrient foraging or improving abiotic stress tolerance. Breeding for new soybean varieties with efficient root systems has tremendous potential in enhancing resource use efficiency and plant adaptation for challenging climates. In this study, root related traits were analyzed in a panel of 260 spring soybean with genome-wide association study (GWAS). Genotyping was done with specific locus amplified fragment sequencing (SLAF-seq), and five GWAS models (GLM, MLM, CMLM, FaST-LMM, and EMMAX) were used for analysis. A total of 179,960 highly consistent SNP markers distributed over the entire genome with an inter-marker distance of 2.36 kb was used for GWAS analysis. Overall, 27 significant SNPs with a phenotypic contribution ranging from 20 to 72% and distributed on chromosomes 2, 6, 8, 9, 13, 16 and 18 were identified and two of them were found to be associated with multiple root-related traits. Based on the linkage disequilibrium (LD) distance of 9.5 kb for the different chromosomes, 11 root and shoot regulating genes were detected based on LD region of a maximum 55-bp and phenotypic contribution greater than 22%. Expression analysis revealed an association between expression levels of those genes and the degree of root branching number. The current study provides new insights into the genetic architecture of soybean roots, and the underlying SNPs/genes could be critical for future breeding of high-efficient root system in soybean.

Assessing Rice Salinity Tolerance: From Phenomics to Association Mapping

Rice is the most salt-sensitive cereal, suffering yield losses above 50% with soil salinity of 6 dS/m. Thus, understanding the mechanisms of rice salinity tolerance is key to address food security. In this chapter, we provide guidelines to assess rice salinity tolerance using a high-throughput phenotyping platform (HTP) with digital imaging at seedling/early tillering stage and suggest improved analysis methods using stress indices. The protocols described here also include computer scripts for users to improve their experimental design, run genome-wide association studies (GWAS), perform multi-testing corrections, and obtain the Manhattan plots, enabling the identification of loci associated with salinity tolerance. Notably, the computer scripts provided here can be used for any stress or GWAS experiment and independently of HTP.

Identifying Anti-Oxidant Biosynthesis Genes in Pearl Millet [Pennisetum glaucum (L.) R. Br.] Using Genome-Wide Association Analysis

Pearl millet [Pennisetum glaucum (L.) R Br.] is an important staple food crop in the semi-arid tropics of Asia and Africa. It is a cereal grain that has the prospect to be used as a substitute for wheat flour for celiac patients. It is an important antioxidant food resource present with a wide range of phenolic compounds that are good sources of natural antioxidants. The present study aimed to identify the total antioxidant content of pearl millet flour and apply it to evaluate the antioxidant activity of its 222 genotypes drawn randomly from the pearl millet inbred germplasm association panel (PMiGAP), a world diversity panel of this crop. The total phenolic content (TPC) significantly correlated with DPPH (1,1-diphenyl-2-picrylhydrazyl) radical scavenging activity (% inhibition), which ranged from 2.32 to 112.45% and ferric-reducing antioxidant power (FRAP) activity ranging from 21.68 to 179.66 (mg ascorbic acid eq./100 g). Genome-wide association studies (GWAS) were conducted using 222 diverse accessions and 67 K SNPs distributed across all the seven pearl millet chromosomes. Approximately, 218 SNPs were found to be strongly associated with DPPH and FRAP activity at high confidence [-log (p) > 3.0-7.4]. Furthermore, flanking regions of significantly associated SNPs were explored for candidate gene harvesting. This identified 18 candidate genes related to antioxidant pathway genes (flavanone 7-O-beta-glycosyltransferase, GDSL esterase/lipase, glutathione S-transferase) residing within or near the association signal that can be selected for further functional characterization. Patterns of genetic variability and the associated genes reported in this study are useful findings, which would need further validation before their utilization in molecular breeding for high antioxidant-containing pearl millet cultivars.