Genome-wide association study of important agronomic traits within a core collection of rice (Oryza sativa L.)

Genome wide association studies on seven yield-related traits of 183 rice varieties in Bangladesh

Rice genetic diversity is regulated by multiple genes and is largely dependent on various environmental factors. Uncovering the genetic variations associated with the diversity in rice populations is the key to breed stable and high yielding rice varieties. We performed genome wide association studies (GWASs) on seven rice yielding traits (grain length, grain width, grain weight, panicle length, leaf length, leaf width, and leaf angle) based on a population of 183 rice landraces of Bangladesh. Our GWASs reveal various chromosomal regions and candidate genes that are associated with different traits in Bangladeshi rice varieties. Noteworthy was the recurrent implication of chromosome 10 in all three grain-shape-related traits (grain length, grain width, and grain weight), indicating its pivotal role in shaping rice grain morphology. Our study also underscores the involvement of transposon gene families across these three traits. For leaf related traits, chromosome 10 was found to harbor regions that are significantly associated with leaf length and leaf width. The results of these association studies support previous findings as well as provide additional insights into the genetic diversity of rice. This is the first known GWAS study on various yield-related traits in the varieties of Oryza sativa available in Bangladesh-the fourth largest rice-producing country. We believe this study will accelerate rice genetics research and breeding stable high-yielding rice in Bangladesh.

Integrated meta-analysis and transcriptomics pinpoint genomic loci and novel candidate genes associated with submergence tolerance in rice

BACKGROUND

Due to rising costs, water shortages, and labour shortages, farmers across the globe now prefer a direct seeding approach. However, submergence stress remains a major bottleneck limiting the success of this approach in rice cultivation. The merger of accumulated rice genetic resources provides an opportunity to detect key genomic loci and candidate genes that influence the flooding tolerance of rice.

RESULTS

In the present study, a whole-genome meta-analysis was conducted on 120 quantitative trait loci (QTL) obtained from 16 independent QTL studies reported from 2004 to 2023. These QTL were confined to 18 meta-QTL (MQTL), and ten MQTL were successfully validated by independent genome-wide association studies from diverse natural populations. The mean confidence interval (CI) of the identified MQTL was 3.44 times narrower than the mean CI of the initial QTL. Moreover, four core MQTL loci with genetic distance less than 2 cM were obtained. By combining differentially expressed genes (DEG) from two transcriptome datasets with 858 candidate genes identified in the core MQTL regions, we found 38 common differentially expressed candidate genes (DECGs). In silico expression analysis of these DECGs led to the identification of 21 genes with high expression in embryo and coleoptile under submerged conditions. These DECGs encode proteins with known functions involved in submergence tolerance including WRKY, F-box, zinc fingers, glycosyltransferase, protein kinase, cytochrome P450, PP2C, hypoxia-responsive family, and DUF domain. By haplotype analysis, the 21 DECGs demonstrated distinct genetic differentiation and substantial genetic distance mainly between indica and japonica subspecies. Further, the MQTL7.1 was successfully validated using flanked marker S2329 on a set of genotypes with phenotypic variation.

CONCLUSION

This study provides a new perspective on understanding the genetic basis of submergence tolerance in rice. The identified MQTL and novel candidate genes lay the foundation for marker-assisted breeding/engineering of flooding-tolerant cultivars conducive to direct seeding.

Deciphering the Genetic Basis of Allelopathy in japonica Rice Cultivated in Temperate Regions Using a Genome-Wide Association Study

Allelopathy has been considered as a natural method of weed control. Despite the nature of allelochemical compounds has been studied, little is known about the genetic basis underlying allelopathy. However, it is known that rice exhibits diverse allelopathic potentials across varieties, and breeding for rice plants exhibiting allelopathic potential conferring an advantage against weeds in paddy fields would be highly desirable. Knowledge of the gene factors and the identification of the genomic regions responsible for allelopathy would facilitate breeding programs. Taking advantage of the existing genetic diversity in rice, particularly in temperate japonica rice, we conducted a comprehensive investigation into the genetic determinants that contribute to rice allelopathy. Employing Genome-Wide Association Study, we identified four Quantitative Trait Loci, with the most promising loci situated on chromosome 2 and 5. Subsequent inspection of the genes located within these QTLs revealed genes associated with the biosynthesis of secondary metabolites such as Phenylalanine Ammonia Lyase (PAL), a key enzyme in the synthesis of phenolic compounds, and two genes coding for R2R3-type MYB transcription factors. The identification of these two QTLs associated to allelopathy in rice provides a useful tool for further exploration and targeted breeding strategies.

Multi-model genome-wide association studies for appearance quality in rice

Improving the quality of the appearance of rice is critical to meet market acceptance. Mining putative quality-related genes has been geared towards the development of effective breeding approaches for rice. In the present study, two SL-GWAS (CMLM and MLM) and three ML-GWAS (FASTmrEMMA, mrMLM, and FASTmrMLM) genome-wide association studies were conducted in a subset of 3K-RGP consisting of 198 rice accessions with 553,831 SNP markers. A total of 594 SNP markers were identified using the mixed linear model method for grain quality traits. Additionally, 70 quantitative trait nucleotides (QTNs) detected by the ML-GWAS models were strongly associated with grain aroma (AR), head rice recovery (HRR, %), and percentage of grains with chalkiness (PGC, %). Finally, 39 QTNs were identified using single- and multi-locus GWAS methods. Among the 39 reliable QTNs, 20 novel QTNs were identified for the above-mentioned three quality-related traits. Based on annotation and previous studies, four functional candidate genes (LOC_Os01g66110, LOC_Os01g66140, LOC_Os07g44910, and LOC_Os02g14120) were found to influence AR, HRR (%), and PGC (%), which could be utilized in rice breeding to improve grain quality traits.

Genetic diversity and genome-wide association study of 13 agronomic traits in 977 Beta vulgaris L. germplasms

BACKGROUND

Sugar beet (Beta vulgaris L.) is an economically essential sugar crop worldwide. Its agronomic traits are highly diverse and phenotypically plastic, influencing taproot yield and quality. The National Beet Medium-term Gene Bank in China maintains more than 1700 beet germplasms with diverse countries of origin. However, it lacks detailed genetic background associated with morphological variability and diversity.

RESULTS

Here, a comprehensive genome-wide association study (GWAS) of 13 agronomic traits was conducted in a panel of 977 sugar beet accessions. Almost all phenotypic traits exhibited wide genetic diversity and high coefficient of variation (CV). A total of 170,750 high-quality single-nucleotide polymorphisms (SNPs) were obtained using the genotyping-by-sequencing (GBS). Neighbour-joining phylogenetic analysis, principal component analysis, population structure and kinship showed no obvious relationships among these genotypes based on subgroups or regional sources. GWAS was carried out using a mixed linear model, and 159 significant associations were detected for these traits. Within the 25 kb linkage disequilibrium decay of the associated markers, NRT1/PTR FAMILY 6.3 (BVRB_5g097760); nudix hydrolase 15 (BVRB_8g182070) and TRANSPORT INHIBITOR RESPONSE 1 (BVRB_8g181550); transcription factor MYB77 (BVRB_2g023500); and ethylene-responsive transcription factor ERF014 (BVRB_1g000090) were predicted to be strongly associated with the taproot traits of root groove depth (RGD); root shape (RS); crown size (CS); and flesh colour (FC), respectively. For the aboveground traits, UDP-glycosyltransferase 79B6 (BVRB_9g223780) and NAC domain-containing protein 7 (BVRB_5g097990); F-box protein At1g10780 (BVRB_6g140760); phosphate transporter PHO1 (BVRB_3g048660); F-box protein CPR1 (BVRB_8g181140); and transcription factor MYB77 (BVRB_2g023500) and alcohol acyltransferase 9 (BVRB_2g023460) might be associated with the hypocotyl colour (HC); plant type (PT); petiole length (PL); cotyledon size (C); and fascicled leaf type (FLT) of sugar beet, respectively. AP-2 complex subunit mu (BVRB_5g106130), trihelix transcription factor ASIL2 (BVRB_2g041790) and late embryogenesis abundant protein 18 (BVRB_5g106150) might be involved in pollen quantity (PQ) variation. The candidate genes extensively participated in hormone response, nitrogen and phosphorus transportation, secondary metabolism, fertilization and embryo maturation.

CONCLUSIONS

The genetic basis of agronomical traits is complicated in heterozygous diploid sugar beet. The putative valuable genes found in this study will help further elucidate the molecular mechanism of each phenotypic trait for beet breeding.

Seed Bio-priming of wheat with a novel bacterial strain to modulate drought stress in Daegu, South Korea

Wheat is one of the major cereal crop grown food worldwide and, therefore, plays has a key role in alleviating the global hunger crisis. The effects of drought stress can reduces crop yields by up to 50% globally. The use of drought-tolerant bacteria for biopriming can improve crop yields by countering the negative effects of drought stress on crop plants. Seed biopriming can reinforce the cellular defense responses to stresses via the stress memory mechanism, that its activates the antioxidant system and induces phytohormone production. In the present study, bacterial strains were isolated from rhizospheric soil taken from around the Artemisia plant at Pohang Beach, located near Daegu, in the South Korea Republic of Korea. Seventy-three isolates were screened for their growth-promoting attributes and biochemical characteristics. Among them, the bacterial strain SH-8 was selected preferred based on its plant growth-promoting bacterial traits, which are as follows: abscisic acid (ABA) concentration = 1.08 ± 0.05 ng/mL, phosphate-solubilizing index = 4.14 ± 0.30, and sucrose production = 0.61 ± 0.13 mg/mL. The novel strain SH-8 demonstrated high tolerance oxidative stress. The antioxidant analysis also showed that SH-8 contained significantly higher levels of catalase (CAT), superoxide dismutase (SOD), and ascorbic peroxidase (APX). The present study also quantified and determined the effects of biopriming wheat (Triticum aestivum) seeds with the novel strain SH-8. SH-8 was highly effective in enhancing the drought tolerance of bioprimed seeds; their drought tolerance and germination potential (GP) were increased by up to 20% and 60%, respectively, compared with those in the control group. The lowest level of impact caused by drought stress and the highest germination potential, seed vigor index (SVI), and germination energy (GE) (90%, 2160, and 80%, respectively), were recorded for seeds bioprimed with with SH-8. These results show that SH-8 enhances drought stress tolerance by up to 20%. Our study suggests that the novel rhizospheric bacterium SH-8 (gene accession number OM535901) is a valuable biostimulant that improves drought stress tolerance in wheat plants and has the potential to be used as a biofertilizer under drought conditions.

Genome-wide association study reveals novel genomic regions governing agronomic and grain quality traits and superior allelic combinations for Basmati rice improvement

BACKGROUND

Basmati is a speciality segment in the rice genepool characterised by explicit grain quality. For the want of suitable populations, genome-wide association study (GWAS) in Basmati rice has not been attempted.

MATERIALS

To address this gap, we have performed a GWAS on a panel of 172 elite Basmati multiparent population comprising of potential restorers and maintainers. Phenotypic data was generated for various agronomic and grain quality traits across seven different environments during two consecutive crop seasons. Based on the observed phenotypic variation, three agronomic traits namely, days to fifty per cent flowering, plant height and panicle length, and three grain quality traits namely, kernel length before cooking, length breadth ratio and kernel length after cooking were subjected to GWAS. Genotyped with 80K SNP array, the population was subjected to principal component analysis to stratify the underlying substructure and subjected to the association analysis using Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK) model.

RESULTS

We identified 32 unique MTAs including 11 robust MTAs for the agronomic traits and 25 unique MTAs including two robust MTAs for the grain quality traits. Six out of 13 robust MTAs were novel. By genome annotation, six candidate genes associated with the robust MTAs were identified. Further analysis of the allelic combinations of the robust MTAs enabled the identification of superior allelic combinations in the population. This information was utilized in selecting 77 elite Basmati rice genotypes from the panel.

CONCLUSION

This is the first ever GWAS study in Basmati rice which could generate valuable information usable for further breeding through marker assisted selection, including enhancing of heterosis.

Genome-Wide Association Study of Zinc Toxicity Tolerance within a Rice Core Collection (Oryza sativa L.)

Zinc (Zn) is an essential micronutrient for rice, but it is toxic at a high concentration, especially in acid soils. It is yet unknown which genes regulate Zn tolerance in rice. In the present study, a genome-wide association study (GWAS) was performed for Zn tolerance in rice at the seedling stage within a rice core collection, named Ting's core collection, which showed extensive phenotypic variations in Zn toxicity with high-density single-nucleotide polymorphisms (SNPs). A total of 7 and 19 quantitative trait loci (QTL) were detected using root elongation (RE) and relative root elongation (RRE) under high Zn toxicity, respectively. Among them, 24 QTL were novel, and qRRE15 was located in the same region where 3 QTL were reported previously. In addition, qRE4 and qRRE9 were identical. Furthermore, we found eight candidate genes that are involved in abiotic and biotic stress, immunity, cell expansion, and phosphate transport in the loci of qRRE8, qRRE9, and qRRE15. Moreover, four candidate genes, i.e., Os01g0200700, Os06g0621900, Os06g0493600, and Os06g0622700, were verified correlating to Zn tolerance in rice by quantitative real time-PCR (qRT-PCR). Taken together, these results provide significant insight into the genetic basis for Zn toxicity tolerance and tolerant germplasm for developing rice tolerance to Zn toxicity and improving rice production in Zn-contaminated soils.

Genetic diversity of wild rice accessions (Oryza rufipogon Griff.) in Guangdong and Hainan Provinces, China, and construction of a wild rice core collection

Oryza rufipogon Griff. is a valuable germplasm resource for rice genetic improvement. However, natural habitat loss has led to the erosion of the genetic diversity of wild rice populations. Genetic diversity analysis of O. rufipogon accessions and development of the core collection are crucial for conserving natural genetic diversity and providing novel traits for rice breeding. In the present study, we developed 1,592 SNPs by multiplex PCR and next-generation sequencing (NGS) technology and used them to genotype 998 O. rufipogon accessions from 14 agroclimatic zones in Guangdong and Hainan Provinces, China. These SNPs were mapped onto 12 chromosomes, and the average MAF value was 0.128 with a minimum of 0.01 and a maximum of 0.499. The O. rufipogon accessions were classified into ten groups. The mean Nei's diversity index and Shannon-Wiener index (I) were 0.187 and 0.308, respectively, in all populations, indicating that O. rufipogon accessions had rich genetic diversity. There were also differences in the genetic diversity of O. rufipogon resources in the 14 regions. Hainan populations possessed higher levels of genetic diversity, whereas the Guangzhou population had lower levels of genetic diversity than did the other populations. Phylogenetic analysis revealed that the genetic relationship among the distribution sites of O. rufipogon was closely related to geographical location. Based on genetic distance, a core collection of 299 accessions captured more than 99% of the genetic variation in the germplasm. This study provides insights into O. rufipogon conservation, and the constructed core collection provides valuable resources for future research and genomics-assisted breeding of rice.

Genome wide association mapping of yield and various desirable agronomic traits in Rice

BACKGROUND

Rice (Oryza sativa L.) is one of the staple foods worldwide. To feed the growing population, the improvement of rice cultivars is important. To make the improvement in the rice breeding program, it is imperative to understand the similarities and differences of the existing rice accessions to find out the genetic diversity. Previous studies demonstrated the existence of abundant elite genes in rice landraces. A genome-wide association study (GWAS) was performed for yield and yield related traits to find the genetic diversity.

DESIGN

Experimental study.

METHODS AND RESULTS

A total of 204 SSRs markers were used among 17 SSRs found to be located on each chromosome in the rice genome. The diversity was analyzed using different genetic characters i.e., the total number of alleles (TNA), polymorphic information content (PIC), and gene diversity by Power markers, and the values for each genetic character per marker ranged from 2 to 9, 0.332 to 0.887 and 0.423 to 0.900 respectively across the whole genome. The results of population structure identified four main groups. MTA identified several markers associated with many agronomically important traits. These results will be very useful for the selection of potential parents, recombinants, and MTAs that govern the improvements and developments of new high yielding rice varieties.

CONCLUSIONS

Analysis of diversity in germplasm is important for the improvement of cultivars in the breeding program. In the present study, the diversity was analyzed with different methods and found that enormous diversity was present in the studied rice germplasm. The structure analysis found the presence of 4 genetic groups in the existing germplasm. A total of 129 marker-trait associations (MTAs) have been found in this study.

Genotyping-by-Sequencing Based Investigation of Population Structure and Genome Wide Association Studies for Seven Agronomically Important Traits in a Set of 346 Oryza rufipogon Accessions

Being one of the most important staple dietary constituents globally, genetic enhancement of cultivated rice for yield, agronomically important traits is of substantial importance. Even though the climatic factors and crop management practices impact complex traits like yield immensely, the contribution of variation by underlying genetic factors surpasses them all. Previous studies have highlighted the importance of utilizing exotic germplasm, landraces in enhancing the diversity of gene pool, leading to better selections and thus superior cultivars. Thus, to fully exploit the potential of progenitor of Asian cultivated rice for productivity related traits, genome wide association study (GWAS) for seven agronomically important traits was conducted on a panel of 346 O. rufipogon accessions using a set of 15,083 high-quality single nucleotide polymorphic markers. The phenotypic data analysis indicated large continuous variation for all the traits under study, with a significant negative correlation observed between grain parameters and agronomic parameters like plant height, culm thickness. The presence of 74.28% admixtures in the panel as revealed by investigating population structure indicated the panel to be very poorly genetically differentiated, with rapid LD decay. The genome-wide association analyses revealed a total of 47 strong MTAs with 19 SNPs located in/close to previously reported QTL/genic regions providing a positive analytic proof for our studies. The allelic differences of significant MTAs were found to be statistically significant at 34 genomic regions. A total of 51 O. rufipogon accessions harboured combination of superior alleles and thus serve as potential candidates for accelerating rice breeding programs. The present study identified 27 novel SNPs to be significantly associated with different traits. Allelic differences between cultivated and wild rice at significant MTAs determined superior alleles to be absent at 12 positions implying substantial scope of improvement by their targeted introgression into cultivars. Introgression of novel significant genomic regions into breeder's pool would broaden the genetic base of cultivated rice, thus making the crop more resilient.

Omics-Facilitated Crop Improvement for Climate Resilience and Superior Nutritive Value

Novel crop improvement approaches, including those that facilitate for the exploitation of crop wild relatives and underutilized species harboring the much-needed natural allelic variation are indispensable if we are to develop climate-smart crops with enhanced abiotic and biotic stress tolerance, higher nutritive value, and superior traits of agronomic importance. Top among these approaches are the "omics" technologies, including genomics, transcriptomics, proteomics, metabolomics, phenomics, and their integration, whose deployment has been vital in revealing several key genes, proteins and metabolic pathways underlying numerous traits of agronomic importance, and aiding marker-assisted breeding in major crop species. Here, citing several relevant examples, we appraise our understanding on the recent developments in omics technologies and how they are driving our quest to breed climate resilient crops. Large-scale genome resequencing, pan-genomes and genome-wide association studies are aiding the identification and analysis of species-level genome variations, whilst RNA-sequencing driven transcriptomics has provided unprecedented opportunities for conducting crop abiotic and biotic stress response studies. Meanwhile, single cell transcriptomics is slowly becoming an indispensable tool for decoding cell-specific stress responses, although several technical and experimental design challenges still need to be resolved. Additionally, the refinement of the conventional techniques and advent of modern, high-resolution proteomics technologies necessitated a gradual shift from the general descriptive studies of plant protein abundances to large scale analysis of protein-metabolite interactions. Especially, metabolomics is currently receiving special attention, owing to the role metabolites play as metabolic intermediates and close links to the phenotypic expression. Further, high throughput phenomics applications are driving the targeting of new research domains such as root system architecture analysis, and exploration of plant root-associated microbes for improved crop health and climate resilience. Overall, coupling these multi-omics technologies to modern plant breeding and genetic engineering methods ensures an all-encompassing approach to developing nutritionally-rich and climate-smart crops whose productivity can sustainably and sufficiently meet the current and future food, nutrition and energy demands.

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.

Molecular and genetic dissection of the USDA rice mini-core collection using high-density SNP markers

Molecular tools and knowledge of crop germplasm are vital for their effective utilization. In this study, we developed 40,866 high-quality and well distributed SNPs for a rice mini-core collection (RMC) developed by the United States Department of Agriculture (USDA). The high-quality SNPs clustered the USDA-RMC into five subpopulations (Ind, indica; Aus, aus; Afr, African rice; TeJ, temperate japonica; TrJ, tropical japonica) and one admixture (Adm). This classification was further confirmed by phylogenetic and principal component analyses. The rice ARO (aromatic) subpopulation of previous studies was re-assigned with Adm and the WD (wild-type) subpopulation was re-defined to the Afr subpopulation because most of its accessions are African cultivated rice. The Aus and Ind subpopulations had a substantially wider genetic variation than the TrJ and TeJ subpopulations. The genetic diversities were much larger between the Ind or Aus subpopulation and the TrJ or TeJ subpopulation than between the Afr subpopulation and the Ind, Aus, TrJ or TeJ subpopulation. Comparative agronomic trait analysis between the subpopulations also supported the genetic structure and variation of the RMC, and suggested the existence of extensive variation in the genes controlling agronomic traits among them. Furthermore, analysis of ancestral membership of the RMC accessions revealed that reproductive barrier or wide incompatibility existed between the Indica and Japonica groups, while gene flow occurred between them. These results provide high-quality SNPs and knowledge of genetic structure and diversity of the USDA-RMC necessary for enhanced rice research and breeding.

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.

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.

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 mapping revealed numerous novel genomic loci for grain nutritional and yield-related traits in rice (Oryza sativa L.) landraces

A core set of 190 rice landraces were used to decipher the genetic structure and to discover the chromosomal regions containing QTLs, affecting the grain micro-nutrients, fatty acids, and yield-related traits by using 148 molecular markers in this study. Landraces were categorized into three sub-groups based on population stratification study and followed by neighbor-joining tree and principal component analysis. Analysis of variance revealed abundant variations among the landraces for studied traits with less influence of environmental factors. Genome Wide Association Studies (GWAS) revealed 22 significant and consistent QTLs through marker trait association (MTAs) for 12 traits based on 2 years and pooled analysis. Out of 22 QTLs, three have been reported earlier while 19 QTLs are novel. Interestingly, 13 QTLs out of 22 were explained more than 10% phenotypic variance. Association of RM1148 and RM205 with Days to 50% flowering was comparable with flowering control genes Ghd8/qDTH8 and qDTH9, respectively. Similarly, Zn content was associated with RM44, which is situated within the QTL qZn8-1. Moreover, significant association of RM25 with oleic acid content was closely positioned with QTL qOle8. Association of RM7434 with grain yield/plant; RM184 with spikelet fertility %; R3M10, R9M42 with hundred seed weight; RM536, RM17467, RM484, RM26063 with Fe content; RM44, RM6839 with Zn content are the major outcomes of this study. In addition, association of R11M23 with days to 50% flowering, panicle length and total spikelets per panicle are explained the possible occurrence of pleiotropism among these traits. Prominent rice landraces viz., Anjani (early maturity); Sihar (extra dwarf); Gangabaru (highest grain yield/plant); Karhani (highest iron content); Byalo-2 (highest zinc content) and Kadamphool (highest oleic acid) were identified through this study. The present study will open many avenues towards utilization of these QTLs and superior landraces in rice breeding for developing nutrition-rich high yielding varieties.

Mining candidate gene for rice aluminum tolerance through genome wide association study and transcriptomic analysis

BACKGROUND

The genetic mechanism of aluminum (Al) tolerance in rice is great complicated. Uncovering genetic mechanism of Al tolerance in rice is the premise for Al tolerance improvement. Mining elite genes within rice landrace is of importance for improvement of Al tolerance in rice.

RESULTS

Genome-wide association study (GWAS) performed in EMMAX for rice Al tolerance was carried out using 150 varieties of Ting's core collection constructed from 2262 Ting's collections with more than 3.8 million SNPs. Within Ting's core collection of clear population structure and kinship relatedness as well as high rate of linkage disequilibrium (LD) decay, 17 genes relating to rice Al tolerance including cloned genes like NRAT1, ART1 and STAR1 were identified in this study. Moreover, 13 new candidate regions with high LD and 69 new candidate genes were detected. Furthermore, 20 of 69 new candidate genes were detected with significant difference between Al treatment and without Al toxicity by transcriptome sequencing. Interestingly, both qRT-PCR and sequence analysis in CDS region demonstrated that the candidate genes in present study might play important roles in rice Al tolerance.

CONCLUSIONS

The present study provided important information for further using these elite genes existing in Ting's core collection for improvement of rice Al tolerance.