Summary: The majority of agronomically important crop traits are quantitative, meaning that they are controlled by multiple genes each with a small effect (quantitative trait loci, QTLs). Mapping and isolation of QTLs is important for efficient crop breeding by marker‐assisted selection (MAS) and for a better understanding of the molecular mechanisms underlying the traits. However, since it requires the development and selection of DNA markers for linkage analysis, QTL analysis has been time‐consuming and labor‐intensive. Here we report the rapid identification of plant QTLs by whole‐genome resequencing of DNAs from two populations each composed of 20–50 individuals showing extreme opposite trait values for a given phenotype in a segregating progeny. We propose to name this approach QTL‐seq as applied to plant species. We applied QTL‐seq to rice recombinant inbred lines and F₂ populations and successfully identified QTLs for important agronomic traits, such as partial resistance to the fungal rice blast disease and seedling vigor. Simulation study showed that QTL‐seq is able to detect QTLs over wide ranges of experimental variables, and the method can be generally applied in population genomics studies to rapidly identify genomic regions that underwent artificial or natural selective sweeps.

To subvert rice (Oryza sativa) host defenses, the devastating ascomycete fungus pathogen Magnaporthe oryzae produces a battery of effector molecules, including some with avirulence (AVR) activity, which are recognized by host resistance (R) proteins resulting in rapid and effective activation of innate immunity. To isolate novel avirulence genes from M. oryzae, we examined DNA polymorphisms of secreted protein genes predicted from the genome sequence of isolate 70-15 and looked for an association with AVR activity. This large-scale study found significantly more presence/absence polymorphisms than nucleotide polymorphisms among 1032 putative secreted protein genes. Nucleotide diversity of M. oryzae among 46 isolates of a worldwide collection was extremely low (θ = 8.2 × 10⁻⁵), suggestive of recent pathogen dispersal. However, no association between DNA polymorphism and AVR was identified. Therefore, we used genome resequencing of Ina168, an M. oryzae isolate that contains nine AVR genes. Remarkably, a total of 1.68 Mb regions, comprising 316 candidate effector genes, were present in Ina168 but absent in the assembled sequence of isolate 70-15. Association analyses of these 316 genes revealed three novel AVR genes, AVR-Pia, AVR-Pii, and AVR-Pik/km/kp, corresponding to five previously known AVR genes, whose products are recognized inside rice cells possessing the cognate R genes. AVR-Pia and AVR-Pii have evolved by gene gain/loss processes, whereas AVR-Pik/km/kp has evolved by nucleotide substitutions and gene gain/loss.

Summary: The Oryza sativa (rice) resistance gene Pia confers resistance to the blast fungus Magnaporthe oryzae carrying the AVR‐Pia avirulence gene. To clone Pia, we employed a multifaceted genomics approach. First, we selected 12 R‐gene analog (RGA) genes encoding nucleotide binding site‐leucine rich repeats (NBS‐LRRs) proteins from a region on chromosome 11 that shows linkage to Pia. By using seven rice accessions, we examined the association between Pia phenotypes and DNA polymorphisms in the 10 genes, which revealed three genes (Os11gRGA3–Os11gRGA5) exhibiting a perfect association with the Pia phenotypes. We also screened ethyl methane sulfonate (EMS)‐treated mutant lines of the rice cultivar ‘Sasanishiki’ harboring Pia, and isolated two mutants that lost the Pia phenotype. DNA sequencing of Os11gRGA3–Os11gRGA5 from the two mutant lines identified independent mutations of major effects in Os11gRGA4. The wild‐type ‘Sasanishiki’ allele of Os11gRGA4 (SasRGA4) complemented Pia function in both mutants, suggesting that SasRGA4 is necessary for Pia function. However, when the rice cultivar ‘Himenomochi’ lacking Pia was transfected with SasRGA4, the Pia phenotype was not recovered. An additional complementation study revealed that the two NBS‐LRR‐type R genes, SasRGA4 and SasRGA5, that are located next to each other and oriented in the opposite direction are necessary for Pia function. A population genetics analysis of SasRGA4 and SasRGA5 suggests that the two genes are under long‐term balancing selection.