SNP-Ratio Mapping (SRM): identifying lethal alleles and mutations in complex genetic backgrounds by next-generation sequencing

Genetic identification of SNP markers and candidate genes associated with sugarcane smut resistance using BSR-Seq

Sugarcane smut caused by Sporisorium scitamineum is one of the most severe fungal diseases worldwide. In this study, a cross was made between a smut-resistant variety YT93-159 and a smut-susceptible variety ROC22, and 312 progenies were obtained. Two bulks of progenies were then constructed, one consisted of 27 highly smut resistant progenies and the other 24 smut susceptible progenies. Total RNAs of the progenies of each bulk, were pooled and subject to bulked segregant RNA-sequence analysis (BSR-Seq). A total of 164.44 Gb clean data containing 2,341,449 SNPs and 64,999 genes were obtained, 7,295 of which were differentially expressed genes (DEGs). These DEGs were mainly enriched in stress-related metabolic pathways, including carbon metabolism, phenylalanine metabolism, plant hormone signal transduction, glutathione metabolism, and plant-pathogen interactions. Besides, 45,946 high-quality, credible SNPs, a 1.27 Mb region at Saccharum spontaneum chromosome Chr5B (68,904,827 to 70,172,982), and 129 candidate genes were identified to be associated with smut resistance. Among them, twenty-four genes, either encoding key enzymes involved in signaling pathways or being transcription factors, were found to be very closely associated with stress resistance. RT-qPCR analysis demonstrated that they played a positive role in smut resistance. Finally, a potential molecular mechanism of sugarcane and S. scitamineum interaction is depicted that activations of MAPK cascade signaling, ROS signaling, Ca²⁺ signaling, and PAL metabolic pathway and initiation of the glyoxalase system jointly promote the resistance to S. scitamineum in sugarcane. This study provides potential SNP markers and candidate gene resources for smut resistance breeding in sugarcane.

The GmSNAP11 Contributes to Resistance to Soybean Cyst Nematode Race 4 in Glycine max

Soybean cyst nematode (SCN) has devastating effects on soybean production, making it crucial to identify genes conferring SCN resistance. Here we employed next-generation sequencing-based bulked segregant analysis (BSA) to discover genomic regions, candidate genes, and diagnostic markers for resistance to SCN race 4 (SCN4) in soybean. Phenotypic analysis revealed highly significant differences among the reactions of 145 recombinant inbred lines (RILs) to SCN4. In combination with euclidean distance (ED) and Δsingle-nucleotide polymorphism (SNP)-index analyses, we identified a genomic region on Gm11 (designated as rhg1-paralog) associated with SCN4 resistance. Overexpression and RNA interference analyzes of the two candidate genes identified in this region (GmPLAC8 and GmSNAP11) revealed that only GmSNAP11 significantly contributes to SCN4 resistance. We developed a diagnostic marker for GmSNAP11. Using this marker, together with previously developed markers for SCN-resistant loci, rhg1 and Rhg4, we evaluated the relationship between genotypes and SCN4 resistance in 145 RILs and 30 soybean accessions. The results showed that all the SCN4-resistant lines harbored all the three loci, however, some lines harboring the three loci were still susceptible to SCN4. This suggests that these three loci are necessary for the resistance to SCN4, but they alone cannot confer full resistance. The GmSNAP11 and the diagnostic markers developed could be used in genomic-assisted breeding to develop soybean varieties with increased resistance to SCN4.

A combinatorial strategy to identify various types of QTLs for quantitative traits using extreme phenotype individuals in an F2 population

Theoretical and applied studies demonstrate the difficulty of detecting extremely over-dominant and small-effect genes for quantitative traits via bulked segregant analysis (BSA) in an F₂ population. To address this issue, we proposed an integrated strategy for mapping various types of quantitative trait loci (QTLs) for quantitative traits via a combination of BSA and whole-genome sequencing. In this strategy, the numbers of read counts of marker alleles in two extreme pools were used to predict the numbers of read counts of marker genotypes. These observed and predicted numbers were used to construct a new statistic, G_w, for detecting quantitative trait genes (QTGs), and the method was named dQTG-seq1. This method was significantly better than existing BSA methods. If the goal was to identify extremely over-dominant and small-effect genes, another reserved DNA/RNA sample from each extreme phenotype F₂ plant was sequenced, and the observed numbers of marker alleles and genotypes were used to calculate G_w to detect QTGs; this method was named dQTG-seq2. In simulated and real rice dataset analyses, dQTG-seq2 could identify many more extremely over-dominant and small-effect genes than BSA and QTL mapping methods. dQTG-seq2 may be extended to other heterogeneous mapping populations. The significance threshold of G_w in this study was determined by permutation experiments. In addition, a handbook for the R software dQTG.seq, which is available at https://cran.r-project.org/web/packages/dQTG.seq/index.html, has been provided in the supplemental materials for the users' convenience. This study provides a new strategy for identifying all types of QTLs for quantitative traits in an F₂ population.

Self-incompatibility requires GPI anchor remodeling by the poppy PGAP1 ortholog HLD1

Glycosylphosphatidylinositol-anchored proteins (GPI-APs) are tethered to the outer leaflet of the plasma membrane where they function as key regulators of a plethora of biological processes in eukaryotes. Self-incompatibility (SI) plays a pivotal role regulating fertilization in higher plants through recognition and rejection of "self" pollen. Here, we used Arabidopsis thaliana lines that were engineered to be self-incompatible by expression of Papaver rhoeas SI determinants for an SI suppressor screen. We identify HLD1/AtPGAP1, an ortholog of the human GPI-inositol deacylase PGAP1, as a critical component required for the SI response. Besides a delay in flowering time, no developmental defects were observed in HLD1/AtPGAP1 knockout plants, but SI was completely abolished. We demonstrate that HLD1/AtPGAP1 functions as a GPI-inositol deacylase and that this GPI-remodeling activity is essential for SI. Using GFP-SKU5 as a representative GPI-AP, we show that the HLD1/AtPGAP1 mutation does not affect GPI-AP production and targeting but affects their cleavage and release from membranes in vivo. Our data not only implicate GPI-APs in SI, providing new directions to investigate SI mechanisms, but also identify a key functional role for GPI-AP remodeling by inositol deacylation in planta.

Broadening the horizon of crop research: a decade of advancements in plant molecular genetics to divulge phenotype governing genes

Advancements in sequencing, genotyping, and computational technologies during the last decade (2011-2020) enabled new forward-genetic approaches, which subdue the impediments of precise gene mapping in varied crops. The modern crop improvement programs rely heavily on two major steps-trait-associated QTL/gene/marker's identification and molecular breeding. Thus, it is vital for basic and translational crop research to identify genomic regions that govern the phenotype of interest. Until the advent of next-generation sequencing, the forward-genetic techniques were laborious and time-consuming. Over the last 10 years, advancements in the area of genome assembly, genotyping, large-scale data analysis, and statistical algorithms have led faster identification of genomic variations regulating the complex agronomic traits and pathogen resistance. In this review, we describe the latest developments in genome sequencing and genotyping along with a comprehensive evaluation of the last 10-year headways in forward-genetic techniques that have shifted the focus of plant research from model plants to diverse crops. We have classified the available molecular genetic methods under bulk-segregant analysis-based (QTL-seq, GradedPool-Seq, QTG-Seq, Exome QTL-seq, and RapMap), target sequence enrichment-based (RenSeq, AgRenSeq, and TACCA), and mutation-based groups (MutMap, NIKS algorithm, MutRenSeq, MutChromSeq), alongside improvements in classical mapping and genome-wide association analyses. Newer methods for outcrossing, heterozygous, and polyploid plant genetics have also been discussed. The use of k-mers has enriched the nature of genetic variants which can be utilized to identify the phenotype-causing genes, independent of reference genomes. We envisage that the recent methods discussed herein will expand the repertoire of useful alleles and help in developing high-yielding and climate-resilient crops.

Fine-mapping and candidate gene analysis of a major locus controlling leaf thickness in rice (Oryza sativa L.)

Leaf thickness is an important trait in rice (Oryza sativa L.). It affects both photosynthesis and sink-resource efficiency. However, compared to leaf length and length width, reports seldom focused on leaf thickness due to the complicated measurement and minor difference. To identify the quantitative trait loci (QTL) and explore the genetic mechanism regulating the natural variation of leaf thickness, we crossed a high leaf thickness variety Aixiuzhan (AXZ) to a thin leaf thickness variety Yangdao No.6 (YD 6) and evaluated 585 F₂ individuals. We further use bulked sergeant analysis with whole-genome resequencing (BSA-seq) to identify five genomic regions, including chromosomes 1, 6, 9, 10, and 12. These regions represented significant allele frequency differentiation between thick and thin leaf thickness among the mixed pool offspring. Moreover, we conducted a linkage mapping using 276 individuals derived from the F₂ population. We fine-mapped and confirmed that chromosome 9 contributed the primary explanation of phenotypic variance. We fine-mapped the candidate regions and confirmed that the chromosome 9 region contributed to flag leaf thickness in rice. We observed the virtual cellular slices and found that the bundle sheath cells in YD 6 flag leaf veins are fewer than AXZ. We analyzed the potential regions on chromosome 9 and narrowed the QTL candidate intervals in the 928-kb region. Candidate genes of this major QTL were listed as potentially controlled leaf thickness. These results provide promising evidence that cloning leaf thickness is associated with yield production in rice.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11032-022-01275-y.

Discovery of genomic regions and candidate genes for grain weight employing next generation sequencing based QTL-seq approach in rice (Oryza sativa L.)

Rice (Oryza sativa L.) yield enhancement is one of the prime objectives of plant breeders. Elucidation of the inheritance of grain weight, a key yield component trait, is of paramount importance for raising the yield thresholds in rice. In the present investigation, we employed Next-Generation Sequencing based QTL-seq approach to identify major genomic regions associated with grain weight using mapping populations derived from a cross between BPT5204 and MTU3626. QTL-seq analysis identified three grain weight quantitative trait loci (QTL) viz., qGW1 (35-40 Mb), qGW7 (10-18 Mb), and qGW8 (2-5 Mb) on chromosomes 1, 7 and 8, respectively and all are found to be novel. Further, qGW8 was confirmed through conventional QTL mapping in F₂, F₃ and BC₁F₂ populations and found to explain the phenotypic variance of 17.88%, 16.70% and 15.00%, respectively, indicating a major QTL for grain weight. Based on previous reports, two candidate genes in the qGW8 QTL were predicted i.e., LOC_Os08g01490 (Cytochrome P450), and LOC_Os08g01680 (WD domain, G-beta repeat domain containing protein) and through in silico analysis they were found to be highly expressed in reproductive organs during different stages of grain development. Here, we have demonstrated that QTL-seq is one of the rapid approaches to uncover novel QTLs controlling complex traits. The candidate genes identified in the present study undoubtedly enhance our understanding of the mechanism and inheritance of the grain weight. These candidate genes can be exploited for yield enhancement after confirmation through complementary studies.

Combining next-generation sequencing and progeny testing for rapid identification of induced recessive and dominant mutations in maize M2 individuals

Molecular identification of mutant alleles responsible for certain phenotypic alterations is a central goal of genetic analyses. In this study we describe a rapid procedure suitable for the identification of induced recessive and dominant mutations applied to two Zea mays mutants expressing a dwarf and a pale green phenotype, respectively, which were obtained through pollen ethyl methanesulfonate (EMS) mutagenesis. First, without prior backcrossing, induced mutations (single nucleotide polymorphisms, SNPs) segregating in a (M₂ ) family derived from a heterozygous (M₁ ) parent were identified using whole-genome shotgun (WGS) sequencing of a small number of (M₂ ) individuals with mutant and wild-type phenotypes. Second, the state of zygosity of the mutation causing the phenotype was determined for each sequenced individual by phenotypic segregation analysis of the self-pollinated (M₃ ) offspring. Finally, we filtered for segregating EMS-induced SNPs whose state of zygosity matched the determined state of zygosity of the mutant locus in each sequenced (M₂ ) individuals. Through this procedure, combining sequencing of individuals and Mendelian inheritance, three and four SNPs in linkage passed our zygosity filter for the homozygous dwarf and heterozygous pale green mutation, respectively. The dwarf mutation was found to be allelic to the an1 locus and caused by an insertion in the largest exon of the AN1 gene. The pale green mutation affected the nuclear W2 gene and was caused by a non-synonymous amino acid exchange in encoded chloroplast DNA polymerase with a predicted deleterious effect. This coincided with lower cpDNA levels in pale green plants.

A high-quality cucumber genome assembly enhances computational comparative genomics

Genetic variation is expressed by the presence of polymorphisms in compared genomes of individuals that can be transferred to next generations. The aim of this work was to reveal genome dynamics by predicting polymorphisms among the genomes of three individuals of the highly inbred B10 cucumber (Cucumis sativus L.) line. In this study, bioinformatic comparative genomics was used to uncover cucumber genome dynamics (also called real-time evolution). We obtained a new genome draft assembly from long single molecule real-time (SMRT) sequencing reads and used short paired-end read data from three individuals to analyse the polymorphisms. Using this approach, we uncovered differentiation aspects in the genomes of the inbred B10 line. The newly assembled genome sequence (B10v3) has the highest contiguity and quality characteristics among the currently available cucumber genome draft sequences. Standard and newly designed approaches were used to predict single nucleotide and structural variants that were unique among the three individual genomes. Some of the variant predictions spanned protein-coding genes and their promoters, and some were in the neighbourhood of annotated interspersed repetitive elements, indicating that the highly inbred homozygous plants remained genetically dynamic. This is the first bioinformatic comparative genomics study of a single highly inbred plant line. For this project, we developed a polymorphism prediction method with optimized precision parameters, which allowed the effective detection of small nucleotide variants (SNVs). This methodology could significantly improve bioinformatic pipelines for comparative genomics and thus has great practical potential in genomic metadata handling.

artMAP: A user-friendly tool for mapping ethyl methanesulfonate-induced mutations in Arabidopsis

Mapping-by-sequencing is a rapid method for identifying both natural as well as induced variations in the genome. However, it requires extensive bioinformatics expertise along with the computational infrastructure to analyze the sequencing data and these requirements have limited its widespread adoption. In the current study, we develop an easy to use tool, artMAP, to discover ethyl methanesulfonate (EMS) induced mutations in the Arabidopsis genome. The artMAP pipeline consists of well-established tools including TrimGalore, BWA, BEDTools, SAMtools, and SnpEff which were integrated in a Docker container. artMAP provides a graphical user interface and can be run on a regular laptop and desktop, thereby limiting the bioinformatics expertise required. artMAP can process input sequencing files generated from single or paired-end sequencing. The results of the analysis are presented in interactive graphs which display the annotation details of each mutation. Due to its ease of use, artMAP makes the identification of EMS-induced mutations in Arabidopsis possible with only a few mouse click. The source code of artMAP is available on Github (https://github.com/RihaLab/artMAP).

Whole Genome Resequencing from Bulked Populations as a Rapid QTL and Gene Identification Method in Rice

Most Quantitative Trait Loci (QTL) and gene isolation approaches, such as positional- or map-based cloning, are time-consuming and low-throughput methods. Understanding and detecting the genetic material that controls a phenotype is a key means to functionally analyzing genes as well as to enhance crop agronomic traits. In this regard, high-throughput technologies have great prospects for changing the paradigms of DNA marker revealing, genotyping, and for discovering crop genetics and genomic study. Bulk segregant analysis, based on whole genome resequencing approaches, permits the rapid isolation of the genes or QTL responsible for the causative mutation of the phenotypes. MutMap, MutMap Gap, MutMap+, modified MutMap, and QTL-seq methods are among those approaches that have been confirmed to be fruitful gene mapping approaches for crop plants, such as rice, irrespective of whether the characters are determined by polygenes. As a result, in the present study we reviewed the progress made by all these methods to identify QTL or genes in rice.

The Protein Phosphatases ATUNIS1 and ATUNIS2 Regulate Cell Wall Integrity in Tip-Growing Cells

Fast tip-growing plant cells such as pollen tubes (PTs) and root hairs (RHs) require a robust coordination between their internal growth machinery and modifications of their extracellular rigid, yet extensible, cell wall (CW). Part of this essential coordination is governed by members of the Catharanthus roseus receptor-like kinase1-like (CrRLK1L) subfamily of RLKs with FERONIA (FER) and its closest homologs, ANXUR1 (ANX1) and ANX2, controlling CW integrity during RH and PT growth, respectively. Recently, Leucine-Rich Repeat Extensin 8 (LRX8) to LRX11 were also shown to be important for CW integrity in PTs. We previously reported an anx1 anx2 suppressor screen in Arabidopsis thaliana that revealed MARIS (MRI) as a positive regulator of both FER- and ANX1/2-dependent CW integrity pathways. Here, we characterize a suppressor that exhibits a weak rescue of the anx1 anx2 PT bursting phenotype and a short RH phenotype. The corresponding suppressor mutation causes a D94N substitution in a Type One Protein Phosphatase we named ATUNIS1 (AUN1). We show that AUN1 and its closest homolog, AUN2, are nucleocytoplasmic negative regulators of tip growth. Moreover, we demonstrate that AUN1^D94N and AUN1^H127A harboring mutations in key amino acids of the conserved catalytic site of phosphoprotein phosphatases function as dominant amorphic variants that repress PT growth. Finally, genetic interaction studies using the hypermorph MRI^R240C and amorph AUN1^D94N dominant variants indicate that LRX8-11 and ANX1/2 function in distinct but converging pathways to fine-tune CW integrity during tip growth.

Exploring DNA variant segregation types in pooled genome sequencing enables effective mapping of weeping trait in Malus

To unlock the power of next generation sequencing-based bulked segregant analysis in allele discovery in out-crossing woody species, and to understand the genetic control of the weeping trait, an F1 population from the cross 'Cheal's Weeping' × 'Evereste' was used to create two genomic DNA pools 'weeping' (17 progeny) and 'standard' (16 progeny). Illumina pair-end (2 × 151 bp) sequencing of the pools to a 27.1× (weeping) and a 30.4× (standard) genome (742.3 Mb) coverage allowed detection of 84562 DNA variants specific to 'weeping', 92148 specific to 'standard', and 173169 common to both pools. A detailed analysis of the DNA variant genotypes in the pools predicted three informative segregation types of variants: (type I) in weeping pool-specific variants, and (type II) and (type III) in variants common to both pools, where the first allele is assumed to be weeping linked and the allele shown in bold is a variant in relation to the reference genome. Conducting variant allele frequency and density-based mappings revealed four genomic regions with a significant association with weeping: a major locus, Weeping (W), on chromosome 13 and others on chromosomes 10 (W2), 16 (W3), and 5 (W4). The results from type I variants were noisier and less certain than those from type II and type III variants, demonstrating that although type I variants are often the first choice, type II and type III variants represent an important source of DNA variants that can be exploited for genetic mapping in out-crossing woody species. Confirmation of the mapping of W and W2, investigation into their genetic interactions, and identification of expressed genes in the W and W2 regions provided insight into the genetic control of weeping and its expressivity in Malus.

Indel-seq: a fast-forward genetics approach for identification of trait-associated putative candidate genomic regions and its application in pigeonpea (Cajanus cajan)

Identification of candidate genomic regions associated with target traits using conventional mapping methods is challenging and time-consuming. In recent years, a number of single nucleotide polymorphism (SNP)-based mapping approaches have been developed and used for identification of candidate/putative genomic regions. However, in the majority of these studies, insertion-deletion (Indel) were largely ignored. For efficient use of Indels in mapping target traits, we propose Indel-seq approach, which is a combination of whole-genome resequencing (WGRS) and bulked segregant analysis (BSA) and relies on the Indel frequencies in extreme bulks. Deployment of Indel-seq approach for identification of candidate genomic regions associated with fusarium wilt (FW) and sterility mosaic disease (SMD) resistance in pigeonpea has identified 16 Indels affecting 26 putative candidate genes. Of these 26 affected putative candidate genes, 24 genes showed effect in the upstream/downstream of the genic region and two genes showed effect in the genes. Validation of these 16 candidate Indels in other FW- and SMD-resistant and FW- and SMD-susceptible genotypes revealed a significant association of five Indels (three for FW and two for SMD resistance). Comparative analysis of Indel-seq with other genetic mapping approaches highlighted the importance of the approach in identification of significant genomic regions associated with target traits. Therefore, the Indel-seq approach can be used for quick and precise identification of candidate genomic regions for any target traits in any crop species.

Forward Genetics by Sequencing EMS Variation-Induced Inbred Lines

In order to leverage novel sequencing techniques for cloning genes in eukaryotic organisms with complex genomes, the false positive rate of variant discovery must be controlled for by experimental design and informatics. We sequenced five lines from three pedigrees of ethyl methanesulfonate (EMS)-mutagenized Sorghum bicolor, including a pedigree segregating a recessive dwarf mutant. Comparing the sequences of the lines, we were able to identify and eliminate error-prone positions. One genomic region contained EMS mutant alleles in dwarfs that were homozygous reference sequences in wild-type siblings and heterozygous in segregating families. This region contained a single nonsynonymous change that cosegregated with dwarfism in a validation population and caused a premature stop codon in the Sorghum ortholog encoding the gibberellic acid (GA) biosynthetic enzyme ent-kaurene oxidase. Application of exogenous GA rescued the mutant phenotype. Our method for mapping did not require outcrossing and introduced no segregation variance. This enables work when line crossing is complicated by life history, permitting gene discovery outside of genetic models. This inverts the historical approach of first using recombination to define a locus and then sequencing genes. Our formally identical approach first sequences all the genes and then seeks cosegregation with the trait. Mutagenized lines lacking obvious phenotypic alterations are available for an extension of this approach: mapping with a known marker set in a line that is phenotypically identical to starting material for EMS mutant generation.

Next-Generation Sequencing from Bulked-Segregant Analysis Accelerates the Simultaneous Identification of Two Qualitative Genes in Soybean

Next-generation sequencing (NGS)-based bulked-segregant analysis (BSA) approaches have been proven successful for rapidly mapping genes in plant species. However, most such methods are based on mutants and usually only one gene controlling the mutant phenotype is identified. In this study, NGS-based BSA was employed to map simultaneously two qualitative genes controlling cotyledon color of seed in soybean. Yellow-cotyledon (YC) and green-cotyledon (GC) bulks from progenies of a biparental population (Zhonghuang 30 × Jiyu 102) were sequenced. The SNP-index of each SNP locus in YC and GC bulks was calculated and two genomic regions on chromosomes 1 and 11 harboring, respectively, loci qCC1 and qCC2 were identified by Δ(SNP-index) analysis. These two BSA-seq-derived loci were further validated with SSR markers and fine-mapped. qCC1 was mapped to a 30.7-kb region containing four annotated genes and qCC2 was mapped to a 67.7-kb region with nine genes. These two regions contained, respectively, genes D1 and D2, which had previously been identified by homology-based cloning as being associated with cotyledon color. Sequence analysis of the NGS data also identified a frameshift deletion in the coding region of D1. These results suggested that BSA-seq could accelerate the mapping of loci controlling qualitative traits, even if a trait is controlled by more than one locus.

Arabinogalactan Proteins as Interactors along the Crosstalk between the Pollen Tube and the Female Tissues

Arabinogalactan proteins (AGPs) have long been considered to be implicated in several steps of the reproductive process of flowering plants. Pollen tube growth along the pistil tissues requires a multiplicity of signaling pathways to be activated and turned off precisely, at crucial timepoints, to guarantee successful fertilization and seed production. In the recent years, an outstanding effort has been made by the plant reproduction scientific community in order to better understand this process. This resulted in the discovery of a fairly substantial number of new players essential for reproduction, as well as their modes of action and interactions. Besides all the indications of AGPs involvement in reproduction, there were no convincing evidences about it. Recently, several studies came out to prove what had long been suggested about this complex family of glycoproteins. AGPs consist of a large family of hydroxyproline-rich proteins, predicted to be anchored to the plasma membrane and extremely rich in sugars. These two last characteristics always made them perfect candidates to be involved in signaling mechanisms, in several plant developmental processes. New findings finally relate AGPs to concrete functions in plant reproduction. In this review, it is intended not only to describe how different molecules and signaling pathways are functioning to achieve fertilization, but also to integrate the recent discoveries about AGPs along this process.

From Mendel's discovery on pea to today's plant genetics and breeding : Commemorating the 150th anniversary of the reading of Mendel's discovery

KEY MESSAGE

This work discusses several selected topics of plant genetics and breeding in relation to the 150th anniversary of the seminal work of Gregor Johann Mendel. In 2015, we celebrated the 150th anniversary of the presentation of the seminal work of Gregor Johann Mendel. While Darwin's theory of evolution was based on differential survival and differential reproductive success, Mendel's theory of heredity relies on equality and stability throughout all stages of the life cycle. Darwin's concepts were continuous variation and "soft" heredity; Mendel espoused discontinuous variation and "hard" heredity. Thus, the combination of Mendelian genetics with Darwin's theory of natural selection was the process that resulted in the modern synthesis of evolutionary biology. Although biology, genetics, and genomics have been revolutionized in recent years, modern genetics will forever rely on simple principles founded on pea breeding using seven single gene characters. Purposeful use of mutants to study gene function is one of the essential tools of modern genetics. Today, over 100 plant species genomes have been sequenced. Mapping populations and their use in segregation of molecular markers and marker-trait association to map and isolate genes, were developed on the basis of Mendel's work. Genome-wide or genomic selection is a recent approach for the development of improved breeding lines. The analysis of complex traits has been enhanced by high-throughput phenotyping and developments in statistical and modeling methods for the analysis of phenotypic data. Introgression of novel alleles from landraces and wild relatives widens genetic diversity and improves traits; transgenic methodologies allow for the introduction of novel genes from diverse sources, and gene editing approaches offer possibilities to manipulate gene in a precise manner.

Next-Generation Sequencing-Based Approaches for Mutation Mapping and Identification in Caenorhabditis elegans

The use of next-generation sequencing (NGS) has revolutionized the way phenotypic traits are assigned to genes. In this review, we describe NGS-based methods for mapping a mutation and identifying its molecular identity, with an emphasis on applications in Caenorhabditis elegans In addition to an overview of the general principles and concepts, we discuss the main methods, provide practical and conceptual pointers, and guide the reader in the types of bioinformatics analyses that are required. Owing to the speed and the plummeting costs of NGS-based methods, mapping and cloning a mutation of interest has become straightforward, quick, and relatively easy. Removing this bottleneck previously associated with forward genetic screens has significantly advanced the use of genetics to probe fundamental biological processes in an unbiased manner.

Exploring Genetic Diversity in Plants Using High-Throughput Sequencing Techniques

Food security has emerged as an urgent concern because of the rising world population. To meet the food demands of the near future, it is required to improve the productivity of various crops, not just of staple food crops. The genetic diversity among plant populations in a given species allows the plants to adapt to various environmental conditions. Such diversity could therefore yield valuable traits that could overcome the food-security challenges. To explore genetic diversity comprehensively and to rapidly identify useful genes and/or allele, advanced high-throughput sequencing techniques, also called next-generation sequencing (NGS) technologies, have been developed. These provide practical solutions to the challenges in crop genomics. Here, we review various sources of genetic diversity in plants, newly developed genetic diversity-mining tools synergized with NGS techniques, and related genetic approaches such as quantitative trait locus analysis and genome-wide association study.

The Cell Wall Arabinose-Deficient Arabidopsis thaliana Mutant murus5 Encodes a Defective Allele of REVERSIBLY GLYCOSYLATED POLYPEPTIDE2

Traditional marker-based mapping and next-generation sequencing was used to determine that the Arabidopsis (Arabidopsis thaliana) low cell wall arabinose mutant murus5 (mur5) encodes a defective allele of REVERSIBLY GLYCOSYLATED POLYPEPTIDE2 (RGP2). Marker analysis of 13 F2 confirmed mutant progeny from a recombinant mapping population gave a rough map position on the upper arm of chromosome 5, and deep sequencing of DNA from these 13 lines gave five candidate genes with G→A (C→T) transitions predicted to result in amino acid changes. Of these five, only insertional mutant alleles of RGP2, a gene that encodes a UDP-arabinose mutase that interconverts UDP-arabinopyranose and UDP-arabinofuranose, exhibited the low cell wall arabinose phenotype. The identities of mur5 and two SALK insertional alleles were confirmed by allelism tests and overexpression of wild-type RGP2 complementary DNA placed under the control of the 35S promoter in the three alleles. The mur5 mutation results in the conversion of cysteine-257 to tyrosine-257 within a conserved hydrophobic cluster predicted to be distal to the active site and essential for protein stability and possible heterodimerization with other isoforms of RGP.

Next-generation sequencing for identification of candidate genes for Fusarium wilt and sterility mosaic disease in pigeonpea (Cajanus cajan)

To map resistance genes for Fusarium wilt (FW) and sterility mosaic disease (SMD) in pigeonpea, sequencing-based bulked segregant analysis (Seq-BSA) was used. Resistant (R) and susceptible (S) bulks from the extreme recombinant inbred lines of ICPL 20096 × ICPL 332 were sequenced. Subsequently, SNP index was calculated between R- and S-bulks with the help of draft genome sequence and reference-guided assembly of ICPL 20096 (resistant parent). Seq-BSA has provided seven candidate SNPs for FW and SMD resistance in pigeonpea. In parallel, four additional genotypes were re-sequenced and their combined analysis with R- and S-bulks has provided a total of 8362 nonsynonymous (ns) SNPs. Of 8362 nsSNPs, 60 were found within the 2-Mb flanking regions of seven candidate SNPs identified through Seq-BSA. Haplotype analysis narrowed down to eight nsSNPs in seven genes. These eight nsSNPs were further validated by re-sequencing 11 genotypes that are resistant and susceptible to FW and SMD. This analysis revealed association of four candidate nsSNPs in four genes with FW resistance and four candidate nsSNPs in three genes with SMD resistance. Further, In silico protein analysis and expression profiling identified two most promising candidate genes namely C.cajan_01839 for SMD resistance and C.cajan_03203 for FW resistance. Identified candidate genomic regions/SNPs will be useful for genomics-assisted breeding in pigeonpea.

A subunit of the oligosaccharyltransferase complex is required for interspecific gametophyte recognition in Arabidopsis

Species-specific gamete recognition is a key premise to ensure reproductive success and the maintenance of species boundaries. During plant pollen tube (PT) reception, gametophyte interactions likely allow the species-specific recognition of signals from the PT (male gametophyte) by the embryo sac (female gametophyte), resulting in PT rupture, sperm release, and double fertilization. This process is impaired in interspecific crosses between Arabidopsis thaliana and related species, leading to PT overgrowth and a failure to deliver the sperm cells. Here we show that ARTUMES (ARU) specifically regulates the recognition of interspecific PTs in A. thaliana. ARU, identified in a genome-wide association study (GWAS), exclusively influences interspecific--but not intraspecific--gametophyte interactions. ARU encodes the OST3/6 subunit of the oligosaccharyltransferase complex conferring protein N-glycosylation. Our results suggest that glycosylation patterns of cell surface proteins may represent an important mechanism of gametophyte recognition and thus speciation.

Receptor-like cytoplasmic kinase MARIS functions downstream of CrRLK1L-dependent signaling during tip growth

Growing plant cells need to rigorously coordinate external signals with internal processes. For instance, the maintenance of cell wall (CW) integrity requires the coordination of CW sensing with CW remodeling and biosynthesis to avoid growth arrest or integrity loss. Despite the involvement of receptor-like kinases (RLKs) of the Catharanthus roseus RLK1-like (CrRLK1L) subfamily and the reactive oxygen species-producing NADPH oxidases, it remains largely unknown how this coordination is achieved. ANXUR1 (ANX1) and ANX2, two redundant members of the CrRLK1L subfamily, are required for tip growth of the pollen tube (PT), and their closest homolog, FERONIA, controls root-hair tip growth. Previously, we showed that ANX1 overexpression mildly inhibits PT growth by oversecretion of CW material, whereas pollen tubes of anx1 anx2 double mutants burst spontaneously after germination. Here, we report the identification of suppressor mutants with improved fertility caused by the rescue of anx1 anx2 pollen tube bursting. Mapping of one these mutants revealed an R240C nonsynonymous substitution in the activation loop of a receptor-like cytoplasmic kinase (RLCK), which we named MARIS (MRI). We show that MRI is a plasma membrane-localized member of the RLCK-VIII subfamily and is preferentially expressed in both PTs and root hairs. Interestingly, mri-knockout mutants display spontaneous PT and root-hair bursting. Moreover, expression of the MRI(R240C) mutant, but not its wild-type form, partially rescues the bursting phenotypes of anx1 anx2 PTs and fer root hairs but strongly inhibits wild-type tip growth. Thus, our findings identify a novel positive component of the CrRLK1L-dependent signaling cascade that coordinates CW integrity and tip growth.

TURAN and EVAN mediate pollen tube reception in Arabidopsis Synergids through protein glycosylation

Pollen tube (PT) reception in flowering plants describes the crosstalk between the male and female gametophytes upon PT arrival at the synergid cells of the ovule. It leads to PT growth arrest, rupture, and sperm cell release, and is thus essential to ensure double fertilization. Here, we describe TURAN (TUN) and EVAN (EVN), two novel members of the PT reception pathway that is mediated by the FERONIA (FER) receptor-like kinase (RLK). Like fer, mutations in these two genes lead to PT overgrowth inside the female gametophyte (FG) without PT rupture. Mapping by next-generation sequencing, cytological analysis of reporter genes, and biochemical assays of glycoproteins in RNAi knockdown mutants revealed both genes to be involved in protein N-glycosylation in the endoplasmic reticulum (ER). TUN encodes a uridine diphosphate (UDP)-glycosyltransferase superfamily protein and EVN a dolichol kinase. In addition to their common role during PT reception in the synergids, both genes have distinct functions in the pollen: whereas EVN is essential for pollen development, TUN is required for PT growth and integrity by affecting the stability of the pollen-specific FER homologs ANXUR1 (ANX1) and ANX2. ANX1- and ANX2-YFP reporters are not expressed in tun pollen grains, but ANX1-YFP is degraded via the ER-associated degradation (ERAD) pathway, likely underlying the anx1/2-like premature PT rupture phenotype of tun mutants. Thus, as in animal sperm–egg interactions, protein glycosylation is essential for the interaction between the female and male gametophytes during PT reception to ensure fertilization and successful reproduction.

Protein glycosylation is essential for gametophyte interactions between the male pollen tube and the female ovule in plants, reminiscent of gamete interactions during fertilization in mammals.

In flowering plants, gametes are produced by the haploid, multicellular male (pollen), and female (embryo sac) gametophytes, which develop within the reproductive organs of the flower. Successful fertilization depends on delivery of the sperm cells to the embryo sac, which is embedded in the ovule, by the pollen tube. Upon arrival of the pollen tube at the opening of the ovule, crosstalk between male and female gametophytes, known as pollen tube reception, ensues; the pollen tube slows or stops its growth, then resumes rapid growth, and finally bursts to release the sperm cells and effect double fertilization. Although several members of the pollen tube reception pathway, including the receptor-like kinase FERONIA, have been identified, the molecular mechanisms underlying this communication process remain unclear. Here, we show that protein N-glycosylation is required for normal pollen tube reception. A mutant screen identified two genes, TURAN and EVAN, which are involved in protein N-glycosylation in the endoplasmic reticulum. Both genes act in the FERONIA-mediated pollen tube reception pathway, which is impaired in these mutants. Thus, in plants, a “dual recognition system,” involving interactions between both protein and glycosyl residues on the surface of male and female gametophytes, appears to be required for successful pollen tube reception, conceptually similar to sperm–egg interactions in mammals, for which N-glycosylation of cell surface proteins also plays an important role.

Whole genome sequences in pulse crops: a global community resource to expedite translational genomics and knowledge-based crop improvement

Unprecedented developments in legume genomics over the last decade have resulted in the acquisition of a wide range of modern genomic resources to underpin genetic improvement of grain legumes. The genome enabled insights direct investigators in various ways that primarily include unearthing novel structural variations, retrieving the lost genetic diversity, introducing novel/exotic alleles from wider gene pools, finely resolving the complex quantitative traits and so forth. To this end, ready availability of cost-efficient and high-density genotyping assays allows genome wide prediction to be increasingly recognized as the key selection criterion in crop breeding. Further, the high-dimensional measurements of agronomically significant phenotypes obtained by using new-generation screening techniques will empower reference based resequencing as well as allele mining and trait mapping methods to comprehensively associate genome diversity with the phenome scale variation. Besides stimulating the forward genetic systems, accessibility to precisely delineated genomic segments reveals novel candidates for reverse genetic techniques like targeted genome editing. The shifting paradigm in plant genomics in turn necessitates optimization of crop breeding strategies to enable the most efficient integration of advanced omics knowledge and tools. We anticipate that the crop improvement schemes will be bolstered remarkably with rational deployment of these genome-guided approaches, ultimately resulting in expanded plant breeding capacities and improved crop performance.

Next-generation sequencing as a tool to quickly identify causative EMS-generated mutations

The advent of next generation sequencing has influenced every aspect of biological research. Many labs are now using whole genome sequencing in Arabidopsis thaliana as a means to quickly identify EMS-generated mutations present in isolated mutants. Following identification of these mutations, examination of T-DNA insertional alleles defective in candidate genes or complementation of the mutant phenotype with a wild type copy of candidate genes can be used to verify which mutation is causative for the phenotype of interest. Here, we discuss the benefits and pitfalls of using this method to identify mutations underlying phenotypes.

SHOREmap v3.0: fast and accurate identification of causal mutations from forward genetic screens

Whole-genome resequencing of pools of recombinant mutant genomes allows direct linking of phenotypic traits to causal mutations. Such analysis, called mapping-by-sequencing, combines classical genetic mapping and next-generation sequencing by relying on selection-induced patterns within genome-wide allele frequency (AF) in pooled genomes. Mapping-by-sequencing can be performed with computational tools such as SHOREmap. Previous versions of SHOREmap, however, did not implement standardized analyses, but were specifically designed for particular experimental settings. Here, we introduce the usage of a novel and advanced implementation of SHOREmap (version 3.0), including several new features like file readers for commonly used file formats, SNP marker selection, and a stable calculation of mapping intervals. SHOREmap can be downloaded at shoremap.org.

Using next-generation sequencing to isolate mutant genes from forward genetic screens

The long-lasting success of forward genetic screens relies on the simple molecular basis of the characterized phenotypes, which are typically caused by mutations in single genes. Mapping the location of causal mutations using genetic crosses has traditionally been a complex, multistep procedure, but next-generation sequencing now allows the rapid identification of causal mutations at single-nucleotide resolution even in complex genetic backgrounds. Recent advances of this mapping-by-sequencing approach include methods that are independent of reference genome sequences, genetic crosses and any kind of linkage information, which make forward genetics amenable for species that have not been considered for forward genetic screens so far.

Short-term exposure to fluconazole induces chromosome loss in Candida albicans: an approach to produce haploid cells

Candida albicans is considered to be an obligate diploid fungus. Here, we describe an approach to isolate aneuploids or haploids induced by the short-term (12-16 h) exposure of diploid reference strains SC5314 and CAI4 to the most commonly used antifungal drug, fluconazole, followed by repeated single-cell separation among small morphologically distinct colonies in the inhibition zone. The isolated strains had altered cell morphology and LOH events in the MTL and other marker alleles of the analyzed loci at 8 chromosomes of C. albicans with decreased DNA content. The present study employed next-generation sequencing (NGS) combined flow cytometry analysis of the DNA content to analyze the haploid, autodiploid, and aneuploid strains that arose from the fluconazole treatment instead of using the conventional single nucleotide polymorphism/comparative genome hybridization (SNP/CGH) method. A multiple-alignment tool was also developed based on sequenced data from NGS to establish haplotype mapping for each chromosome of the selected strains. These findings revealed that C. albicans experiences 'concerted chromosome loss' to form strains with homozygous alleles and that it even has a haploid status after short-term exposure to fluconazole. Additionally, we developed a new platform to analyze chromosome copy number using NGS.

A novel viable allele of Arabidopsis CULLIN1 identified in a screen for superroot2 suppressors by next generation sequencing-assisted mapping

Map-based cloning (MBC) is the conventional approach for linking phenotypes to genotypes, and has been successfully used to identify causal mutations in diverse organisms. Next-generation sequencing (NGS) technologies offer unprecedented possibilities to sequence the entire genomes of organisms, thereby in principle enabling direct identification of causal mutations without mapping. However, although mapping-by-sequencing has proven to be a cost effective alternative to classical MBC in particular situations, methods based solely on NGS still have limitations and need to be refined. Aiming to identify the causal mutations in suppressors of Arabidopsis thaliana superroot2 phenotype, generated by ethyl methane sulfonate (EMS) treatment, we combined NGS and classical mapping, to rapidly identify the point mutations and restrict the number of testable candidates by defining the chromosomal intervals containing the causal mutations, respectively. The NGS-assisted mapping approach we describe here facilitates unbiased identification of virtually any causal EMS-generated mutation by overlapping the identification (deep sequencing) and validation (mapping) steps. To exemplify the useful marriage of the two approaches we discuss the strategy used to identify a new viable recessive allele of the Arabidopsis CULLIN1 gene in the non-reference Wassilewskija (Ws-4) accession.

Mapping-by-sequencing accelerates forward genetics in barley

Mapping-by-sequencing has emerged as a powerful technique for genetic mapping in several plant and animal species. As this resequencing-based method requires a reference genome, its application to complex plant genomes with incomplete and fragmented sequence resources remains challenging. We perform exome sequencing of phenotypic bulks of a mapping population of barley segregating for a mutant phenotype that increases the rate of leaf initiation. Read depth analysis identifies a candidate gene, which is confirmed by the analysis of independent mutant alleles. Our method illustrates how the genomic resources of barley together with exome resequencing can underpin mapping-by-sequencing.

Mapping-by-sequencing accelerates forward genetics in barley

Mapping-by-sequencing has emerged as a powerful technique for genetic mapping in several plant and animal species. As this resequencing-based method requires a reference genome, its application to complex plant genomes with incomplete and fragmented sequence resources remains challenging. We perform exome sequencing of phenotypic bulks of a mapping population of barley segregating for a mutant phenotype that increases the rate of leaf initiation. Read depth analysis identifies a candidate gene, which is confirmed by the analysis of independent mutant alleles. Our method illustrates how the genomic resources of barley together with exome resequencing can underpin mapping-by-sequencing.

Growth Control: A Saga of Cell Walls, ROS, and Peptide Receptors

Despite an increasingly detailed understanding of endogenous and environmental growth-controlling signals and their signaling networks, little is known on how these networks are integrated with the cell expansion machinery. Members of the CrRLK1L family control cell wall properties and cell expansion in a variety of developmental and environmental contexts. Two recent reports provide exciting new insights into the mode of action of these RLKs. One study shows that one family member, FERONIA (FER), is required for the production of hydroxyl radicals in the female gametophyte, which causes pollen tube rupture and sperm cell release during fertilization. Another study shows that FER is a receptor for a signaling peptide (Rapid Alkalinization Factor 1 [RALF1]) that triggers cell wall alkalinization and growth arrest, possibly through the inhibition of plasma membrane H⁺-ATPase activity. RALF1 belongs to a large gene family, with a wide range of expression patterns. Other CrRLK1L family members therefore may also be receptors for RALF-like peptides. These findings have important implications for our understanding of the control of cell wall integrity during growth and raise new intriguing questions.

Genetic resources for quantitative trait analysis: novelty and efficiency in design from an Arabidopsis perspective

The use of genetic resources for the analysis of quantitative traits finds its roots in crop breeding but has seen a rejuvenation in Arabidopsis thaliana thanks to specific tools and genomic approaches. Although widely used in numerous crop and natural species, many approaches were first developed in this reference plant. We will discuss the scientific background and historical use of mapping populations in Arabidopsis and highlight the technological innovations that drove the development of novel strategies. We will especially lay emphasis on the methodologies used to generate the diverse population types and designate possible applications. Finally we highlight some of the most recent developments in generating genetic mapping resources and suggest specific usage for these novel tools and concepts.

Forward genetics by genome sequencing reveals that rapid cyanide release deters insect herbivory of Sorghum bicolor

Whole genome sequencing has allowed rapid progress in the application of forward genetics in model species. In this study, we demonstrated an application of next-generation sequencing for forward genetics in a complex crop genome. We sequenced an ethyl methanesulfonate-induced mutant of Sorghum bicolor defective in hydrogen cyanide release and identified the causal mutation. A workflow identified the causal polymorphism relative to the reference BTx623 genome by integrating data from single nucleotide polymorphism identification, prior information about candidate gene(s) implicated in cyanogenesis, mutation spectra, and polymorphisms likely to affect phenotypic changes. A point mutation resulting in a premature stop codon in the coding sequence of dhurrinase2, which encodes a protein involved in the dhurrin catabolic pathway, was responsible for the acyanogenic phenotype. Cyanogenic glucosides are not cyanogenic compounds but their cyanohydrins derivatives do release cyanide. The mutant accumulated the glucoside, dhurrin, but failed to efficiently release cyanide upon tissue disruption. Thus, we tested the effects of cyanide release on insect herbivory in a genetic background in which accumulation of cyanogenic glucoside is unchanged. Insect preference choice experiments and herbivory measurements demonstrate a deterrent effect of cyanide release capacity, even in the presence of wild-type levels of cyanogenic glucoside accumulation. Our gene cloning method substantiates the value of (1) a sequenced genome, (2) a strongly penetrant and easily measurable phenotype, and (3) a workflow to pinpoint a causal mutation in crop genomes and accelerate in the discovery of gene function in the postgenomic era.

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.

User guide for mapping-by-sequencing in Arabidopsis

Mapping-by-sequencing combines genetic mapping with whole-genome sequencing in order to accelerate mutant identification. However, application of mapping-by-sequencing requires decisions on various practical settings on the experimental design that are not intuitively answered. Following an experimentally determined recombination landscape of Arabidopsis and next generation sequencing-specific biases, we simulated more than 400,000 mapping-by-sequencing experiments. This allowed us to evaluate a broad range of different types of experiments and to develop general rules for mapping-by-sequencing in Arabidopsis. Most importantly, this informs about the properties of different crossing scenarios, the number of recombinants and sequencing depth needed for successful mapping experiments.

CrRLK1L receptor-like kinases: not just another brick in the wall

In plants, receptor-like kinases regulate many processes during reproductive and vegetative development. The Arabidopsis subfamily of Catharanthus roseus RLK1-like kinases (CrRLK1Ls) comprises 17 members with a putative extracellular carbohydrate-binding malectin-like domain. Only little is known about the functions of these proteins, although mutant analyses revealed a role during cell elongation, polarized growth, and fertilization. However, the molecular nature of the underlying signal transduction cascades remains largely unknown. CrRLK1L proteins are also involved in biotic and abiotic stress responses. It is likely that carbohydrate-rich ligands transmit a signal, which could originate from cell wall components, an arriving pollen tube, or a pathogen attack. Thus, post-translational modifications could be crucial for CrRLK1L signal transduction and ligand binding.