A SIMPLE Pipeline for Mapping Point Mutations

MAPtools: command-line tools for mapping-by-sequencing and QTL-Seq analysis and visualization

BACKGROUND

Classical mutagenesis is a powerful tool that has allowed researchers to elucidate the molecular and genetic basis of a plethora of processes in many model species. The integration of these methods with modern massively parallel sequencing techniques, initially in model species but currently also in many crop species, is accelerating the identification of genes underlying a wide range of traits of agronomic interest.

RESULTS

We have developed MAPtools, an open-source Python3 application designed specifically for the analysis of genomic data from bulked segregant analysis experiments, including mapping-by-sequencing (MBS) and quantitative trait locus sequencing (QTL-seq) experiments. We have extensively tested MAPtools using datasets published in recent literature.

CONCLUSIONS

MAPtools gives users the flexibility to customize their bioinformatics pipeline with various commands for calculating allele count-based statistics, generating plots to pinpoint candidate regions, and annotating the effects of SNP and indel mutations. While extensively tested with plants, the program is versatile and applicable to any species for which a mapping population can be generated and a sequenced genome is available.

AVAILABILITY AND IMPLEMENTATION

MAPtools is available under GPL v3.0 license and documented as a Python3 package at https://github.com/hcandela/MAPtools .

Leaf dissection and margin serration are independently regulated by two regulators converging on the CUC2-auxin module in strawberry

The remarkable diversity of leaf forms allows plants to adapt to their living environment. In general, leaf diversity is shaped by leaf complexity (compound or simple) and leaf margin pattern (entire, serrated, or lobed). Prior studies in multiple species have uncovered a conserved module of CUC2-auxin that regulates both leaf complexity and margin serration. How this module is regulated in different species to contribute to the species-specific leaf form is unclear. Furthermore, the mechanistic connection between leaf complexity and leaf serration regulation is not well studied. Strawberry has trifoliate compound leaves with serrations at the margin. In the wild strawberry Fragaria vesca, a mutant named salad was isolated that showed deeper leaf serrations but normal leaf complexity. SALAD encodes a single-Myb domain protein and is expressed at the leaf margin. Genetic analysis showed that cuc2a is epistatic to salad, indicating that SALAD normally limits leaf serration depth by repressing CUC2a expression. When both Arabidopsis homologs of SALAD were knocked out, deeper serrations were observed in Arabidopsis rosette leaves, supporting a conserved function of SALAD in leaf serration regulation. We incorporated the analysis of a third strawberry mutant simple leaf 1 (sl1) with reduced leaf complexity but normal leaf serration. We showed that SL1 and SALAD independently regulate CUC2a at different stages of leaf development to, respectively, regulate leaf complexity and leaf serration. Our results provide a clear and simple mechanism of how leaf complexity and leaf serration are coordinately as well as independently regulated to achieve diverse leaf forms.

An auxin-mediated ultradian rhythm positively influences root regeneration via EAR1/EUR1 in Arabidopsis

Ultradian rhythms have been proved to be critical for diverse biological processes. However, comprehensive understanding of the short-period rhythms remains limited. Here, we discover that leaf excision triggers a gene expression rhythm with ~3-h periodicity, named as the excision ultradian rhythm (UR), which is regulated by the plant hormone auxin. Promoter-luciferase analyses showed that the spatiotemporal patterns of the excision UR were positively associated with de novo root regeneration (DNRR), a post-embryonic developmental process. Transcriptomic analysis indicated more than 4,000 genes including DNRR-associated genes were reprogramed toward ultradian oscillation. Genetic studies showed that EXCISION ULTRADIAN RHYTHM 1 (EUR1) encoding ENHANCER OF ABSCISIC ACID CO-RECEPTOR1 (EAR1), an abscisic acid signaling regulator, was required to generate the excision ultradian rhythm and enhance root regeneration. The eur1 mutant exhibited the absence of auxin-induced excision UR generation and partial failure during rescuing root regeneration. Our results demonstrate a link between the excision UR and adventitious root formation via EAR1/EUR1, implying an additional regulatory layer in plant regeneration.

The multi-BRCT domain protein DDRM2 promotes the recruitment of RAD51 to DNA damage sites to facilitate homologous recombination

DNA double-strand breaks (DSBs) are the most toxic form of DNA damage in cells. Homologous recombination (HR) is an error-free repair mechanism for DSBs as well as a basis for gene targeting using genome-editing techniques. Despite the importance of HR, the HR mechanism in plants is poorly understood. Through genetic screens for DNA damage response mutants (DDRMs), we find that the Arabidopsis ddrm2 mutant is hypersensitive to DSB-inducing reagents. DDRM2 encodes a protein with four BRCA1 C-terminal (BRCT) domains and is highly conserved in plants including bryophytes, the earliest land plant lineage. The plant-specific transcription factor SOG1 binds to the promoter of DDRM2 and activates its expression. In consistence, the expression of DDRM2 is induced by DSBs in a SOG1-dependent manner. In support, genetic analysis suggests that DDRM2 functions downstream of SOG1. Similar to the sog1 mutant, the ddrm2 mutant shows dramatically reduced HR efficiency. Mechanistically, DDRM2 interacts with the core HR protein RAD51 and is required for the recruitment of RAD51 to DSB sites. Our study reveals that SOG1-DDRM2-RAD51 is a novel module for HR, providing a potential target for improving the efficiency of gene targeting.

MicroRNA156 conditions auxin sensitivity to enable growth plasticity in response to environmental changes in Arabidopsis

MicroRNAs (miRNAs) play diverse roles in plant development, but whether and how miRNAs participate in thermomorphogenesis remain ambiguous. Here we show that HYPONASTIC LEAVES 1 (HYL1)-a key component of miRNA biogenesis-acts downstream of the thermal regulator PHYTOCHROME INTERACTING FACTOR 4 in the temperature-dependent plasticity of hypocotyl growth in Arabidopsis. A hyl1-2 suppressor screen identified a dominant dicer-like1 allele that rescues hyl1-2's defects in miRNA biogenesis and thermoresponsive hypocotyl elongation. Genome-wide miRNA and transcriptome analysis revealed microRNA156 (miR156) and its target SQUAMOSA PROMOTER-BINDING-PROTEIN-LIKE 9 (SPL9) to be critical regulators of thermomorphogenesis. Surprisingly, perturbation of the miR156/SPL9 module disengages seedling responsiveness to warm temperatures by impeding auxin sensitivity. Moreover, miR156-dependent auxin sensitivity also operates in the shade avoidance response at lower temperatures. Thus, these results unveil the miR156/SPL9 module as a previously uncharacterized genetic circuit that enables plant growth plasticity in response to environmental temperature and light changes.

Cyclin A participates in the TSO1-MYB3R1 regulatory module to maintain shoot meristem size and fertility in Arabidopsis

The stem cell pools at the shoot apex and root tip give rise to all the above- and below-ground tissues of a plant. Previous studies in Arabidopsis identified a TSO1-MYB3R1 transcriptional module that controls the number and size of the stem cell pools at the shoot apex and root tip. As TSO1 and MYB3R1 are homologous to components of an animal cell cycle regulatory complex, DREAM, Arabidopsis mutants of TSO1 and MYB3R1 provide valuable tools for investigations into the link between cell cycle regulation and stem cell maintenance in plants. In this study, an Arabidopsis cyclin A gene, CYCA3;4, was identified as a member of the TSO1-MYB3R1 regulatory module and cyca3;4 mutations suppressed the tso1-1 mutant phenotype specifically in the shoot. The work reveals how the TSO1-MYB3R1 module is integrated with the cell cycle machinery to control cell division at the shoot meristem.

Versatile mapping-by-sequencing with Easymap v.2

Mapping-by-sequencing combines Next Generation Sequencing (NGS) with classical genetic mapping by linkage analysis to establish gene-to-phenotype relationships. Although numerous tools have been developed to analyze NGS datasets, only a few are available for mapping-by-sequencing. One such tool is Easymap, a versatile, easy-to-use package that performs automated mapping of point mutations and large DNA insertions. Here, we describe Easymap v.2, which also maps small insertion/deletions (InDels), and includes workflows to perform QTL-seq and variant density mapping analyses. Each mapping workflow can accommodate different experimental designs, including outcrossing and backcrossing, F₂, M₂, and M₃ mapping populations, chemically induced mutation and natural variant mapping, input files containing single-end or paired-end reads of genomic or complementary DNA sequences, and alternative control sample files in FASTQ and VCF formats. Easymap v.2 can also be used as a variant analyzer in the absence of a mapping algorithm and includes a multi-threading option.

CORK1, A LRR-Malectin Receptor Kinase, Is Required for Cellooligomer-Induced Responses in Arabidopsis thaliana

Cell wall integrity (CWI) maintenance is central for plant cells. Mechanical and chemical distortions, pH changes, and breakdown products of cell wall polysaccharides activate plasma membrane-localized receptors and induce appropriate downstream responses. Microbial interactions alter or destroy the structure of the plant cell wall, connecting CWI maintenance to immune responses. Cellulose is the major polysaccharide in the primary and secondary cell wall. Its breakdown generates short-chain cellooligomers that induce Ca²⁺-dependent CWI responses. We show that these responses require the malectin domain-containing CELLOOLIGOMER-RECEPTOR KINASE 1 (CORK1) in Arabidopsis and are preferentially activated by cellotriose (CT). CORK1 is required for cellooligomer-induced cytoplasmic Ca²⁺ elevation, reactive oxygen species (ROS) production, mitogen-associated protein kinase (MAPK) activation, cellulose synthase phosphorylation, and the regulation of CWI-related genes, including those involved in biosynthesis of cell wall material, secondary metabolites and tryptophan. Phosphoproteome analyses identified early targets involved in signaling, cellulose synthesis, the endoplasmic reticulum/Golgi secretory pathway, cell wall repair and immune responses. Two conserved phenylalanine residues in the malectin domain are crucial for CORK1 function. We propose that CORK1 is required for CWI and immune responses activated by cellulose breakdown products.

Bulk segregation analysis in the NGS era: a review of its teenage years

Bulk segregation analysis (BSA) utilizes a strategy of pooling individuals with extreme phenotypes to conduct economical and rapidly linked marker screening or quantitative trait locus (QTL) mapping. With the development of next-generation sequencing (NGS) technology in the past 10 years, BSA methods and technical systems have been gradually developed and improved. At the same time, the ever-decreasing costs of sequencing accelerate NGS-based BSA application in different species, including eukaryotic yeast, grain crops, economic crops, horticultural crops, trees, aquatic animals, and insects. This paper provides a landscape of BSA methods and reviews the BSA development process in the past decade, including the sequencing method for BSA, different populations, different mapping algorithms, associated region threshold determination, and factors affecting BSA mapping. Finally, we summarize related strategies in QTL fine mapping combining BSA.

Mapping-by-Sequencing of Point and Insertional Mutations with Easymap

Random mutagenesis followed by screening for phenotypes of interest is a widely used strategy for genetic dissection of biological pathways; however, identifying the causal gene traditionally required time-consuming mapping approaches based on iterative linkage analysis. Mapping-by-sequencing accelerates this process, efficiently linking the phenotype of a mutant to a narrow candidate genomic region, using next-generation sequencing (NGS) data from a mapping population segregating for the mutant phenotype. To enable researchers at any bioinformatics skill level to conduct mapping-by-sequencing, we developed the Easymap mapping software. In this protocol we break down the steps involved in mapping-by-sequencing. First, we describe different ways of obtaining a mapping population and the steps used to generate NGS data. Next, we show how to analyze the NGS data using Easymap and how to obtain a list of candidate mutations, along with comprehensive information for assessing the potential causality of each candidate. Thus, this protocol enables the user to conduct mapping-by-sequencing using Easymap, facilitating the identification of causal loci for a mutant phenotype of interest.

Connecting research and teaching introductory cell and molecular biology using an Arabidopsis mutant screen

A complex research project was translated into a Course-based Undergraduate Research Experience (CURE), which was implemented in sections of an introductory Cell and Molecular Biology laboratory course. The research laboratory generated an engineered plant line producing a growth-inhibiting, lipid-derived plant hormone and mutagenized this line. Students in the CURE cultured the mutagenized plant population and selected and characterized suppressor mutants. They learned to observe phenotypes related to the biosynthesis and perception of the plant hormone and explored the genetic and biochemical basis of these phenotypes. As the students studied the relevant genetic, molecular and biochemical concepts during this CURE, they were able to translate this knowledge into practice and develop scientific arguments. This CURE was a successful collaboration between the teaching lab and the research lab. It benefited both parties as the students had a real-life, deep learning experience in scientific methodology, while the research lab gathered data and materials for further studies.

Random mutagenesis in vegetatively propagated crops: opportunities, challenges and genome editing prospects

In order to meet the growing human food and nutrition demand a perpetual process of crop improvement is idealized. It has seen changing trends and varying concepts throughout human history; from simple selection to complex gene-editing. Among these techniques, random mutagenesis has been shown to be a promising technology to achieve desirable genetic gain with less time and minimal efforts. Over the decade, several hundred varieties have been released through random mutagenesis, but the production is falling behind the demand. Several food crops like banana, potato, cassava, sweet potato, apple, citrus, and others are vegetatively propagated. Since such crops are not propagated through seed, genetic improvement through classical breeding is impractical for them. Besides, in the case of polyploids, accomplishment of allelic homozygosity requires a considerable land area, extensive fieldwork with huge manpower, and hefty funding for an extended period of time. Apart from induction, mapping of induced genes to facilitate the knowledge of biological processes has been performed only in a few selected facultative vegetative crops like banana and cassava which can form a segregating population. During the last few decades, there has been a shift in the techniques used for crop improvement. With the introduction of the robust technologies like meganucleases, zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and clustered regularly interspaced short palindromic repeats (CRISPR) more and more crops are being subjected to gene editing. However, more work needs to be done in case of vegetatively propagated crops.

Loss of an ABC transporter in Arabidopsis thaliana confers hypersensitivity to the anti-cancer drug bleomycin

Bleomycin (BLM) is used as an anti-cancer drug clinically. However, some cancer cells are resistant to BLM, which limits the usage of BLM in chemotherapy. But the underlying mechanism of such resistance is poorly understood. Here we show that the ATP binding cassette (ABC) transporter ABCC3 is required for the BLM-resistance in Arabidopsis. In a genetic screen for ddrm (DNA damage response mutants), we found that loss of ABCC3 confers the hypersensitivity to BLM. In contrast, overexpression of ABCC3 enhances the resistance to BLM. We further found that the expression of ABCC3 is induced by BLM, which is dependent on the protein kinase ATM and the transcription factor SOG1, two master regulators of DNA damage response. Our study revealed that the ABC transporter contributes to BLM-resistance, indicating that the combination of ABC transporter inhibitors and BLM may enhance the efficacy of BLM in cancer therapy.

Anno genominis XX: 20 years of Arabidopsis genomics

Twenty years ago, the Arabidopsis thaliana genome sequence was published. This was an important moment as it was the first sequenced plant genome and explicitly brought plant science into the genomics era. At the time, this was not only an outstanding technological achievement, but it was characterized by a superb global collaboration. The Arabidopsis genome was the seed for plant genomic research. Here, we review the development of numerous resources based on the genome that have enabled discoveries across plant species, which has enhanced our understanding of how plants function and interact with their environments.

Tuning self-renewal in the Arabidopsis stomatal lineage by hormone and nutrient regulation of asymmetric cell division

Asymmetric and self-renewing divisions build and pattern tissues. In the Arabidopsis stomatal lineage, asymmetric cell divisions, guided by polarly localized cortical proteins, generate most cells on the leaf surface. Systemic and environmental signals modify tissue development, but the mechanisms by which plants incorporate such cues to regulate asymmetric divisions are elusive. In a screen for modulators of cell polarity, we identified CONSTITUTIVE TRIPLE RESPONSE1, a negative regulator of ethylene signaling. We subsequently revealed antagonistic impacts of ethylene and glucose signaling on the self-renewing capacity of stomatal lineage stem cells. Quantitative analysis of cell polarity and fate dynamics showed that developmental information may be encoded in both the spatial and temporal asymmetries of polarity proteins. These results provide a framework for a mechanistic understanding of how nutritional status and environmental factors tune stem-cell behavior in the stomatal lineage, ultimately enabling flexibility in leaf size and cell-type composition.

A novel WEE1 pathway for replication stress responses

DNA replication stress poses a severe threat to genome stability and is a hallmark of cancer as well as a target for cancer therapy. It is well known that the evolutionarily conserved protein kinase WEE1 regulates replication stress responses by directly phosphorylating and inhibiting the major cell cycle driver CDKs in many organisms. Here, we report a novel WEE1 pathway. We found that Arabidopsis WEE1 directly interacts with and phosphorylates the E3 ubiquitin ligase FBL17 that promotes the degradation of CDK inhibitors. The phosphorylated FBL17 is further polyubiquitinated and degraded, thereby leading to the accumulation of CDK inhibitors and the inhibition of CDKs. In strong support for this model, either loss of function of FBL17 or overexpression of CDK inhibitors suppresses the hypersensitivity of the wee1 mutant to replication stress. Intriguingly, human WEE1 also phosphorylates and destabilizes the FBL17 equivalent protein SKP2, indicating that this is a conserved mechanism. This study reveals that the WEE1-FBL17/SKP2-CKIs-CDKs axis is a molecular framework for replication stress responses, which may have clinical implications because the WEE1 inhibitor AZD1775 is currently in phase II clinical trial as an anticancer drug.

Easymap: A User-Friendly Software Package for Rapid Mapping-by-Sequencing of Point Mutations and Large Insertions

Mapping-by-sequencing strategies combine next-generation sequencing (NGS) with classical linkage analysis, allowing rapid identification of the causal mutations of the phenotypes exhibited by mutants isolated in a genetic screen. Computer programs that analyze NGS data obtained from a mapping population of individuals derived from a mutant of interest to identify a causal mutation are available; however, the installation and usage of such programs requires bioinformatic skills, modifying or combining pieces of existing software, or purchasing licenses. To ease this process, we developed Easymap, an open-source program that simplifies the data analysis workflows from raw NGS reads to candidate mutations. Easymap can perform bulked segregant mapping of point mutations induced by ethyl methanesulfonate (EMS) with DNA-seq or RNA-seq datasets, as well as tagged-sequence mapping for large insertions, such as transposons or T-DNAs. The mapping analyses implemented in Easymap have been validated with experimental and simulated datasets from different plant and animal model species. Easymap was designed to be accessible to all users regardless of their bioinformatics skills by implementing a user-friendly graphical interface, a simple universal installation script, and detailed mapping reports, including informative images and complementary data for assessment of the mapping results. Easymap is available at http://genetics.edu.umh.es/resources/easymap; its Quickstart Installation Guide details the recommended procedure for installation.

Recent Advances in Carbon and Nitrogen Metabolism in C3 Plants

C and N are the most important essential elements constituting organic compounds in plants. The shoots and roots depend on each other by exchanging C and N through the xylem and phloem transport systems. Complex mechanisms regulate C and N metabolism to optimize plant growth, agricultural crop production, and maintenance of the agroecosystem. In this paper, we cover the recent advances in understanding C and N metabolism, regulation, and transport in plants, as well as their underlying molecular mechanisms. Special emphasis is given to the mechanisms of starch metabolism in plastids and the changes in responses to environmental stress that were previously overlooked, since these changes provide an essential store of C that fuels plant metabolism and growth. We present general insights into the system biology approaches that have expanded our understanding of core biological questions related to C and N metabolism. Finally, this review synthesizes recent advances in our understanding of the trade-off concept that links C and N status to the plant's response to microorganisms.

Cell wall remodeling and vesicle trafficking mediate the root clock in Arabidopsis

In Arabidopsis thaliana, lateral roots initiate in a process preceded by periodic gene expression known as the root clock. We identified the vesicle-trafficking regulator GNOM and its suppressor, ADENOSINE PHOSPHATE RIBOSYLATION FACTOR GTPase ACTIVATION PROTEIN DOMAIN3, as root clock regulators. GNOM is required for the proper distribution of pectin, a mediator of intercellular adhesion, whereas the pectin esterification state is essential for a functional root clock. In sites of lateral root primordia emergence, both esterified and de-esterified pectin variants are differentially distributed. Using a reverse-genetics approach, we show that genes controlling pectin esterification regulate the root clock and lateral root initiation. These results indicate that the balance between esterified and de-esterified pectin states is essential for proper root clock function and the subsequent initiation of lateral root primordia.

Mutagenomics: A Rapid, High-Throughput Method to Identify Causative Mutations from a Genetic Screen

Genetic screens are powerful tools to dissect complex biological processes, but a rate-limiting step is often the cloning of targeted genes. Here, we present a strategy, "mutagenomics," to identify causal mutations from a screen in a high throughput fashion in the absence of backcrossing. Mutagenomics is initiated by sequencing the genomes of the mutants identified, which are then subjected to a three-stage pipeline. The first stage identifies sequence changes in genes previously linked to the targeted pathway. The second stage uses heuristics derived from a simulation strategy to identify genes that are represented by multiple independent alleles more often than expected by chance. The third stage identifies candidate genes for the remaining lines by sequencing multiple lines of common descent. Our simulations indicate that sequencing as few as three to four sibling lines generally results in fewer than five candidate genes. We applied mutagenomics to a screen for Arabidopsis (Arabidopsis thaliana) mutants involved in the response to the phytohormone cytokinin. Mutagenomics identified likely causative genes for many of the mutant lines analyzed from this screen, including 13 alleles of the gene encoding the ARABIDOPSIS HIS KINASE4 cytokinin receptor. The screen also identified 1-AMINOCYCLOPROPANE-1-CARBOXYLATE (ACC) SYNTHASE7, an ACC synthase homolog involved in ethylene biosynthesis, and ELONGATED HYPOCOTYL5 (HY5), a master transcriptional regulator of photomorphogenesis. HY5 was found to mediate a subset of the transcriptional response to cytokinin. Mutagenomics has the potential to accelerate the pace and utility of genetic screens in Arabidopsis.

A quick protocol for the identification and characterization of early growth mutants in tomato

Root system architecture (RSA) manipulation may improve water and nutrient capture by plants under normal and extreme climate conditions. With the aim of initiating the genetic dissection of RSA in tomato, we established a defined ontology that allowed the curated annotation of the observed phenotypes on 12 traits at four consecutive growth stages. In addition, we established a quick approach for the molecular identification of the mutations associated with the trait-of-interest by using a whole-genome sequencing approach that does not require the building of an additional mapping population. As a proof-of-concept, we screened 4543 seedlings from 300 tomato M₃ lines (Solanum lycopersicum L. cv. Micro-Tom) generated by chemical mutagenesis with ethyl methanesulfonate. We studied the growth and early development of both the root system (primary and lateral roots) and the aerial part of the seedlings as well as the wound-induced adventitious roots emerging from the hypocotyl. We identified 659 individuals (belonging to 203 M₃ lines) whose early seedling and RSA phenotypes differed from those of their reference background. We confirmed the genetic segregation of the mutant phenotypes affecting primary root length, seedling viability and early RSA in 31 M₄ families derived from 15 M₃ lines selected in our screen. Finally, we identified a missense mutation in the SlCESA3 gene causing a seedling-lethal phenotype with short roots. Our results validated the experimental approach used for the identification of tomato mutants during early growth, which will allow the molecular identification of the genes involved.

An ethyl methanesulfonate-induced neutral mutant-bridging method efficiently identifies spontaneously mutated genes in rice

Spontaneous mutants are mainly obtained from tissue culture or natural occurrences in plants. The traditional strategy for identifying spontaneously mutated genes is to continuously backcross these mutants to another variety and develop a near-isogenic F₂ population for map-based cloning or bulked segregant analysis. However, this strategy is time-consuming. Here, we have developed a new method to efficiently accelerate the identification process. The chemical mutagen ethyl methanesulfonate was first used to treat the wild type of the spontaneous mutants to induce thousands of neutral mutations. An induced individual without any statistically significant phenotypic changes which was compared with the wild type was chosen as the neutral mutant. The spontaneous mutant was then crossed with the neutral mutant to develop a pseudo-near-isogenic F₂ population in which only the induced neutral mutations and the causal mutation were segregated in the genome. This population ensures that the variation of the mutated trait is controlled only by the spontaneously mutated gene. Finally, after sequencing the neutral mutant and the mutant-type DNA pool of the F₂ population the spontaneous mutation will be identified quickly by bioinformatics analysis. Using this method, two spontaneously mutated genes were identified successfully. Therefore, the neutral mutant-bridging method efficiently identifies spontaneously mutated genes in rice, and its value in other plants is discussed.

Next Generation Sequencing Based Forward Genetic Approaches for Identification and Mapping of Causal Mutations in Crop Plants: A Comprehensive Review

The recent advancements in forward genetics have expanded the applications of mutation techniques in advanced genetics and genomics, ahead of direct use in breeding programs. The advent of next-generation sequencing (NGS) has enabled easy identification and mapping of causal mutations within a short period and at relatively low cost. Identifying the genetic mutations and genes that underlie phenotypic changes is essential for understanding a wide variety of biological functions. To accelerate the mutation mapping for crop improvement, several high-throughput and novel NGS based forward genetic approaches have been developed and applied in various crops. These techniques are highly efficient in crop plants, as it is relatively easy to grow and screen thousands of individuals. These approaches have improved the resolution in quantitative trait loci (QTL) position/point mutations and assisted in determining the functional causative variations in genes. To be successful in the interpretation of NGS data, bioinformatics computational methods are critical elements in delivering accurate assembly, alignment, and variant detection. Numerous bioinformatics tools/pipelines have been developed for such analysis. This article intends to review the recent advances in NGS based forward genetic approaches to identify and map the causal mutations in the crop genomes. The article also highlights the available bioinformatics tools/pipelines for reducing the complexity of NGS data and delivering the concluding outcomes.

Next-generation forward genetic screens: using simulated data to improve the design of mapping-by-sequencing experiments in Arabidopsis

Forward genetic screens have successfully identified many genes and continue to be powerful tools for dissecting biological processes in Arabidopsis and other model species. Next-generation sequencing technologies have revolutionized the time-consuming process of identifying the mutations that cause a phenotype of interest. However, due to the cost of such mapping-by-sequencing experiments, special attention should be paid to experimental design and technical decisions so that the read data allows to map the desired mutation. Here, we simulated different mapping-by-sequencing scenarios. We first evaluated which short-read technology was best suited for analyzing gene-rich genomic regions in Arabidopsis and determined the minimum sequencing depth required to confidently call single nucleotide variants. We also designed ways to discriminate mutagenesis-induced mutations from background Single Nucleotide Polymorphisms in mutants isolated in Arabidopsis non-reference lines. In addition, we simulated bulked segregant mapping populations for identifying point mutations and monitored how the size of the mapping population and the sequencing depth affect mapping precision. Finally, we provide the computational basis of a protocol that we already used to map T-DNA insertions with paired-end Illumina-like reads, using very low sequencing depths and pooling several mutants together; this approach can also be used with single-end reads as well as to map any other insertional mutagen. All these simulations proved useful for designing experiments that allowed us to map several mutations in Arabidopsis.

Identification of Arabidopsis Mutants with Altered Freezing Tolerance

Low temperature is an important determinant in the configuration of natural plant communities and defines the range of distribution and growth of important crops. Some plants, including Arabidopsis thaliana, have evolved sophisticated adaptive mechanisms to tolerate freezing temperatures. Central to this adaptation is the process of cold acclimation. By means of this process, many plants from temperate regions are able to develop or increase their freezing tolerance in response to low, nonfreezing temperatures. The identification and characterization of factors involved in freezing tolerance is crucial to understand the molecular mechanisms underlying the cold acclimation response and has a potential interest to improve crop tolerance to freezing temperatures. Many genes implicated in cold acclimation have been identified in numerous plant species by using molecular approaches followed by reverse genetic analysis. Remarkably, however, direct genetic analyses have not been conveniently exploited in their capacity for identifying genes with pivotal roles in that adaptive response. In this chapter, we describe a protocol for evaluating the freezing tolerance of both nonacclimated and cold acclimated Arabidopsis plants. This protocol allows for the accurate and simple screening of mutant collections for the identification of novel factors involved in freezing tolerance and cold acclimation.

Genome-Wide Identification of Splicing Quantitative Trait Loci (sQTLs) in Diverse Ecotypes of Arabidopsis thaliana

Alternative splicing (AS) of pre-mRNAs contributes to transcriptome diversity and enables plants to generate different protein isoforms from a single gene and/or fine-tune gene expression during different development stages and environmental changes. Although AS is pervasive, the genetic basis for differential isoform usage in plants is still emerging. In this study, we performed genome-wide analysis in 666 geographically distributed diverse ecotypes of Arabidopsis thaliana to identify genomic regions [splicing quantitative trait loci (sQTLs)] that may regulate differential AS. These ecotypes belong to different microclimatic conditions and are part of the relict and non-relict populations. Although sQTLs were spread across the genome, we observed enrichment for trans-sQTL (trans-sQTLs hotspots) on chromosome one. Furthermore, we identified several sQTL (911) that co-localized with trait-linked single nucleotide polymorphisms (SNP) identified in the Arabidopsis genome-wide association studies (AraGWAS). Many sQTLs were enriched among circadian clock, flowering, and stress-responsive genes, suggesting a role for differential isoform usage in regulating these important processes in diverse ecotypes of Arabidopsis. In conclusion, the current study provides a deep insight into SNPs affecting isoform ratios/genes and facilitates a better mechanistic understanding of trait-associated SNPs in GWAS studies. To the best of our knowledge, this is the first report of sQTL analysis in a large set of Arabidopsis ecotypes and can be used as a reference to perform sQTL analysis in the Brassicaceae family. Since whole genome and transcriptome datasets are available for these diverse ecotypes, it could serve as a powerful resource for the biological interpretation of trait-associated loci, splice isoform ratios, and their phenotypic consequences to help produce more resilient and high yield crop varieties.

Expanding Avenue of Fast Neutron Mediated Mutagenesis for Crop Improvement

: Fast neutron (FN) radiation mediated mutagenesis is a unique approach among the several induced mutagenesis methods being used in plant science in terms of impacted mutations. The FN mutagenesis usually creates deletions from few bases to several million bases (Mb). A library of random deletion generated using FN mutagenesis lines can provide indispensable resources for the reverse genetic approaches. In this review, information from several efforts made using FN mutagenesis has been compiled to understand the type of induced mutations, frequency, and genetic stability. Concerns regarding the utilization of FN mutagenesis technique for a plant with different level of ploidy and genome complexity are discussed. We have highlighted the utility of next-generation sequencing techniques that can be efficiently utilized for the characterization of mutant lines as well as for the mapping of causal mutations. Pros and cons of mapping by mutation (MutMap), mutant chromosome sequencing (MutChromSeq), exon capture, whole genome sequencing, MutRen-Seq, and different tilling approaches that can be used for the detection of FN-induced mutation has also been discussed. Genomic resources developed using the FN mutagenesis have been catalogued wooing to meaningful utilization of the available resources. The information provided here will be helpful for the efficient exploration for the crop improvement programs and for better understanding of genetic regulations.

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

Candidate gene identification of existing or induced mutations with pipelines applicable to large genomes

Bulked segregant analysis (BSA) is used to identify existing or induced variants that are linked to phenotypes. Although it is widely used in Arabidopsis and rice, it remains challenging for crops with large genomes, such as maize. Moreover, analysis of huge data sets can present a bottleneck linking phenotypes to their molecular basis, especially for geneticists without programming experience. Here, we identified two genes of maize defective kernel mutants with newly developed analysis pipelines that require no programing skills and should be applicable to any large genome. In the 1970s, Neuffer and Sheridan generated a chemically induced defective kernel (dek) mutant collection with the potential to uncover critical genes for seed development. To locate such mutations, the dek phenotypes were introgressed into two inbred lines to take advantage of maize haplotype variations and their sequenced genomes. We generated two pipelines that take fastq files derived from next-generation (nextGen) paired-end DNA and cDNA sequencing as input, call on several well established and freely available genomic analysis tools to call SNPs and INDELs, and generate lists of the most likely causal mutations together with variant index plots to locate the mutation to a specific sequence position on a chromosome. The pipelines were validated with a known strawberry mutation before cloning the dek mutants, thereby enabling phenotypic analysis of large genomes by next-generation sequencing.

Synteny-based mapping of causal point mutations relevant to sand rice (Agriophyllum squarrosum) trichomeless1 mutant by RNA-sequencing

Sand rice (Agriophyllum squarrosum), a diploid Amaranthaceae species within the order Caryophyllalles, has potential as a crop in response to concerns about ongoing climate change and future food security. Modifying the weedy traits, such as dense trichomes, is important for the domestication of sand rice. In this study, an ethyl methanesulfonate (EMS) induced Agriophyllum squarrosum trichomeless mutant (astcl1) was isolated and genetic analysis revealed that this mutant was controlled by a single recessive locus. RNA sequencing was employed to analyze sequence variations between the mutant and wild-type individuals based on allele frequencies. Synteny-based mapping against two closely related and sequenced species finally delimited the causal mutations into an approximately 18.97 mega base pair (Mb) interval on the top portion of the chromosome 6 (Bv6) of sugar beet (Beta vulgaris) and two intervals (5.56 and 2.69 Mb) on the chromosomes 14 (CqA14) and 6 (CqB06) of quinoa (Chenopodium quinoa). These two quinoa intervals were located in the orthologous chromosomal regions against sugar beet Bv6. Differential expression analysis revealed that the astcl1 mutation only affects the expression of a small subset of genes. Most significantly, 17 lipid transport and metabolism related genes, such as GDSL esterases and very-long-chain 3-oxoacyl-CoA reductase 1, and two R2R3 MYB genes, MYB39 and RAX3, were down-regulated in astcl1 mutants. These results imply that the Astcl1 protein coordinately regulates trichome initiation and cuticle biosynthesis in sand rice.