Genome-wide association mapping for yield and other agronomic traits in an elite breeding population of tropical rice (Oryza sativa)

Evaluation of deep learning for predicting rice traits using structural and single-nucleotide genomic variants

BACKGROUND

Structural genomic variants (SVs) are prevalent in plant genomes and have played an important role in evolution and domestication, as they constitute a significant source of genomic and phenotypic variability. Nevertheless, most methods in quantitative genetics focusing on crop improvement, such as genomic prediction, consider only Single Nucleotide Polymorphisms (SNPs). Deep Learning (DL) is a promising strategy for genomic prediction, but its performance using SVs and SNPs as genetic markers remains unknown.

RESULTS

We used rice to investigate whether combining SVs and SNPs can result in better trait prediction over SNPs alone and examine the potential advantage of Deep Learning (DL) networks over Bayesian Linear models. Specifically, the performances of BayesC (considering additive effects) and a Bayesian Reproducible Kernel Hilbert space (RKHS) regression (considering both additive and non-additive effects) were compared to those of two different DL architectures, the Multilayer Perceptron, and the Convolution Neural Network, to explore their prediction ability by using various marker input strategies. We found that exploiting structural and nucleotide variation slightly improved prediction ability on complex traits in 87% of the cases. DL models outperformed Bayesian models in 75% of the studied cases, considering the four traits and the two validation strategies used. Finally, DL systematically improved prediction ability of binary traits against the Bayesian models.

CONCLUSIONS

Our study reveals that the use of structural genomic variants can improve trait prediction in rice, independently of the methodology used. Also, our results suggest that Deep Learning (DL) networks can perform better than Bayesian models in the prediction of binary traits, and in quantitative traits when the training and target sets are not closely related. This highlights the potential of DL to enhance crop improvement in specific scenarios and the importance to consider SVs in addition to SNPs in genomic selection.

TrG2P: A transfer-learning-based tool integrating multi-trait data for accurate prediction of crop yield

Yield prediction is the primary goal of genomic selection (GS)-assisted crop breeding. Because yield is a complex quantitative trait, making predictions from genotypic data is challenging. Transfer learning can produce an effective model for a target task by leveraging knowledge from a different, but related, source domain and is considered a great potential method for improving yield prediction by integrating multi-trait data. However, it has not previously been applied to genotype-to-phenotype prediction owing to the lack of an efficient implementation framework. We therefore developed TrG2P, a transfer-learning-based framework. TrG2P first employs convolutional neural networks (CNN) to train models using non-yield-trait phenotypic and genotypic data, thus obtaining pre-trained models. Subsequently, the convolutional layer parameters from these pre-trained models are transferred to the yield prediction task, and the fully connected layers are retrained, thus obtaining fine-tuned models. Finally, the convolutional layer and the first fully connected layer of the fine-tuned models are fused, and the last fully connected layer is trained to enhance prediction performance. We applied TrG2P to five sets of genotypic and phenotypic data from maize (Zea mays), rice (Oryza sativa), and wheat (Triticum aestivum) and compared its model precision to that of seven other popular GS tools: ridge regression best linear unbiased prediction (rrBLUP), random forest, support vector regression, light gradient boosting machine (LightGBM), CNN, DeepGS, and deep neural network for genomic prediction (DNNGP). TrG2P improved the accuracy of yield prediction by 39.9%, 6.8%, and 1.8% in rice, maize, and wheat, respectively, compared with predictions generated by the best-performing comparison model. Our work therefore demonstrates that transfer learning is an effective strategy for improving yield prediction by integrating information from non-yield-trait data. We attribute its enhanced prediction accuracy to the valuable information available from traits associated with yield and to training dataset augmentation. The Python implementation of TrG2P is available at https://github.com/lijinlong1991/TrG2P. The web-based tool is available at http://trg2p.ebreed.cn:81.

GWAS analysis reveals the genetic basis of blast resistance associated with heading date in rice

Rice blast is a destructive fungal disease affecting rice plants at various growth stages, significantly threatening global yield stability. Development of resistant rice cultivars stands as a practical means of disease control. Generally, association mapping with a diversity panel powerfully identifies new alleles controlling trait of interest. On the other hand, utilization of a breeding panel has its advantage that can be directly applied in a breeding program. In this study, we conducted a genome-wide association study (GWAS) for blast resistance using 296 commercial rice cultivars with low population structure but large phenotypic diversity. We attempt to answer the genetic basis behind rice blast resistance among early maturing cultivars by subdividing the population based on its Heading date 1 (Hd1) functionality. Subpopulation-specific GWAS using the mixed linear model (MLM) based on blast nursery screening conducted in three years revealed a total of 26 significant signals, including three nucleotide-binding site leucine-rich repeat (NBS-LRR) genes (Os06g0286500, Os06g0286700, and Os06g0287500) located at Piz locus on chromosome 6, and one at the Pi-ta locus (Os12g0281300) on chromosome 12. Haplotype analysis revealed blast resistance associated with Piz locus was exclusively specific to Type 14 hd1 among japonica rice. Our findings provide valuable insights for breeding blast resistant rice and highlight the applicability of our elite cultivar panel to detect superior alleles associated with important agronomic traits.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

Genome-wide association studies for phenological and agronomic traits in mungbean (Vigna radiata L. Wilczek)

Mungbean (Vigna radiata L. Wilczek) is one of the important warm-season food legumes, contributing substantially to nutritional security and environmental sustainability. The genetic complexity of yield-associated agronomic traits in mungbean is not well understood. To dissect the genetic basis of phenological and agronomic traits, we evaluated 153 diverse mungbean genotypes for two phenological (days to heading and days to maturity) and eight agronomic traits (leaf nitrogen status using SPAD, plant height, number of primary branches, pod length, number of pods per plant, seeds per pod, 100-seed weight, and yield per plant) under two environmental conditions. A wide array of phenotypic variability was apparent among the studied genotypes for all the studied traits. The broad sense of heritability of traits ranged from 0.31 to 0.95 and 0.21 to 0.94 at the Delhi and Ludhiana locations, respectively. A total of 55,634 genome-wide single nucleotide polymorphisms (SNPs) were obtained by the genotyping-by-sequencing method, of which 15,926 SNPs were retained for genome-wide association studies (GWAS). GWAS with Bayesian information and linkage-disequilibrium iteratively nested keyway (BLINK) model identified 50 SNPs significantly associated with phenological and agronomic traits. In total, 12 SNPs were found to be significantly associated with phenological traits across environments, explaining 7%-18.5% of phenotypic variability, and 38 SNPs were significantly associated with agronomic traits, explaining 4.7%-27.6% of the phenotypic variability. The maximum number of SNPs (15) were located on chromosome 1, followed by seven SNPs each on chromosomes 2 and 8. The BLAST search identified 19 putative candidate genes that were involved in light signaling, nitrogen responses, phosphorus (P) transport and remobilization, photosynthesis, respiration, metabolic pathways, and regulating growth and development. Digital expression analysis of 19 genes revealed significantly higher expression of 12 genes, viz. VRADI01G08170, VRADI11G09170, VRADI02G00450, VRADI01G00700, VRADI07G14240, VRADI03G06030, VRADI02G14230, VRADI08G01540, VRADI09G02590, VRADI08G00110, VRADI02G14240, and VRADI02G00430 in the roots, cotyledons, seeds, leaves, shoot apical meristems, and flowers. The identified SNPs and putative candidate genes provide valuable genetic information for fostering genomic studies and marker-assisted breeding programs that improve yield and agronomic traits in mungbean.

Environmental Sustainability in Selected OPEC Countries: Do the Influence of FDI and ICT Matter?

Considering the undeniable role and importance of the environment in people's lives, the present study is designed to investigate the combined effect of information and communication technology (ICT) and foreign direct investment (FDI) on achieving environmental sustainability. Since the increasing emission of carbon in society and its destructive environmental effects on social economic aspects and even political tensions have become a challenge, the main question of the research is what strategies have governments, especially oil exporting countries, used in the past to reduce the level they have discovered pollution and what policies do they want to follow in the future? Among the policies undertaken by the OPEC oil exporting countries, has the action for foreign direct investment (FDI) and the development of information and communication technology (ICT) been effective in preventing harmful environmental effects? For this purpose, data on renewable energy consumption, the intensity of use of information and communication technology, foreign direct investment (FDI), and urbanization have been used as explanatory variables, and carbon dioxide (CO2) emission as a dependent variable. The target countries selected are oil exporting countries (OPEC) for the period 2000 to 2020, and the analysis method used is panel VAR. The results showed that creating a shock in FDI, labor force, urban population, and renewable energy consumption decreases CO2 while creating a shock in Gross capital formation increases CO2. The impact of shock of ICT on CO2 is also insignificant and can be ignored. The results of variance analysis also showed that urban population, labor force, and FDI variables have the largest contribution in explaining the behavior of CO2; therefore, it is necessary to pay attention to FDI and try to increase the attraction of foreign direct investment to reduce CO2 in OPEC countries. JEL: C23, F43, F64.

Genome-Wide Association Mapping for Yield and Yield-Related Traits in Rice (Oryza Sativa L.) Using SNPs Markers

Rice (Oryza sativa L.) is a staple food for more than 50% of the world's population. Rice cultivar improvement is critical in order to feed the world's growing population. Improving yield is one of the main aims of rice breeders. However, yield is a complex quantitative trait controlled by many genes. The presence of genetic diversity is the key factor to improve the yield hence, the presence of diversity in any germplasm is important for yield improvement. In the current study, the rice germplasm was collected from Pakistan and the United States of America and a panel of 100 diverse genotypes was utilized to identify important yield and yield-related traits. For this, a genome-wide association study (GWAS) was performed to identify the genetic loci related to yield. The GWAS on the diverse germplasm will lead to the identification of new genes which can be utilized in the breeding program for improvement of yield. For this reason, firstly, the germplasm was phenotypically evaluated in two growing seasons for yield and yield-related traits. The analysis of variance results showed significant differences among traits which showed the presence of diversity in the current germplasm. Secondly, the germplasm was also genotypically evaluated using 10K SNP. Genetic structure analysis showed the presence of four groups which showed that enough genetic diversity was present in the rice germplasm to be used for association mapping analysis. The results of GWAS identified 201 significant marker trait associations (MTAs. 16 MTAs were identified for plant height, 49 for days to flowering, three for days to maturity, four for tillers per plant, four for panicle length, eight for grains per panicle, 20 unfilled grains per panicle, 81 for seed setting %, four for thousand-grain weight, five for yield per plot and seven for yield per hectare. Apart from this, some pleiotropic loci were also identified. The results showed that panicle length (PL) and thousand-grain weight (TGW) were controlled by a pleiotropic locus OsGRb23906 on chromosome 1 at 10,116,371 cM. The loci OsGRb25803 and OsGRb15974 on chromosomes 4 and 8 at the position of 14,321,111 cM and 6,205,816 cM respectively, showed pleiotropic effects for seed setting % (SS) and unfilled grain per panicle (UG/P). A locus OsGRb09180 on chromosome 4 at 19,850,601 cM was significantly linked with SS and yield/ha. Furthermore, gene annotation was performed, and results indicated that the 190 candidate genes or QTLs that closely linked with studied traits. These candidate genes and novel significant markers could be useful in marker-assisted gene selection and QTL pyramiding to improve rice yield and the selection of potential parents, recombinants and MTAs which could be used in rice breeding programs to develop high-yielding rice varieties for sustainable food security.

Intelligent reprogramming of wheat for enhancement of fungal and nematode disease resistance using advanced molecular techniques

Wheat (Triticum aestivum L.) diseases are major factors responsible for substantial yield losses worldwide, which affect global food security. For a long time, plant breeders have been struggling to improve wheat resistance against major diseases by selection and conventional breeding techniques. Therefore, this review was conducted to shed light on various gaps in the available literature and to reveal the most promising criteria for disease resistance in wheat. However, novel techniques for molecular breeding in the past few decades have been very fruitful for developing broad-spectrum disease resistance and other important traits in wheat. Many types of molecular markers such as SCAR, RAPD, SSR, SSLP, RFLP, SNP, and DArT, etc., have been reported for resistance against wheat pathogens. This article summarizes various insightful molecular markers involved in wheat improvement for resistance to major diseases through diverse breeding programs. Moreover, this review highlights the applications of marker assisted selection (MAS), quantitative trait loci (QTL), genome wide association studies (GWAS) and the CRISPR/Cas-9 system for developing disease resistance against most important wheat diseases. We also reviewed all reported mapped QTLs for bunts, rusts, smuts, and nematode diseases of wheat. Furthermore, we have also proposed how the CRISPR/Cas-9 system and GWAS can assist breeders in the future for the genetic improvement of wheat. If these molecular approaches are used successfully in the future, they can be a significant step toward expanding food production in wheat crops.

High-resolution quantitative trait locus mapping for rice grain quality traits using genotyping by sequencing

Rice is a major food crop that sustains approximately half of the world population. Recent worldwide improvements in the standard of living have increased the demand for high-quality rice. Accurate identification of quantitative trait loci (QTLs) for rice grain quality traits will facilitate rice quality breeding and improvement. In the present study, we performed high-resolution QTL mapping for rice grain quality traits using a genotyping-by-sequencing approach. An F₂ population derived from a cross between an elite japonica variety, Koshihikari, and an indica variety, Nona Bokra, was used to construct a high-density genetic map. A total of 3,830 single nucleotide polymorphism markers were mapped to 12 linkage groups spanning a total length of 2,456.4 cM, with an average genetic distance of 0.82 cM. Seven grain quality traits-the percentage of whole grain, percentage of head rice, percentage of area of head rice, transparency, percentage of chalky rice, percentage of chalkiness area, and degree of chalkiness-of the F₂ population were investigated. In total, 15 QTLs with logarithm of the odds (LOD) scores >4 were identified, which mapped to chromosomes 6, 7, and 9. These loci include four QTLs for transparency, four for percentage of chalky rice, four for percentage of chalkiness area, and three for degree of chalkiness, accounting for 0.01%-61.64% of the total phenotypic variation. Of these QTLs, only one overlapped with previously reported QTLs, and the others were novel. By comparing the major QTL regions in the rice genome, several key candidate genes reported to play crucial roles in grain quality traits were identified. These findings will expedite the fine mapping of these QTLs and QTL pyramiding, which will facilitate the genetic improvement of rice grain quality.

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

BACKGROUND

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

MATERIALS

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

RESULTS

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

CONCLUSION

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

HGD: an integrated homologous gene database across multiple species

Homology is fundamental to infer genes' evolutionary processes and relationships with shared ancestry. Existing homolog gene resources vary in terms of inferring methods, homologous relationship and identifiers, posing inevitable difficulties for choosing and mapping homology results from one to another. Here, we present HGD (Homologous Gene Database, https://ngdc.cncb.ac.cn/hgd), a comprehensive homologs resource integrating multi-species, multi-resources and multi-omics, as a complement to existing resources providing public and one-stop data service. Currently, HGD houses a total of 112 383 644 homologous pairs for 37 species, including 19 animals, 16 plants and 2 microorganisms. Meanwhile, HGD integrates various annotations from public resources, including 16 909 homologs with traits, 276 670 homologs with variants, 398 573 homologs with expression and 536 852 homologs with gene ontology (GO) annotations. HGD provides a wide range of omics gene function annotations to help users gain a deeper understanding of gene function.

Transposable element polymorphisms improve prediction of complex agronomic traits in rice

Transposon insertion polymorphisms can improve prediction of complex agronomic traits in rice compared to using SNPs only, especially when accessions to be predicted are less related to the training set. Transposon insertion polymorphisms (TIPs) are significant sources of genetic variation. Previous work has shown that TIPs can improve detection of causative loci on agronomic traits in rice. Here, we quantify the fraction of variance explained by single nucleotide polymorphisms (SNPs) compared to TIPs, and we explore whether TIPs can improve prediction of traits when compared to using only SNPs. We used eleven traits of agronomic relevance from by five different rice population groups (Aus, Indica, Aromatic, Japonica, and Admixed), 738 accessions in total. We assess prediction by applying data split validation in two scenarios. In the within-population scenario, we predicted performance of improved Indica varieties using the rest of Indica accessions. In the across population scenario, we predicted all Aromatic and Admixed accessions using the rest of populations. In each scenario, Bayes C and a Bayesian reproducible kernel Hilbert space regression were compared. We find that TIPs can explain an important fraction of total genetic variance and that they also improve genomic prediction. In the across population prediction scenario, TIPs outperformed SNPs in nine out of the eleven traits analyzed. In some traits like leaf senescence or grain width, using TIPs increased predictive correlation by 30-50%. Our results evidence, for the first time, that TIPs genotyping can improve prediction on complex agronomic traits in rice, especially when accessions to be predicted are less related to training accessions.

Whole-genome analysis of hard winter wheat germplasm identifies genomic regions associated with spike and kernel traits

Genetic dissection of yield component traits including spike and kernel characteristics is essential for the continuous improvement in wheat yield. Genome-wide association studies (GWAS) have been frequently used to identify genetic determinants for spike and kernel-related traits in wheat, though none have been employed in hard winter wheat (HWW) which represents a major class in US wheat acreage. Further, most of these studies relied on assembled diversity panels instead of adapted breeding lines, limiting the transferability of results to practical wheat breeding. Here we assembled a population of advanced/elite breeding lines and well-adapted cultivars and evaluated over four environments for phenotypic analysis of spike and kernel traits. GWAS identified 17 significant multi-environment marker-trait associations (MTAs) for various traits, representing 12 putative quantitative trait loci (QTLs), with five QTLs affecting multiple traits. Four of these QTLs mapped on three chromosomes 1A, 5B, and 7A for spike length, number of spikelets per spike (NSPS), and kernel length are likely novel. Further, a highly significant QTL was detected on chromosome 7AS that has not been previously associated with NSPS and putative candidate genes were identified in this region. The allelic frequencies of important quantitative trait nucleotides (QTNs) were deduced in a larger set of 1,124 accessions which revealed the importance of identified MTAs in the US HWW breeding programs. The results from this study could be directly used by the breeders to select the lines with favorable alleles for making crosses, and reported markers will facilitate marker-assisted selection of stable QTLs for yield components in wheat breeding.

Evaluation of nine statistics to identify QTLs in bulk segregant analysis using next generation sequencing approaches

Background

Bulk segregant analysis (BSA) combined with next generation sequencing is a powerful tool to identify quantitative trait loci (QTL). The impact of the size of the study population and the percentage of extreme genotypes analysed have already been assessed. But a good comparison of statistical approaches designed to identify QTL regions using next generation sequencing (NGS) technologies for BSA is still lacking.

Results

We developed an R code to simulate QTLs in bulks of F2 contrasted lines. We simulated a range of recombination rates based on estimations using different crop species. The simulations were used to benchmark the ability of statistical methods identify the exact location of true QTLs. A single QTL led to a shift in allele frequency across a large fraction of the chromosome for plant species with low recombination rate. The smoothed version of all statistics performed best notably the smoothed Euclidean distance-based statistics was always found to be more accurate in identifying the location of QTLs. We propose a simulation approach to build confidence interval statistics for the detection of QTLs.

Conclusion

We highlight the statistical methods best suited for BSA studies using NGS technologies in crops even when recombination rate is low. We also provide simulation codes to build confidence intervals and to assess the impact of recombination for application to other studies. This computational study will help select NGS-based BSA statistics that are useful to the broad scientific community.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08718-y.

Functional Allele Validation by Gene Editing to Leverage the Wealth of Genetic Resources for Crop Improvement

Advances in molecular technologies over the past few decades, such as high-throughput DNA marker genotyping, have provided more powerful plant breeding approaches, including marker-assisted selection and genomic selection. At the same time, massive investments in plant genetics and genomics, led by whole genome sequencing, have led to greater knowledge of genes and genetic pathways across plant genomes. However, there remains a gap between approaches focused on forward genetics, which start with a phenotype to map a mutant locus or QTL with the goal of cloning the causal gene, and approaches using reverse genetics, which start with large-scale sequence data and work back to the gene function. The recent establishment of efficient CRISPR-Cas-based gene editing promises to bridge this gap and provide a rapid method to functionally validate genes and alleles identified through studies of natural variation. CRISPR-Cas techniques can be used to knock out single or multiple genes, precisely modify genes through base and prime editing, and replace alleles. Moreover, technologies such as protoplast isolation, in planta transformation, and the use of developmental regulatory genes promise to enable high-throughput gene editing to accelerate crop improvement.

Genomic Regions and Floral Traits Contributing to Low Temperature Tolerance at Young Microspore Stage in a Rice (Oryza sativa L.) Recombinant Inbred Line Population of Sherpa/IRAT109

Aerobic rice production (AP) consumes less water compared to flooded systems. Developing genotypes and identifying genomic regions associated with low temperature (LT) tolerance at the young microspore stage (YMS) is imperative for AP, particularly for temperate regions. Using a recombinant inbred line population derived from the Australian LT tolerant variety Sherpa, experiments were conducted to map and dissect quantitative trait loci (QTL) associated with spikelet sterility (SS) after exposure to LT and to investigate floral traits contributing to the development of lower SS. Significant genotypic variation for SS was observed in the population after exposure to LT at YMS. Three genomic regions associated with SS, qYMCT3, qYMCT4, and qYMCT8.1 were identified in chromosomes 3, 4, and 8 respectively, using multiple QTL models explaining 22.4% of the genotypic variation. Introgression of the favorable allele from qYMCT3 was estimated to reduce SS by up to 15.4%. A co-locating genomic region with qYMCT3, qDTHW3.1 was identified as the major QTL affecting days to heading and explained as much as 44.7% of the genotypic variation. Whole-genome sequence and bioinformatic analyses demonstrated OsMADS50 as the candidate gene for qYMCT3/qDTHW3.1 and to our knowledge, this was the first attempt in connecting the role of OsMADS50 in both LT and flowering in rice. Differential sets selected for extreme SS showed LT tolerant genotype group produced higher total pollen per spikelet resulting in a higher number of dehisced anthers and pollen on stigma and eventually, lower SS than THE sensitive group. The relationship between these key floral traits with SS was induced only after exposure to LT and was not observed in warm ideal temperature conditions. Identification of elite germplasm with favorable QTL allele and combinations, gene cloning, and pyramiding with additional high-value QTL for key traits should empower breeders to develop AP adapted genotypes for temperate growing regions, and ultimately produce climate-resilient rice.

Utilization Strategies of Two Environment Phenotypes in Genomic Prediction

Multiple environment phenotypes may be utilized to implement genomic prediction in plant breeding, while it is unclear about optimal utilization strategies according to its different availability. It is necessary to assess the utilization strategies of genomic prediction models based on different availability of multiple environment phenotypes. Here, we compared the prediction accuracy of three genomic prediction models (genomic prediction model (genomic best linear unbiased prediction (GBLUP), genomic best linear unbiased prediction (GFBLUP), and multi-trait genomic best linear unbiased prediction (mtGBLUP)) which leveraged diverse information from multiple environment phenotypes using a rice dataset containing 19 agronomic traits in two disparate seasons. We found that the prediction accuracy of genomic prediction models considering multiple environment phenotypes (GFBLUP and mtGBLUP) was better than the classical genomic prediction model (GBLUP model). The deviation of prediction accuracy of between GBLUP and mtGBLUP or GFBLUP was associated with the phenotypic correlation. In summary, the genomic prediction models considering multiple environment phenotypes (GFBLUP and mtGBLUP) demonstrated better prediction accuracy. In addition, we could utilize different genomic prediction strategies according to different availability of multiple environment phenotypes.

DNA-Based Tools to Certify Authenticity of Rice Varieties—An Overview

Rice (Oryza sativa L.) is one of the most cultivated and consumed crops worldwide. It is mainly produced in Asia but, due to its large genetic pool, it has expanded to several ecosystems, latitudes and climatic conditions. Europe is a rice producing region, especially in the Mediterranean countries, that grow mostly typical japonica varieties. The European consumer interest in rice has increased over the last decades towards more exotic types, often more expensive (e.g., aromatic rice) and Europe is a net importer of this commodity. This has increased food fraud opportunities in the rice supply chain, which may deliver mixtures with lower quality rice, a problem that is now global. The development of tools to clearly identify undesirable mixtures thus became urgent. Among the various tools available, DNA-based markers are considered particularly reliable and stable for discrimination of rice varieties. This review covers aspects ranging from rice diversity and fraud issues to the DNA-based methods used to distinguish varieties and detect unwanted mixtures. Although not exhaustive, the review covers the diversity of strategies and ongoing improvements already tested, highlighting important advantages and disadvantages in terms of costs, reliability, labor-effort and potential scalability for routine fraud detection.

Genetic diversity and population structure of rice (Oryza sativa L.) landraces from Kerala, India analyzed through genotyping-by-sequencing

Researchers stand at the vanguard of advancement and application of next-generation sequencing technology for developing dominant strategies for the sustainable management of genetically diverse crops. We attempt to fill the existing research lacuna in the molecular characterization of potent rice landraces in Kerala. Genotyping-by-sequencing (GBS) was performed on 96 Kerala rice accessions to identify single-nucleotide polymorphisms (SNPs), to examine the genetic diversity, population structure, and to delineate linkage disequilibrium (LD) pattern. GBS identified 5856 high-quality SNPs. The structure analysis indicated three subpopulations with the highest probability for population clustering with significant genetic differentiation, confirmed by principal component analysis. The genome-wide LD decay distance was 772 kb, at which the r² dropped to half its maximum value. The analysis of genetic properties of the identified SNP panel with an average polymorphism information content (PIC) value of 0.22 and a minor allele frequency (MAF) > 0.1 unveiled their efficacy in genome-wide association studies (GWAS). High F_ST (0.266) and low Nm (0.692) portray a strong genetic differentiation among the rice landraces, complementing the genetic structuring observed in the studied population. Slow LD decay in the rice landraces reflects their self-pollinating behavior and the indirect selection of desired traits by domestication. Moreover, the high LD entails only a minimum number of SNP markers for detecting marker-trait association. The diverse germplasm utilized in this study can be further utilized to disclose genetic variants associated with phenotypic traits and define signatures of selection via GWAS and selective sweep, respectively.

Genome-wide association mapping and genomic prediction of yield-related traits and starch pasting properties in cassava

GWAS identified eight yield-related, peak starch type of waxy and wild-type starch and 21 starch pasting property-related traits (QTLs). Prediction ability of eight GS models resulted in low to high predictability, depending on trait, heritability, and genetic architecture. Cassava is both a food and an industrial crop in Africa, South America, and Asia, but knowledge of the genes that control yield and starch pasting properties remains limited. We carried out a genome-wide association study to clarify the molecular mechanisms underlying these traits and to explore marker-based breeding approaches. We estimated the predictive ability of genomic selection (GS) using parametric, semi-parametric, and nonparametric GS models with a panel of 276 cassava genotypes from Thai Tapioca Development Institute, International Center for Tropical Agriculture, International Institute of Tropical Agriculture, and other breeding programs. The cassava panel was genotyped via genotyping-by-sequencing, and 89,934 single-nucleotide polymorphism (SNP) markers were identified. A total of 31 SNPs associated with yield, starch type, and starch properties traits were detected by the fixed and random model circulating probability unification (FarmCPU), Bayesian-information and linkage-disequilibrium iteratively nested keyway and compressed mixed linear model, respectively. GS models were developed, and forward predictabilities using all the prediction methods resulted in values of - 0.001-0.71 for the four yield-related traits and 0.33-0.82 for the seven starch pasting property traits. This study provides additional insight into the genetic architecture of these important traits for the development of markers that could be used in cassava breeding programs.

Genetic and phenotypic characterization of rice grain quality traits to define research strategies for improving rice milling, appearance, and cooking qualities in Latin America and the Caribbean

Rice (Oryza sativa L.)grain quality is a set of complex interrelated traits that include grain milling, appearance, cooking, and edible properties. As consumer preferences in Latin America and the Caribbean evolve, determining what traits best capture regional grain quality preferences is fundamental for breeding and cultivar release. In this study, a genome-wide association study (GWAS), marker-assisted selection (MAS), and genomic selection (GS) were evaluated to help guide the development of new breeding strategies for rice grain quality improvement. For this purpose, 284 rice lines representing over 20 yr of breeding in Latin America and the Caribbean were genotyped and phenotyped for 10 different traits including grain milling, appearance, cooking, and edible quality traits. Genetic correlations among the 10 traits ranged from -0.83 to 0.85. A GWAS identified 19 significant marker/trait combinations associated with eight grain quality traits. Four functional markers, three located in the Waxy and one in the starch synthase IIa genes, were significantly associated with six grain-quality traits. These markers individually explained 51-75% of the phenotypic variance depending on the trait, clearly indicating their potential utility for MAS. Cross-validation studies to evaluate predictive abilities of four different GS models for each of the 10 quality traits were conducted and predictive abilities ranged from 0.3 to 0.72. Overall, the machine learning model random forest had the highest predictive abilities and was especially effective for traits where large effect quantitative trait loci were identified. This study provides the foundation for deploying effective molecular breeding strategies for grain quality in Latin American rice breeding programs.

Advances and Challenges for QTL Analysis and GWAS in the Plant-Breeding of High-Yielding: A Focus on Rapeseed

Yield is one of the most important agronomic traits for the breeding of rapeseed (Brassica napus L), but its genetic dissection for the formation of high yield remains enigmatic, given the rapid population growth. In the present review, we review the discovery of major loci underlying important agronomic traits and the recent advancement in the selection of complex traits. Further, we discuss the benchmark summary of high-throughput techniques for the high-resolution genetic breeding of rapeseed. Biparental linkage analysis and association mapping have become powerful strategies to comprehend the genetic architecture of complex agronomic traits in crops. The generation of improved crop varieties, especially rapeseed, is greatly urged to enhance yield productivity. In this sense, the whole-genome sequencing of rapeseed has become achievable to clone and identify quantitative trait loci (QTLs). Moreover, the generation of high-throughput sequencing and genotyping techniques has significantly enhanced the precision of QTL mapping and genome-wide association study (GWAS) methodologies. Furthermore, this study demonstrates the first attempt to identify novel QTLs of yield-related traits, specifically focusing on ovule number per pod (ON). We also highlight the recent breakthrough concerning single-locus-GWAS (SL-GWAS) and multi-locus GWAS (ML-GWAS), which aim to enhance the potential and robust control of GWAS for improved complex traits.

A high-resolution genome-wide association study of the grain ionome and agronomic traits in rice Oryza sativa subsp. indica

This study presents a comprehensive study of the genetic bases controlling variation in the rice ionome employing genome-wide association studies (GWAS) with a diverse panel of indica accessions, each genotyped with 5.2 million markers. GWAS was performed for twelve elements including B, Ca, Co, Cu, Fe, K, Mg, Mn, Mo, Na, P, and Zn and four agronomic traits including days to 50% flowering, grain yield, plant height and thousand grain weight. GWAS identified 128 loci associated with the grain elements and 57 associated with the agronomic traits. There were sixteen co-localization regions containing QTL for two or more traits. Fourteen grain element quantitative trait loci were stable across growing environments, which can be strong candidates to be used in marker-assisted selection to improve the concentrations of nutritive elements in rice grain. Potential candidate genes were revealed including OsNAS3 linked to the locus that controls the variation of Zn and Co concentrations. The effects of starch synthesis and grain filling on multiple grain elements were elucidated through the likely involvement of OsSUS1 and OsGSSB1 genes. Overall, our study provides crucial insights into the genetic basis of ionomic variations in rice and will facilitate improvement in breeding for trace mineral content.

Genome-Wide Association Mapping of Mixed Linkage (1,3;1,4)-β-Glucan and Starch Contents in Rice Whole Grain

The glucan content of rice is a key factor defining its nutritional and economic value. Starch and its derivatives have many industrial applications such as in fuel and material production. Non-starch glucans such as (1,3;1,4)-β-D-glucan (mixed-linkage β-glucan, MLG) have many benefits in human health, including lowering cholesterol, boosting the immune system, and modulating the gut microbiome. In this study, the genetic variability of MLG and starch contents were analyzed in rice (Oryza sativa L.) whole grain, by performing a new quantitative analysis of the polysaccharide content of rice grains. The 197 rice accessions investigated had an average MLG content of 252 μg/mg, which was negatively correlated with the grain starch content. A new genome-wide association study revealed seven significant quantitative trait loci (QTLs) associated with the MLG content and two QTLs associated with the starch content in rice whole grain. Novel genes associated with the MLG content were a hexose transporter and anthocyanidin 5,3-O-glucosyltransferase. Also, the novel gene associated with the starch content was a nodulin-like domain. The data pave the way for a better understanding of the genes involved in determining both MLG and starch contents in rice grains and should facilitate future plant breeding programs.

Marker-Assisted Introgression and Stacking of Major QTLs Controlling Grain Number (Gn1a) and Number of Primary Branching (WFP) to NERICA Cultivars

The era of the green revolution has significantly improved rice yield productivity. However, with the growing population and decreasing arable land, rice scientists must find new ways to improve rice productivity. Although hundreds of rice yield-related QTLs were already mapped and some of them were cloned, only a few were utilized for actual systematic introgression breeding programs. In this study, the major yield QTLs Grain Number 1a (Gn1a) and Wealthy Farmer’s Panicle (WFP) were introgressed and stacked in selected NERICA cultivars by marker-assisted backcross breeding (MABB). The DNA markers RM3360, RM3452, and RM5493 were used for foreground selection. At BC₃F₄ and BC₃F₅ generation, a combination of marker-assisted selection and phenotypic evaluation were carried out to select lines with target alleles and traits. Further, genotyping-by-sequencing (GBS) was conducted to validate the introgression and determine the recurrent parent genome recovery (RPGR) of the selected lines. The Gn1a and/or WFP introgression lines showed significantly higher numbers of spikelets per panicle and primary branching compared to the recurrent parents. In addition, lines with Gn1a and/or WFP alleles were comparatively similar to the recurrent parents (RP) in most yield-related traits. This study demonstrates the success of utilizing yield QTLs and marker-assisted selection to develop and improve rice cultivars.

Assessing Rice Salinity Tolerance: From Phenomics to Association Mapping

Rice is the most salt-sensitive cereal, suffering yield losses above 50% with soil salinity of 6 dS/m. Thus, understanding the mechanisms of rice salinity tolerance is key to address food security. In this chapter, we provide guidelines to assess rice salinity tolerance using a high-throughput phenotyping platform (HTP) with digital imaging at seedling/early tillering stage and suggest improved analysis methods using stress indices. The protocols described here also include computer scripts for users to improve their experimental design, run genome-wide association studies (GWAS), perform multi-testing corrections, and obtain the Manhattan plots, enabling the identification of loci associated with salinity tolerance. Notably, the computer scripts provided here can be used for any stress or GWAS experiment and independently of HTP.

Genome-wide association studies of fruit quality traits in jujube germplasm collections using genotyping-by-sequencing

Chinese jujube (Ziziphus jujuba Mill.) is an important fruit crop and harbors many highly diverse traits of potential economic importance. Fruit size, stone size, and fruit cracking have an important influence on the commercial value of jujube. This study is the first to conduct a genome-wide association study (GWAS) on 180 accessions of jujube and focuses on locating single-nucleotide polymorphisms (SNPs) associated with nine important fruit quality traits. Genotyping was performed using genotyping-by-sequencing and 4651 high-quality SNPs were identified. A genetic diversity analysis revealed the presence of three distinct groups, and rapid linkage disequilibrium decay was observed in this jujube population. Using a mixed linear model, a total of 45 significant SNP-trait associations were detected, among which 33 SNPs had associations with fruit size-related traits, nine were associated with stone size-related traits, and three with fruit cracking-related traits. In total, 21 candidate genes involved in cell expansion, abiotic stress responses, hormone signaling, and growth development were identified from the genome sequences of jujube. These results are useful as basic data for GWAS of other jujube traits, and these significant SNP loci and candidate genes should aid marker-assisted breeding and genomic selection of improved jujube cultivars.

PANOMICS meets germplasm

Genotyping-by-sequencing has enabled approaches for genomic selection to improve yield, stress resistance and nutritional value. More and more resource studies are emerging providing 1000 and more genotypes and millions of SNPs for one species covering a hitherto inaccessible intraspecific genetic variation. The larger the databases are growing, the better statistical approaches for genomic selection will be available. However, there are clear limitations on the statistical but also on the biological part. Intraspecific genetic variation is able to explain a high proportion of the phenotypes, but a large part of phenotypic plasticity also stems from environmentally driven transcriptional, post-transcriptional, translational, post-translational, epigenetic and metabolic regulation. Moreover, regulation of the same gene can have different phenotypic outputs in different environments. Consequently, to explain and understand environment-dependent phenotypic plasticity based on the available genotype variation we have to integrate the analysis of further molecular levels reflecting the complete information flow from the gene to metabolism to phenotype. Interestingly, metabolomics platforms are already more cost-effective than NGS platforms and are decisive for the prediction of nutritional value or stress resistance. Here, we propose three fundamental pillars for future breeding strategies in the framework of Green Systems Biology: (i) combining genome selection with environment-dependent PANOMICS analysis and deep learning to improve prediction accuracy for marker-dependent trait performance; (ii) PANOMICS resolution at subtissue, cellular and subcellular level provides information about fundamental functions of selected markers; (iii) combining PANOMICS with genome editing and speed breeding tools to accelerate and enhance large-scale functional validation of trait-specific precision breeding.

PANOMICS meets germplasm

Genotyping-by-sequencing has enabled approaches for genomic selection to improve yield, stress resistance and nutritional value. More and more resource studies are emerging providing 1000 and more genotypes and millions of SNPs for one species covering a hitherto inaccessible intraspecific genetic variation. The larger the databases are growing, the better statistical approaches for genomic selection will be available. However, there are clear limitations on the statistical but also on the biological part. Intraspecific genetic variation is able to explain a high proportion of the phenotypes, but a large part of phenotypic plasticity also stems from environmentally driven transcriptional, post-transcriptional, translational, post-translational, epigenetic and metabolic regulation. Moreover, regulation of the same gene can have different phenotypic outputs in different environments. Consequently, to explain and understand environment-dependent phenotypic plasticity based on the available genotype variation we have to integrate the analysis of further molecular levels reflecting the complete information flow from the gene to metabolism to phenotype. Interestingly, metabolomics platforms are already more cost-effective than NGS platforms and are decisive for the prediction of nutritional value or stress resistance. Here, we propose three fundamental pillars for future breeding strategies in the framework of Green Systems Biology: (i) combining genome selection with environment-dependent PANOMICS analysis and deep learning to improve prediction accuracy for marker-dependent trait performance; (ii) PANOMICS resolution at subtissue, cellular and subcellular level provides information about fundamental functions of selected markers; (iii) combining PANOMICS with genome editing and speed breeding tools to accelerate and enhance large-scale functional validation of trait-specific precision breeding.

Role of NGS and SNP genotyping methods in sugarcane improvement programs

Sugarcane (Saccharum spp.) is one of the most economically significant crops because of its high sucrose content and it is a promising biomass feedstock for biofuel production. Sugarcane genome sequencing and analysis is a difficult task due to its heterozygosity and polyploidy. Long sequence read technologies, PacBio Single-Molecule Real-Time (SMRT) sequencing, the Illumina TruSeq, and the Oxford Nanopore sequencing could solve the problem of genome assembly. On the applications side, next generation sequencing (NGS) technologies played a major role in the discovery of single nucleotide polymorphism (SNP) and the development of low to high throughput genotyping platforms. The two mainstream high throughput genotyping platforms are the SNP microarray and genotyping by sequencing (GBS). This paper reviews the NGS in sugarcane genomics, genotyping methodologies, and the choice of these methods. Array-based SNP genotyping is robust, provides consistent SNPs, and relatively easier downstream data analysis. The GBS method identifies large scale SNPs across the germplasm. A combination of targeted GBS and array-based genotyping methods should be used to increase the accuracy of genomic selection and marker-assisted breeding.

Recommendations for Choosing the Genotyping Method and Best Practices for Quality Control in Crop Genome-Wide Association Studies

High-throughput genotyping boosts genome-wide association studies (GWAS) in crop species, leading to the identification of single-nucleotide polymorphisms (SNPs) associated with economically important traits. Choosing a cost-effective genotyping method for crop GWAS requires careful examination of several aspects, namely, the purpose and the scale of the study, crop-specific genomic features, and technical and economic matters associated with each genotyping option. Once genotypic data have been obtained, quality control (QC) procedures must be applied to avoid bias and false signals in genotype–phenotype association tests. QC for human GWAS has been extensively reviewed; however, QC for crop GWAS may require different actions, depending on the GWAS population type. Here, we review most popular genotyping methods based on next-generation sequencing (NGS) and array hybridization, and report observations that should guide the investigator in the choice of the genotyping method for crop GWAS. We provide recommendations to perform QC in crop species, and deliver an overview of bioinformatics tools that can be used to accomplish all needed tasks. Overall, this work aims to provide guidelines to harmonize those procedures leading to SNP datasets ready for crop GWAS.

Choosing the optimal population for a genome-wide association study: A simulation of whole-genome sequences from rice

A genome-wide association study (GWAS) needs to have a suitable population. The factors that affect a GWAS (e.g. population structure, sample size, and sequence analysis and field testing costs) need to be considered. Mixed populations containing subpopulations of different genetic backgrounds may be suitable populations. We conducted simulation experiments to see if a population with high genetic diversity, such as a diversity panel, should be added to a target population, especially when the target population harbors small genetic diversity. The target population was 112 accessions of Oryza sativa L. subsp. japonica, mainly developed in Japan. We combined the target population with three populations that had higher genetic diversity. These were 100 indica accessions, 100 japonica accessions, and 100 accessions with various genetic backgrounds. The results showed that the GWAS's power with a mixed population was generally higher than with a separate population. Also, the optimal GWAS populations varied depending on the fixation index (F_ST ) of the quantitative trait nucleotides (QTNs) and the polymorphism of QTNs in each population. When a QTN was polymorphic in a target population, a target population combined with a higher diversity population improved the QTN's detection power. By investigating F_ST and the expected heterozygosity (H_e ) as factors influencing the detection power, we showed that single nucleotide polymorphisms with high F_ST or low H_e are less likely to be detected by GWAS with mixed populations. Sequenced or genotyped germplasm collections can improve the GWAS's detection power by using a subset of the collections with a target population.

Genome-Wide Associations and Transcriptional Profiling Reveal ROS Regulation as One Underlying Mechanism of Sheath Blight Resistance in Rice

Rice sheath blight, caused by the necrotrophic fungus Rhizoctonia solani Kühn, continues to be an important and challenging rice disease worldwide. Here, we used genome-wide association studies over a high-density rice array to facilitate the identification of potential novel genes and quantitative trait loci related to sheath blight resistance. We identified multiple regions that significantly associated with independent disease components in chromosomes 1, 4, and 11 under controlled condition. In particular, we investigated qLN11²⁸, a quantitative trait locus enriched with defense-related genes that reduce disease lesions in a near-isogenic line. RNA profiling of the line carrying qLN11²⁸ showed a number of differentially expressed genes related to the reactive oxygen species (ROS)-redox pathway. Histochemical staining revealed less ROS accumulation on the resistant line, suggesting efficient ROS deregulation that delays pathogen colonization. The detection of genomic regions controlling multiple mechanisms of resistance to sheath blight will provide tools to design effective breeding interventions in rice.

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.

Multiple Small-Effect Alleles of Indica Origin Enhance High Iron-Associated Stress Tolerance in Rice Under Field Conditions in West Africa

Understanding the genetics of field-based tolerance to high iron-associated (HIA) stress in rice can accelerate the development of new varieties with enhanced yield performance in West African lowland ecosystems. To date, few field-based studies have been undertaken to rigorously evaluate rice yield performance under HIA stress conditions. In this study, two NERICA × O. sativa bi-parental rice populations and one O.sativa diversity panel consisting of 296 rice accessions were evaluated for grain yield and leaf bronzing symptoms over multiple years in four West African HIA stress and control sites. Mapping of these traits identified a large number of QTLs and single nucleotide polymorphisms (SNPs) associated with stress tolerance in the field. Favorable alleles associated with tolerance to high levels of iron in anaerobic rice soils were rare and almost exclusively derived from the indica subpopulation, including the most favorable alleles identified in NERICA varieties. These findings highlight the complex genetic architecture underlying rice response to HIA stress and suggest that a recurrent selection program focusing on an expanded indica genepool could be productively used in combination with genomic selection to increase the efficiency of selection in breeding programs designed to enhance tolerance to this prevalent abiotic stress in West Africa.

Exploring genetic architecture of grain yield and quality traits in a 16-way indica by japonica rice MAGIC global population

Identification of Quantitative Trait Loci (QTL) has been a challenge for complex traits due to the use of populations with narrow genetic base. Most of QTL mapping studies were carried out from crosses made within the subspecies, either indica × indica or japonica × japonica. In this study we report advantages of using Multi-parent Advanced Generation Inter-Crosses global population, derived from a combination of eight indica and eight japonica elite parents, in QTL discovery for yield and grain quality traits. Genome-wide association study and interval mapping identified 38 and 34 QTLs whereas Bayesian networking detected 60 QTLs with 22 marker-marker associations, 32 trait-trait associations and 65 marker-trait associations. Notably, nine known QTLs/genes qPH₁/OsGA20ox2, qDF₃/OsMADS50, PL, QDg1, qGW-5b, grb7-2, qGL₃/GS3, Amy6/Wx gene and OsNAS3 were consistently identified by all approaches for nine traits whereas qDF₃/OsMADS50 was co-located for both yield and days-to-flowering traits on chromosome 3. Moreover, we identified a number of candidate QTLs in either one or two analyses but further validations will be needed. The results indicate that this new population has enabled identifications of significant QTLs and interactions for 16 traits through multiple approaches. Pyramided recombinant inbred lines provide a valuable source for integration into future breeding programs.

Association genetics of bunch weight and its component traits in East African highland banana (Musa spp. AAA group)

The major quantitative trait loci associated with bunch weight and its component traits in the East African highland banana-breeding population are located on chromosome 3. Bunch weight increase is one of the major objectives of banana improvement programs, but little is known about the loci controlling bunch weight and its component traits. Here we report for the first time some genomic loci associated with bunch weight and its component traits in banana as revealed through a genome-wide association study. A banana-breeding population of 307 genotypes varying in ploidy was phenotyped in three locations under different environmental conditions, and data were collected on bunch weight, number of hands and fruits; fruit length and circumference; and diameter of both fruit and pulp for three crop cycles. The population was genotyped with genotyping by sequencing and 27,178 single nucleotide polymorphisms (SNPs) were generated. The association between SNPs and the best linear unbiased predictors of traits was performed with TASSEL v5 using a mixed linear model accounting for population structure and kinship. Using Bonferroni correction, false discovery rate, and long-range linkage disequilibrium (LD), 25 genomic loci were identified with significant SNPs and most were localized on chromosome 3. Most SNPs were located in genes encoding uncharacterized and hypothetical proteins, but some mapped to transcription factors and genes involved in cell cycle regulation. Inter-chromosomal LD of SNPs was present in the population, but none of the SNPs were significantly associated with the traits. The clustering of significant SNPs on chromosome 3 supported our hypothesis that fruit filling in this population was under control of a few quantitative trait loci with major effects.

Reference Genome Anchoring of High-Density Markers for Association Mapping and Genomic Prediction in European Winter Wheat

In this study, we anchored genotyping-by-sequencing data to the International Wheat Genome Sequencing Consortium Reference Sequence v1.0 assembly to generate over 40,000 high quality single nucleotide polymorphism markers on a panel of 376 elite European winter wheat varieties released between 1946 and 2007. We compared association mapping and genomic prediction accuracy for a range of productivity traits with previous results based on lower density dominant DArT markers. The results demonstrate that the availability of RefSeq v1.0 supports higher precision trait mapping and provides the density of markers required to obtain accurate predictions of traits controlled by multiple small effect loci, including grain yield.

Genome-wide association study for leaf area, rachis length and total dry weight in oil palm (Eleaeisguineensis) using genotyping by sequencing

The marker-trait association for complex traits using genotyping by sequencing (GBS) method is being widely spread in plants. The study aimed to identify significant single nucleotide polymorphism (SNP) associations for rachis length (RL), leaf area (LA) and total dry weight (TrDW) in oil palm among diverse African germplasm. The Illumina NextSeq platform has been used for SNP genotyping and retained 4031 fully informative SNPs after applying the filter criterion. These 4031 SNPs were used for genome wide association study for the above three traits. The LD decay rates of the African germplasm using GBS data of SNP is observed to be 25 Kb at 0.45 of average pair wise correlation coefficient (r²). Association mapping led to the identification of seven significant associations for three traits using MLM approach at a P value of ≤ 0.001. Three associations were identified for total dry weight, two each for leaf area index and rachis length. The qtlLA1 was found to be highly significant at a P value of 7.39E^-05 (18.4% phenotypic variance) which is located on chromosome 4. Two QTLs (qtlLA2 and qtlRL1) were located on chromosome 1, which explained 11.9% and 12.4% of phenotypic variance respectively. Three QTLs for total dry weight were located on chromosome 2, 14 and 16, all-together explained 40% phenotypic variance. The results showed that the SNP-trait associations identified in the present study could be used in selection of elite oil palm germplasm for higher yields.

1k-RiCA (1K-Rice Custom Amplicon) a novel genotyping amplicon-based SNP assay for genetics and breeding applications in rice

BACKGROUND

While a multitude of genotyping platforms have been developed for rice, the majority of them have not been optimized for breeding where cost, turnaround time, throughput and ease of use, relative to density and informativeness are critical parameters of their utility. With that in mind we report the development of the 1K-Rice Custom Amplicon, or 1k-RiCA, a robust custom sequencing-based amplicon panel of ~ 1000-SNPs that are uniformly distributed across the rice genome, designed to be highly informative within indica rice breeding pools, and tailored for genomic prediction in elite indica rice breeding programs.

RESULTS

Empirical validation tests performed on the 1k-RiCA showed average marker call rates of 95% with marker repeatability and concordance rates of 99%. These technical properties were not affected when two common DNA extraction protocols were used. The average distance between SNPs in the 1k-RiCA was 1.5 cM, similar to the theoretical distance which would be expected between 1,000 uniformly distributed markers across the rice genome. The average minor allele frequencies on a panel of indica lines was 0.36 and polymorphic SNPs estimated on pairwise comparisons between indica by indica accessions and indica by japonica accessions were on average 430 and 450 respectively. The specific design parameters of the 1k-RiCA allow for a detailed view of genetic relationships and unambiguous molecular IDs within indica accessions and good cost vs. marker-density balance for genomic prediction applications in elite indica germplasm. Predictive abilities of Genomic Selection models for flowering time, grain yield, and plant height were on average 0.71, 0.36, and 0.65 respectively based on cross-validation analysis. Furthermore the inclusion of important trait markers associated with 11 different genes and QTL adds value to parental selection in crossing schemes and marker-assisted selection in forward breeding applications.

CONCLUSIONS

This study validated the marker quality and robustness of the 1k-RiCA genotypic platform for genotyping populations derived from indica rice subpopulation for genetic and breeding purposes including MAS and genomic selection. The 1k-RiCA has proven to be an alternative cost-effective genotyping system for breeding applications.

Molecular characterization of bacterial leaf streak resistance in hard winter wheat

Bacterial leaf streak (BLS) caused by Xanthomonas campestris pv. translucens is one of the major bacterial diseases threatening wheat production in the United States Northern Great Plains (NGP) region. It is a sporadic but widespread wheat disease that can cause significant loss in grain yield and quality. Identification and characterization of genomic regions in wheat that confer resistance to BLS will help track resistance genes/QTLs in future wheat breeding. In this study, we evaluated a hard winter wheat association mapping panel (HWWAMP) containing 299 hard winter wheat lines from the US hard winter wheat growing region for their reactions to BLS. We observed a range of BLS responses among the lines, importantly, we identified ten genotypes that showed a resistant reaction both in greenhouse and field evaluation. -Genome-wide association analysis with 15,990 SNPs was conducted using an exponentially compressed mixed linear model. Five genomic regions (p < 0.001) that regulate the resistance to BLS were identified on chromosomes 1AL, 1BS, 3AL, 4AL, and 7AS. The QTLs Q.bls.sdsu-1AL, Q.bls.sdsu-1BS, Q.bls.sdsu-3AL, Q.bls.sdsu-4AL, and Q.bls.sdsu-7AS explain a total of 42% of the variation. In silico analysis of sequences in the candidate regions on chromosomes 1AL, 1BS, 3AL, 4AL, and 7AS identified 10, 25, 22, eight, and nine genes, respectively with known plant defense-related functions. Comparative analysis with rice showed two syntenic regions in rice that harbor genes for bacterial leaf streak resistance. The ten BLS resistant genotypes and SNP markers linked to the QTLs identified in our study could facilitate breeding for BLS resistance in winter wheat.

Shovelomics for phenotyping root architectural traits of rapeseed/canola (Brassica napus L.) and genome-wide association mapping

Root system in plants plays an important role in mining moisture and nutrients from the soil and is positively correlated to yield in many crops including rapeseed/canola (Brassica napus L.). Substantial phenotypic diversity in root architectural traits among the B. napus growth types leads to a scope of root system improvement in breeding populations. In this study, 216 diverse genotypes were phenotyped for five different root architectural traits following shovelomics approach in the field condition during 2015 and 2016. A single nucleotide polymorphism (SNP) marker panel consisting of 30,262 SNPs was used to conduct genome-wide association study to detect marker/trait association. A total of 31 significant marker loci were identified at 0.01 percentile tail P value cutoff for different root traits. Six marker loci for soil-level taproot diameter (R₁Dia), six loci for belowground taproot diameter (R₂Dia), seven loci for number of primary root branches (PRB), eight loci for root angle, and eight loci for root score (RS) were detected in this study. Several markers associated with root diameters R₁Dia and R₂Dia were also associated with PRB and RS. Significant phenotypic correlation between these traits was observed in both environments. Therefore, taproot diameter appears to be a major determinant of the canola root system architecture and can be used as proxy for other root traits. Fifteen candidate genes related to root traits and root development were detected within 100 kbp upstream and downstream of different significant markers. The identified markers associated with different root architectural traits can be considered for marker-assisted selection for root traits in canola in future.

Edit at will: Genotype independent plant transformation in the era of advanced genomics and genome editing

The combination of advanced genomics, genome editing and plant transformation biology presents a powerful platform for basic plant research and crop improvement. Together these advances provide the tools to identify genes as targets for direct editing as single base pair changes, deletions, insertions and site specific homologous recombination. Recent breakthrough technologies using morphogenic regulators in plant transformation creates the ability to introduce reagents specific toward their identified targets and recover stably transformed and/or edited plants which are genotype independent. These technologies enable the possibility to alter a trait in any variety, without genetic disruption which would require subsequent extensive breeding, but rather to deliver the same variety with one trait changed. Regulatory issues regarding this technology will predicate how broadly these technologies will be implemented. In addition, education will play a crucial role for positive public acceptance. Taken together these technologies comprise a platform for advanced breeding which is an imperative for future world food security.

Distinguishing CPT gene family members and vetting the sequence structure of a putative rubber synthesizing variant in Hevea brasiliensis

cis-Prenyltransferases (cis-PTs) constitute a large family of enzymes conserved during evolution and present in all domains of life. cis-PTs catalyze the cis-1,4-polymerization of isoprene units to generate isoprenoids with carbon skeletons varying from C10 (neryl pyrophosphate) to C > 10,000 (natural rubber). Though the previously reported CPTs in Hevea are designated based on sequence variations, their classification was done mostly by phylogenetic analysis using a mixture of partial as well as full length sequences often excluding the UTRs. In this context an attempt was made to reclassify the CPTs strictly based on their sequence similarity and distinguish the members putatively associated with rubber biosynthesis from the others. Extensive computational analysis was carried out on CPT sequences obtained from public resources and whole genome assemblies of Hevea. Based on the results from BLAST analysis, multiple sequence alignments of protein, nucleotide and untranslated regions, open reading frame analysis, gene prediction analysis and sequence length variations, we conclude that there exists mainly three CPTs namely RubCPT1, RubCPT2 and RubCPT3 putatively associated with rubber biosynthesis in Hevea brasiliensis. The rest were categorised as variants of dehydrodolichyl diphosphate synthase (DHDDS) involved in the synthesis of dolichols having short chain isoprenoids. Analysis of the sequence structure of the most highly expressed RubCPT1 in latex revealed the allele richness and diversity of this important variant prevailing in the popular rubber clones. Haplotypes consisting of SNPs with high degree of heterozygosity were also identified. Segregation and linkage disequilibrium analysis confirmed that recombination is the major contributor towards the generation of allelic diversity rather than point mutations. Alternatively, gene expression analysis indicated the possibility of association between specific haplotypes and RubCPT1 expression in Hevea clones which may have downstream impact up to the level of rubber production. The conclusions from this study may pave way for the identification and better understanding of CPTs directly involved with natural rubber biosynthesis in Hevea and the SNP data generated may aid in the development of molecular markers putatively associated with yield in rubber.

Validation of Yield Component Traits Identified by Genome-Wide Association Mapping in a tropical japonica × tropical japonica Rice Biparental Mapping Population

The Rice Diversity Panel 1 (RDP1) was developed for genome-wide association (GWA) studies to explore five rice ( L.) subpopulations (, , , , and ). The RDP1 was evaluated for over 30 traits, including agronomic, panicle architecture, seed, and disease traits and genotyped with 700,000 single nucleotide polymorphisms (SNPs). Most rice grown in the southern United States is and thus the diversity in this subpopulation is interesting to U.S. breeders. Among the RDP1 accessions, 'Estrela' and 'NSFTV199' are both phenotypically and genotypically diverse, thus making them excellent parents for a biparental mapping population. The objectives were to (i) ascertain the GWA QTLs from the RDP1 GWA studies that overlapped with the QTLs uncovered in an Estrela × NSFTV199 recombinant inbred line (RIL) population evaluated for 15 yield traits, and (ii) identify known or novel genes potentially controlling specific yield component traits. The 256 RILs were genotyped with 132 simple sequence repeat markers and 70 QTLs were found. Perl scripts were developed for automatic identification of the underlying candidate genes in the GWA QTL regions. Approximately 100 GWA QTLs overlapped with 41 Estrela × NSFTV199 QTL (RIL QTL) regions and 47 known genes were identified. Two seed trait RIL QTLs with overlapping GWA QTLs were not associated with a known gene. Segregating SNPs in the overlapping GWA QTLs for RIL QTLs with high values will be evaluated as potential DNA markers useful to breeding programs for the associated yield trait.

Criteria for evaluating molecular markers: Comprehensive quality metrics to improve marker-assisted selection

Despite strong interest over many years, the usage of quantitative trait loci in plant breeding has often failed to live up to expectations. A key weak point in the utilisation of QTLs is the “quality” of markers used during marker-assisted selection (MAS): unreliable markers result in variable outcomes, leading to a perception that MAS products fail to achieve reliable improvement. Most reports of markers used for MAS focus on markers derived from the mapping population. There are very few studies that examine the reliability of these markers in other genetic backgrounds, and critically, no metrics exist to describe and quantify this reliability. To improve the MAS process, this work proposes five core metrics that fully describe the reliability of a marker. These metrics give a comprehensive and quantitative measure of the ability of a marker to correctly classify germplasm as QTL[+]/[–], particularly against a background of high allelic diversity. Markers that score well on these metrics will have far higher reliability in breeding, and deficiencies in specific metrics give information on circumstances under which a marker may not be reliable. The metrics are applicable across different marker types and platforms, allowing an objective comparison of the performance of different markers irrespective of the platform. Evaluating markers using these metrics demonstrates that trait-specific markers consistently out-perform markers designed for other purposes. These metrics also provide a superb set of criteria for designing superior marker systems for a target QTL, enabling the selection of an optimal marker set before committing to design.