Population improvement via recurrent selection drives genetic gain in upland rice breeding.

One of the main challenges of breeding programs is to identify superior genotypes from a large number of candidates. By gradually increasing the frequency of favorable alleles in the breeding population, recurrent selection improves the population mean for target traits, increasing the chance to identify promising genotypes. In rice, population improvement through recurrent selection has been used very little to date, except in Latin America. At Embrapa (Brazilian Agricultural Research Corporation), the upland rice breeding program is conducted in two phases: population improvement followed by product development. In this study, the CNA6 population, evaluated over five cycles (3 to 7) of selection, including 20 field trials, was used to assess the realized genetic gain. A high rate of genetic gain was observed for grain yield, at 215 kg.ha-1 per cycle or 67.8 kg.ha-1 per year (3.08%). The CNA6 population outperformed the controls only for the last cycle, with a yield difference of 1128 kg.ha-1. An analysis of the product development pipeline, based on 29 advanced yield trials with lines derived from cycles 3 to 6, showed that lines derived from the CNA6 population had high grain yield, but did not outperform the controls. These results demonstrate that the application of recurrent selection to a breeding population with sufficient genetic variability can result in significant genetic gains for quantitative traits, such as grain yield. The integration of this strategy into a two-phase breeding program also makes it possible to increase quantitative traits while selecting for other traits of interest.

The pulse of the tree is under genetic control: eucalyptus as a case study.

The pulse of the tree (diurnal cycle of stem radius fluctuations) has been widely studied as a way of analyzing tree responses to the environment, including the phenotypic plasticity of tree-water relationships in particular. However, the genetic basis of this daily phenotype and its interplay with the environment remain largely unexplored. We characterized the genetic and environmental determinants of this response, by monitoring daily stem radius fluctuation (dSRF) on 210 trees from a Eucalyptus urophylla × E. grandis full-sib family over 2 years. The dSRF signal was broken down into hydraulic capacitance, assessed as the daily amplitude of shrinkage (DA), and net growth, estimated as the change in maximum radius between two consecutive days (ΔR). The environmental determinants of these two traits were clearly different: DA was positively correlated with atmospheric variables relating to water demand, while ΔR was associated with soil water content. The heritability for these two traits ranged from low to moderate over time, revealing a time-dependent or environment-dependent complex genetic determinism. We identified 686 and 384 daily quantitative trait loci (QTL) representing 32 and 31 QTL regions for DA and ΔR, respectively. The identification of gene networks underlying the 27 major genomics regions for both traits generated additional hypotheses concerning the biological mechanisms involved in response to water demand and supply. This study highlights that environmentally induced changes in daily stem radius fluctuation are genetically controlled in trees and suggests that these daily responses integrated over time shape the genetic architecture of mature traits.

The genetics of exapted resistance to two exotic pathogens in pedunculate oak.

Exotic pathogens cause severe damage in natural populations in the absence of coevolutionary dynamics with their hosts. However, some resistance to such pathogens may occur in naive populations. The objective of this study was to investigate the genetics of this so-called 'exapted' resistance to two pathogens of Asian origin (Erysiphe alphitoides and Phytophthora cinnamomi) in European oak. Host-pathogen compatibility was assessed by recording infection success and pathogen growth in a full-sib family of Quercus robur under controlled and natural conditions. Two high-resolution genetic maps anchored on the reference genome were used to study the genetic architecture of resistance and to identify positional candidate genes. Two genomic regions, each containing six strong and stable quantitative trait loci (QTLs) accounting for 12-19% of the phenotypic variation, were mainly associated with E. alphitoides infection. Candidate genes, especially genes encoding receptor-like-kinases and galactinol synthases, were identified in these regions. The three QTLs associated with P. cinnamomi infection did not colocate with QTLs found for E. alphitoides. These findings provide evidence that exapted resistance to E. alphitoides and P. cinnamomi is present in Q. robur and suggest that the underlying molecular mechanisms involve genes encoding proteins with extracellular signaling functions.

Performance of genomic prediction within and across generations in maritime pine.

BACKGROUND: Genomic selection (GS) is a promising approach for decreasing breeding cycle length in forest trees. Assessment of progeny performance and of the prediction accuracy of GS models over generations is therefore a key issue. RESULTS: A reference population of maritime pine (Pinus pinaster) with an estimated effective inbreeding population size (status number) of 25 was first selected with simulated data. This reference population (n = 818) covered three generations (G0, G1 and G2) and was genotyped with 4436 single-nucleotide polymorphism (SNP) markers. We evaluated the effects on prediction accuracy of both the relatedness between the calibration and validation sets and validation on the basis of progeny performance. Pedigree-based (best linear unbiased prediction, ABLUP) and marker-based (genomic BLUP and Bayesian LASSO) models were used to predict breeding values for three different traits: circumference, height and stem straightness. On average, the ABLUP model outperformed genomic prediction models, with a maximum difference in prediction accuracies of 0.12, depending on the trait and the validation method. A mean difference in prediction accuracy of 0.17 was found between validation methods differing in terms of relatedness. Including the progenitors in the calibration set reduced this difference in prediction accuracy to 0.03. When only genotypes from the G0 and G1 generations were used in the calibration set and genotypes from G2 were used in the validation set (progeny validation), prediction accuracies ranged from 0.70 to 0.85. CONCLUSIONS: This study suggests that the training of prediction models on parental populations can predict the genetic merit of the progeny with high accuracy: an encouraging result for the implementation of GS in the maritime pine breeding program.

Quantitative Proteomic and Phosphoproteomic Approaches for Deciphering the Signaling Pathway for Tension Wood Formation in Poplar.

Trees adjust their growth following forced changes in orientation to re-establish a vertical position. In angiosperms, this adjustment involves the differential regulation of vascular cambial activity between the lower (opposite wood) and upper (tension wood) sides of the leaning stem. We investigated the molecular mechanisms leading to the formation of differential wood types through a quantitative proteomic and phosphoproteomic analysis on poplar subjected to a gravitropic stimulus. We identified and quantified 675 phosphopeptides, corresponding to 468 phosphoproteins, and 3 763 nonphosphorylated peptides, corresponding to 1 155 proteins, in the differentiating xylem of straight-growing trees (control) and trees subjected to a gravitational stimulus during 8 weeks. About 1% of the peptides were specific to a wood type (straight, opposite, or tension wood). Proteins quantified in more than one type of wood were more numerous: a mixed linear model showed 389 phosphopeptides and 556 proteins to differ in abundance between tension wood and opposite wood. Twenty-one percent of the phosphoproteins identified here were described in their phosphorylated form for the first time. Our analyses revealed remarkable developmental molecular plasticity, with wood type-specific phosphorylation events, and highlighted the involvement of different proteins in the biosynthesis of cell wall components during the formation of the three types of wood.

The oak gene expression atlas: insights into Fagaceae genome evolution and the discovery of genes regulated during bud dormancy release.

BACKGROUND: Many northern-hemisphere forests are dominated by oaks. These species extend over diverse environmental conditions and are thus interesting models for studies of plant adaptation and speciation. The genomic toolbox is an important asset for exploring the functional variation associated with natural selection. RESULTS: The assembly of previously available and newly developed long and short sequence reads for two sympatric oak species, Quercus robur and Quercus petraea, generated a comprehensive catalog of transcripts for oak. The functional annotation of 91 k contigs demonstrated the presence of a large proportion of plant genes in this unigene set. Comparisons with SwissProt accessions and five plant gene models revealed orthologous relationships, making it possible to decipher the evolution of the oak genome. In particular, it was possible to align 9.5 thousand oak coding sequences with the equivalent sequences on peach chromosomes. Finally, RNA-seq data shed new light on the gene networks underlying vegetative bud dormancy release, a key stage in development allowing plants to adapt their phenology to the environment. CONCLUSION: In addition to providing a vast array of expressed genes, this study generated essential information about oak genome evolution and the regulation of genes associated with vegetative bud phenology, an important adaptive traits in trees. This resource contributes to the annotation of the oak genome sequence and will provide support for forward genetics approaches aiming to link genotypes with adaptive phenotypes.

Genetic architecture of carbon isotope composition and growth in Eucalyptus across multiple environments.

In the context of climate change, the water-use efficiency (WUE) of highly productive tree varieties, such as eucalypts, has become a major issue for breeding programmes. This study set out to dissect the genetic architecture of carbon isotope composition (δ(13) C), a proxy of WUE, across several environments. A family of Eucalyptus urophylla × E. grandis was planted in three trials and phenotyped for δ(13) C and growth traits. High-resolution genetic maps enabled us to target genomic regions underlying δ(13) C quantitative trait loci (QTLs) on the E. grandis genome. Of the 15 QTLs identified for δ(13) C, nine were stable across the environments and three displayed significant QTL-by-environment interaction, suggesting medium to high genetic determinism for this trait. Only one colocalization was found between growth and δ(13) C. Gene ontology (GO) term enrichment analysis suggested candidate genes related to foliar δ(13) C, including two involved in the regulation of stomatal movements. This study provides the first report of the genetic architecture of δ(13) C and its relation to growth in Eucalyptus. The low correlations found between the two traits at phenotypic and genetic levels suggest the possibility of improving the WUE of Eucalyptus varieties without having an impact on breeding for growth.

High-resolution genetic maps of Eucalyptus improve Eucalyptus grandis genome assembly.

Genetic maps are key tools in genetic research as they constitute the framework for many applications, such as quantitative trait locus analysis, and support the assembly of genome sequences. The resequencing of the two parents of a cross between Eucalyptus urophylla and Eucalyptus grandis was used to design a single nucleotide polymorphism (SNP) array of 6000 markers evenly distributed along the E. grandis genome. The genotyping of 1025 offspring enabled the construction of two high-resolution genetic maps containing 1832 and 1773 markers with an average marker interval of 0.45 and 0.5 cM for E. grandis and E. urophylla, respectively. The comparison between genetic maps and the reference genome highlighted 85% of collinear regions. A total of 43 noncollinear regions and 13 nonsynthetic regions were detected and corrected in the new genome assembly. This improved version contains 4943 scaffolds totalling 691.3 Mb of which 88.6% were captured by the 11 chromosomes. The mapping data were also used to investigate the effect of population size and number of markers on linkage mapping accuracy. This study provides the most reliable linkage maps for Eucalyptus and version 2.0 of the E. grandis genome.

Genomic selection for tolerance to aluminum toxicity in a synthetic population of upland rice.

Over half of the world's arable land is acidic, which constrains cereal production. In South America, different rice-growing regions (Cerrado in Brazil and Llanos in Colombia and Venezuela) are particularly affected due to high aluminum toxicity levels. For this reason, efforts have been made to breed for tolerance to aluminum toxicity using synthetic populations. The breeding program of CIAT-CIRAD is a good example of the use of recurrent selection to increase productivity for the Llanos in Colombia. In this study, we evaluated the performance of genomic prediction models to optimize the breeding scheme by hastening the development of an improved synthetic population and elite lines. We characterized 334 families at the S0:4 generation in two conditions. One condition was the control, managed with liming, while the other had high aluminum toxicity. Four traits were considered: days to flowering (FL), plant height (PH), grain yield (YLD), and zinc concentration in the polished grain (ZN). The population presented a high tolerance to aluminum toxicity, with more than 72% of the families showing a higher yield under aluminum conditions. The performance of the families under the aluminum toxicity condition was predicted using four different models: a single-environment model and three multi-environment models. The multi-environment models differed in the way they integrated genotype-by-environment interactions. The best predictive abilities were achieved using multi-environment models: 0.67 for FL, 0.60 for PH, 0.53 for YLD, and 0.65 for ZN. The gain of multi-environment over single-environment models ranged from 71% for YLD to 430% for FL. The selection of the best-performing families based on multi-trait indices, including the four traits mentioned above, facilitated the identification of suitable families for recombination. This information will be used to develop a new cycle of recurrent selection through genomic selection.

Plasticity of primary and secondary growth dynamics in Eucalyptus hybrids: a quantitative genetics and QTL mapping perspective.

BACKGROUND: The genetic basis of growth traits has been widely studied in forest trees. Quantitative trait locus (QTL) studies have highlighted the presence of both stable and unstable genomic regions accounting for biomass production with respect to tree age and genetic background, but results remain scarce regarding the interplay between QTLs and the environment. In this study, our main objective was to dissect the genetic architecture of the growth trajectory with emphasis on genotype x environment interaction by measuring primary and secondary growth covering intervals connected with environmental variations. RESULTS: Three different trials with the same family of Eucalyptus urophylla x E. grandis hybrids (with different genotypes) were planted in the Republic of Congo, corresponding to two QTL mapping experiments and one clonal test. Height and radial growths were monitored at regular intervals from the seedling stage to five years old. The correlation between growth increments and an aridity index revealed that growth before two years old (r = 0.5; 0.69) was more responsive to changes in water availability than late growth (r = 0.39; 0.42) for both height and circumference. We found a regular increase in heritability with time for cumulative growth for both height [0.06 - 0.33] and circumference [0.06 - 0.38]. Heritabilities for incremental growth were more heterogeneous over time even if ranges of variation were similar (height [0-0.31]; circumference [0.19 to 0.48]). Within the trials, QTL analysis revealed collocations between primary and secondary growth QTLs as well as between early growth increments and final growth QTLs. Between trials, few common QTLs were detected highlighting a strong environmental effect on the genetic architecture of growth, validated by significant QTL x E interactions. CONCLUSION: These results suggest that early growth responses to water availability determine the genetic architecture of total growth at the mature stage and highlight the importance of considering growth as a composite trait (such as yields for annual plants) for a better understanding of its genetic bases.

Realized Genetic Gain in Rice: Achievements from Breeding Programs.

Genetic improvement is crucial for ensuring food security globally. Indeed, plant breeding has contributed significantly to increasing the productivity of major crops, including rice, over the last century. Evaluating the efficiency of breeding strategies necessitates a quantification of this progress. One approach involves assessing the genetic gain achieved through breeding programs based on quantitative traits. This study aims to provide a theoretical understanding of genetic gain, summarize the major results of genetic gain studies in rice breeding, and suggest ways of improving breeding program strategies and future studies on genetic gain. To achieve this, we present the concept of genetic gain and the essential aspects of its estimation. We also provide an extensive literature review of genetic gain studies in rice (Oryza sativa L.) breeding programs to understand the advances made to date. We reviewed 29 studies conducted between 1999 and 2023, covering different regions, traits, periods, and estimation methods. The genetic gain for grain yield, in particular, showed significant variation, ranging from 1.5 to 167.6 kg/ha/year, with a mean value of 36.3 kg/ha/year. This translated into a rate of genetic gain for grain yield ranging from 0.1% to over 3.0%. The impact of multi-trait selection on grain yield was clarified by studies that reported genetic gains for other traits, such as plant height, days to flowering, and grain quality. These findings reveal that while breeding programs have achieved significant gains, further improvements are necessary to meet the growing demand for rice. We also highlight the limitations of these studies, which hinder accurate estimations of genetic gain. In conclusion, we offer suggestions for improving the estimation of genetic gain based on quantitative genetic principles and computer simulations to optimize rice breeding strategies.

Genomic selection for salinity tolerance in japonica rice.

Improving plant performance in salinity-prone conditions is a significant challenge in breeding programs. Genomic selection is currently integrated into many plant breeding programs as a tool for increasing selection intensity and precision for complex traits and for reducing breeding cycle length. A rice reference panel (RP) of 241 Oryza sativa L. japonica accessions genotyped with 20,255 SNPs grown in control and mild salinity stress conditions was evaluated at the vegetative stage for eight morphological traits and ion mass fractions (Na and K). Weak to strong genotype-by-condition interactions were found for the traits considered. Cross-validation showed that the predictive ability of genomic prediction methods ranged from 0.25 to 0.64 for multi-environment models with morphological traits and from 0.05 to 0.40 for indices of stress response and ion mass fractions. The performances of a breeding population (BP) comprising 393 japonica accessions were predicted with models trained on the RP. For validation of the predictive performances of the models, a subset of 41 accessions was selected from the BP and phenotyped under the same experimental conditions as the RP. The predictive abilities estimated on this subset ranged from 0.00 to 0.66 for the multi-environment models, depending on the traits, and were strongly correlated with the predictive abilities on cross-validation in the RP in salt condition (r = 0.69). We show here that genomic selection is efficient for predicting the salt stress tolerance of breeding lines. Genomic selection could improve the efficiency of rice breeding strategies for salinity-prone environments.

Optimization of Multi-Generation Multi-location Genomic Prediction Models for Recurrent Genomic Selection in an Upland Rice Population.

Genomic selection is a worthy breeding method to improve genetic gain in recurrent selection breeding schemes. The integration of multi-generation and multi-location information could significantly improve genomic prediction models in the context of shuttle breeding. The Cirad-CIAT upland rice breeding program applies recurrent genomic selection and seeks to optimize the scheme to increase genetic gain while reducing phenotyping efforts. We used a synthetic population (PCT27) of which S0 plants were all genotyped and advanced by selfing and bulk seed harvest to the S0:2, S0:3, and S0:4 generations. The PCT27 was then divided into two sets. The S0:2 and S0:3 progenies for PCT27A and the S0:4 progenies for PCT27B were phenotyped in two locations: Santa Rosa the target selection location, within the upland rice growing area, and Palmira, the surrogate location, far from the upland rice growing area but easier for experimentation. While the calibration used either one of the two sets phenotyped in one or two locations, the validation population was only the PCT27B phenotyped in Santa Rosa. Five scenarios of genomic prediction and 24 models were performed and compared. Training the prediction model with the PCT27B phenotyped in Santa Rosa resulted in predictive abilities ranging from 0.19 for grain zinc concentration to 0.30 for grain yield. Expanding the training set with the inclusion of the PCT27A resulted in greater predictive abilities for all traits but grain yield, with increases from 5% for plant height to 61% for grain zinc concentration. Models with the PCT27B phenotyped in two locations resulted in higher prediction accuracy when the models assumed no genotype-by-environment (G × E) interaction for flowering (0.38) and grain zinc concentration (0.27). For plant height, the model assuming a single G × E variance provided higher accuracy (0.28). The gain in predictive ability for grain yield was the greatest (0.25) when environment-specific variance deviation effect for G × E was considered. While the best scenario was specific to each trait, the results indicated that the gain in predictive ability provided by the multi-location and multi-generation calibration was low. Yet, this approach could lead to increased selection intensity, acceleration of the breeding cycle, and a sizable economic advantage for the program.

Multi-environment Genomic Selection in Rice Elite Breeding Lines.

BACKGROUND: Assessing the performance of elite lines in target environments is essential for breeding programs to select the most relevant genotypes. One of the main complexities in this task resides in accounting for the genotype by environment interactions. Genomic prediction models that integrate information from multi-environment trials and environmental covariates can be efficient tools in this context. The objective of this study was to assess the predictive ability of different genomic prediction models to optimize the use of multi-environment information. We used 111 elite breeding lines representing the diversity of the international rice research institute breeding program for irrigated ecosystems. The lines were evaluated for three traits (days to flowering, plant height, and grain yield) in 15 environments in Asia and Africa and genotyped with 882 SNP markers. We evaluated the efficiency of genomic prediction to predict untested environments using seven multi-environment models and three cross-validation scenarios. RESULTS: The elite lines were found to belong to the indica group and more specifically the indica-1B subgroup which gathered improved material originating from the Green Revolution. Phenotypic correlations between environments were high for days to flowering and plant height (33% and 54% of pairwise correlation greater than 0.5) but low for grain yield (lower than 0.2 in most cases). Clustering analyses based on environmental covariates separated Asia's and Africa's environments into different clusters or subclusters. The predictive abilities ranged from 0.06 to 0.79 for days to flowering, 0.25-0.88 for plant height, and - 0.29-0.62 for grain yield. We found that models integrating genotype-by-environment interaction effects did not perform significantly better than models integrating only main effects (genotypes and environment or environmental covariates). The different cross-validation scenarios showed that, in most cases, the use of all available environments gave better results than a subset. CONCLUSION: Multi-environment genomic prediction models with main effects were sufficient for accurate phenotypic prediction of elite lines in targeted environments. These results will help refine the testing strategy to update the genomic prediction models to improve predictive ability.

Genomic Prediction: Progress and Perspectives for Rice Improvement.

Genomic prediction can be a powerful tool to achieve greater rates of genetic gain for quantitative traits if thoroughly integrated into a breeding strategy. In rice as in other crops, the interest in genomic prediction is very strong with a number of studies addressing multiple aspects of its use, ranging from the more conceptual to the more practical. In this chapter, we review the literature on rice (Oryza sativa) and summarize important considerations for the integration of genomic prediction in breeding programs. The irrigated breeding program at the International Rice Research Institute is used as a concrete example on which we provide data and R scripts to reproduce the analysis but also to highlight practical challenges regarding the use of predictions. The adage "To someone with a hammer, everything looks like a nail" describes a common psychological pitfall that sometimes plagues the integration and application of new technologies to a discipline. We have designed this chapter to help rice breeders avoid that pitfall and appreciate the benefits and limitations of applying genomic prediction, as it is not always the best approach nor the first step to increasing the rate of genetic gain in every context.

Genetic Trends Estimation in IRRIs Rice Drought Breeding Program and Identification of High Yielding Drought-Tolerant Lines.

Estimating genetic trends using historical data is an important parameter to check the success of the breeding programs. The estimated genetic trends can act as a guideline to target the appropriate breeding strategies and optimize the breeding program for improved genetic gains. In this study, 17 years of historical data from IRRI's rice drought breeding program was used to estimate the genetic trends and assess the breeding program's success. We also identified top-performing lines based on grain yield breeding values as an elite panel for implementing future population improvement-based breeding schemes. A two-stage approach of pedigree-based mixed model analysis was used to analyze the data and extract the breeding values and estimate the genetic trends for grain yield under non-stress, drought, and in combined data of non-stress and drought. Lower grain yield values were observed in all the drought trials. Heritability for grain yield estimates ranged between 0.20 and 0.94 under the drought trials and 0.43-0.83 under non-stress trials. Under non-stress conditions, the genetic gain of 0.21% (10.22 kg/ha/year) for genotypes and 0.17% (7.90 kg/ha/year) for checks was observed. The genetic trend under drought conditions exhibited a positive trend with the genetic gain of 0.13% (2.29 kg/ha/year) for genotypes and 0.55% (9.52 kg/ha/year) for checks. For combined analysis showed a genetic gain of 0.27% (8.32 kg/ha/year) for genotypes and 0.60% (13.69 kg/ha/year) for checks was observed. For elite panel selection, 200 promising lines were selected based on higher breeding values for grain yield and prediction accuracy of > 0.40. The breeding values of the 200 genotypes formulating the core panel ranged between 2366.17 and 4622.59 (kg/ha). A positive genetic rate was observed under all the three conditions; however, the rate of increase was lower than the required rate of 1.5% genetic gain. We propose a recurrent selection breeding strategy within the elite population with the integration of modern tools and technologies to boost the genetic gains in IRRI's drought breeding program. The elite breeding panel identified in this study forms an easily available and highly enriched genetic resource for future recurrent selection programs to boost the genetic gains.

Impact of early genomic prediction for recurrent selection in an upland rice synthetic population.

Population breeding through recurrent selection is based on the repetition of evaluation and recombination among best-selected individuals. In this type of breeding strategy, early evaluation of selection candidates combined with genomic prediction could substantially shorten the breeding cycle length, thus increasing the rate of genetic gain. The objective of this study was to optimize early genomic prediction in an upland rice (Oryza sativa L.) synthetic population improved through recurrent selection via shuttle breeding in two sites. To this end, we used genomic prediction on 334 S0 genotypes evaluated with early generation progeny testing (S0:2 and S0:3) across two sites. Four traits were measured (plant height, days to flowering, grain yield, and grain zinc concentration) and the predictive ability was assessed for the target site. For days to flowering and plant height, which correlate well among sites (0.51-0.62), an increase of up to 0.4 in predictive ability was observed when the model was trained using the two sites. For grain zinc concentration, adding the phenotype of the predicted lines in the nontarget site to the model improved the predictive ability (0.51 with two-site and 0.31 with single-site model), whereas for grain yield the gain was less (0.42 with two-site and 0.35 with single-site calibration). Through these results, we found a good opportunity to optimize the genomic recurrent selection scheme and maximize the use of resources by performing early progeny testing in two sites for traits with best expression and/or relevance in each specific environment.

Identification of an Elite Core Panel as a Key Breeding Resource to Accelerate the Rate of Genetic Improvement for Irrigated Rice.

Rice genetic improvement is a key component of achieving and maintaining food security in Asia and Africa in the face of growing populations and climate change. In this effort, the International Rice Research Institute (IRRI) continues to play a critical role in creating and disseminating rice varieties with higher productivity. Due to increasing demand for rice, especially in Africa, there is a strong need to accelerate the rate of genetic improvement for grain yield. In an effort to identify and characterize the elite breeding pool of IRRI's irrigated rice breeding program, we analyzed 102 historical yield trials conducted in the Philippines during the period 2012-2016 and representing 15,286 breeding lines (including released varieties). A mixed model approach based on the pedigree relationship matrix was used to estimate breeding values for grain yield, which ranged from 2.12 to 6.27 t·ha-1. The rate of genetic gain for grain yield was estimated at 8.75 kg·ha-1 year-1 (0.23%) for crosses made in the period from 1964 to 2014. Reducing the data to only IRRI released varieties, the rate doubled to 17.36 kg·ha-1 year-1 (0.46%). Regressed against breeding cycle the rate of gain for grain yield was 185 kg·ha-1 cycle-1 (4.95%). We selected 72 top performing lines based on breeding values for grain yield to create an elite core panel (ECP) representing the genetic diversity in the breeding program with the highest heritable yield values from which new products can be derived. The ECP closely aligns with the indica 1B sub-group of Oryza sativa that includes most modern varieties for irrigated systems. Agronomic performance of the ECP under multiple environments in Asia and Africa confirmed its high yield potential. We found that the rate of genetic gain for grain yield found in this study was limited primarily by long cycle times and the direct introduction of non-improved material into the elite pool. Consequently, the current breeding scheme for irrigated rice at IRRI is based on rapid recurrent selection among highly elite lines. In this context, the ECP constitutes an important resource for IRRI and NAREs breeders to carefully characterize and manage that elite diversity.

Genome-Wide Association Study for Yield and Yield Related Traits under Reproductive Stage Drought in a Diverse indica-aus Rice Panel.

BACKGROUND: Reproductive-stage drought stress is a major impediment to rice production in rainfed areas. Conventional and marker-assisted breeding strategies for developing drought-tolerant rice varieties are being optimized by mining and exploiting adaptive traits, genetic diversity; identifying the alleles, and understanding their interactions with genetic backgrounds for their increased contribution to drought tolerance. Field experiments were conducted in this study to identify marker-trait associations (MTAs) involved in response to yield under reproductive-stage (RS) drought. A diverse set of 280 indica-aus accessions was phenotyped for ten agronomic traits including yield and yield-related traits under normal irrigated condition and under two managed reproductive-stage drought environments. The accessions were genotyped with 215,250 single nucleotide polymorphism markers. RESULTS: The study identified a total of 219 significant MTAs for 10 traits and candidate gene analysis within a 200 kb window centred from GWAS identified SNP peaks detected these MTAs within/ in close proximity to 38 genes, 4 earlier reported major grain yield QTLs and 6 novel QTLs for 7 traits out of the 10. The significant MTAs were mainly located on chromosomes 1, 2, 5, 6, 9, 11 and 12 and the percent phenotypic variance captured for these traits ranged from 5 to 88%. The significant positive correlation of grain yield with yield-related and other agronomic traits except for flowering time, observed under different environments point towards their contribution in improving rice yield under drought. Seven promising accessions were identified for use in future genomics-assisted breeding programs targeting grain yield improvement under drought. CONCLUSION: These results provide a promising insight into the complex genetic architecture of grain yield under reproductive-stage drought in different environments. Validation of major genomic regions reported in the study will enable their effectiveness to develop drought-tolerant varieties following marker-assisted selection as well as to identify genes and understanding the associated physiological mechanisms.

Fine Scale Genomic Signals of Admixture and Alien Introgression among Asian Rice Landraces.

Modern rice cultivars are adapted to a range of environmental conditions and human preferences. At the root of this diversity is a marked genetic structure, owing to multiple foundation events. Admixture and recurrent introgression from wild sources have played upon this base to produce the myriad adaptations existing today. Genome-wide studies bring support to this idea, but understanding the history and nature of particular genetic adaptations requires the identification of specific patterns of genetic exchange. In this study, we explore the patterns of haplotype similarity along the genomes of a subset of rice cultivars available in the 3,000 Rice Genomes data set. We begin by establishing a custom method of classification based on a combination of dimensionality reduction and kernel density estimation. Through simulations, the behavior of this classifier is studied under scenarios of varying genetic divergence, admixture, and alien introgression. Finally, the method is applied to local haplotypes along the genome of a Core set of Asian Landraces. Taking the Japonica, Indica, and cAus groups as references, we find evidence of reciprocal introgressions covering 2.6% of reference genomes on average. Structured signals of introgression among reference accessions are discussed. We extend the analysis to elucidate the genetic structure of the group circum-Basmati: we delimit regions of Japonica, cAus, and Indica origin, as well as regions outlier to these groups (13% on average). Finally, the approach used highlights regions of partial to complete loss of structure that can be attributed to selective pressures during domestication.