Genome-wide association mapping of leaf mass traits in a Vietnamese rice landrace panel

Multi-model genome-wide association studies of leaf anatomical traits and vein architecture in rice

Introduction

The anatomy of rice leaves is closely related to photosynthesis and grain yield. Therefore, exploring insight into the quantitative trait loci (QTLs) and alleles related to rice flag leaf anatomical and vein traits is vital for rice improvement.

Methods

Here, we aimed to explore the genetic architecture of eight flag leaf traits using one single-locus model; mixed-linear model (MLM), and two multi-locus models; fixed and random model circulating probability unification (FarmCPU) and Bayesian information and linkage disequilibrium iteratively nested keyway (BLINK). We performed multi-model GWAS using 329 rice accessions of RDP1 with 700K single-nucleotide polymorphisms (SNPs) markers.

Results

The phenotypic correlation results indicated that rice flag leaf thickness was strongly correlated with leaf mesophyll cells layer (ML) and thickness of both major and minor veins. All three models were able to identify several significant loci associated with the traits. MLM identified three non-synonymous SNPs near NARROW LEAF 1 (NAL1) in association with ML and the distance between minor veins (IVD) traits.

Discussion

Several numbers of significant SNPs associated with known gene function in leaf development and yield traits were detected by multi-model GWAS performed in this study. Our findings indicate that flag leaf traits could be improved via molecular breeding and can be one of the targets in high-yield rice development.

Genomic regions and candidate genes selected during the breeding of rice in Vietnam

Vietnam harnesses a rich diversity of rice landraces adapted to a range of conditions, which constitute a largely untapped source of diversity for the continuous improvement of cultivars. We previously identified a strong population structure in Vietnamese rice, which is captured in five Indica and four Japonica subpopulations, including an outlying Indica‐5 group. Here, we leveraged that strong differentiation and 672 native rice genomes to identify genomic regions and genes putatively selected during the breeding of rice in Vietnam. We identified significant distorted patterns in allele frequency (XP‐CLR) and population differentiation scores (F_ST) resulting from differential selective pressures between native subpopulations, and later annotated them with QTLs previously identified by GWAS in the same panel. We particularly focussed on the outlying Indica‐5 subpopulation because of its likely novelty and differential evolution, where we annotated 52 selected regions, which represented 8.1% of the rice genome. We annotated the 4576 genes in these regions and selected 65 candidate genes as promising breeding targets, several of which harboured alleles with nonsynonymous substitutions. Our results highlight genomic differences between traditional Vietnamese landraces, which are likely the product of adaption to multiple environmental conditions and regional culinary preferences in a very diverse country. We also verified the applicability of this genome scanning approach to identify potential regions harbouring novel loci and alleles to breed a new generation of sustainable and resilient rice.

Genome-wide association study for lignocellulosic compounds and fermentable sugar in rice straw

Cellulose and lignin are the two main components of secondary plant cell walls with substantial impact on stalk in the field and on straw during industrial processing. The amount of fermentable sugar that can be accessed is another important parameter affecting various industrial applications. In the present study, genetic variability of rice (Oryza sativa L.) genotypes for cellulose, lignin, and fermentable sugars contents was analyzed in rice straw. A genome-wide association study of 33,484 single nucleotide polymorphisms (SNPs) with a minor allele frequency (MAF) >0.05 was performed. The genome-wide association study identified seven, three, and three genomic regions to be significantly associated with cellulose, lignin, and fermentable sugar contents, respectively. Candidate genes in the associated genomic regions were enzymes mainly involved in cell wall metabolism. Novel SNP markers associated with cellulose were tagged to GH16, peroxidase, GT6, GT8, and CSLD2. For lignin content, Villin protein, OsWAK1/50/52/53, and GH16 were identified. For fermentable sugar content, UTP-glucose-1-phosphate uridylyltransferase, BRASSINOSTEROID INSENSITIVE 1, and receptor-like protein kinase 5 were found. The results of this study should improve our understanding of the genetic basis of the factors that might be involved in biosynthesis, turnover, and modification of major cell wall components and saccharides in rice straw.

QTL mapping and candidate gene mining of flag leaf size traits in Japonica rice based on linkage mapping and genome-wide association study

BACKGROUND

As one of the most important factors of the japonica rice plant, leaf shape affects the photosynthesis and carbohydrate accumulation directly. Mining and using new leaf shape related genes/QTLs can further enrich the theory of molecular breeding and accelerate the breeding process of japonica rice.

METHODS

In the present study, 2 RILs and a natural population with 295 japonica rice varieties were used to map QTLs for flag leaf length (FL), flag leaf width (FW) and flag leaf area (FLA) by linkage analysis and genome-wide association study (GWAS) throughout 2 years.

RESULTS

A total of 64 QTLs were detected by 2 ways, and pleiotropic QTLs qFL2 (Chr2_33,332,579) and qFL10 (Chr10_10,107,835; Chr10_10,230,100) consisted of overlapping QTLs mapped by linkage analysis and GWAS throughout the 2 years were identified.

CONCLUSIONS

The candidate genes LOC_Os02g54254, LOC_Os02g54550, LOC_Os10g20160, LOC_Os10g20240, LOC_Os10g20260 were obtained, filtered by linkage disequilibrium (LD), and haplotype analysis. LOC_Os10g20160 (SD-RLK-45) showed outstanding characteristics in quantitative real-time PCR (qRT-PCR) analysis in leaf development period, belongs to S-domain receptor-like protein kinases gene and probably to be a main gene regulating flag leaf width of japonica rice. The results of this study provide valuable resources for mining the main genes/QTLs of japonica rice leaf development and molecular breeding of japonica rice ideal leaf shape.

Importance of Landraces in Cereal Breeding for Stress Tolerance

The renewed focus on cereal landraces is a response to some negative consequences of modern agriculture and conventional breeding which led to a reduction of genetic diversity. Cereal landraces are still cultivated on marginal lands due to their adaptability to unfavourable conditions, constituting an important source of genetic diversity usable in modern plant breeding to improve the adaptation to abiotic or biotic stresses, yield performance and quality traits in limiting environments. Traditional agricultural production systems have played an important role in the evolution and conservation of wide variability in gene pools within species. Today, on-farm and ex situ conservation in gene bank collections, together with data sharing among researchers and breeders, will greatly benefit cereal improvement. Many efforts are usually made to collect, organize and phenotypically and genotypically analyse cereal landrace collections, which also utilize genomic approaches. Their use in breeding programs based on genomic selection, and the discovery of beneficial untapped QTL/genes/alleles which could be introgressed into modern varieties by MAS, pyramiding or biotechnological tools, increase the potential for their better deployment and exploitation in breeding for a more sustainable agricultural production, particularly enhancing adaptation and productivity in stress-prone environments to cope with current climate changes.

Resequencing of 672 Native Rice Accessions to Explore Genetic Diversity and Trait Associations in Vietnam

BACKGROUND

Vietnam possesses a vast diversity of rice landraces due to its geographical situation, latitudinal range, and a variety of ecosystems. This genetic diversity constitutes a highly valuable resource at a time when the highest rice production areas in the low-lying Mekong and Red River Deltas are enduring increasing threats from climate changes, particularly in rainfall and temperature patterns.

RESULTS

We analysed 672 Vietnamese rice genomes, 616 newly sequenced, that encompass the range of rice varieties grown in the diverse ecosystems found throughout Vietnam. We described four Japonica and five Indica subpopulations within Vietnam likely adapted to the region of origin. We compared the population structure and genetic diversity of these Vietnamese rice genomes to the 3000 genomes of Asian cultivated rice. The named Indica-5 (I5) subpopulation was expanded in Vietnam and contained lowland Indica accessions, which had very low shared ancestry with accessions from any other subpopulation and were previously overlooked as admixtures. We scored phenotypic measurements for nineteen traits and identified 453 unique genotype-phenotype significant associations comprising twenty-one QTLs (quantitative trait loci). The strongest associations were observed for grain size traits, while weaker associations were observed for a range of characteristics, including panicle length, heading date and leaf width.

CONCLUSIONS

We showed how the rice diversity within Vietnam relates to the wider Asian rice diversity by using a number of approaches to provide a clear picture of the novel diversity present within Vietnam, mainly around the Indica-5 subpopulation. Our results highlight differences in genome composition and trait associations among traditional Vietnamese rice accessions, which are likely the product of adaption to multiple environmental conditions and regional preferences in a very diverse country. Our results highlighted traits and their associated genomic regions that are a potential source of novel loci and alleles to breed a new generation of low input sustainable and climate resilient rice.

Genome-Wide Association Mapping of Salinity Tolerance at the Seedling Stage in a Panel of Vietnamese Landraces Reveals New Valuable QTLs for Salinity Stress Tolerance Breeding in Rice

Rice tolerance to salinity stress involves diverse and complementary mechanisms, such as the regulation of genome expression, activation of specific ion-transport systems to manage excess sodium at the cell or plant level, and anatomical changes that avoid sodium penetration into the inner tissues of the plant. These complementary mechanisms can act synergistically to improve salinity tolerance in the plant, which is then interesting in breeding programs to pyramidize complementary QTLs (quantitative trait loci), to improve salinity stress tolerance of the plant at different developmental stages and in different environments. This approach presupposes the identification of salinity tolerance QTLs associated with different mechanisms involved in salinity tolerance, which requires the greatest possible genetic diversity to be explored. To contribute to this goal, we screened an original panel of 179 Vietnamese rice landraces genotyped with 21,623 SNP markers for salinity stress tolerance under 100 mM NaCl treatment, at the seedling stage, with the aim of identifying new QTLs involved in the salinity stress tolerance via a genome-wide association study (GWAS). Nine salinity tolerance-related traits, including the salt injury score, chlorophyll and water content, and K⁺ and Na⁺ contents were measured in leaves. GWAS analysis allowed the identification of 26 QTLs. Interestingly, ten of them were associated with several different traits, which indicates that these QTLs act pleiotropically to control the different levels of plant responses to salinity stress. Twenty-one identified QTLs colocalized with known QTLs. Several genes within these QTLs have functions related to salinity stress tolerance and are mainly involved in gene regulation, signal transduction or hormone signaling. Our study provides promising QTLs for breeding programs to enhance salinity tolerance and identifies candidate genes that should be further functionally studied to better understand salinity tolerance mechanisms in rice.

Deep placement of nitrogen fertilizer improves yield, nitrogen use efficiency and economic returns of transplanted fine rice

Rice (Oryza sativa L.) feeds to two-third of the global population by serving as staple food. It is the main export commodity of several countries; thus, contributes towards foreign exchange earnings. Unfortunately, average global rice yield is far below than its genetic potential. Low nitrogen (N) use efficiency (NUE) is among the major reasons for low average yield. Current study evaluated the impact of nitrogen fertilizer application methods (conventional and deep placement) on growth, yield-related traits, chlorophyll contents, photosynthesis rate, agronomic N-use efficiency (ANUE), partial factors productivity of applied N (PFP) and economic returns of two different transplanted rice varieties (Basmati-515 and Super-Basmati). Fertilizer application methods significantly affected allometry, yield-related traits, chlorophyll contents, photosynthesis rate, ANUE, PFP and economic returns. Deep placement of N-fertilizer (DPNF) observed better allometric traits, high chlorophyll contents, photosynthesis rate, ANUE, PFP, yield attributes and economic returns compared to conventional application of N-fertilizer (CANF). Similarly, Basmati-515 had better allometric and yield-related traits, chlorophyll contents, photosynthesis rate, ANUE, PFP and economic returns than Super-Basmati. Regarding interactions among N-fertilizer application methods and rice varieties, Basmati-515 with DPNF resulted in higher chlorophyll contents, photosynthesis rate, ANUE, PFP, allometric and yield related traits and economic returns than CANF. The lowest values of these traits were observed for Super-Basmati with no application of N-fertilizer. Both varieties had better yield and economic returns with DPNF compared to CANF. It is concluded that DPNF improved yield, ANUE and economic returns; therefore, should be opted to improve productivity of transplanted fine rice. Nonetheless, lower nitrogen doses need to be tested for DPNF to infer whether it could lower N use in rice crop.

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.

Genome-wide association and epistatic interactions of flowering time in soybean cultivar

Genome-wide association studies (GWAS) have enabled the discovery of candidate markers that play significant roles in various complex traits in plants. Recently, with increased interest in the search for candidate markers, studies on epistatic interactions between single nucleotide polymorphism (SNP) markers have also increased, thus enabling the identification of more candidate markers along with GWAS on single-variant-additive-effect. Here, we focused on the identification of candidate markers associated with flowering time in soybean (Glycine max). A large population of 2,662 cultivated soybean accessions was genotyped using the 180k Axiom® SoyaSNP array, and the genomic architecture of these accessions was investigated to confirm the population structure. Then, GWAS was conducted to evaluate the association between SNP markers and flowering time. A total of 93 significant SNP markers were detected within 59 significant genes, including E1 and E3, which are the main determinants of flowering time. Based on the GWAS results, multilocus epistatic interactions were examined between the significant and non-significant SNP markers. Two significant and 16 non-significant SNP markers were discovered as candidate markers affecting flowering time via interactions with each other. These 18 candidate SNP markers mapped to 18 candidate genes including E1 and E3, and the 18 candidate genes were involved in six major flowering pathways. Although further biological validation is needed, our results provide additional information on the existing flowering time markers and present another option to marker-assisted breeding programs for regulating flowering time of soybean.