Genome-wide association studies for agronomical traits in winter rice accessions of Assam

New insights into QTNs and potential candidate genes governing rice yield via a multi-model genome-wide association study

BACKGROUND

Rice (Oryza sativa L.) is one of the globally important staple food crops, and yield-related traits are prerequisites for improved breeding efficiency in rice. Here, we used six different genome-wide association study (GWAS) models for 198 accessions, with 553,229 single nucleotide markers (SNPs) to identify the quantitative trait nucleotides (QTNs) and candidate genes (CGs) governing rice yield.

RESULTS

Amongst the 73 different QTNs in total, 24 were co-localized with already reported QTLs or loci in previous mapping studies. We obtained fifteen significant QTNs, pathway analysis revealed 10 potential candidates within 100kb of these QTNs that are predicted to govern plant height, days to flowering, and plot yield in rice. Based on their superior allelic information in 20 elite and 6 inferior genotypes, we found a higher percentage of superior alleles in the elite genotypes in comparison to inferior genotypes. Further, we implemented expression analysis and enrichment analysis enabling the identification of 73 candidate genes and 25 homologues of Arabidopsis, 19 of which might regulate rice yield traits. Of these candidate genes, 40 CGs were found to be enriched in 60 GO terms of the studied traits for instance, positive regulator metabolic process (GO:0010929), intracellular part (GO:0031090), and nucleic acid binding (GO:0090079). Haplotype and phenotypic variation analysis confirmed that LOC_OS09G15770, LOC_OS02G36710 and LOC_OS02G17520 are key candidates associated with rice yield.

CONCLUSIONS

Overall, we foresee that the QTNs, putative candidates elucidated in the study could summarize the polygenic regulatory networks controlling rice yield and be useful for breeding high-yielding varieties.

Grain iron and zinc content is independent of anthocyanin accumulation in pigmented rice genotypes of Northeast region of India

The traditional rice genotypes of Assam are considered to have biological value due to the presence of several bioactive compounds like flavonoids, polyphenols, and anthocyanins, which have antioxidant, anti-cancer, anti-diabetic, and anti-aging properties. The pigmented genotypes are considered to have high iron (Fe) content. However, the effect of Fe and Zinc (Zn) accumulation on anthocyanin content is yet to be studied in pigmented rice of Assam. We studied the Fe, Zn, and anthocyanin content in grains of 204 traditional rice of Assam, which are traditionally preferred for their nutraceutical properties. We performed phenotypic and biochemical compositional analyses of 204 genotypes to identify those having high Fe, Zn, and anthocyanin. We also carried out the differential expression of a few selected Fe and Zn transporter genes along with the expression of anthocyanin biosynthesis genes. Interestingly, all pigmented rice genotypes contained a higher amount of phenolic compound than the non-pigmented form of rice. We found the highest (32.73 g) seed yield per plant for genotype Jengoni followed by Kajoli chokuwa and Khau Pakhi 1. We also listed 30 genotypes having high levels of Fe and Zn content. The genotype Jengoni accumulated the highest (186.9 μg g-1) Fe, while the highest Zn (119.9 μg g-1) content was measured in genotype Bora (Nagaon), The levels of Ferritin 2 gene expression were found to be significantly higher in Bora (Nagaon) (> 2-fold). For Zn accumulation, the genotype DRR Dhan-45, which was released as a high Zn content variety, showed significant up-regulation of the ZIP4 gene at booting (> 7-fold), post-anthesis (7.8-fold) and grain filling (> 5-fold) stages followed by Bora (Nagaon) (> 3-fold) at post-anthesis. Anthocyanidin synthase gene, Flavanone 3-dioxygenase 1-like (FDO1), and Chalcone-flavanone isomerase-like genes were up-regulated in highly pigmented genotype Bora (Nagaon) followed by Jengoni. Based on our data there was no significant correlation between iron and zinc content on the accumulation of anthocyanin. This challenges the present perception of the higher nutritive value in terms of the micronutrient content of the colored rice of Assam. This is the first report on the detailed characterization of traditional rice genotypes inclusive of phenotypic, biochemical, nutritional, and molecular attributes, which would be useful for designing the breeding program to improve Fe, Zn, or anthocyanin content in rice.

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.

Mapping of QTLs associated with yield and related traits under reproductive stage drought stress in rice using SNP linkage map

BACKGROUND

Drought stress is a major constraint for rice production worldwide. Reproductive stage drought stress (RSDS) leads to heavy yield losses in rice. The prospecting of new donor cultivars for identification and introgression of QTLs of major effect (Quantitative trait locus) for drought tolerance is crucial for the development of drought-resilient rice varieties.

METHODS AND RESULTS

Our study aimed to map QTLs associated with yield and its related traits under RSDS conditions. A saturated linkage map was constructed using 3417 GBS (Genotyping by sequencing) derived SNP (Single nucleotide polymorphism) markers spanning 1924.136 cM map length with an average marker density of 0.56 cM, in the F3 mapping population raised via cross made between the traditional ahu rice cultivar, Koniahu (drought tolerant) and a high-yielding variety, Disang (drought susceptible). Using the Inclusive composite interval mapping approach, 35 genomic regions governing yield and related traits were identified in pooled data from 198 F3 and F4 segregating lines evaluated for two consecutive seasons under both RSDS and irrigated control conditions. Of the 35 QTLs, 23 QTLs were identified under RSDS with LOD (Logarithm of odds) values ranging between 2.50 and 7.83 and PVE (phenotypic variance explained) values of 2.95-12.42%. Two major QTLs were found to be linked to plant height (qPH1.29) and number of filled grains per panicle (qNOG5.12) under RSDS. Five putative QTLs for grain yield namely, qGY2.00, qGY5.05, qGY6.16, qGY9.19, and qGY10.20 were identified within drought conditions. Fourteen QTL regions having ≤ 10 Mb QTL interval size were further analysed for candidate gene identification and a total of 4146 genes were detected out of these 2263 (54.63%) genes were annotated to at least one gene ontology (GO) term.

CONCLUSION

Several QTLs associated with grain yield and yield components and putative candidate genes were identified. The putative QTLs and candidate genes identified could be employed to augment drought resilience in rice after further validation through MAS strategies.

Genetic analysis for detection of genes associated to drought tolerance in rice accessions belonging to north east India

INTRODUCTION

The North East (NE) India is rich in biodiversity and also considered as the secondary centre for origin of rice. The NE rice accessions was characterized previously using genetic markers and morphological traits. Simultaneously, genome-wide association studies (GWAS) reveal significant marker-trait associations for the drought tolerance traits.

METHODS AND RESULTS

The genetic diversity and population structure of 296 NE rice accessions were studied using 96,712 single nucleotide polymorphism (SNP) markers distributed across 12 chromosomes. The accessions were clustered into two major sub-groups (SG). A total of 91 accessions were assembled as SG1 and 114 accessions as SG2, while the remaining 91 were admixture genotypes. A total of 200 genotypes belonging to different groups were phenotyped for yield component traits under drought and control conditions. The GWAS was performed to identify significant marker-trait associations (MTAs). Consequently, 47 MTAs were detected under drought, exhibiting 0.02-9.95% of phenotypic variance (P.V.). Whereas 58 MTAs were discovered under control conditions, showing a 0.01-9.74% contribution to the phenotype. Through in-silico mining of QTLs, 2999 genes were identified. Among these; only 22 genes were directly associated with stress response.

CONCLUSION

These QTLs/genes may be deployed for marker-assisted pyramiding to improve drought tolerance in popular drought susceptible rice varieties.

Molecular mapping of drought-responsive QTLs during the reproductive stage of rice using a GBS (genotyping-by-sequencing) based SNP linkage map

BACKGROUND

In rice, drought stress at reproductive stage drastically reduces yield, which in turn hampers farmer's efforts towards crop production. The majority of the rice varieties have resistance genes against several abiotic and biotic stresses. Therefore, the traditional landraces were studied to identify QTLs/candidate genes associated with drought tolerance.

METHODS AND RESULTS

A high-density SNP-based genetic map was constructed using a Genotyping-by-sequencing (GBS) approach. The recombinant inbred lines (RILs) derived from crossing 'Banglami × Ranjit' were used for QTL analysis. A total map length of 1306.424 cM was constructed, which had an average inter-marker distance of 0.281 cM. The phenotypic evaluation of F₆ and F₇ RILs were performed under drought stress and control conditions. A total of 42 QTLs were identified under drought stress and control conditions for yield component traits explaining 1.95-13.36% of the total phenotypic variance (PVE). Among these, 19 QTLs were identified under drought stress conditions, whereas 23 QTLs were located under control conditions. A total of 4 QTLs explained a PVE ≥ 10% which are considered as the major QTLs. Moreover, bioinformatics analysis revealed the presence of 6 candidate genes, which showed differential expression under drought and control conditions.

CONCLUSION

These QTLs/genes may be deployed for marker-assisted pyramiding to improve drought tolerance in the existing rice varieties.

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.

Identification and characterization of genes for drought tolerance in upland rice cultivar 'Banglami' of North East India

INTRODUCTION

Rice is a major crop in Assam, North East (NE) India. The rice accessions belonging to NE India possess unique traits of breeder's interest, i.e., tolerant to biotic and abiotic stresses. In the present research programme, the stress responsive genes were identified within the QTLs associated with drought tolerance. The differential expression profiling of genes were performed under drought stress and control conditions. Thus, the 'candidate genes' associated with drought tolerance were recognised and may be deployed in a breeding programme.

METHODS AND RESULTS

A drought-tolerant traditional rice cultivar, Banglami, was crossed with a high-yielding, drought-susceptible variety, Ranjit. The mapping population (F₄) was raised through the single seed descent (SSD) method and used in QTL analysis. Under drought stress, a total of 4752 genes were identified through in-silico mining of QTLs. Among these, only 21 genes primarily associated with the stress response. The maximum of four stress-responsive genes were located within the QTLs, qNOG12.1 and qGY1.1. However, under control conditions, 2088 genes were identified, out of which, only 15 were categorised as the major stress responsive genes. The functional characterization of genes recognized 24 different types of proteins. Among these, peroxidase and heat shock proteins (Hsp) are the principal proteins encoded during stress. In addition to that, OsbZIP23, inorganic pyrophosphatase, universal stress protein, serine threonine kinase, NADPH oxidoreductase, and proteins belonging to the ABC1 family were also produced during stress condition. The differential expression profiling showed a profound expression pattern of three candidate genes under drought stress condition, i.e., OsI_32199 (Ascorbate peroxidase), OsI_37694 (Universal stress protein) and OsI_32167 (Heat shock protein 81 - 1).

CONCLUSION

The novel candidate genes identified for drought tolerance, may be used in the breeding programme for the development of 'climate smart rice varieties'.

GWAS to spot candidate genes associated with grain quality traits in diverse rice accessions of North East India

INTRODUCTION

North East (NE) India is the second centre for the origin of rice and is enriched with a diverse collection of traditional rice accessions. These genotypes possess unique traits of breeding interest and are rich in grain nutritional and cooking qualities. Therefore, quantitative trait loci (QTLs)/genes associated with the various quality traits may be identified through genome-wide association studies (GWAS) and used in crop improvement programmes.

METHODS AND RESULTS

A pool of 526 unique rice accessions from Assam, North East (NE) India were characterized by using 9 grain-quality traits and grouped into 16 clusters. Among these, the highest number of 156 (29.65%) genotypes belongs to diverse phenotypic classes; Sali, Lahi, and Chokuwa were grouped into cluster 6. The first three principal components showed 54.76% of morphological variability with Eigenvalue >1. Genome-wide association studies (GWAS) was performed in 103 rice accessions using 42,446 SNP markers. A total of 11 significant marker-trait associations were detected for 5 grain-quality traits, explaining 0.22-8.86% of phenotypic variation (PV). In-silico mining of QTLs detected 'candidate genes' associated with the quality traits.

CONCLUSIONS

The phenotypic diversity among the 526 rice accessions of NE India was studied using grain quality traits and grouped into 16 significantly different clusters. The QTLs, or candidate genes identified for various grain quality traits, may be used in breeding programmes for the development of improved rice varieties.

Identification of Novel Quantitative Trait Nucleotides and Candidate Genes for Bacterial Wilt Resistance in Tobacco (Nicotiana tabacum L.) Using Genotyping-by-Sequencing and Multi-Locus Genome-Wide Association Studies

Tobacco bacterial wilt (TBW) is a devastating soil-borne disease threatening the yield and quality of tobacco. However, its genetic foundations are not fully understood. In this study, we identified 126,602 high-quality single-nucleotide polymorphisms (SNPs) in 94 tobacco accessions using genotyping-by-sequencing (GBS) and a 94.56 KB linkage disequilibrium (LD) decay rate for candidate gene selection. The population structure analysis revealed two subpopulations with 37 and 57 tobacco accessions. Four multi-locus genome-wide association study (ML-GWAS) approaches identified 142 quantitative trait nucleotides (QTNs) in E1-E4 and the best linear unbiased prediction (BLUP), explaining 0.49-22.52% phenotypic variance. Of these, 38 novel stable QTNs were identified across at least two environments/methods, and their alleles showed significant TBW-DI differences. The number of superior alleles associated with TBW resistance for each accession ranged from 4 to 24; eight accessions had more than 18 superior alleles. Based on TBW-resistant alleles, the five best cross combinations were predicted, including MC133 × Ruyuan No. 1 and CO258 × ROX28. We identified 52 candidate genes around 38 QTNs related to TBW resistance based on homologous functional annotation and KEGG enrichment analysis, e.g., CYCD3;2, BSK1, Nitab4.5_0000641g0050, Nitab4.5_0000929g0030. To the best of our knowledge, this is the first comprehensive study to identify QTNs, superior alleles, and their candidate genes for breeding TBW-resistant tobacco varieties. The results provide further insight into the genetic architecture, marker-assisted selection, and functional genomics of TBW resistance, improving future breeding efforts to increase crop productivity.