Identification of Fusarium head blight sources of resistance and associated QTLs in historical and modern Canadian spring wheat.

Fusarium head blight (FHB) is one the most globally destructive fungal diseases in wheat and other small grains, causing a reduction in grain yield by 10-70%. The present study was conducted in a panel of historical and modern Canadian spring wheat (Triticum aestivum L.) varieties and lines to identify new sources of FHB resistance and map associated quantitative trait loci (QTLs). We evaluated 249 varieties and lines for reaction to disease incidence, severity, and visual rating index (VRI) in seven environments by artificially spraying a mixture of four Fusarium graminearum isolates. A subset of 198 them were genotyped with the Wheat 90K iSelect single nucleotide polymorphisms (SNPs) array. Genome-wide association mapping performed on the overall best linear unbiased estimators (BLUE) computed from all seven environments and the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map of 26,449 polymorphic SNPs out of the 90K identified sixteen FHB resistance QTLs that individually accounted for 5.7-10.2% of the phenotypic variance. The positions of two of the FHB resistance QTLs overlapped with plant height and flowering time QTLs. Four of the QTLs (QFhb.dms-3B.1, QFhb.dms-5A.5, QFhb.dms-5A.7, and QFhb.dms-6A.4) were simultaneously associated with disease incidence, severity, and VRI, which accounted for 27.0-33.2% of the total phenotypic variance in the combined environments. Three of the QTLs (QFhb.dms-2A.2, QFhb.dms-2D.2, and QFhb.dms-5B.8) were associated with both incidence and VRI and accounted for 20.5-22.1% of the total phenotypic variance. In comparison with the VRI of the checks, we identified four highly resistant and thirty-three moderately resistant lines and varieties. The new FHB sources of resistance and the physical map of the associated QTLs would provide wheat breeders valuable information towards their efforts in developing improved varieties in western Canada.

Identification of Disease Resistance Parents and Genome-Wide Association Mapping of Resistance in Spring Wheat.

The likelihood of success in developing modern cultivars depend on multiple factors, including the identification of suitable parents to initiate new crosses, and characterizations of genomic regions associated with target traits. The objectives of the present study were to (a) determine the best economic weights of four major wheat diseases (leaf spot, common bunt, leaf rust, and stripe rust) and grain yield for multi-trait restrictive linear phenotypic selection index (RLPSI), (b) select the top 10% cultivars and lines (hereafter referred as genotypes) with better resistance to combinations of the four diseases and acceptable grain yield as potential parents, and (c) map genomic regions associated with resistance to each disease using genome-wide association study (GWAS). A diversity panel of 196 spring wheat genotypes was evaluated for their reaction to stripe rust at eight environments, leaf rust at four environments, leaf spot at three environments, common bunt at two environments, and grain yield at five environments. The panel was genotyped with the Wheat 90K SNP array and a few KASP SNPs of which we used 23,342 markers for statistical analyses. The RLPSI analysis performed by restricting the expected genetic gain for yield displayed significant (p < 0.05) differences among the 3125 economic weights. Using the best four economic weights, a subset of 22 of the 196 genotypes were selected as potential parents with resistance to the four diseases and acceptable grain yield. GWAS identified 37 genomic regions, which included 12 for common bunt, 13 for leaf rust, 5 for stripe rust, and 7 for leaf spot. Each genomic region explained from 6.6 to 16.9% and together accounted for 39.4% of the stripe rust, 49.1% of the leaf spot, 94.0% of the leaf rust, and 97.9% of the common bunt phenotypic variance combined across all environments. Results from this study provide valuable information for wheat breeders selecting parental combinations for new crosses to develop improved germplasm with enhanced resistance to the four diseases as well as the physical positions of genomic regions that confer resistance, which facilitates direct comparisons for independent mapping studies in the future.

Identification of Spring Wheat with Superior Agronomic Performance under Contrasting Nitrogen Managements Using Linear Phenotypic Selection Indices.

Both the Linear Phenotypic Selection Index (LPSI) and the Restrictive Linear Phenotypic Selection Index (RLPSI) have been widely used to select parents and progenies, but the effect of economic weights on the selection parameters (the expected genetic gain, response to selection, and the correlation between the indices and genetic merits) have not been investigated in detail. Here, we (i) assessed combinations of 2304 economic weights using four traits (maturity, plant height, grain yield and grain protein content) recorded under four organically (low nitrogen) and five conventionally (high nitrogen) managed environments, (ii) compared single-trait and multi-trait selection indices (LPSI vs. RLPSI by imposing restrictions to the expected genetic gain of either yield or grain protein content), and (iii) selected a subset of about 10% spring wheat cultivars that performed very well under organic and/or conventional management systems. The multi-trait selection indices, with and without imposing restrictions, were superior to single trait selection. However, the selection parameters differed quite a lot depending on the economic weights, which suggests the need for optimizing the weights. Twenty-two of the 196 cultivars that showed superior performance under organic and/or conventional management systems were consistently selected using all five of the selected economic weights, and at least two of the selection scenarios. The selected cultivars belonged to the Canada Western Red Spring (16 cultivars), the Canada Northern Hard Red (3), and the Canada Prairie Spring Red (3), and required 83-93 days to maturity, were 72-100 cm tall, and produced from 4.0 to 6.2 t ha-1 grain yield with 14.6-17.7% GPC. The selected cultivars would be highly useful, not only as potential trait donors for breeding under an organic management system, but also for other studies, including nitrogen use efficiency.

Genomic Prediction Accuracy of Stripe Rust in Six Spring Wheat Populations by Modeling Genotype by Environment Interaction.

Some previous studies have assessed the predictive ability of genome-wide selection on stripe (yellow) rust resistance in wheat, but the effect of genotype by environment interaction (GEI) in prediction accuracies has not been well studied in diverse genetic backgrounds. Here, we compared the predictive ability of a model based on phenotypic data only (M1), the main effect of phenotype and molecular markers (M2), and a model that incorporated GEI (M3) using three cross-validations (CV1, CV2, and CV0) scenarios of interest to breeders in six spring wheat populations. Each population was evaluated at three to eight field nurseries and genotyped with either the DArTseq technology or the wheat 90K single nucleotide polymorphism arrays, of which a subset of 1,058- 23,795 polymorphic markers were used for the analyses. In the CV1 scenario, the mean prediction accuracies of the M1, M2, and M3 models across the six populations varied from -0.11 to -0.07, from 0.22 to 0.49, and from 0.19 to 0.48, respectively. Mean accuracies obtained using the M3 model in the CV1 scenario were significantly greater than the M2 model in two populations, the same in three populations, and smaller in one population. In both the CV2 and CV0 scenarios, the mean prediction accuracies of the three models varied from 0.53 to 0.84 and were not significantly different in all populations, except the Attila/CDC Go in the CV2, where the M3 model gave greater accuracy than both the M1 and M2 models. Overall, the M3 model increased prediction accuracies in some populations by up to 12.4% and decreased accuracy in others by up to 17.4%, demonstrating inconsistent results among genetic backgrounds that require considering each population separately. This is the first comprehensive genome-wide prediction study that investigated details of the effect of GEI on stripe rust resistance across diverse spring wheat populations.

Comparison of single-trait and multi-trait genomic predictions on agronomic and disease resistance traits in spring wheat.

KEY MESSAGE: This study performed comprehensive analyses on the predictive abilities of single-trait and two multi-trait models in three populations. Our results demonstrated the superiority of multi-traits over single-trait models across seven agronomic and four to seven disease resistance traits of different genetic architecture. The predictive ability of multi-trait and single-trait prediction models has not been investigated on diverse traits evaluated under organic and conventional management systems. Here, we compared the predictive abilities of 25% of a testing set that has not been evaluated for a single trait (ST), not evaluated for multi-traits (MT1), and evaluated for some traits but not others (MT2) in three spring wheat populations genotyped either with the wheat 90K single nucleotide polymorphisms array or DArTseq. Analyses were performed on seven agronomic traits evaluated under conventional and organic management systems, four to seven disease resistance traits, and all agronomic and disease resistance traits simultaneously. The average prediction accuracies of the ST, MT1, and MT2 models varied from 0.03 to 0.78 (mean 0.41), from 0.05 to 0.82 (mean 0.47), and from 0.05 to 0.92 (mean 0.67), respectively. The predictive ability of the MT2 model was significantly greater than the ST model in all traits and populations except common bunt with the MT1 model being intermediate between them. The MT2 model increased prediction accuracies over the ST and MT1 models in all traits by 9.0-82.4% (mean 37.3%) and 2.9-82.5% (mean 25.7%), respectively, except common bunt that showed up to 7.7% smaller accuracies in two populations. A joint analysis of all agronomic and disease resistance traits further improved accuracies within the MT1 and MT2 models on average by 21.4% and 17.4%, respectively, as compared to either the agronomic or disease resistance traits, demonstrating the high potential of the multi-traits models in improving prediction accuracies.

Genomic Predictions for Common Bunt, FHB, Stripe Rust, Leaf Rust, and Leaf Spotting Resistance in Spring Wheat.

Some studies have investigated the potential of genomic selection (GS) on stripe rust, leaf rust, Fusarium head blight (FHB), and leaf spot in wheat, but none of them have assessed the effect of the reaction norm model that incorporated GE interactions. In addition, the prediction accuracy on common bunt has not previously been studied. Here, we investigated within-population prediction accuracies using the baseline M1 model and two reaction norm models (M2 and M3) with three random cross-validation (CV1, CV2, and CV0) schemes. Three Canadian spring wheat populations were evaluated in up to eight field environments and genotyped with 3158, 5732, and 23,795 polymorphic markers. The M3 model that incorporated GE interactions reduced residual variance by an average of 10.2% as compared with the main effect M2 model and increased prediction accuracies on average by 2-6%. In some traits, the M3 model increased prediction accuracies up to 54% as compared with the M2 model. The average prediction accuracies of the M3 model with CV1, CV2, and CV0 schemes varied from 0.02 to 0.48, from 0.25 to 0.84, and from 0.14 to 0.87, respectively. In both CV2 and CV0 schemes, stripe rust in all three populations, common bunt and leaf rust in two populations, as well as FHB severity, FHB index, and leaf spot in one population had high to very high (0.54-0.87) prediction accuracies. This is the first comprehensive genomic selection study on five major diseases in spring wheat.

Genome-based prediction of agronomic traits in spring wheat under conventional and organic management systems.

KEY MESSAGE: Using phenotype data of three spring wheat populations evaluated at 6-15 environments under two management systems, we found moderate to very high prediction accuracies across seven traits. The phenotype data collected under an organic management system effectively predicted the performance of lines in the conventional management and vice versa. There is growing interest in developing wheat cultivars specifically for organic agriculture, but we are not aware of the effect of organic management on the predictive ability of genomic selection (GS). Here, we evaluated within populations prediction accuracies of four GS models, four combinations of training and testing sets, three reaction norm models, and three random cross-validations (CV) schemes in three populations phenotyped under organic and conventional management systems. Our study was based on a total of 578 recombinant inbred lines and varieties from three spring wheat populations, which were evaluated for seven traits at 3-9 conventionally and 3-6 organically managed field environments and genotyped either with the wheat 90 K SNP array or DArTseq. We predicted the management systems (CV0M) or environments (CV0), a subset of lines that have been evaluated in either management (CV2M) or some environments (CV2), and the performance of newly developed lines in either management (CV1M) or environments (CV1). The average prediction accuracies of the model that incorporated genotype × environment interactions with CV0 and CV2 schemes varied from 0.69 to 0.97. In the CV1 and CV1M schemes, prediction accuracies ranged from - 0.12 to 0.77 depending on the reaction norm models, the traits, and populations. In most cases, grain protein showed the highest prediction accuracies. The phenotype data collected under the organic management effectively predicted the performance of lines under conventional management and vice versa. This is the first comprehensive GS study that investigated the effect of the organic management system in wheat.

Genetic diversity and selective sweeps in historical and modern Canadian spring wheat cultivars using the 90K SNP array.

Previous molecular characterization studies conducted in Canadian wheat cultivars shed some light on the impact of plant breeding on genetic diversity, but the number of varieties and markers used was small. Here, we used 28,798 markers of the wheat 90K single nucleotide polymorphisms to (a) assess the extent of genetic diversity, relationship, population structure, and divergence among 174 historical and modern Canadian spring wheat varieties registered from 1905 to 2018 and 22 unregistered lines (hereinafter referred to as cultivars), and (b) identify genomic regions that had undergone selection. About 91% of the pairs of cultivars differed by 20-40% of the scored alleles, but only 7% of the pairs had kinship coefficients of < 0.250, suggesting the presence of a high proportion of redundancy in allelic composition. Although the 196 cultivars represented eight wheat classes, our results from phylogenetic, principal component, and the model-based population structure analyses revealed three groups, with no clear structure among most wheat classes, breeding programs, and breeding periods. FST statistics computed among different categorical variables showed little genetic differentiation (< 0.05) among breeding periods and breeding programs, but a diverse level of genetic differentiation among wheat classes and predicted groups. Diversity indices were the highest and lowest among cultivars registered from 1970 to 1980 and from 2011 to 2018, respectively. Using two outlier detection methods, we identified from 524 to 2314 SNPs and 41 selective sweeps of which some are close to genes with known phenotype, including plant height, photoperiodism, vernalization, gluten strength, and disease resistance.

Physical mapping of QTL associated with agronomic and end-use quality traits in spring wheat under conventional and organic management systems.

KEY MESSAGE: Using phenotypic data of four biparental spring wheat populations evaluated at multiple environments under two management systems, we discovered 152 QTL and 22 QTL hotspots, of which two QTL accounted for up to 37% and 58% of the phenotypic variance, consistently detected in all environments, and fell within genomic regions harboring known genes. Identification of the physical positions of quantitative trait loci (QTL) would be highly useful for developing functional markers and comparing QTL results across multiple independent studies. The objectives of the present study were to map and characterize QTL associated with nine agronomic and end-use quality traits (tillering ability, plant height, lodging, grain yield, grain protein content, thousand kernel weight, test weight, sedimentation volume, and falling number) in hard red spring wheat recombinant inbred lines (RILs) using the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map. We evaluated a total of 698 RILs from four populations derived from crosses involving seven parents at 3-8 conventionally (high N) and organically (low N) managed field environments. Using the phenotypic data combined across all environments per management, and the physical map between 1058 and 6526 markers per population, we identified 152 QTL associated with the nine traits, of which 29 had moderate and 2 with major effects. Forty-nine of the 152 QTL mapped across 22 QTL hotspot regions with each region coincident to 2-6 traits. Some of the QTL hotspots were physically located close to known genes. QSv.dms-1A and QPht.dms-4B.1 individually explained up to 37% and 58% of the variation in sedimentation volume and plant height, respectively, and had very large LOD scores that varied from 19.0 to 35.7 and from 16.7 to 55.9, respectively. We consistently detected both QTL in the combined and all individual environments, laying solid ground for further characterization and possibly for cloning.

Analysis of unigenes involved in lateral root development in Bupleurum chinense and B. scorzonerifolium.

MAIN CONCLUSION: We identified IAA13 negatively associated with lateral root number by comparing the differential expressed genes between Bupleurum chinense and B. scorzonerifolium. Dried roots of the genus Bupleurum L. are used as a herbal medicine for diseases in Asia. Bupleurum chinense has a greater number of lateral roots than B. scorzonerifolium, but the genetic mechanisms for such differences are largely unknown. We (a) compared the transcriptome profiles of the two species and (b) identified a subset of candidate genes involved in auxin signal transduction and explored their functions in lateral root development. By isoform sequencing (Iso-Seq) analyses of the whole plant, more unigenes were found in B. scorzonerifolium (118,868) than in B. chinense (93,485). Given the overarching role of indole-3-acetic acid (IAA) as one of the major regulators of lateral root development, we identified 539 unigenes associated with auxin signal transduction. Fourteen and 44 unigenes in the pathway were differentially expressed in B. chinense and B. scorzonerifolium, respectively, and 3 unigenes (LAX2, LAX4, and IAA13) were expressed in both species. The number of lateral root primordia increased after exogenous auxin application at 8 h and 12 h in B. scorzonerifolium and B. chinense, respectively. Since overexpression of IAA13 in Arabidopsis reduced the number of lateral roots, we hypothesized that IAA13 is involved in the reduction of the number of lateral roots in B. scorzonerifolium.

Physical Mapping of QTL in Four Spring Wheat Populations under Conventional and Organic Management Systems. I. Earliness.

In previous studies, we reported quantitative trait loci (QTL) associated with the heading, flowering, and maturity time in four hard red spring wheat recombinant inbred line (RIL) populations but the results are scattered in population-specific genetic maps, which is challenging to exploit efficiently in breeding. Here, we mapped and characterized QTL associated with these three earliness traits using the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map. Our data consisted of (i) 6526 single nucleotide polymorphisms (SNPs) and two traits evaluated at five conventionally managed environments in the 'Cutler' × 'AC Barrie' population; (ii) 3158 SNPs and two traits evaluated across three organic and seven conventional managements in the 'Attila' × 'CDC Go' population; (iii) 5731 SilicoDArT and SNP markers and the three traits evaluated at four conventional and organic management systems in the 'Peace' × 'Carberry' population; and (iv) 1058 SNPs and two traits evaluated across two conventionally and organically managed environments in the 'Peace' × 'CDC Stanley' population. Using composite interval mapping, the phenotypic data across all environments, and the IWGSC RefSeq v2.0 physical maps, we identified a total of 44 QTL associated with days to heading (11), flowering (10), and maturity (23). Fifteen of the 44 QTL were common to both conventional and organic management systems, and the remaining QTL were specific to either the conventional (21) or organic (8) management systems. Some QTL harbor known genes, including the Vrn-A1, Vrn-B1, Rht-A1, and Rht-B1 that regulate photoperiodism, flowering time, and plant height in wheat, which lays a solid basis for cloning and further characterization.

Comparisons of sampling methods for assessing intra- and inter-accession genetic diversity in three rice species using genotyping by sequencing.

To minimize the cost of sample preparation and genotyping, most genebank genomics studies in self-pollinating species are conducted on a single individual to represent an accession, which may be heterogeneous with larger than expected intra-accession genetic variation. Here, we compared various population genetics parameters among six DNA (leaf) sampling methods on 90 accessions representing a wild species (O. barthii), cultivated and landraces (O. glaberrima, O. sativa), and improved varieties derived through interspecific hybridizations. A total of 1,527 DNA samples were genotyped with 46,818 polymorphic single nucleotide polymorphisms (SNPs) using DArTseq. Various statistical analyses were performed on eleven datasets corresponding to 5 plants per accession individually and in a bulk (two sets), 10 plants individually and in a bulk (two sets), all 15 plants individually (one set), and a randomly sampled individual repeated six times (six sets). Overall, we arrived at broadly similar conclusions across 11 datasets in terms of SNP polymorphism, heterozygosity/heterogeneity, diversity indices, concordance among genetic dissimilarity matrices, population structure, and genetic differentiation; there were, however, a few discrepancies between some pairs of datasets. Detailed results of each sampling method, the concordance in their outputs, and the technical and cost implications of each method were discussed.

Molecular diversity and selective sweeps in maize inbred lines adapted to African highlands.

Little is known on maize germplasm adapted to the African highland agro-ecologies. In this study, we analyzed high-density genotyping by sequencing (GBS) data of 298 African highland adapted maize inbred lines to (i) assess the extent of genetic purity, genetic relatedness, and population structure, and (ii) identify genomic regions that have undergone selection (selective sweeps) in response to adaptation to highland environments. Nearly 91% of the pairs of inbred lines differed by 30-36% of the scored alleles, but only 32% of the pairs of the inbred lines had relative kinship coefficient <0.050, which suggests the presence of substantial redundancy in allelic composition that may be due to repeated use of fewer genetic backgrounds (source germplasm) during line development. Results from different genetic relatedness and population structure analyses revealed three different groups, which generally agrees with pedigree information and breeding history, but less so by heterotic groups and endosperm modification. We identified 944 single nucleotide polymorphic (SNP) markers that fell within 22 selective sweeps that harbored 265 protein-coding candidate genes of which some of the candidate genes had known functions. Details of the candidate genes with known functions and differences in nucleotide diversity among groups predicted based on multivariate methods have been discussed.

Comparisons of molecular diversity indices, selective sweeps and population structure of African rice with its wild progenitor and Asian rice.

KEY MESSAGE: The extent of molecular diversity parameters across three rice species was compared using large germplasm collection genotyped with genomewide SNPs and SNPs that fell within selective sweep regions. Previous studies conducted on limited number of accessions have reported very low genetic variation in African rice (Oryza glaberrima Steud.) as compared to its wild progenitor (O. barthii A. Chev.) and to Asian rice (O. sativa L.). Here, we characterized a large collection of African rice and compared its molecular diversity indices and population structure with the two other species using genomewide single nucleotide polymorphisms (SNPs) and SNPs that mapped within selective sweeps. A total of 3245 samples representing African rice (2358), Asian rice (772) and O. barthii (115) were genotyped with 26,073 physically mapped SNPs. Using all SNPs, the level of marker polymorphism, average genetic distance and nucleotide diversity in African rice accounted for 59.1%, 63.2% and 37.1% of that of O. barthii, respectively. SNP polymorphism and overall nucleotide diversity of the African rice accounted for 20.1-32.1 and 16.3-37.3% of that of the Asian rice, respectively. We identified 780 SNPs that fell within 37 candidate selective sweeps in African rice, which were distributed across all 12 rice chromosomes. Nucleotide diversity of the African rice estimated from the 780 SNPs was 8.3 × 10-4, which is not only 20-fold smaller than the value estimated from all genomewide SNPs (π = 1.6 × 10-2), but also accounted for just 4.1%, 0.9% and 2.1% of that of O. barthii, lowland Asian rice and upland Asian rice, respectively. The genotype data generated for a large collection of rice accessions conserved at the AfricaRice genebank will be highly useful for the global rice community and promote germplasm use.

Development of species diagnostic SNP markers for quality control genotyping in four rice (Oryza L.) species.

Species misclassification (misidentification) and handling errors have been frequently reported in various plant species conserved at diverse gene banks, which could restrict use of germplasm for correct purpose. The objectives of the present study were to (i) determine the extent of genotyping error (reproducibility) on DArTseq-based single-nucleotide polymorphisms (SNPs); (ii) determine the proportion of misclassified accessions across 3134 samples representing three African rice species complex (Oryza glaberrima, O. barthii, and O. longistaminata) and an Asian rice (O. sativa), which are conserved at the AfricaRice gene bank; and (iii) develop species- and sub-species (ecotype)-specific diagnostic SNP markers for rapid and low-cost quality control (QC) analysis. Genotyping error estimated from 15 accessions, each replicated from 2 to 16 times, varied from 0.2 to 3.1%, with an overall average of 0.8%. Using a total of 3134 accessions genotyped with 31,739 SNPs, the proportion of misclassified samples was 3.1% (97 of the 3134 accessions). Excluding the 97 misclassified accessions, we identified a total of 332 diagnostic SNPs that clearly discriminated the three indigenous African species complex from Asian rice (156 SNPs), O. longistaminata accessions from both O. barthii and O. glaberrima (131 SNPs), and O. sativa spp. indica from O. sativa spp. japonica (45 SNPs). Using chromosomal position, minor allele frequency, and polymorphic information content as selection criteria, we recommended a subset of 24 to 36 of the 332 diagnostic SNPs for routine QC genotyping, which would be highly useful in determining the genetic identity of each species and correct human errors during routine gene bank operations.

Discovery and validation of genomic regions associated with resistance to maize lethal necrosis in four biparental populations.

In sub-Saharan Africa, maize is the key determinant of food security for smallholder farmers. The sudden outbreak of maize lethal necrosis (MLN) disease is seriously threatening the maize production in the region. Understanding the genetic basis of MLN resistance is crucial. In this study, we used four biparental populations applied linkage mapping and joint linkage mapping approaches to identify and validate the MLN resistance-associated genomic regions. All populations were genotyped with low to high density markers and phenotyped in multiple environments against MLN under artificial inoculation. Phenotypic variation for MLN resistance was significant and heritability was moderate to high in all four populations for both early and late stages of disease infection. Linkage mapping revealed three major quantitative trait loci (QTL) on chromosomes 3, 6, and 9 that were consistently detected in at least two of the four populations. Phenotypic variance explained by a single QTL in each population ranged from 3.9% in population 1 to 43.8% in population 2. Joint linkage association mapping across three populations with three biometric models together revealed 16 and 10 main effect QTL for MLN-early and MLN-late, respectively. The QTL identified on chromosomes 3, 5, 6, and 9 were consistent with the QTL identified by linkage mapping. Ridge regression best linear unbiased prediction with five-fold cross-validation revealed high accuracy for prediction across populations for both MLN-early and MLN-late. Overall, the study discovered and validated the presence of major effect QTL on chromosomes 3, 6, and 9 which can be potential candidates for marker-assisted breeding to improve the MLN resistance.

Assessment of Genetic Variation and Population Structure of Diverse Rice Genotypes Adapted to Lowland and Upland Ecologies in Africa Using SNPs.

Using interspecific crosses involving Oryza glaberrima Steud. as donor and O. sativa L. as recurrent parents, rice breeders at the Africa Rice Center developed several 'New Rice for Africa (NERICA)' improved varieties. A smaller number of interspecific and intraspecific varieties have also been released as 'Advanced Rice for Africa (ARICA)'. The objective of the present study was to investigate the genetic variation, relatedness, and population structure of 330 widely used rice genotypes in Africa using DArTseq-based single nucleotide polymorphisms (SNPs). A sample of 11 ARICAs, 85 NERICAs, 62 O. sativa spp. japonica, and 172 O. sativa spp. indica genotypes were genotyped with 27,560 SNPs using diversity array technology (DArT)-based sequencing (DArTseq) platform. Nearly 66% of the SNPs were polymorphic, of which 15,020 SNPs were mapped to the 12 rice chromosomes. Genetic distance between pairs of genotypes that belong to indica, japonica, ARICA, and NERICA varied from 0.016 to 0.623, from 0.020 to 0.692, from 0.075 to 0.763, and from 0.014 to 0.644, respectively. The proportion of pairs of genotypes with genetic distance > 0.400 was the largest within NERICAs (35.1% of the pairs) followed by ARICAs (18.2%), japonica (17.4%), and indica (5.6%). We found one pair of japonica, 11 pairs of indica, and 35 pairs of NERICA genotypes differing by <2% of the total scored alleles, which was due to 26 pairs of genotypes with identical pedigrees. Cluster analysis, principal component analysis, and the model-based population structure analysis all revealed two distinct groups corresponding to the lowland (primarily indica and lowland NERICAs) and upland (japonica and upland NERICAs) growing ecologies. Most of the interspecific lowland NERICAs formed a sub-group, likely caused by differences in the O. glaberrima genome as compared with the indica genotypes. Analysis of molecular variance revealed very great genetic differentiation (FST = 0.688) between the lowland and upland ecologies, and 31.2% of variation attributable to differences within cluster groups. About 8% (1,197 of 15,020) of the 15,020 SNPs were significantly (P < 0.05) different between the lowland and upland ecologies and formed contrasting haplotypes that could clearly discriminate lowland from upland genotypes. This is the first study using high density markers that characterized NERICA and ARICA varieties in comparison with indica and japonica varieties widely used in Africa, which could aid rice breeders on parent selection for developing new improved rice germplasm.

Genetic Variation and Population Structure of Oryza glaberrima and Development of a Mini-Core Collection Using DArTseq.

The sequence variation present in accessions conserved in genebanks can best be used in plant improvement when it is properly characterized and published. Using low cost and high density single nucleotide polymorphism (SNP) assays, the genetic diversity, population structure, and relatedness between pairs of accessions can be quickly assessed. This information is relevant for different purposes, including creating core and mini-core sets that represent the maximum possible genetic variation contained in the whole collection. Here, we studied the genetic variation and population structure of 2,179 Oryza glaberrima Steud. accessions conserved at the AfricaRice genebank using 27,560 DArTseq-based SNPs. Only 14% (3,834 of 27,560) of the SNPs were polymorphic across the 2,179 accessions, which is much lower than diversity reported in other Oryza species. Genetic distance between pairs of accessions varied from 0.005 to 0.306, with 1.5% of the pairs nearly identical, 8.0% of the pairs similar, 78.1% of the pairs moderately distant, and 12.4% of the pairs very distant. The number of redundant accessions that contribute little or no new genetic variation to the O. glaberrima collection was very low. Using the maximum length sub-tree method, we propose a subset of 1,330 and 350 accessions to represent a core and mini-core collection, respectively. The core and mini-core sets accounted for ~61 and 16%, respectively, of the whole collection, and captured 97-99% of the SNP polymorphism and nearly all allele and genotype frequencies observed in the whole O. glaberrima collection available at the AfricaRice genebank. Cluster, principal component and model-based population structure analyses all divided the 2,179 accessions into five groups, based roughly on country of origin but less so on ecology. The first, third and fourth groups consisted of accessions primarily from Liberia, Nigeria, and Mali, respectively; the second group consisted primarily of accessions from Togo and Nigeria; and the fifth and smallest group was a mixture of accessions from multiple countries. Analysis of molecular variance showed between 10.8 and 28.9% of the variation among groups with the remaining 71.1-89.2% attributable to differences within groups.

Effect of Trait Heritability, Training Population Size and Marker Density on Genomic Prediction Accuracy Estimation in 22 bi-parental Tropical Maize Populations.

Genomic selection is being used increasingly in plant breeding to accelerate genetic gain per unit time. One of the most important applications of genomic selection in maize breeding is to predict and select the best un-phenotyped lines in bi-parental populations based on genomic estimated breeding values. In the present study, 22 bi-parental tropical maize populations genotyped with low density SNPs were used to evaluate the genomic prediction accuracy (rMG ) of the six trait-environment combinations under various levels of training population size (TPS) and marker density (MD), and assess the effect of trait heritability (h2 ), TPS and MD on rMG estimation. Our results showed that: (1) moderate rMG values were obtained for different trait-environment combinations, when 50% of the total genotypes was used as training population and ~200 SNPs were used for prediction; (2) rMG increased with an increase in h2 , TPS and MD, both correlation and variance analyses showed that h2 is the most important factor and MD is the least important factor on rMG estimation for most of the trait-environment combinations; (3) predictions between pairwise half-sib populations showed that the rMG values for all the six trait-environment combinations were centered around zero, 49% predictions had rMG values above zero; (4) the trend observed in rMG differed with the trend observed in rMG /h, and h is the square root of heritability of the predicted trait, it indicated that both rMG and rMG /h values should be presented in GS study to show the accuracy of genomic selection and the relative accuracy of genomic selection compared with phenotypic selection, respectively. This study provides useful information to maize breeders to design genomic selection workflow in their breeding programs.

Genetic variation and population structure of maize inbred lines adapted to the mid-altitude sub-humid maize agro-ecology of Ethiopia using single nucleotide polymorphic (SNP) markers.

BACKGROUND: Molecular characterization is important for efficient utilization of germplasm and development of improved varieties. In the present study, we investigated the genetic purity, relatedness and population structure of 265 maize inbred lines from the Ethiopian Institute of Agricultural Research (EIAR), the International Maize and Wheat Improvement Centre (CIMMYT) and the International Institute of Tropical Agriculture (IITA) using 220,878 single nucleotide polymorphic (SNP) markers obtained using genotyping by sequencing (GBS). RESULTS: Only 22% of the inbred lines were considered pure with <5% heterogeneity, while the remaining 78% of the inbred lines had a heterogeneity ranging from 5.1 to 31.5%. Pairwise genetic distances among the 265 inbred lines varied from 0.011 to 0.345, with 89% of the pairs falling between 0.301 and 0.345. Only <1% of the pairs had a genetic distance lower than 0.200, which included 14 pairs of sister lines that were nearly identical. Relative kinship analysis showed that the kinship coefficients for 59% of the pairs of lines was close to zero, which agrees with the genetic distance estimates. Principal coordinate analysis, discriminant analysis of principal components (DAPC) and the model-based population structure analysis consistently suggested the presence of three groups, which generally agreed with pedigree information (genetic background). Although not distinct enough, the SNP markers showed some level of separation between the two CIMMYT heterotic groups A and B established based on pedigree and combining ability information. CONCLUSIONS: The high level of heterogeneity detected in most of the inbred lines suggested the requirement for purification or further inbreeding except those deliberately maintained at early inbreeding level. The genetic distance and relative kinship analysis clearly indicated the uniqueness of most of the inbred lines in the maize germplasm available for breeders in the mid-altitude maize breeding program of Ethiopia. Results from the present study facilitate the maize breeding work in Ethiopia and germplasm exchange among breeding programs in Africa. We suggest the incorporation of high density molecular marker information in future heterotic group assignments.