Genome-wide association mapping for root traits in a panel of rice accessions from Vietnam

Dissecting the genetic architecture of root-related traits in a grafted wild Vitis berlandieri population for grapevine rootstock breeding

In woody perennial plants, quantitative genetics and association studies remain scarce for root-related traits, due to the time required to obtain mature plants and the complexity of phenotyping. In grapevine, a grafted cultivated plant, most of the rootstocks used are hybrids between American Vitis species (V. rupestris, V. riparia, and V. berlandieri). In this study, we used a wild population of an American Vitis species (V. berlandieri) to analyze the genetic architecture of the root-related traits of rootstocks in a grafted context. We studied a population consisting of 211 genotypes, with one to five replicates each (n = 846 individuals), plus four commercial rootstocks as control genotypes (110R, 5BB, Börner, and SO4). After two independent years of experimentation, the best linear unbiased estimates method revealed root-related traits with a moderate-to-high heritability (0.36-0.82) and coefficient of genetic variation (0.15-0.45). A genome-wide association study was performed with the BLINK model, leading to the detection of 11 QTL associated with four root-related traits (one QTL was associated with the total number of roots, four were associated with the number of small roots (< 1 mm in diameter), two were associated with the number of medium-sized roots (1 mm < diameter < 2 mm), and four were associated with mean diameter) accounting for up to 25.1% of the variance. Three genotypes were found to have better root-related trait performances than the commercial rootstocks and therefore constitute possible new candidates for use in grapevine rootstock breeding programs.

Effects of irrigation on root growth and development of soybean: A 3-year sandy field experiment

Increasing the water use efficiency of crops is an important agricultural goal closely related to the root system -the primary plant organ for water and nutrient acquisition. In an attempt to evaluate the response of root growth and development of soybean to water supply levels, 200 genotypes were grown in a sandy field for 3 years under irrigated and non-irrigated conditions, and 14 root traits together with shoot fresh weight and plant height were investigated. Three-way ANOVA revealed a significant effect of treatments and years on growth of plants, accounting for more than 80% of the total variability. The response of roots to irrigation was consistent over the years as most root traits were improved by irrigation. However, the actual values varied between years because the growth of plants was largely affected by the field microclimatic conditions (i.e., temperature, sunshine duration, and precipitation). Therefore, the best linear unbiased prediction values for each trait were calculated using the original data. Principal component analysis showed that most traits contributed to principal component (PC) 1, whereas average diameter, the ratio of thin and medium thickness root length to total root length contributed to PC2. Subsequently, we focused on selecting genotypes that exhibited significant improvements in root traits under irrigation than under non-irrigated conditions using the increment (I-index) and relative increment (RI-index) indices calculated for all traits. Finally, we screened for genotypes with high stability and root growth over the 3 years using the multi-trait selection index (MTSI).Six genotypes namely, GmJMC130, GmWMC178, GmJMC092, GmJMC068, GmWMC075, and GmJMC081 from the top 10% of genotypes scoring MTSI less than the selection threshold of 7.04 and 4.11 under irrigated and non-irrigated conditions, respectively, were selected. The selected genotypes have great potential for breeding cultivars with improved water usage abilities, meeting the goal of water-saving agriculture.

Response of Southeast Asian rice root architecture and anatomy phenotypes to drought stress

Drought stress in Southeast Asia greatly affects rice production, and the rice root system plays a substantial role in avoiding drought stress. In this study, we examined the phenotypic and genetic correlations among root anatomical, morphological, and agronomic phenotypes over multiple field seasons. A set of >200 rice accessions from Southeast Asia (a subset of the 3000 Rice Genomes Project) was characterized with the aim to identify root morphological and anatomical phenotypes related to productivity under drought stress. Drought stress resulted in slight increases in the basal metaxylem and stele diameter of nodal roots. Although few direct correlations between root phenotypes and grain yield were identified, biomass was consistently positively correlated with crown root number and negatively correlated with stele diameter. The accessions with highest grain yield were characterized by higher crown root numbers and median metaxylem diameter and smaller stele diameter. Genome-wide association study (GWAS) revealed 162 and 210 significant SNPs associated with root phenotypes in the two seasons which resulted in identification of 59 candidate genes related to root development. The gene OsRSL3 was found in a QTL region for median metaxylem diameter. Four SNPs in OsRSL3 were found that caused amino acid changes and significantly associated with the root phenotype. Based on the haplotype analysis for median metaxylem diameter, the rice accessions studied were classified into five allele combinations in order to identify the most favorable haplotypes. The candidate genes and favorable haplotypes provide information useful for the genetic improvement of root phenotypes under drought stress.

CROWN ROOTLESS1 binds DNA with a relaxed specificity and activates OsROP and OsbHLH044 genes involved in crown root formation in rice

In cereals, the root system is mainly composed of post-embryonic shoot-borne roots, named crown roots. The CROWN ROOTLESS1 (CRL1) transcription factor, belonging to the ASYMMETRIC LEAVES2-LIKE/LATERAL ORGAN BOUNDARIES DOMAIN (ASL/LBD) family, is a key regulator of crown root initiation in rice (Oryza sativa). Here, we show that CRL1 can bind, both in vitro and in vivo, not only the LBD-box, a DNA sequence recognized by several ASL/LBD transcription factors, but also another not previously identified DNA motif that was named CRL1-box. Using rice protoplast transient transactivation assays and a set of previously identified CRL1-regulated genes, we confirm that CRL1 transactivates these genes if they possess at least a CRL1-box or an LBD-box in their promoters. In planta, ChIP-qPCR experiments targeting two of these genes that include both a CRL1- and an LBD-box in their promoter show that CRL1 binds preferentially to the LBD-box in these promoter contexts. CRISPR/Cas9-targeted mutation of these two CRL1-regulated genes, which encode a plant Rho GTPase (OsROP) and a basic helix-loop-helix transcription factor (OsbHLH044), show that both promote crown root development. Finally, we show that OsbHLH044 represses a regulatory module, uncovering how CRL1 regulates specific processes during crown root formation.

Genome-Wide Association Analyses to Identify SNPs Related to Drought Tolerance

Drought stress is a serious agronomic problem resulting in significant yield losses globally. Breeding cultivars with drought tolerance is an important strategy that can be used to address this problem. Drought tolerance, however, is a complex multigenic trait, making advancements with conventional breeding approaches very challenging. This emphasizes the importance of dissecting the genetics of this trait and the identification and cloning of genes responsible for drought tolerance. With the rapid development of sequencing technologies and analytic methodologies, genome-wide association study (GWAS) has become an important tool for detecting natural variations underlying complex traits in crops. Identified loci can serve as targets for genomic selection or precise editing that enables the molecular design of new cultivars. This chapter describes the pipeline of statistical methods used in GWAS analysis, and covers field design, quality control, population structure control, association tests, and visualization of data. GWAS methodology used to dissect the genetic basis of drought tolerance is presented, and perspectives for optimizing the design and analysis of GWAS are discussed. The provided information serves as a valuable resource for researchers interested in GWAS technology.

Association mapping of drought tolerance and agronomic traits in rice (Oryza sativa L.) landraces

BACKGROUND

Asian cultivars were predominantly represented in global rice panel selected for sequencing and to identify novel alleles for drought tolerance. Diverse genetic resources adapted to Indian subcontinent were not represented much in spite harboring useful alleles that could improve agronomic traits, stress resilience and productivity. These rice accessions are valuable genetic resource in developing rice varieties suited to different rice ecosystem that experiences varying drought stress level, and at different crop stages. A core collection of rice germplasm adapted to Southwestern Indian peninsular genotyped using SSR markers and characterized by contrasting water regimes to associate genomic regions for physiological, root traits and yield related traits. Genotyping-By-Sequencing of selected accessions within the diverse panel revealed haplotype variation in genic content within genomic regions mapped for physiological, morphological and root traits.

RESULTS

Diverse rice panel (99 accessions) were evaluated in field and measurements on plant physiological, root traits and yield related traits were made over five different seasons experiencing varying drought stress intensity at different crop stages. Traits like chlorophyll stability index, leaf rolling, days to 50% flowering, chlorophyll content, root volume and root biomass were identified as best predictors of grain yield under stress. Association mapping revealed genetic variation among accessions and revealed 14 genomic targets associated with different physiological, root and plant production traits. Certain accessions were found to have beneficial allele to improve traits, plant height, root length and spikelet fertility, that contribute to the grain yield under stress. Genomic characterization of eleven accessions revealed haplotype variation within key genomic targets on chromosomes 1, 4, 6 and 11 for potential use as molecular markers to combine drought avoidance and tolerance traits. Genes mined within the genomic QTL intervals identified were prioritized based on tissue specific expression level in publicly available rice transcriptome data.

CONCLUSION

The genetic and genomic resources identified will enable combining traits with agronomic value to optimize yield under stress and hasten trait introgression into elite cultivars. Alleles associated with plant height, specific leaf area, root length from PTB8 and spikelet fertility and grain weight from PTB26 can be harnessed in future rice breeding program.

Nitrate Regulates Maize Root Transcriptome through Nitric Oxide Dependent and Independent Mechanisms

Maize root responds to nitrate by modulating its development through the coordinated action of many interacting players. Nitric oxide is produced in primary root early after the nitrate provision, thus inducing root elongation. In this study, RNA sequencing was applied to discover the main molecular signatures distinguishing the response of maize root to nitrate according to their dependency on, or independency of, nitric oxide, thus discriminating the signaling pathways regulated by nitrate through nitric oxide from those regulated by nitrate itself of by further downstream factors. A set of subsequent detailed functional annotation tools (Gene Ontology enrichment, MapMan, KEGG reconstruction pathway, transcription factors detection) were used to gain further information and the lateral root density was measured both in the presence of nitrate and in the presence of nitrate plus cPTIO, a specific NO scavenger, and compared to that observed for N-depleted roots. Our results led us to identify six clusters of transcripts according to their responsiveness to nitric oxide and to their regulation by nitrate provision. In general, shared and specific features for the six clusters were identified, allowing us to determine the overall root response to nitrate according to its dependency on nitric oxide.

Identification of Markers for Root Traits Related to Drought Tolerance Using Traditional Rice Germplasm

Drought is one of the important constraints affecting rice productivity worldwide. The vigorous shoot and deep root system help to improve drought resistance. In present era, genome-wide association study (GWAS) is the preferred method for mapping of QTLs for complex traits such as root and drought tolerance traits. In the present study, 114 rice genotypes were evaluated for various root and shoot traits under water stress conditions. All genotypes showed a significant amount of variation for various root and shoot traits. Correlation analysis revealed that high dry shoot weight and fresh shoot weight is associated with root length, root volume, fresh root weight and dry root weight. A total of 11 significant marker-trait associations were detected for various root, shoot and drought tolerance traits with the coefficient of determination (R²) ranging from 18.99 to 53.41%. Marker RM252 and RM212 showed association with three root traits which suggests their scope for improvement of root system. In the present study, a novel QTL was detected for root length associated with RM127, explaining 19.30% of variation. The marker alleles with increasing phenotypic effects for root and drought-tolerant traits can be exploited for improvement of root and drought tolerance traits using marker-assisted selection.

Genetic basis and network underlying synergistic roots and shoots biomass accumulation revealed by genome-wide association studies in rice

Genetic basis and network studies underlying synergistic biomass accumulation of roots and shoots (SBA) are conducive for rational design of high-biomass rice breeding. In this study, association signals for root weight, shoot weight, and the ratio of root-to-shoot mass (R/S) were identified using 666 rice accessions by genome-wide association study, together with their sub-traits, root length, root thickness and shoot length. Most association signals for root weight and shoot weight did not show association with their sub-traits. Based on the results, we proposed a top-to-bottom model for SBA, i.e. root weight, shoot weight and R/S were determined by their highest priority in contributing to biomass in the regulatory pathway, followed by a lower priority pathway for their sub-traits. Owing to 37 enriched clusters with more than two association signals identified, the relationship among the six traits could be also involved in linkage and pleiotropy. Furthermore, a discrimination of pleiotropy and LD at sequencing level using the known gene OsPTR9 for root weight, R/S and root length was provided. The results of given moderate correlation between traits and their corresponding sub-traits, and moderate additive effects between a trait and the accumulation of excellent alleles corresponding to its sub-traits supported a bottom-to-top regulation model for SBA. This model depicted each lowest-order trait (root length, root thickness and shoot length) was determined by its own regulation loci, and competition among different traits, as well as the pleiotropy and LD. All above ensure the coordinated development of each trait and the accumulation of the total biomass, although the predominant genetic basis of SBA is still indistinguishable. The presentation of the above two models and evidence of this study shed light on dissecting the genetic architecture of SBA.

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.

Assessment of root phenotypes in mungbean mini-core collection (MMC) from the World Vegetable Center (AVRDC) Taiwan

Mungbean (Vigna radiata L.) is an important food grain legume, but its production capacity is threatened by global warming, which can intensify plant stress and limit future production. Identifying new variation of key root traits in mungbean will provide the basis for breeding lines with effective root characteristics for improved water uptake to mitigate heat and drought stress. The AVRDC mungbean mini core collection consisting of 296 genotypes was screened under modified semi-hydroponic screening conditions to determine the variation for fourteen root-related traits. The AVRDC mungbean mini core collection displayed wide variations for the primary root length, total surface area, and total root length, and based on agglomerative hierarchical clustering eight homogeneous groups displaying different root traits could be identified. Germplasm with potentially favorable root traits has been identified for further studies to identify the donor genotypes for breeding cultivars with enhanced adaptation to water-deficit stress and other stress conditions.

Discovery of new genetic determinants of morphological plasticity in rice roots and shoots under phosphate starvation using GWAS

Phosphorus is an essential nutrient for plants that is often in short supply. In rice (Oryza sativa L.), inorganic phosphate (P_i) deficiency leads to various physiological disorders that consequently affect plant productivity. In this study, a large-scale phenotyping experiment using 160 Vietnamese rice landraces was performed under greenhouse conditions, by employing an alpha lattice design with three replicates, to identify quantitative trait loci (QTLs) associated with plant growth inhibition caused by P_i deficiency. Rice plantlets were grown for six weeks in the PVC sand column (16 cm diameter × 80 cm height) supplied with P_i-deficient medium (10 μM P) or full-P_i Yoshida medium (320 μM P). The effects of P_i deficiency on the number of crown roots, root length, shoot length, root weight, shoot weight and total weight were studied. From 36 significant markers identified using a genome-wide association study, 21 QTLs associated with plant growth inhibition under P_i starvation were defined. In total, 158 candidate genes co-located with the defined QTLs were identified. Interestingly, one QTL (qRST9.14) was associated with all three weight-traits. The co-located gene GLYCEROPHOSPHODIESTER PHOSPHODIESTERASE 13 was found to be potentially involved in P_i transport. Understanding the molecular mechanisms of Pi-starvation responses, and identifying the potential QTLs responsible for low-P_i stress tolerance, will provide valuable information for developing new varieties tolerant of low-P_i conditions.

Transcriptome profiling of laser-captured crown root primordia reveals new pathways activated during early stages of crown root formation in rice

Crown roots constitute the main part of the rice root system. Several key genes involved in crown root initiation and development have been identified by functional genomics approaches. Nevertheless, these approaches are impaired by functional redundancy and mutant lethality. To overcome these limitations, organ targeted transcriptome analysis can help to identify genes involved in crown root formation and early development. In this study, we generated an atlas of genes expressed in developing crown root primordia in comparison with adjacent stem cortical tissue at three different developmental stages before emergence, using laser capture microdissection. We identified 3975 genes differentially expressed in crown root primordia. About 30% of them were expressed at the three developmental stages, whereas 10.5%, 19.5% and 12.8% were specifically expressed at the early, intermediate and late stages, respectively. Sorting them by functional ontology highlighted an active transcriptional switch during the process of crown root primordia formation. Cross-analysis with other rice root development-related datasets revealed genes encoding transcription factors, chromatin remodeling factors, peptide growth factors, and cell wall remodeling enzymes that are likely to play a key role during crown root primordia formation. This atlas constitutes an open primary data resource for further studies on the regulation of crown root initiation and development.

Harnessing High-throughput Phenotyping and Genotyping for Enhanced Drought Tolerance in Crop Plants

Development of drought-tolerant cultivars is one of the challenging tasks for the plant breeders due to its complex inheritance and polygenic regulation. Evaluating genetic material for drought tolerance is a complex process due to its spatiotemporal interactions with environmental factors. The conventional breeding approaches are costly, lengthy, and inefficient to achieve the expected gain in drought tolerance. In this regard, genomics-assisted breeding (GAB) offers promise to develop cultivars with improved drought tolerance in a more efficient, quicker, and cost-effective manner. The success of GAB depends upon the precision in marker-trait association and estimation of genomic estimated breeding values (GEBVs), which mostly depends on coverage and precision of genotyping and phenotyping. A wide gap between the discovery and practical use of quantitative trait loci (QTL) for crop improvement has been observed for many important agronomical traits. Such a limitation could be due to the low accuracy in QTL detection, mainly resulting from low marker density and manually collected phenotypes of complex agronomic traits. Increasing marker density using the high-throughput genotyping (HTG), and accurate and precise phenotyping using high-throughput digital phenotyping (HTP) platforms can improve the precision and power of QTL detection. Therefore, both HTG and HTP can enhance the practical utility of GAB along with a faster characterization of germplasm and breeding material. In the present review, we discussed how the recent innovations in HTG and HTP would assist in the breeding of improved drought-tolerant varieties. We have also discussed strategies, tools, and analytical advances made on the HTG and HTP along with their pros and cons.

Low nitrogen induces root elongation via auxin-induced acid growth and auxin-regulated target of rapamycin (TOR) pathway in maize

Under low nitrogen (N) supply, an important adaption of the maize root system is to promote the root elongation so as to increase N uptake from a larger soil space. The underlying physiological mechanism is largely unknown. In the present study, two maize inbred lines (Ye478 and Wu312) were used to study the possible involvement of the auxin and target of rapamycin (TOR) pathway in low-N-induced root elongation. Compared to Wu312, primary root elongation of Ye478 was more sensitive to low nitrate supply. Correspondingly, more auxin was accumulated in the root tip, and more protons were secreted, increasing the acidity of the apoplast space. On the other hand, low-N-induced root elongation was greatly reduced when shoot-to-root auxin transport was inhibited by applying N-1-naphthylphthalamic acid (NPA) at the plant base or by pruning the top leaf where auxin is mostly synthesized. Furthermore, exogenous application of TOR inhibitor also eliminated the response of root elongation under low N. The content of TOR kinase and the expression of TOR pathway-related genes were significantly changed when shoot-to-root auxin transport was reduced by NPA treatment. Taken together, it is concluded that low-N stress increases shoot-to-root auxin transport which enhances root elongation via auxin-dependent acid growth and the auxin-regulated TOR pathway in maize.

Phenotyping Root Systems in a Set of Japonica Rice Accessions: Can Structural Traits Predict the Response to Drought?

BACKGROUND

The root system plays a major role in plant growth and development and root system architecture is reported to be the main trait related to plant adaptation to drought. However, phenotyping root systems in situ is not suited to high-throughput methods, leading to the development of non-destructive methods for evaluations in more or less controlled root environments. This study used a root phenotyping platform with a panel of 20 japonica rice accessions in order to: (i) assess their genetic diversity for a set of structural and morphological root traits and classify the different types; (ii) analyze the plastic response of their root system to a water deficit at reproductive phase and (iii) explore the ability of the platform for high-throughput phenotyping of root structure and morphology.

RESULTS

High variability for the studied root traits was found in the reduced set of accessions. Using eight selected traits under irrigated conditions, five root clusters were found that differed in root thickness, branching index and the pattern of fine and thick root distribution along the profile. When water deficit occurred at reproductive phase, some accessions significantly reduced root growth compared to the irrigated treatment, while others stimulated it. It was found that root cluster, as defined under irrigated conditions, could not predict the plastic response of roots under drought.

CONCLUSIONS

This study revealed the possibility of reconstructing the structure of root systems from scanned images. It was thus possible to significantly class root systems according to simple structural traits, opening up the way for using such a platform for medium to high-throughput phenotyping. The study also highlighted the uncoupling between root structures under non-limiting water conditions and their response to drought.

Exploiting natural variation in root system architecture via genome-wide association studies

Root growth and development has become an important research topic for breeders and researchers based on a growing need to adapt plants to changing and more demanding environmental conditions worldwide. Over the last few years, genome-wide association studies (GWASs) became an important tool to identify the link between traits in the field and their genetic background. Here we give an overview of the current literature concerning GWASs performed on root system architecture (RSA) in plants. We summarize which root traits and approaches have been used for GWAS, mentioning their respective success rate towards a successful gene discovery. Furthermore, we zoom in on the current technical hurdles in root phenotyping and GWAS, and discuss future possibilities in this field of research.

The physiological mechanism underlying root elongation in response to nitrogen deficiency in crop plants

In response to low nitrogen stress, multiple hormones together with nitric oxide signaling pathways work synergistically and antagonistically in crop root elongation. Changing root morphology allows plants to adapt to soil nutrient availability. Nitrogen is the most important essential nutrient for plant growth. An important adaptive strategy for crops responding to nitrogen deficiency is root elongation, thereby accessing increased soil space and nitrogen resources. Multiple signaling pathways are involved in this regulatory network, working together to fine-tune root elongation in response to soil nitrogen availability. Based on existing research, we propose a model to explain how different signaling pathways interact to regulate root elongation in response to low nitrogen stress. In response to a low shoot nitrogen status signal, auxin transport from the shoot to the root increases. High auxin levels in the root tip stimulate the production of nitric oxide, which promotes the synthesis of strigolactones to accelerate cell division. In this process, cytokinin, ethylene, and abscisic acid play an antagonistic role, while brassinosteroids and auxin play a synergistic role in regulating root elongation. Further study is required to identify the QTLs, genes, and favorable alleles which control the root elongation response to low nitrogen stress in crops.

Unraveling the Genetic Elements Involved in Shoot and Root Growth Regulation by Jasmonate in Rice Using a Genome-Wide Association Study

BACKGROUND

Due to their sessile life style, plant survival is dependent on the ability to build up fast and highly adapted responses to environmental stresses by modulating defense response and organ growth. The phytohormone jasmonate plays an essential role in regulating these plant responses to stress.

RESULTS

To assess variation of plant growth responses and identify genetic determinants associated to JA treatment, we conducted a genome-wide association study (GWAS) using an original panel of Vietnamese rice accessions. The phenotyping results showed a high natural genetic variability of the 155 tested rice accessions in response to JA for shoot and root growth. The level of growth inhibition by JA is different according to the rice varieties tested. We conducted genome-wide association study and identified 28 significant associations for root length (RTL), shoot length (SHL), root weight (RTW), shoot weight (SHW) and total weight (TTW) in response to JA treatment. Three common QTLs were found for RTL, RTW and SHL. Among a list of 560 candidate genes found to co-locate with the QTLs, a transcriptome analysis from public database for the JA response allows us to identify 232 regulated genes including several JA-responsive transcription factors known to play a role in stress response.

CONCLUSION

Our genome-wide association study shows that common and specific genetic elements are associated with inhibition of shoot and root growth under JA treatment suggesting the involvement of a complex JA-dependent genetic control of rice growth inhibition at the whole plant level. Besides, numerous candidate genes associated to stress and JA response are co-located with the association loci, providing useful information for future studies on genetics and breeding to optimize the growth-defense trade-off in rice.

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

Leaf traits are often strongly correlated with yield, which poses a major challenge in rice breeding. In the present study, using a panel of Vietnamese rice landraces genotyped with 21,623 single-nucleotide polymorphism markers, a genome-wide association study (GWAS) was conducted for several leaf traits during the vegetative stage. Vietnamese landraces are often poorly represented in panels used for GWAS, even though they are adapted to contrasting agrosystems and can contain original, valuable genetic determinants. A panel of 180 rice varieties was grown in pots for four weeks with three replicates under nethouse conditions. Different leaf traits were measured on the second fully expanded leaf of the main tiller, which often plays a major role in determining the photosynthetic capacity of the plant. The leaf fresh weight, turgid weight and dry weight were measured; then, from these measurements, the relative tissue weight and leaf dry matter percentage were computed. The leaf dry matter percentage can be considered a proxy for the photosynthetic efficiency per unit leaf area, which contributes to yield. By a GWAS, thirteen QTLs associated with these leaf traits were identified. Eleven QTLs were identified for fresh weight, eleven for turgid weight, one for dry weight, one for relative tissue weight and one for leaf dry matter percentage. Eleven QTLs presented associations with several traits, suggesting that these traits share common genetic determinants, while one QTL was specific to leaf dry matter percentage and one QTL was specific to relative tissue weight. Interestingly, some of these QTLs colocalize with leaf- or yield-related QTLs previously identified using other material. Several genes within these QTLs with a known function in leaf development or physiology are reviewed.

Next-Generation Sequencing Accelerates Crop Gene Discovery

The identification and isolation of genes underlying quantitative trait loci (QTLs) associated with agronomic traits in crops have been recently accelerated thanks to next-generation sequencing (NGS)-based technologies combined with plant genetics. With NGS, various revisited genetic approaches, which benefited from higher marker density, have been elaborated. These approaches improved resolution in QTL position and assisted in determining functional causative variations in genes. Examples of QTLs/genes associated with agronomic traits in crops and identified using different strategies based on whole-genome sequencing (WGS)/whole-genome resequencing (WGR) or RNA-seq are presented and discussed in this review. More specifically, we summarize and illustrate how NGS boosted bulk-segregant analysis (BSA), expression profiling, and the construction of polymorphism databases to facilitate the detection of QTLs and causative genes.

Genome-wide Association Study of a Panel of Vietnamese Rice Landraces Reveals New QTLs for Tolerance to Water Deficit During the Vegetative Phase

BACKGROUND

Drought tolerance is a major challenge in breeding rice for unfavorable environments. In this study, we used a panel of 180 Vietnamese rice landraces genotyped with 21,623 single-nucleotide polymorphism markers to perform a genome-wide association study (GWAS) for different drought response and recovery traits during the vegetative stage. These landraces originate from different geographical locations and are adapted to different agrosystems characterized by contrasted water regimes. Vietnamese landraces are often underrepresented in international panels used for GWAS, but they can contain original genetic determinants related to drought resistance.

RESULTS

The panel of 180 rice varieties was phenotyped under greenhouse conditions for several drought-related traits in an experimental design with 3 replicates. Plants were grown in pots for 4 weeks and drought-stressed by stopping irrigation for an additional 4 weeks. Drought sensitivity scores and leaf relative water content were measured throughout the drought stress. The recovery capacity was measured 2 weeks after plant rewatering. Several QTLs associated with these drought tolerance traits were identified by GWAS using a mixed model with control of structure and kinship. The number of detected QTLs consisted of 14 for leaf relative water content, 9 for slope of relative water content, 12 for drought sensitivity score, 3 for recovery ability and 1 for relative crop growth rate. This set of 39 QTLs actually corresponded to a total of 17 different QTLs because 9 were simultaneously associated with two or more traits, which indicates that these common loci may have pleiotropic effects on drought-related traits. No QTL was found in association with the same traits in both the indica and japonica subpanels. The possible candidate genes underlying the quantitative trait loci are reviewed.

CONCLUSIONS

Some of the identified QTLs contain promising candidate genes with a function related to drought tolerance by osmotic stress adjustment.

A genome-wide association study using a Vietnamese landrace panel of rice (Oryza sativa) reveals new QTLs controlling panicle morphological traits

CONTEXT

Yield improvement is an important issue for rice breeding. Panicle architecture is one of the key components of rice yield and exhibits a large diversity. To identify the morphological and genetic determinants of panicle architecture, we performed a detailed phenotypic analysis and a genome-wide association study (GWAS) using an original panel of Vietnamese landraces.

RESULTS

Using a newly developed image analysis tool, morphological traits of the panicles were scored over two years: rachis length; primary, secondary and tertiary branch number; average length of primary and secondary branches; average length of internode on rachis and primary branch. We observed a high contribution of spikelet number and secondary branch number per panicle to the overall phenotypic diversity in the dataset. Twenty-nine stable QTLs associated with seven traits were detected through GWAS over the two years. Some of these QTLs were associated with genes already implicated in panicle development. Importantly, the present study revealed the existence of new QTLs associated with the spikelet number, secondary branch number and primary branch number traits.

CONCLUSIONS

Our phenotypic analysis of panicle architecture variation suggests that with the panel of samples used, morphological diversity depends largely on the balance between indeterminate vs. determinate axillary meristem fate on primary branches, supporting the notion of differences in axillary meristem fate between rachis and primary branches. Our genome-wide association study led to the identification of numerous genomic sites covering all the traits studied and will be of interest for breeding programs aimed at improving yield. The new QTLs detected in this study provide a basis for the identification of new genes controlling panicle development and yield in rice.

Genetic analysis of roots and shoots in rice seedling by association mapping

The vigorous shoots and roots help to improve drought resistance and post-transplanting recovery in rice seedlings (Oryza sativa L.). Hundreds of loci related to root system have been identified recently, but little research has been done on shoot traits, and the relationship between roots and shoots development is also still unclear. The objective of this study was to identify associated loci for roots and shoots in rice seedlings as well as to screen pleiotropic QTLs involved in coordinated development of roots and shoots. Using mini core collection of 273 cultivated rice accessions and 280 simple-sequence repeat markers, we investigated six traits [root length (RL), root thickness (RT), root weight (RW), shoot length (SL), shoot weight (SW) and ratio of root-to-shoot mass] in seedlings. Study was performed in hydroponic medium and genetic analysis was performed by association mapping using general linear model (GLM) with population structure (Q) and mixed linear model (MLM) involving Q and familial relatedness (K). Two subgroups indica and japonica showed significant differences in RT, RW and SW. Maximum correlation was observed between RW and SW. Using GLM 65 QTLs for root and 43 QTLs associated with shoot traits were detected. Among them, seven QTLs were present between RL and RW and five common QTLs were detected between SL and SW with high phenotypic variation effects (PVEs). Two key pleiotropic QTLs were also identified involved in collaborative development of roots and shoots in rice seedlings. Importantly, 17 and 10 QTLs were identified for root and shoot traits respectively in both studies of GLM and MLM. More common QTLs with high PVEs between root and shoot traits suggested that longitudinal growth (RL and SL) played an important role in accumulation of biomass (RW and SW). Considering the obvious phenotypic differences and fewer common QTLs between indica and japonica, we suggested that there could be different mechanisms of seedling development between both subpopulations. Key pleiotropic QTLs and QTLs identified for root and shoot traits in both studies of GLM and MLM could be preferentially used in marker-assisted breeding for strong rice seedling.

Identification of candidate genes for gelatinization temperature, gel consistency and pericarp color by GWAS in rice based on SLAF-sequencing

Rice is an important cereal in the world. The study of the genetic basis of important agronomic traits in rice landraces and identification of genes will facilitate the breed improvement. Gelatinization temperature (GT), gel consistency (GC) and pericarp color (PC) are important indices of rice cooking and eating quality evaluation and potential nutritional importance, which attract wide attentions in the application of genetic and breeding. To dissect the genetic basis of GT, GC and PC, a total of 419 rice landraces core germplasm collections consisting of 330 indica lines, 78 japonica lines and 11 uncertain varieties were planted, collected, then GT, GC, PC were measured for two years, and sequenced using specific-locus amplified fragment sequencing (SLAF-seq) technology. In this study, 261,385,070 clean reads and 56,768 polymorphic SLAF tags were obtained, which a total of 211,818 single nucleotide polymorphisms (SNPs) were discovered. With 208,993 SNPs meeting the criterion of minor allele frequency (MAF) > 0.05 and integrity> 0.5, the phylogenetic tree and population structure analysis were performed for all 419 rice landraces, and the whole panel mainly separated into six subpopulations based on population structure analysis. Genome-wide association study (GWAS) was carried out for the whole panel, indica subpanel and japonica subpanel with subset SNPs respectively. One quantitative trait locus (QTL) on chromosome 6 for GT was detected in the whole panel and indica subpanel, and one QTL associated with GC was located on chromosome 6 in the whole panel and indica subpanel. For the PC trait, 8 QTLs were detected in the whole panel on chromosome 1, 3, 4, 7, 8, 10 and 11, and 7 QTLs in the indica subpanel on chromosome 3, 4, 7, 8, 10 and 11. For the three traits, no QTL was detected in japonica subpanel, probably because of the polymorphism repartition between the subpanel, or small population size of japonica subpanel. This paper provides new gene resources and insights into the molecular mechanisms of important agricultural trait of rice phenotypic variation and genetic improvement of rice quality variety breeding.

Tolerance to mild salinity stress in japonica rice: A genome-wide association mapping study highlights calcium signaling and metabolism genes

Salinity tolerance is an important quality for European rice grown in river deltas. We evaluated the salinity tolerance of a panel of 235 temperate japonica rice accessions genotyped with 30,000 SNP markers. The panel was exposed to mild salt stress (50 mM NaCl; conductivity of 6 dS m^-1) at the seedling stage. Eight different root and shoot growth parameters were measured for both the control and stressed treatments. The Na⁺ and K⁺ mass fractions of the stressed plants were measured using atomic absorption spectroscopy. The salt treatment affected plant growth, particularly the shoot parameters. The panel showed a wide range of Na⁺/K⁺ ratio and the temperate accessions were distributed over an increasing axis, from the most resistant to the most susceptible checks. We conducted a genome-wide association study on indices of stress response and ion mass fractions in the leaves using a classical mixed model controlling structure and kinship. A total of 27 QTLs validated by sub-sampling were identified. For indices of stress responses, we also used another model that focused on marker × treatment interactions and detected 50 QTLs, three of which were also identified using the classical method. We compared the positions of the significant QTLs to those of approximately 300 genes that play a role in rice salt tolerance. The positions of several QTLs were close to those of genes involved in calcium signaling and metabolism, while other QTLs were close to those of kinases. These results reveal the salinity tolerance of accessions with a temperate japonica background. Although the detected QTLs must be confirmed by other approaches, the number of associations linked to candidate genes involved in calcium-mediated ion homeostasis highlights pathways to explore in priority to understand the salinity tolerance of temperate rice.

Favorable Marker Alleles for Panicle Exsertion Length in Rice (Oryza sativa L.) Mined by Association Mapping and the RSTEP-LRT Method

The panicle exsertion length (PEL) in rice (Oryza sativa L.) is an important trait for hybrid seed production. We investigated the PEL in a chromosome segment substitution line (CSSL) population consisting of 66 lines and a natural population composed of 540 varieties. In the CSSL population, a total of seven QTLs for PEL were detected across two environments. The percentage of phenotypic variance explained (PVE) ranged from 10.22 to 50.18%, and the additive effect ranged from -1.77 to 6.47 cm. Among the seven QTLs, qPEL10.2 had the largest PVE, 44.05 and 50.18%, with an additive effect of 5.91 and 6.47 cm in 2015 and in 2016, respectively. In the natural population, 13 SSR marker loci were detected that were associated with PEL in all four environments, with the PVE ranging from 1.20 to 6.26%. Among the 13 loci, 7 were novel. The RM5746-170 bp allele had the largest phenotypic effect (5.11 cm), and the typical carrier variety was Qiaobinghuang. An RM5620-RM6100 region harboring the EUI2 locus on chromosome 10 was detected in both populations. The sequencing results showed that the accessions with a shorter PEL contained the A base, while the accessions with a longer PEL contained the G base at the 1,475 bp location of the EUI2 gene.

Genome-wide association mapping for root cone angle in rice

BACKGROUND

Plant root systems play a major role in anchoring and in water and nutrient uptake from the soil. The root cone angle is an important parameter of the root system architecture because, combined with root depth, it helps to determine the volume of soil explored by the plant. Two genes, DRO1 and SOR1, and several QTLs for root cone angle have been discovered in the last 5 years.

RESULTS

To find other QTLs linked to root cone angle, a genome-wide association mapping study was conducted on two panels of 162 indica and 169 japonica rice accessions genotyped with two sets of SNP markers (genotyping-by-sequencing set with approximately 16,000 markers and high-density-rice-array set with approximately 300,000 markers). The root cone angle of all accessions was measured using a screen protractor on images taken after 1 month of plant growth in the Rhizoscope phenotyping system. The distribution of the root cone angle in the indica panel was Gaussian, but several accessions of the japonica panel (all the bulus from Indonesia and three temperate japonicas from Nepal or India) appeared as outliers with a very wide root cone angle. The data were submitted to association mapping using a mixed model with control of structure and kinship. A total of 15 QTLs for the indica panel and 40 QTLs for the japonica panel were detected. Genes underlying these QTLs (+/-50 kb from the significant markers) were analyzed. We focused our analysis on auxin-related genes, kinases, and genes involved in root developmental processes and identified 8 particularly interesting genes.

CONCLUSIONS

The present study identifies new sources of wide root cone angle in rice, proposes ways to bypass some drawbacks of association mapping to further understand the genetics of the trait and identifies candidate genes deserving further investigation.

Genetic Control of Plasticity in Root Morphology and Anatomy of Rice in Response to Water Deficit

Elucidating the genetic control of rooting behavior under water-deficit stress is essential to breed climate-robust rice (Oryza sativa) cultivars. Using a diverse panel of 274 indica genotypes grown under control and water-deficit conditions during vegetative growth, we phenotyped 35 traits, mostly related to root morphology and anatomy, involving 45,000 root-scanning images and nearly 25,000 cross sections from the root-shoot junction. The phenotypic plasticity of these traits was quantified as the relative change in trait value under water-deficit compared with control conditions. We then carried out a genome-wide association analysis on these traits and their plasticity, using 45,608 high-quality single-nucleotide polymorphisms. One hundred four significant loci were detected for these traits under control conditions, 106 were detected under water-deficit stress, and 76 were detected for trait plasticity. We predicted 296 (control), 284 (water-deficit stress), and 233 (plasticity) a priori candidate genes within linkage disequilibrium blocks for these loci. We identified key a priori candidate genes regulating root growth and development and relevant alleles that, upon validation, can help improve rice adaptation to water-deficit stress.

Genetic variation and association mapping of waterlogging tolerance in chrysanthemum

Forty-five molecular markers were detected significantly associated with chrysanthemum' waterlogging tolerance, and four favorable parental lines were identified as potential donors for improving waterlogging tolerance in chrysanthemum. The productivity of chrysanthemum is downgraded by waterlogging soils, which has driven a search for germplasm showing an enhanced level of waterlogging tolerance (WT). As yet little is known regarding the mode of inheritance of WT in chrysanthemum. The study set out to characterize the extent of genetic variation for WT represented in a collection of one hundred chrysanthemum accessions by testing them under both greenhouse and field conditions. A membership function value of waterlogging (MFVW), which integrated a wilting index, a chlorosis score and the proportion of dead leaf in waterlogged plants, was used as a measure of WT. The variation for MFVW among plants grown in the greenhouse (two experiments) was generally higher than that generated in field-grown (one experiment) plants. The MFVW broad sense heritability was 0.82, and the phenotypic coefficient of variation (31.8 %) was larger than the genetic one (28.8 %). Association mapping (AM) identified 45 markers related to WT: 25 by applying the general linear model (GLM) + principal component (PC) model, 16 by applying the mixed linear model (MLM), 31 by applying the MLM + Q matrix model and 12 by applying the MLM + PC model. Of the associated markers, eight and two were predictive in two and three experiments within all models, respectively; the proportion of the phenotypic variance explained by the eight associations ranged from 6.3 to 16.4 %. On the basis of their harboring all four of the leading markers E2M16-2, SSR150-6, E19M16-1 and E10M10-12, the varieties 'Nannong Xuefeng', 'Qx097', 'Nannong Xunzhang' and 'Finch' were identified as potential donors for future improvement of WT in chrysanthemum.