WGNAM: whole-genome nested association mapping

Deciphering the genetic basis of agronomic, yield, and nutritional traits in rice (Oryza sativa L.) using a saturated GBS-based SNP linkage map

Developing high-yielding rice varieties that possess favorable agronomic characteristics and enhanced grain Zn content is crucial in ensuring food security and addressing nutritional needs. This research employed ICIM, IM, and multi-parent population QTL mapping methods to identify important genetic regions associated with traits such as DF, PH, NT, NP, PL, YLD, TGW, GL, GW, Zn, and Fe. Two populations of recombinant inbred lines consisting of 373 lines were phenotyped for agronomic, yield and grain micronutrient traits for three seasons at IRRI, and genotyped by sequencing. Most of the traits demonstrated moderate to high broad-sense heritability. There was a positive relationship between Zn and Fe contents. The principal components and correlation results revealed a significant negative association between YLD and Zn/Fe. ICIM identified 81 QTLs, while IM detected 36 QTLs across populations. The multi-parent population analysis detected 27 QTLs with six of them consistently detected across seasons. We shortlisted eight candidate genes associated with yield QTLs, 19 genes with QTLs for agronomic traits, and 26 genes with Zn and Fe QTLs. Notable candidate genes included CL4 and d35 for YLD, dh1 for DF, OsIRX10, HDT702, sd1 for PH, OsD27 for NP, whereas WFP and OsIPI1 were associated with PL, OsRSR1 and OsMTP1 were associated to TGW. The OsNAS1, OsRZFP34, OsHMP5, OsMTP7, OsC3H33, and OsHMA1 were associated with Fe and Zn QTLs. We identified promising RILs with acceptable yield potential and high grain Zn content from each population. The major effect QTLs, genes and high Zn RILs identified in our study are useful for efficient Zn biofortification of rice.

Characterization of adaptation mechanisms in sorghum using a multireference back-cross nested association mapping design and envirotyping

Identifying the genetic factors impacting the adaptation of crops to environmental conditions is of key interest for conservation and selection purposes. It can be achieved using population genomics, and evolutionary or quantitative genetics. Here we present a sorghum multireference back-cross nested association mapping population composed of 3,901 lines produced by crossing 24 diverse parents to 3 elite parents from West and Central Africa-back-cross nested association mapping. The population was phenotyped in environments characterized by differences in photoperiod, rainfall pattern, temperature levels, and soil fertility. To integrate the multiparental and multi-environmental dimension of our data we proposed a new approach for quantitative trait loci (QTL) detection and parental effect estimation. We extended our model to estimate QTL effect sensitivity to environmental covariates, which facilitated the integration of envirotyping data. Our models allowed spatial projections of the QTL effects in agro-ecologies of interest. We utilized this strategy to analyze the genetic architecture of flowering time and plant height, which represents key adaptation mechanisms in environments like West Africa. Our results allowed a better characterization of well-known genomic regions influencing flowering time concerning their response to photoperiod with Ma6 and Ma1 being photoperiod-sensitive and the region of possible candidate gene Elf3 being photoperiod-insensitive. We also accessed a better understanding of plant height genetic determinism with the combined effects of phenology-dependent (Ma6) and independent (qHT7.1 and Dw3) genomic regions. Therefore, we argue that the West and Central Africa-back-cross nested association mapping and the presented analytical approach constitute unique resources to better understand adaptation in sorghum with direct application to develop climate-smart varieties.

Multi-parental fungal mapping population study to detect genomic regions associated with Pyrenophora teres f. teres virulence

In recent years multi-parental mapping populations (MPPs) have been widely adopted in many crops to detect quantitative trait loci (QTLs) as this method can compensate for the limitations of QTL analyses using bi-parental mapping populations. Here we report the first multi-parental nested association mapping (MP-NAM) population study used to detect genomic regions associated with host-pathogenic interactions. MP-NAM QTL analyses were conducted on 399 Pyrenophora teres f. teres individuals using biallelic, cross-specific and parental QTL effect models. A bi-parental QTL mapping study was also conducted to compare the power of QTL detection between bi-parental and MP-NAM populations. Using MP-NAM with 399 individuals detected a maximum of eight QTLs with a single QTL effect model whilst only a maximum of five QTLs were detected with an individual bi-parental mapping population of 100 individuals. When reducing the number of isolates in the MP-NAM to 200 individuals the number of QTLs detected remained the same for the MP-NAM population. This study confirms that MPPs such as MP-NAM populations can be successfully used in detecting QTLs in haploid fungal pathogens and that the power of QTL detection with MPPs is greater than with bi-parental mapping populations.

Fine mapping of a novel QTL DM9.1 conferring downy mildew resistance in melon

Downy mildew (DM) is a major foliar disease globally causing great economic loss in melon production. Utilizing disease-resistant cultivars is the most efficient approach for disease control, while discovery of disease-resistant genes is crucial for the success of DM-resistant breeding. To address this problem, two F₂ populations were constructed using the DM-resistant accession PI 442177 in this study, and QTLs conferring DM resistance were mapped using linkage map and QTL-seq analysis, respectively. A high-density genetic map with the length of 1096.7 cM and density of 0.7 cM was generated by using the genotyping-by-sequencing data of a F₂ population. A major QTL DM9.1 with the phenotypic variance explained proportion of 24.3-37.7% was consistently detected at the early, middle, and late growth stages using the genetic map. QTL-seq analyses on the two F₂ populations also validated the presence of DM9.1. Kompetitive Allele-Specific PCR (KASP) assay was further carried out to fine map DM9.1 into 1.0 Mb interval. A KASP marker co-segregating with DM9.1 was successfully developed. These results not only provided valuable information for DM-resistant gene cloning, but also offered useful markers for melon DM-resistant breeding programs.

A multi-locus linear mixed model methodology for detecting small-effect QTLs for quantitative traits in MAGIC, NAM, and ROAM populations

Although multi-parent populations (MPPs) integrate the advantages of linkage and association mapping populations in the genetic dissection of complex traits and especially combine genetic analysis with crop breeding, it is difficult to detect small-effect quantitative trait loci (QTL) for complex traits in multiparent advanced generation intercross (MAGIC), nested association mapping (NAM), and random-open-parent association mapping (ROAM) populations. To address this issue, here we proposed a multi-locus linear mixed model method, namely mppQTL, to detect QTLs, especially small-effect QTLs, in these MPPs. The new method includes two steps. The first is genome-wide scanning based on a single-locus linear mixed model; the P-values are obtained from likelihood-ratio test, the peaks of negative logarithm P-value curve are selected by group-lasso, and all the selected peaks are regarded as potential QTLs. In the second step, all the potential QTLs are placed on a multi-locus linear mixed model, all the effects are estimated using expectation-maximization empirical Bayes algorithm, and all the non-zero effect vectors are further evaluated via likelihood-ratio test for significant QTLs. In Monte Carlo simulation studies, the new method has higher power in QTL detection, lower false positive rate, lower mean absolute deviation for QTL position estimate, and lower mean squared error for the estimate of QTL size (r²) than existing methods because the new method increases the power of detecting small-effect QTLs. In real dataset analysis, the new method (19) identified five more known genes than the existing three methods (14). This study provides an effective method for detecting small-effect QTLs in any MPPs.