Dissecting the forces that shape crops: an interview with Ryohei Terauchi

Ryohei Terauchi is a Professor at Kyoto University and a Group Leader at the Iwate Biotechnology Research Center, Japan, studying the evolution of crops and their pathogens. In this interview, Ryohei describes his research interests, how the revolution in sequencing technology helped improve our understanding of orphan crops, and who are the scientists that inspire him

RIL-StEp: epistasis analysis of rice recombinant inbred lines (RILs) reveals candidate interacting genes that control seed hull color and leaf chlorophyll content

Characterizing epistatic gene interactions is fundamental for understanding the genetic architecture of complex traits. However, due to the large number of potential gene combinations, detecting epistatic gene interactions is computationally demanding. A simple, easy-to-perform method for sensitive detection of epistasis is required. Due to their homozygous nature, use of recombinant inbred lines excludes the dominance effect of alleles and interactions involving heterozygous genotypes, thereby allowing detection of epistasis in a simple and interpretable model. Here, we present an approach called RIL-StEp (recombinant inbred lines stepwise epistasis detection) to detect epistasis using single-nucleotide polymorphisms in the genome. We applied the method to reveal epistasis affecting rice (Oryza sativa) seed hull color and leaf chlorophyll content and successfully identified pairs of genomic regions that presumably control these phenotypes. This method has the potential to improve our understanding of the genetic architecture of various traits of crops and other organisms

The Power of QTL Mapping with RILs

QTL (quantitative trait loci) mapping is commonly used to identify genetic regions responsible to important phenotype variation. A common strategy of QTL mapping is to use recombinant inbred lines (RILs), which are usually established by several generations of inbreeding of an F1 population (usually up to F6 or F7 populations). As this inbreeding process involves a large amount of labor, we are particularly interested in the effect of the number of inbreeding generations on the power of QTL mapping; a part of the labor could be saved if a smaller number of inbreeding provides sufficient power. By using simulations, we investigated the performance of QTL mapping with recombinant inbred lines (RILs). As expected, we found that the power of F4 population could be almost comparable to that of F6 and F7 populations. A potential problem in using F4 population is that a large proportion of RILs are heterozygotes. We here introduced a new method to partly relax this problem. The performance of this method was verified by simulations with a wide range of parameters including the size of the segregation population, recombination rate, genome size and the density of markers. We found our method works better than the commonly used standard method especially when there are a number of heterozygous markers. Our results imply that in most cases, QTL mapping does not necessarily require RILs at F6 or F7 generations; rather, F4 (or even F3) populations would be almost as useful as F6 or F7 populations. Because the cost to establish a number of RILs for many generations is enormous, this finding will cause a reduction in the cost of QTL mapping, thereby accelerating gene mapping in many species

Gill transcriptome response to changes in environmental calcium in the green spotted puffer fish

Abstract Background Calcium ion is tightly regulated in body fluids and for euryhaline fish, which are exposed to rapid changes in environmental [Ca2+], homeostasis is especially challenging. The gill is the main organ of active calcium uptake and therefore plays a crucial role in the maintenance of calcium ion homeostasis. To study the molecular basis of the short-term responses to changing calcium availability, the whole gill transcriptome obtained by Super Serial Analysis of Gene Expression (SuperSAGE) of the euryhaline teleost green spotted puffer fish, Tetraodon nigroviridis, exposed to water with altered [Ca2+] was analysed. Results Transfer of T. nigroviridis from 10 ppt water salinity containing 2.9 mM Ca2+ to high (10 mM Ca2+ ) and low (0.01 mM Ca2+) calcium water of similar salinity for 2-12 h resulted in 1,339 differentially expressed SuperSAGE tags (26-bp transcript identifiers) in gills. Of these 869 tags (65%) were mapped to T. nigroviridis cDNAs or genomic DNA and 497 (57%) were assigned to known proteins. Thirteen percent of the genes matched multiple tags indicating alternative RNA transcripts. The main enriched gene ontology groups belong to Ca2+ signaling/homeostasis but also muscle contraction, cytoskeleton, energy production/homeostasis and tissue remodeling. K-means clustering identified co-expressed transcripts with distinct patterns in response to water [Ca2+] and exposure time. Conclusions The generated transcript expression patterns provide a framework of novel water calcium-responsive genes in the gill during the initial response after transfer to different [Ca2+]. This molecular response entails initial perception of alterations, activation of signaling networks and effectors and suggests active remodeling of cytoskeletal proteins during the initial acclimation process. Genes related to energy production and energy homeostasis are also up-regulated, probably reflecting the increased energetic needs of the acclimation response. This study is the first genome-wide transcriptome analysis of fish gills and is an important resource for future research on the short-term mechanisms involved in the gill acclimation responses to environmental Ca2+ changes and osmoregulation.Peer Reviewe

Whole-genome analysis of recombinant inbred rice lines reveals a quantitative trait locus on chromosome 3 with genotype-by-environment interaction effects

Elucidating genotype-by-environment interactions is fundamental for understanding the interplay between genetic and environmental factors that shape complex traits in crops. Genotype-by-environment interactions are of practical importance, as they determine the performance of cultivars grown in different environments, prompting the need for an efficient approach for evaluating genotype-by-environment interactions. Here, we describe a method for genotype-by-environment detection that involves comparing linear mixed models. This method successfully detected genotype-by-environment interactions in rice (Oryza sativa) recombinant inbred lines grown at 3 locations. We identified a quantitative trait locus (QTL) on chromosome 3 that was associated with heading date, grain number, and leaf length. The effect of this QTL on plant growth–related traits varied with environmental conditions, indicating the presence of genotype-by-environment interactions. Therefore, our method enables a powerful genotype-by-environment detection pipeline that should facilitate the production of high-yielding crops in a given environment

Cross-reactivity of a rice NLR immune receptor to distinct effectors from the rice blast pathogen Magnaporthe oryzae provides partial disease resistance

Unconventional integrated domains in plant intracellular immune receptors of the nucleotide-binding leucine-rich repeat (NLRs) type can directly bind translocated effector proteins from pathogens and thereby initiate an immune response. The rice (Oryza sativa) immune receptor pairs Pik-1/Pik-2 and RGA5/RGA4 both use integrated heavy metal-associated (HMA) domains to bind the effectors AVR-Pik and AVR-Pia, respectively, from the rice blast fungal pathogen Magnaporthe oryzae These effectors both belong to the MAX effector family and share a core structural fold, despite being divergent in sequence. How integrated domains in NLRs maintain specificity of effector recognition, even of structurally similar effectors, has implications for understanding plant immune receptor evolution and function. Here, using plant cell death and pathogenicity assays and protein-protein interaction analyses, we show that the rice NLR pair Pikp-1/Pikp-2 triggers an immune response leading to partial disease resistance toward the "mis-matched" effector AVR-Pia in planta and that the Pikp-HMA domain binds AVR-Pia in vitro We observed that the HMA domain from another Pik-1 allele, Pikm, cannot bind AVR-Pia, and it does not trigger a plant response. The crystal structure of Pikp-HMA bound to AVR-Pia at 1.9 Å resolution revealed a binding interface different from those formed with AVR-Pik effectors, suggesting plasticity in integrated domain-effector interactions. The results of our work indicate that a single NLR immune receptor can bait multiple pathogen effectors via an integrated domain, insights that may enable engineering plant immune receptors with extended disease resistance profiles

Protein engineering expands the effector recognition profile of a rice NLR immune receptor

Plant nucleotide binding, leucine-rich repeat (NLR) receptors detect pathogen effectors and initiate an immune response. Since their discovery, NLRs have been the focus of protein engineering to improve disease resistance. However, this approach has proven challenging, in part due to their narrow response specificity. Previously, we revealed the structural basis of pathogen recognition by the integrated heavy metal associated (HMA) domain of the rice NLR Pikp (Maqbool et al., 2015). Here, we used structure-guided engineering to expand the response profile of Pikp to variants of the rice blast pathogen effector AVR-Pik. A mutation located within an effector-binding interface of the integrated Pikp-HMA domain increased the binding affinity for AVR-Pik variants in vitro and in vivo. This translates to an expanded cell-death response to AVR-Pik variants previously unrecognized by Pikp in planta. The structures of the engineered Pikp-HMA in complex with AVR-Pik variants revealed the mechanism of expanded recognition. These results provide a proof-of-concept that protein engineering can improve the utility of plant NLR receptors where direct interaction between effectors and NLRs is established, particularly where this interaction occurs via integrated domains

Genome Analysis Revives a Forgotten Hybrid Crop Edo-dokoro in the Genus Dioscorea

忘れられた作物「えどどころ」の起原 --ゲノム解析が明らかにする青森県三八上北地域に残る栽培イモの歴史--. 京都大学プレスリリース. 2022-08-10.A rhizomatous Dioscorea crop “Edo-dokoro” was described in old records of Japan, but its botanical identify has not been characterized. We found that Edo-dokoro is still produced by four farmers in Tohoku-machi of Aomori Prefecture, Japan. Rhizomes of Edo-dokoro are a delicacy to the local people and are sold in the markets. Morphological characters of Edo-dokoro suggest its hybrid origin between the two species, D. tokoro and D. tenuipes. Genome analysis revealed that Edo-dokoro is likely originated by hybridization of a male D. tokoro to a female D. tenuipes, followed by a backcross with a male plant of D. tokoro. Edo-dokoro is a typical minor crop possibly maintained for more than 300 years but now almost forgotten from the public. We hypothesize that there are many such uncharacterized genetic heritages passed over generations by small scale farmers that await serious scientific investigation for future use and improvement by using modern genomics information

High-performance pipeline for MutMap and QTL-seq

[Summary] Bulked segregant analysis implemented in MutMap and QTL-seq is a powerful and efficient method to identify loci contributing to important phenotypic traits. However, the previous pipelines were not user-friendly to install and run. Here, we describe new pipelines for MutMap and QTL-seq. These updated pipelines are approximately 5–8 times faster than the previous pipeline, are easier for novice users to use, and can be easily installed through bioconda with all dependencies. [Availability] The new pipelines of MutMap and QTL-seq are written in Python and can be installed via bioconda. The source code and manuals are available online (MutMap: https://github.com/YuSugihara/MutMap, QTL-seq: https://github.com/YuSugihara/QTL-seq)