Genetic Diversity and Elite Allele Mining for Grain Traits in Rice (Oryza sativa L.) by Association Mapping

Genetic Diversity and Association Mapping of Grain-Size Traits in Rice Landraces from the Honghe Hani Rice Terraces System in Yunnan Province

The Honghe Hani Rice Terraces System (HHRTS) of Yunnan Province is an important agricultural and cultural heritage landscape. Until now, a large number of local rice landraces have been planted. Mining excellent genes contained in these landraces provides a reference for variety improvement and new variety breeding. In this study, 96 rice landraces collected from the Hani terraces were planted in Honghe Mengzi, Yunnan Province, in 2013, 2014, 2015, and 2021, and five major grain traits were measured and analyzed. The genomic variation of 96 rice landraces was scanned by 201 simple sequence repeat (SSR) markers. The genetic diversity, population structure, and genetic relationships of the natural population were analyzed. The mixed linear model (MLM) method of the TASSEL software was used to analyze the associations between markers and traits. A total of 936 alleles were amplified by 201 pairs of SSR primers. The average number of observed alleles (Na), the effective number of alleles (Ne), Shannon's information index (I), heterozygosity (H), and the polymorphism information content (PIC) per marker were 4.66, 2.71, 1.08, 0.15, and 0.55, respectively. Ninety-six landraces were divided into two groups by population structure, clustering, and principal component analysis, and indica rice was the main group. The coefficients of variation of the five traits ranged from 6.80 to 15.24%, and their broad heritabilities were more than 70%. In addition, there were positive correlations among the same grain traits between different years. Through MLM analysis, 2, 36, 7, 7, and 4 SSR markers were significantly associated with grain length (GL), grain width (GW), grain thickness (GT), grain length-width ratio (LWR), and thousand-grain weight (TGW), respectively. The explanation rates of phenotypic variation were 16.31 (RM449, Chr. 1)-23.51% (RM316, Chr. 9), 10.84 (RM523, Chr. 3; RM161/RM305, Chr. 5)-43.01% (RM5496, Chr. 1), 11.98 (RM161/RM305, Chr. 5)-24.72% (RM275, Chr. 6), 12.68 (RM126, Chr. 8)-36.96% (RM5496, Chr. 1), and 17.65 (RM4499, Chr. 2)-26.32% (RM25, Chr. 8), respectively. The associated markers were distributed on 12 chromosomes of the genome.

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.

Using aquatic animals as partners to increase yield and maintain soil nitrogen in the paddy ecosystems

Whether species coculture can overcome the shortcomings of crop monoculture requires additional study. Here, we show how aquatic animals (i.e. carp, crabs, and softshell turtles) benefit paddy ecosystems when cocultured with rice. Three separate field experiments and three separate mesocosm experiments were conducted. Each experiment included a rice monoculture (RM) treatment and a rice-aquatic animal (RA) coculture treatment; RA included feed addition for aquatic animals. In the field experiments, rice yield was higher with RA than with RM, and RA also produced aquatic animal yields that averaged 0.52–2.57 t ha^-1. Compared to their corresponding RMs, the three RAs had significantly higher apparent nitrogen (N)-use efficiency and lower weed infestation, while soil N contents were stable over time. Dietary reconstruction analysis based on ¹³C and ¹⁵N showed that 16.0–50.2% of aquatic animal foods were from naturally occurring organisms in the rice fields. Stable-isotope-labeling (¹³C) in the field experiments indicated that the organic matter decomposition rate was greater with RA than with RM. Isotope ¹⁵N labeling in the mesocosm experiments indicated that rice used 13.0–35.1% of the aquatic animal feed-N. All these results suggest that rice-aquatic animal coculture increases food production, increases N-use efficiency, and maintains soil N content by reducing weeds and promoting decomposition and complementary N use. Our study supports the view that adding species to monocultures may enhance agroecosystem functions.

Monoculture, where only one type of crop is grown to the exclusion of any other organism, is a pillar of modern agriculture. Yet this narrow focus disregards how complex inter-species interactions can increase crop yield and biodiversity while decreasing the need for fertilizers or pesticides. For example, many farmers across Asia introduce carps, crabs, turtles or other freshwater grazers into their rice paddies. This coculture approach yields promising results but remains poorly understood. In particular, it is unclear how these animals’ behaviours and biological processes benefit the ecosystem.

To examine these questions, Guo, Zhao et al. conducted three separate four-year field experiments; they compared rice plots inhabited by either carp, mitten crabs or Chinese softshell turtles with fields where these organisms were not present.

With animals, the rice paddies had less weeds, better crop yields and steady levels of nitrogen (a natural fertiliser) in their soil. These ecosystems could breakdown organic matter faster, use it better and had a reduced need for added fertilizer. While animal feed was provided in the areas that were studied, carp, crabs and turtles obtained up to half their food from the field itself, eating weeds, algae and pests and therefore reducing competition for the crops.

This work helps to understand the importance of species interactions, showing that diversifying monocultures may boost yields and make agriculture more sustainable.

Genetic architecture of Al3+ toxicity tolerance in rice F2:3 populations determined through QTL mapping

Aluminum (Al³⁺) toxicity is one of the factors limiting crop production in acidic soils. Identifying quantitative trait loci (QTLs)/genes for tolerance to Al³⁺ toxicity at seed germination can aid the development of new tolerant cultivars. The segregating population derived from Pak Basmati (Indica) × Pokkali (Indica) was used for mapping QTLs linked with tolerance to Al³⁺ toxicity ranging from 0 to 20 mM at pH 4 ± 0.2 at germination. The favorable alleles for all new QTLs were analyzed based on germination traits, i.e., final germination percentage (FG%), germination energy (GE), germination speed (GS), germination index (GI), mean germination time (MGT), germination value (GV), germination velocity (GVe), peak value of germination (GPV), and germination capacity (GC), and growth traits, such as root length (RL), shoot length (SL), total dry biomass (TDB) and germination vigor index (GVI). The phenotypic evolution showed transgressive variations. For genome-wide mapping, 90 polymorphic SSRs with 4 gene-specific markers and Win QTL Cart were used for QTL analysis. In all, 35 QTLs for germination and 11 QTLs for seedling growth were detected in distinct chromosomal regions by composite interval mapping (CIM), and multiple interval mapping (MIM) confirmed the pleiotropy at region RM128 on chromosome 1. Based on our genetic mapping studies, the genes/QTLs underlying tolerance to Al³⁺ toxicity could differ for both the germination and seedling stages in segregated populations. The QTLs identified in this study could be a source of new alleles for improving tolerance to Al³⁺ toxicity in rice.

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.

Identification of QTLs for high grain yield and component traits in new plant types of rice

A panel of 60 genotypes comprising New Plant Types (NPTs) along with indica, tropical and temperate japonica genotypes was phenotypically evaluated for four seasons in irrigated situation for grain yield per se and component traits. Twenty NPT genotypes were found promising with an average grain yield varying from 5.45 to 8.8 t/ha. A total of 85 SSR markers were used in the study to identify QTLs associated with grain yield per se and related traits. Sixty-six (77.65%) markers were found to be polymorphic. The PIC values varied from 0.516 to 0.92 with an average of 0.704. A moderate level of genetic diversity (0.39) was detected among genotypes. Variation to the tune of 8% within genotypes, 68% among the genotypes within the population and 24% among the populations were observed (AMOVA). This information may help in identification of potential parents for development of transgressive segregants with very high yield. The association analysis using GLM and MLM models led to the identification of 30 and 10 SSR markers associated with 70 and 16 QTLs, respectively. Thirty novel QTLs linked with 16 SSRs were identified to be associated with eleven traits, namely tiller number (qTL-6.1, qTL-11.1, qTL-4.1), panicle length (qPL-1.1, qPL-5.1, qPL-7.1, qPL-8.1), flag leaf length (qFLL-8.1, qFLL-9.1), flag leaf width (qFLW-6.2, qFLW-5.1, qFLW-8.1, qFLW-7.1), total no. of grains (qTG-2.2, qTG-a7.1), thousand-grain weight (qTGW-a1.1, qTGW-a9.2, qTGW-5.1, qTGW-8.1), fertile grains (qFG-7.1), seed length-breadth ratio (qSlb-3.1), plant height (qPHT-6.1, qPHT-9.1), days to 50% flowering (qFD-1.1) and grain yield per se (qYLD-5.1, qYLD-6.1a, qYLD-11.1).Some of the SSRs were co-localized with more than two traits. The highest co-localization was identified with RM5709 linked to nine traits, followed by RM297 with five traits. Similarly, RM5575, RM204, RM168, RM112, RM26499 and RM22899 were also recorded to be co-localized with more than one trait and could be rated as important for marker-assisted backcross breeding programs, for pyramiding of these QTLs for important yield traits, to produce new-generation rice for prospective increment in yield potentiality and breaking yield ceiling.

Genome-Wide Association Study Reveals Novel Marker-Trait Associations (MTAs) Governing the Localization of Fe and Zn in the Rice Grain

Micronutrient malnutrition due to Fe and Zn, affects around two billion people globally particularly in the developing countries. More than 90% of the Asian population is dependent on rice-based diets, which is low in these micronutrients. In the present study, a set of 192 Indian rice germplasm accessions, grown at two locations, were evaluated for Fe and Zn in brown rice (BR) and milled rice (MR). A significant variation was observed in the rice germplasm for these micronutrients. The grain Fe concentration was in the range of 6.2–23.1 ppm in BR and 0.8–12.3 ppm in MR, while grain Zn concentration was found to be in the range of 11.0–47.0 ppm and 8.2–40.8 ppm in the BR and MR, respectively. Grain Fe exhibited maximum loss upon milling with a mean retention of 24.9% in MR, while Zn showed a greater mean retention of 74.2% in MR. A genome-wide association study (GWAS) was carried out implementing the FarmCPU model to control the population structure and kinship, and resulted in the identification of 29 marker-trait associations (MTAs) with significant associations for traits viz. FeBR (6 MTAs), FeMR (7 MTAs), ZnBR (11 MTAs), and ZnMR (5 MTAs), which could explain the phenotypic variance from 2.1 to as high as 53.3%. The MTAs governing the correlated traits showed co-localization, signifying the possibility of their simultaneous improvement. The robust MTAs identified in the study could be valuable resource for enhancing Fe and Zn concentration in the rice grain and addressing the problem of Fe and Zn malnutrition among rice consumers.

Variability Assessment for Root and Drought Tolerance Traits and Genetic Diversity Analysis of Rice Germplasm using SSR Markers

The studies on genetic variation, diversity and population structure of rice germplasm of North East India could be an important step for improvements of abiotic and biotic stress tolerance in rice. Genetic diversity and genetic relatedness among 114 rice genotypes of North East India were assessed using genotypic data of 65 SSR markers and phenotypic data. The phenotypic diversity analysis showed the considerable variation across genotypes for root, shoot and drought tolerance traits. The principal component analysis (PCA) revealed the fresh shoot weight, root volume, dry shoot weight, fresh root weight and drought score as a major contributor to diversity. Genotyping of 114 rice genotypes using 65 SSR markers detected 147 alleles with the average polymorphic information content (PIC) value of 0.51. Population structure analysis using the Bayesian clustering model approach, distance-based neighbor-joining cluster and principal coordinate analysis using genotypic data grouped the accession into three sub-populations. Population structure analysis revealed that rice accession was moderately structured based on FST value estimates. Analysis of molecular variance (AMOVA) and pairwise FST values showed significant differentiation among all the pairs of sub-population ranging from 0.152 to 0.222 suggesting that all the three subpopulations were significantly different from each other. AMOVA revealed that most of the variation in rice accession mainly occurred among individuals. The present study suggests that diverse germplasm of NE India could be used for the improvement of root and drought tolerance in rice breeding programmes.

InDel Marker Based Estimation of Multi-Gene Allele Contribution and Genetic Variations for Grain Size and Weight in Rice (Oryza sativa L.)

The market success of any rice cultivar is exceedingly dependent on its grain appearance, as well as its grain yield, which define its demand by consumers as well as growers. The present study was undertaken to explore the contribution of nine major genes, qPE9~1, GW2, SLG7, GW5, GS3, GS7, GW8, GS5, and GS2, in regulating four size and weight related traits, i.e., grain length (GL), grain width (GW), grain thickness (GT), and thousand grain weight (TGW) in 204 diverse rice germplasms using Insertion/Deletion (InDel) markers. The studied germplasm displayed wide-ranging variability in the four studied traits. Except for three genes, all six genes showed considerable association with these traits with varying strengths. Whole germplasm of 204 genotypes could be categorized into three major clusters with different grain sizes and weights that could be utilized in rice breeding programs where grain appearance and weight are under consideration. The study revealed that TGW was 24.9% influenced by GL, 37.4% influenced by GW, and 49.1% influenced by GT. Hence, assuming the trend of trait selection, i.e., GT > GW > GL, for improving TGW in the rice yield enhancement programs. The InDel markers successfully identified a total of 38 alleles, out of which 27 alleles were major and were found in more than 20 genotypes. GL was associated with four genes (GS3, GS7, GW8, and GS2). GT was also found to be regulated by four different genes (GS3, GS7, GW8, and GS2) out of the nine studied genes. GW was found to be under the control of three studied genes (GW5, GW8, and GS2), whereas TGW was found to be under the influence of four genes (SLG7, GW5, GW8, and GS5) in the germplasm under study. The Unweighted Pair Group Method with Arithmetic means (UPGMA) tree based on the studied InDel marker loci segregated the whole germplasm into three distinct clusters with dissimilar grain sizes and weights. A two-dimensional scatter plot constructed using Principal Coordinate Analysis (PCoA) based on InDel markers further separated the 204 rice germplasms into four sub-populations with prominent demarcations of extra-long, long, medium, and short grain type germplasms that can be utilized in breeding programs accordingly. The present study could help rice breeders to select a suitable InDel marker and in formulation of breeding strategies for improving grain appearance, as well as weight, to develop rice varieties to compete international market demands with higher yield returns. This study also confirms the efficient application of InDel markers in studying diverse types of rice germplasm, allelic frequencies, multiple-gene allele contributions, marker-trait associations, and genetic variations that can be explored further.

Genetic Mapping and Validation of Quantitative Trait Loci (QTL) for the Grain Appearance and Quality Traits in Rice (Oryza sativa L.) by Using Recombinant Inbred Line (RIL) Population

Rice grain shape and nutritional quality traits have high economic value for commercial production of rice and largely determine the market price, besides influencing the global food demand for high-quality rice. In order to understand the genetic components of grain appearance traits in paddy, brown, and head rice, 15 traits were evaluated by using 157 recombinant inbred lines (RILs) derived from a cross between two Iranian rice cultivars Ali-Kazemi (A) and Kadous (K). A significant variation was observed and showed transgressive segregation among the RILs. Correlations between the visual appearances of grain traits were studied. A linkage map with 65 polymorphic SSR markers was constructed, which covered 1517.32 cM of the rice genome. A total of seven QTLs were identified on four chromosomes, 1, 6, 9, and 12, associated with four traits, which are explained by the total phenotypic variation of 44.27% and LOD score of 32.77 in 2014 and 2015, respectively. Among these, four QTLs for two traits were consistently flanked by RM23904 and RM24432 on chromosome 9. Single QTL for head grain length (HGL) expressed in both the years on chromosomes 1 and 9. A major QTL for seed weight was detected on chromosome 9, which explained 10.18% of the phenotypic variation. The additive effect of all the QTLs was positively contributed by Ali-Kazemi allele, except one QTL on chromosome 6 (qBGL_6) that showed a negative additive effect being contributed by the Kadous allele. The study also validated the identified QTLs with the polymorphic SSR markers that were previously reported. Novel QTLs were identified on chromosomes 6 and 9, and many of the polymorphic markers were found to be associated with milling processing of grain quality, cooking, and nutraceutical properties of rice by extensive literature and database analysis. Therefore, these validated QTLs and marker information could be utilized in the marker-assisted selection to improve grain appearance and nutritional grain quality traits in rice.

The first identification of genomic loci in plants associated with resistance to galling insects: a case study in Eucalyptus L'Hér. (Myrtaceae)

Genomic loci related with resistance to gall-inducing insects have not been identified in any plants. Here, association mapping was used to identify molecular markers for resistance to the gall wasp Leptocybe invasa in two Eucalyptus species. A total of 86 simple sequence repeats (SSR) markers were screened out from 839 SSRs and used for association mapping in E. grandis. By applying the mixed linear model, seven markers were identified to be associated significantly (P ≤ 0.05) with the gall wasp resistance in E. grandis, including two validated with a correction of permutation test (P ≤ 0.008). The proportion of the variance in resistance explained by a significant marker ranged from 3.3% to 37.8%. Four out of the seven significant associations in E. grandis were verified and also validated (P ≤ 0.073 in a permutation test) in E. tereticornis, with the variation explained ranging from 24.3% to 48.5%. Favourable alleles with positive effect were also mined from the significant markers in both species. These results provide insight into the genetic control of gall wasp resistance in plants and have great potential for marker-assisted selection for resistance to L. invasa in the important tree genus Eucalyptus.

The first identification of genomic loci in plants associated with resistance to galling insects: a case study in Eucalyptus L'Hér. (Myrtaceae)

Genomic loci related with resistance to gall-inducing insects have not been identified in any plants. Here, association mapping was used to identify molecular markers for resistance to the gall wasp Leptocybe invasa in two Eucalyptus species. A total of 86 simple sequence repeats (SSR) markers were screened out from 839 SSRs and used for association mapping in E. grandis. By applying the mixed linear model, seven markers were identified to be associated significantly (P ≤ 0.05) with the gall wasp resistance in E. grandis, including two validated with a correction of permutation test (P ≤ 0.008). The proportion of the variance in resistance explained by a significant marker ranged from 3.3% to 37.8%. Four out of the seven significant associations in E. grandis were verified and also validated (P ≤ 0.073 in a permutation test) in E. tereticornis, with the variation explained ranging from 24.3% to 48.5%. Favourable alleles with positive effect were also mined from the significant markers in both species. These results provide insight into the genetic control of gall wasp resistance in plants and have great potential for marker-assisted selection for resistance to L. invasa in the important tree genus Eucalyptus.

Identification of putative markers linked to grain plumpness in rice (Oryza sativa L.) via association mapping

BACKGROUND

Poor grain plumpness (GP) is one of the main constraints to reaching the yield potential of hybrid rice.

RESULTS

In this study, the GP of 177 rice varieties was investigated in three locations across 2 years. By combining the genotype data of 261 simple sequence repeat (SSR) markers, association mapping was conducted to identify the marker-GP association loci. Among 31 marker-GP association loci detected in two or more environments and determined using general linear model (GLM) analyses, seven association loci were also detected using mixed linear model (MLM) analyses. The seven common loci detected by the two analytical methods were located on chromosomes 2, 3 (2), 7, 8 and 12 (2) and explained 7.24~22.28% of the variance. Of these 7 association loci, five markers linked to GP were newly detected: RM5340 on Chr2, RM5480 and RM148 on Chr3, RM1235 on Chr8, and RM5479 on Chr12.

CONCLUSIONS

Five marker-GP association loci were newly detected using both the GLM and MLM analytical methods. Elite allele RM505-170 bp had the highest average phenotypic effect on increasing the GP, and the typical carrier variety was 'Maozitou'. Based on the distribution of the elite alleles among the carrier varieties, the top 10 parental combinations for improving the GP in rice via cross-breeding were predicted.

New Candidate Genes Affecting Rice Grain Appearance and Milling Quality Detected by Genome-Wide and Gene-Based Association Analyses

Appearance and milling quality are two crucial properties of rice grains affecting its market acceptability. Understanding the genetic base of rice grain quality could considerably improve the high quality breeding. Here, we carried out an association analysis to identify QTL affecting nine rice grain appearance and milling quality traits using a diverse panel of 258 accessions selected from 3K Rice Genome Project and evaluated in two environments Sanya and Shenzhen. Genome-wide association analyses using 22,488 high quality SNPs identified 72 QTL affecting the nine traits. Combined gene-based association and haplotype analyses plus functional annotation allowed us to shortlist 19 candidate genes for seven important QTL regions affecting the grain quality traits, including two cloned genes (GS3 and TUD), two fine mapped QTL (qGRL7.1 and qPGWC7) and three newly identified QTL (qGL3.4, qGW1.1, and qGW10.2). The most likely candidate gene(s) for each important QTL were also discussed. This research demonstrated the superior power to shortlist candidate genes affecting complex phenotypes by the strategy of combined GWAS, gene-based association and haplotype analyses. The identified candidate genes provided valuable sources for future functional characterization and genetic improvement of rice appearance and milling quality.

Rice grain nutritional traits and their enhancement using relevant genes and QTLs through advanced approaches

BACKGROUND

Rice breeding program needs to focus on development of nutrient dense rice for value addition and helping in reducing malnutrition. Mineral and vitamin deficiency related problems are common in the majority of the population and more specific to developing countries as their staple food is rice.

RESULTS

Genes and QTLs are recently known for the nutritional quality of rice. By comprehensive literature survey and public domain database, we provided a critical review on nutritional aspects like grain protein and amino acid content, vitamins and minerals, glycemic index value, phenolic and flavonoid compounds, phytic acid, zinc and iron content along with QTLs linked to these traits. In addition, achievements through transgenic and advanced genomic approaches have been discussed. The information available on genes and/or QTLs involved in enhancement of micronutrient element and amino acids are summarized with graphical representation.

CONCLUSION

Compatible QTLs/genes may be combined together to design a desirable genotype with superior in multiple grain quality traits. The comprehensive review will be helpful to develop nutrient dense rice cultivars by integrating molecular markers and transgenic assisted breeding approaches with classical breeding.