V-primer: software for the efficient design of genome-wide InDel and SNP markers from multi-sample variant call format (VCF) genotyping data.

DNA markers are indispensable tools in genetics and genomics research as well as in crop breeding, particularly for marker-assisted selection. Recent advances in next-generation sequencing technology have made it easier to obtain genome sequences for various crop species, enabling the large-scale identification of DNA polymorphisms among varieties, which in turn has made DNA marker design more accessible. However, existing primer design software is not suitable for designing many types of genome-wide DNA markers from next-generation sequencing data. Here, we describe the development of V-primer, high-throughput software for designing insertion/deletion, cleaved amplified polymorphic sequence, and single-nucleotide polymorphism (SNP) markers. We validated the applicability of these markers in different crops. In addition, we performed multiplex PCR targeted amplicon sequencing using SNP markers designed with V-primer. Our results demonstrate that V-primer facilitates the efficient and accurate design of primers and is thus a useful tool for genetics, genomics, and crop breeding. V-primer is freely available at https://github.com/ncod3/vprimer.

Rice apoplastic CBM1-interacting protein counters blast pathogen invasion by binding conserved carbohydrate binding module 1 motif of fungal proteins.

When infecting plants, fungal pathogens secrete cell wall-degrading enzymes (CWDEs) that break down cellulose and hemicellulose, the primary components of plant cell walls. Some fungal CWDEs contain a unique domain, named the carbohydrate binding module (CBM), that facilitates their access to polysaccharides. However, little is known about how plants counteract pathogen degradation of their cell walls. Here, we show that the rice cysteine-rich repeat secretion protein OsRMC binds to and inhibits xylanase MoCel10A of the blast fungus pathogen Magnaporthe oryzae, interfering with its access to the rice cell wall and degradation of rice xylan. We found binding of OsRMC to various CBM1-containing enzymes, suggesting that it has a general role in inhibiting the action of CBM1. OsRMC is localized to the apoplast, and its expression is strongly induced in leaves infected with M. oryzae. Remarkably, knockdown and overexpression of OsRMC reduced and enhanced rice defense against M. oryzae, respectively, demonstrating that inhibition of CBM1-containing fungal enzymes by OsRMC is crucial for rice defense. We also identified additional CBM-interacting proteins (CBMIPs) from Arabidopsis thaliana and Setaria italica, indicating that a wide range of plants counteract pathogens through this mechanism.

A jacalin-like lectin domain-containing protein of Sclerospora graminicola acts as an apoplastic virulence effector in plant-oomycete interactions.

The plant extracellular space, including the apoplast and plasma membrane, is the initial site of plant-pathogen interactions. Pathogens deliver numerous secreted proteins, called effectors, into this region to suppress plant immunity and establish infection. Downy mildew caused by the oomycete pathogen Sclerospora graminicola (Sg) is an economically important disease of Poaceae crops including foxtail millet (Setaria italica). We previously reported the genome sequence of Sg and showed that the jacalin-related lectin (JRL) gene family has significantly expanded in this lineage. However, the biological functions of JRL proteins remained unknown. Here, we show that JRL from Sg (SgJRL) functions as an apoplastic virulence effector. We identified eight SgJRLs by protein mass spectrometry analysis of extracellular fluid from Sg-inoculated foxtail millet leaves. SgJRLs consist of a jacalin-like lectin domain and an N-terminal putative secretion signal; SgJRL expression is induced by Sg infection. Heterologous expression of three SgJRLs with N-terminal secretion signal peptides in Nicotiana benthamiana enhanced the virulence of the pathogen Phytophthora palmivora inoculated onto the same leaves. Of the three SgJRLs, SG06536 fused with green fluorescent protein (GFP) localized to the apoplastic space in N. benthamiana leaves. INF1-mediated induction of defence-related genes was suppressed by co-expression of SG06536-GFP. These findings suggest that JRLs are novel apoplastic effectors that contribute to pathogenicity by suppressing plant defence responses.

Sat-BSA: an NGS-based method using local de novo assembly of long reads for rapid identification of genomic structural variations associated with agronomic traits.

Advances in next generation sequencing (NGS)-based methodologies have accelerated the identifications of simple genetic variants such as point mutations and small insertions/deletions (InDels). Structural variants (SVs) including large InDels and rearrangements provide vital sources of genetic diversity for plant breeding. However, their analysis remains a challenge due to their complex nature. Consequently, novel NGS-based approaches are needed to rapidly and accurately identify SVs. Here, we present an NGS-based bulked-segregant analysis (BSA) technique called Sat-BSA (SVs associated with traits) for identifying SVs controlling traits of interest in crops. Sat-BSA targets allele frequencies at all SNP positions to first identify candidate genomic regions associated with a trait, which is then reconstructed by long reads-based local de novo assembly. Finally, the association between SVs, RNA-seq-based gene expression patterns and trait is evaluated for multiple cultivars to narrow down the candidate genes. We applied Sat-BSA to segregating F2 progeny obtained from crosses between turnip cultivars with different tuber colors and successfully isolated two genes harboring SVs that are responsible for tuber phenotypes. The current study demonstrates the utility of Sat-BSA for the identification of SVs associated with traits of interest in species with large and heterozygous genomes.

Lumi-Map, a Real-Time Luciferase Bioluminescence Screen of Mutants Combined with MutMap, Reveals Arabidopsis Genes Involved in PAMP-Triggered Immunity.

Plants recognize pathogen-associated molecular patterns (PAMPs) to activate PAMP-triggered immunity (PTI). However, our knowledge of PTI signaling remains limited. In this report, we introduce Lumi-Map, a high-throughput platform for identifying causative single-nucleotide polymorphisms (SNPs) for studying PTI signaling components. In Lumi-Map, a transgenic reporter plant line is produced that contains a firefly luciferase (LUC) gene driven by a defense gene promoter, which generates luminescence upon PAMP treatment. The line is mutagenized and the mutants with altered luminescence patterns are screened by a high-throughput real-time bioluminescence monitoring system. Selected mutants are subjected to MutMap analysis, a whole-genome sequencing-based method of rapid mutation identification, to identify the causative SNP responsible for the luminescence pattern change. We generated nine transgenic Arabidopsis reporter lines expressing the LUC gene fused to multiple promoter sequences of defense-related genes. These lines generate luminescence upon activation of FLAGELLIN-SENSING 2 (FLS2) by flg22, a PAMP derived from bacterial flagellin. We selected the WRKY29-promoter reporter line to identify mutants in the signaling pathway downstream of FLS2. After screening 24,000 ethylmethanesulfonate-induced mutants of the reporter line, we isolated 22 mutants with altered WRKY29 expression upon flg22 treatment (abbreviated as awf mutants). Although five flg22-insensitive awf mutants harbored mutations in FLS2 itself, Lumi-Map revealed three genes not previously associated with PTI. Lumi-Map has the potential to identify novel PAMPs and their receptors as well as signaling components downstream of the receptors.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Rice Exo70 interacts with a fungal effector, AVR-Pii, and is required for AVR-Pii-triggered immunity.

Vesicle trafficking including the exocytosis pathway is intimately associated with host immunity against pathogens. However, we still have insufficient knowledge about how it contributes to immunity, and how pathogen factors affect it. In this study, we explore host factors that interact with the Magnaporthe oryzae effector AVR-Pii. Gel filtration chromatography and co-immunoprecipitation assays identified a 150 kDa complex of proteins in the soluble fraction comprising AVR-Pii and OsExo70-F2 and OsExo70-F3, two rice Exo70 proteins presumably involved in exocytosis. Simultaneous knockdown of OsExo70-F2 and F3 totally abrogated Pii immune receptor-dependent resistance, but had no effect on Pia- and Pik-dependent resistance. Knockdown levels of OsExo70-F3 but not OsExo70-F2 correlated with reduction of Pii function, suggesting that OsExo70-F3 is specifically involved in Pii-dependent resistance. Under our current experimental conditions, over-expression of AVR-Pii or knockdown of OsExo70-F2 and -F3 genes in rice did not affect the virulence of compatible isolates of M. oryzae. AVR-Pii interaction with OsExo70-F3 appears to play a crucial role in immunity triggered by Pii, suggesting a role for OsExo70 as a decoy or helper in Pii/AVR-Pii interactions.

WIPPER: an accurate and efficient field phenotyping platform for large-scale applications.

More accurate, rapid, and easy phenotyping tools are required to match the recent advances in high-throughput genotyping for accelerating breeding and genetic analysis. The conventional data recording in field notebooks and then inputting data to computers for further analysis is inefficient, time-consuming, laborious, and prone to human error. Here, we report WIPPER (for Wireless Plant Phenotyper), a new phenotyping platform that combines field phenotyping and data recording with the aid of Bluetooth communication, thus saving time and labor not only for field data recoding but also for inputting data to computers. Additionally, it eliminates the risk of human error associated with phenotyping and inputting data. We applied WIPPER to 100 individuals of a rice recombinant inbred line (RIL) for measuring leaf width and relative chlorophyll content (SPAD value), and were able to record an accurate data in a significantly reduced time compared with the conventional method of data collection. We are currently using WIPPER for routine management of rice germplasm including recording and documenting information on phenotypic data, seeds, and DNA for their accelerated utilization in crop breeding.

The tillering phenotype of the rice plastid terminal oxidase (PTOX) loss-of-function mutant is associated with strigolactone deficiency.

The significance of plastid terminal oxidase (PTOX) in phytoene desaturation and chloroplast function has been demonstrated using PTOX-deficient mutants, particularly in Arabidopsis. However, studies on its role in monocots are lacking. Here, we report cloning and characterization of the rice (Oryza sativa) PTOX1 gene. Using Ecotype Targeting Induced Local Lesions IN Genomes (EcoTILLING) and TILLING as forward genetic tools, we identified the causative mutation of an EMS mutant characterized by excessive tillering, semi-dwarfism and leaf variegation that corresponded to the PTOX1 gene. The tillering and semi-dwarf phenotypes of the ptox1 mutant are similar to phenotypes of known strigolactone (SL)-related rice mutants, and both phenotypic traits could be rescued by application of the synthetic SL GR24. The ptox1 mutant accumulated phytoene in white leaf sectors with a corresponding deficiency in ß-carotene, consistent with the expected function of PTOX1 in promoting phytoene desaturase activity. There was also no accumulation of the carotenoid-derived SL ent-2'-epi-5-deoxystrigol in root exudates. Elevated concentrations of auxin were detected in the mutant, supporting previous observations that SL interaction with auxin is important in shoot branching control. Our results demonstrate that PTOX1 is required for both carotenoid and SL synthesis resulting in SL-deficient phenotypes in rice.

MutMap+: genetic mapping and mutant identification without crossing in rice.

Advances in genome sequencing technologies have enabled researchers and breeders to rapidly associate phenotypic variation to genome sequence differences. We recently took advantage of next-generation sequencing technology to develop MutMap, a method that allows rapid identification of causal nucleotide changes of rice mutants by whole genome resequencing of pooled DNA of mutant F2 progeny derived from crosses made between candidate mutants and the parental line. Here we describe MutMap+, a versatile extension of MutMap, that identifies causal mutations by comparing SNP frequencies of bulked DNA of mutant and wild-type progeny of M3 generation derived from selfing of an M2 heterozygous individual. Notably, MutMap+ does not necessitate artificial crossing between mutants and the wild-type parental line. This method is therefore suitable for identifying mutations that cause early development lethality, sterility, or generally hamper crossing. Furthermore, MutMap+ is potentially useful for gene isolation in crops that are recalcitrant to artificial crosses.

MutMap-Gap: whole-genome resequencing of mutant F2 progeny bulk combined with de novo assembly of gap regions identifies the rice blast resistance gene Pii.

Next-generation sequencing allows the identification of mutations responsible for mutant phenotypes by whole-genome resequencing and alignment to a reference genome. However, when the resequenced cultivar/line displays significant structural variation from the reference genome, mutations in the genome regions missing from the reference (gaps) cannot be identified by simple alignment. Here we report on a method called 'MutMap-Gap', which involves delineating a candidate region harboring a mutation of interest using the recently reported MutMap method, followed by de novo assembly, alignment, and identification of the mutation within genome gaps. We applied MutMap-Gap to isolate the blast resistant gene Pii from the rice cv Hitomebore using mutant lines that have lost Pii function. MutMap-Gap should prove useful for cloning genes that exhibit significant structural variations such as disease resistance genes of the nucleotide-binding site-leucine rich repeat (NBS-LRR) class.

QTL-seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations.

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 F2 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.

Discovery of nucleotide polymorphisms in the Musa gene pool by Ecotilling.

Musa (banana and plantain) is an important genus for the global export market and in local markets where it provides staple food for approximately 400 million people. Hybridization and polyploidization of several (sub)species, combined with vegetative propagation and human selection have produced a complex genetic history. We describe the application of the Ecotilling method for the discovery and characterization of nucleotide polymorphisms in diploid and polyploid accessions of Musa. We discovered over 800 novel alleles in 80 accessions. Sequencing and band evaluation shows Ecotilling to be a robust and accurate platform for the discovery of polymorphisms in homologous and homeologous gene targets. In the process of validating the method, we identified two single nucleotide polymorphisms that may be deleterious for the function of a gene putatively important for phototropism. Evaluation of heterozygous polymorphism and haplotype blocks revealed a high level of nucleotide diversity in Musa accessions. We further applied a strategy for the simultaneous discovery of heterozygous and homozygous polymorphisms in diploid accessions to rapidly evaluate nucleotide diversity in accessions of the same genome type. This strategy can be used to develop hypotheses for inheritance patterns of nucleotide polymorphisms within and between genome types. We conclude that Ecotilling is suitable for diversity studies in Musa, that it can be considered for functional genomics studies and as tool in selecting germplasm for traditional and mutation breeding approaches.