Mapping mutations in plant genomes with the user-friendly web application CandiSNP

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 .

Spliceosomal complex components are critical for adjusting the C:N balance during high-light acclimation

Plant acclimation to an ever-changing environment is decisive for growth, reproduction, and survival. Light availability limits biomass production on both ends of the intensity spectrum. Therefore, the adjustment of plant metabolism is central to high-light (HL) acclimation, and the accumulation of photoprotective anthocyanins is commonly observed. However, mechanisms and factors regulating the HL acclimation response are less clear. Two Arabidopsis mutants of spliceosome components exhibiting a pronounced anthocyanin overaccumulation in HL were isolated from a forward genetic screen for new factors crucial for plant acclimation. Time-resolved physiological, transcriptome, and metabolome analysis revealed a vital function of the spliceosome components for rapidly adjusting gene expression and metabolism. Deficiency of INCREASED LEVEL OF POLYPLOIDY1 (ILP1), NTC-RELATED PROTEIN1 (NTR1), and PLEIOTROPIC REGULATORY LOCUS1 (PRL1) resulted in a marked overaccumulation of carbohydrates and strongly diminished amino acid biosynthesis in HL. While not generally limited in N-assimilation, ilp1, ntr1, and prl1 showed higher glutamate levels and reduced amino acid biosynthesis in HL. The comprehensive analysis reveals a function of the spliceosome components in the conditional regulation of the carbon:nitrogen balance and the accumulation of anthocyanins during HL acclimation. The importance of gene expression, metabolic regulation, and re-direction of carbon towards anthocyanin biosynthesis for HL acclimation are discussed.

A pgr5 suppressor screen uncovers two distinct suppression mechanisms and links cytochrome b6f complex stability to PGR5

PROTON GRADIENT REGULATION5 (PGR5) is thought to promote cyclic electron flow, and its deficiency impairs photosynthetic control and increases photosensitivity of photosystem (PS) I, leading to seedling lethality under fluctuating light (FL). By screening for Arabidopsis (Arabidopsis thaliana) suppressor mutations that rescue the seedling lethality of pgr5 plants under FL, we identified a portfolio of mutations in 12 different genes. These mutations affect either PSII function, cytochrome b6f (cyt b6f) assembly, plastocyanin (PC) accumulation, the CHLOROPLAST FRUCTOSE-1,6-BISPHOSPHATASE1 (cFBP1), or its negative regulator ATYPICAL CYS HIS-RICH THIOREDOXIN2 (ACHT2). The characterization of the mutants indicates that the recovery of viability can in most cases be explained by the restoration of PSI donor side limitation, which is caused by reduced electron flow to PSI due to defects in PSII, cyt b6f, or PC. Inactivation of cFBP1 or its negative regulator ACHT2 results in increased levels of the NADH dehydrogenase-like complex. This increased activity may be responsible for suppressing the pgr5 phenotype under FL conditions. Plants that lack both PGR5 and DE-ETIOLATION-INDUCED PROTEIN1 (DEIP1)/NEW TINY ALBINO1 (NTA1), previously thought to be essential for cyt b6f assembly, are viable and accumulate cyt b6f. We suggest that PGR5 can have a negative effect on the cyt b6f complex and that DEIP1/NTA1 can ameliorate this negative effect.

Arabidopsis translation initiation factor binding protein CBE1 negatively regulates accumulation of the NADPH oxidase respiratory burst oxidase homolog D

Cell surface pattern recognition receptors sense invading pathogens by binding microbial or endogenous elicitors to activate plant immunity. These responses are under tight control to avoid excessive or untimely activation of cellular responses, which may otherwise be detrimental to host cells. How this fine-tuning is accomplished is an area of active study. We previously described a suppressor screen that identified Arabidopsis thaliana mutants with regained immune signaling in the immunodeficient genetic background bak1-5, which we named modifier of bak1-5 (mob) mutants. Here, we report that bak1-5 mob7 mutant restores elicitor-induced signaling. Using a combination of map-based cloning and whole-genome resequencing, we identified MOB7 as conserved binding of eIF4E1 (CBE1), a plant-specific protein that interacts with the highly conserved eukaryotic translation initiation factor eIF4E1. Our data demonstrate that CBE1 regulates the accumulation of respiratory burst oxidase homolog D, the NADPH oxidase responsible for elicitor-induced apoplastic reactive oxygen species production. Furthermore, several mRNA decapping and translation initiation factors colocalize with CBE1 and similarly regulate immune signaling. This study thus identifies a novel regulator of immune signaling and provides new insights into reactive oxygen species regulation, potentially through translational control, during plant stress responses.

Commonalities and specialties in photosynthetic functions of PROTON GRADIENT REGULATION5 variants in Arabidopsis

The PROTON GRADIENT REGULATION5 (PGR5) protein is required for trans-thylakoid proton gradient formation and acclimation to fluctuating light (FL). PGR5 functionally interacts with two other thylakoid proteins, PGR5-like 1 (PGRL1) and 2 (PGRL2); however, the molecular details of these interactions are largely unknown. In the Arabidopsis (Arabidopsis thaliana) pgr5-1 mutant, the PGR5G130S protein accumulates in only small amounts. In this work, we generated a knockout allele of PGR5 (pgr5-Cas) using CRISPR-Cas9 technology. Like pgr5-1, pgr5-Cas is seedling-lethal under FL, but photosynthesis and particularly cyclic electron flow, as well as chlorophyll content, are less severely affected in both pgr5-Cas and pgrl1ab (which lacks PGRL1 and PGR5) than in pgr5-1. These differences are associated with changes in the levels of 260 proteins, including components of the Calvin-Benson cycle, photosystems II and I, and the NDH complex, in pgr5-1 relative to the wild type (WT), pgr5-Cas, and pgrl1ab. Some of the differences between pgr5-1 and the other mutant lines could be tentatively assigned to second-site mutations in the pgr5-1 line, identified by whole-genome sequencing. However, others, particularly the more pronounced photosynthetic defects and PGRL1 depletion (compared to pgr5-Cas), are clearly due to specific negative effects of the amino-acid substitution in PGR5G130S, as demonstrated by complementation analysis. Moreover, pgr5-1 and pgr5-Cas plants are less tolerant to long-term exposure to high light than pgrl1ab plants. These results imply that, in addition to the previously reported necessity of PGRL1 for optimal PGR5 function, PGR5 is required alongside PGRL1 to avoid harmful effects on plant performance.

Harnessing the potential of bulk segregant analysis sequencing and its related approaches in crop breeding

Most plant traits are governed by polygenes including both major and minor genes. Linkage mapping and positional cloning have contributed greatly to mapping genomic loci controlling important traits in crop species. However, they are low-throughput, time-consuming, and have low resolution due to which their efficiency in crop breeding is reduced. In this regard, the bulk segregant analysis sequencing (BSA-seq) and its related approaches, viz., quantitative trait locus (QTL)-seq, bulk segregant RNA-Seq (BSR)-seq, and MutMap, have emerged as efficient methods to identify the genomic loci/QTLs controlling specific traits at high resolution, accuracy, reduced time span, and in a high-throughput manner. These approaches combine BSA with next-generation sequencing (NGS) and enable the rapid identification of genetic loci for qualitative and quantitative assessments. Many previous studies have shown the successful identification of the genetic loci for different plant traits using BSA-seq and its related approaches, as discussed in the text with details. However, the efficiency and accuracy of the BSA-seq depend upon factors like sequencing depth and coverage, which enhance the sequencing cost. Recently, the rapid reduction in the cost of NGS together with the expected cost reduction of third-generation sequencing in the future has further increased the accuracy and commercial applicability of these approaches in crop improvement programs. This review article provides an overview of BSA-seq and its related approaches in crop breeding together with their merits and challenges in trait mapping.

Loss of a pyridoxal-phosphate phosphatase rescues Arabidopsis lacking an endoplasmic reticulum ATP carrier

The endoplasmic reticulum (ER)-located ATP/ADP-antiporter (ER-ANT1) occurs specifically in vascular plants. Structurally different transporters mediate energy provision to the ER, but the cellular function of ER-ANT1 is still unknown. Arabidopsis (Arabidopsis thaliana) mutants lacking ER-ANT1 (er-ant1 plants) exhibit a photorespiratory phenotype accompanied by high glycine levels and stunted growth, pointing to an inhibition of glycine decarboxylase (GDC). To reveal whether it is possible to suppress this marked phenotype, we exploited the power of a forward genetic screen. Absence of a so far uncharacterized member of the HaloAcid Dehalogenase (HAD)-like hydrolase family strongly suppressed the dwarf phenotype of er-ant1 plants. Localization studies suggested that the corresponding protein locates to chloroplasts, and activity assays showed that the enzyme dephosphorylates, with high substrate affinity, the B6 vitamer pyridoxal 5'-phosphate (PLP). Additional physiological experiments identified imbalances in vitamin B6 homeostasis in er-ant1 mutants. Our data suggest that impaired chloroplast metabolism, but not decreased GDC activity, causes the er-ant1 mutant dwarf phenotype. We present a hypothesis, setting transport of PLP by ER-ANT1 and chloroplastic PLP dephosphorylation in the cellular context. With the identification of this HAD-type PLP phosphatase, we also provide insight into B6 vitamer homeostasis.

The RNA-binding protein RBP45D of Arabidopsis promotes transgene silencing and flowering time

RNA-directed DNA methylation (RdDM) helps to defend plants against invasive nucleic acids. In the canonical form of RdDM, 24-nt small interfering RNAs (siRNAs) are produced by DICER-LIKE 3 (DCL3). The siRNAs are loaded onto ARGONAUTE (AGO) proteins leading ultimately to de novo DNA methylation. Here, we introduce the Arabidopsis thaliana prors1 (LUC) transgenic system, in which 24-nt siRNAs are generated to silence the promoter-LUC construct. A forward genetic screen performed with this system identified, besides known components of RdDM (NRPD2A, RDR2, AGO4 and AGO6), the RNA-binding protein RBP45D. RBP45D is involved in CHH (where H is A, C or T) DNA methylation, and maintains siRNA production originating from the LUC transgene. RBP45D is localized to the nucleus, where it is associated with small nuclear RNAs (snRNAs) and small nucleolar RNAs (snoRNAs). RNA-Seq analysis showed that in CRISPR/Cas-mediated rbp-ko lines FLOWERING LOCUS C (FLC) mRNA levels are upregulated and several loci differentially spliced, among them FLM. In consequence, loss of RBP45D delays flowering, presumably mediated by the release of FLC levels and/or alternative splicing of FLM. Moreover, because levels and processing of transcripts of known RdDM genes are not altered in rbp-ko lines, RBP45D should have a more direct function in transgene silencing, probably independent of the canonical RdDM pathway. We suggest that RBP45D facilitates siRNA production by stabilizing either the precursor RNA or the slicer protein. Alternatively, RBP45D could be involved in chromatin modifications, participate in retention of Pol IV transcripts and/or in Pol V-dependent lncRNA retention in chromatin to enable their scaffold function.

A novel allele of the Arabidopsis thaliana MACPF protein CAD1 results in deregulated immune signaling

Immune recognition in plants is governed by two major classes of receptors: pattern recognition receptors (PRRs) and nucleotide-binding leucine-rich repeat receptors (NLRs). Located at the cell surface, PRRs bind extracellular ligands originating from microbes (indicative of "non-self") or damaged plant cells (indicative of "infected-self"), and trigger signaling cascades to protect against infection. Located intracellularly, NLRs sense pathogen-induced physiological changes and trigger localized cell death and systemic resistance. Immune responses are under tight regulation in order to maintain homeostasis and promote plant health. In a forward-genetic screen to identify regulators of PRR-mediated immune signaling, we identified a novel allele of the membrane-attack complex and perforin (MACPF)-motif containing protein CONSTITUTIVE ACTIVE DEFENSE 1 (CAD1) resulting from a missense mutation in a conserved N-terminal cysteine. We show that cad1-5 mutants display deregulated immune signaling and symptoms of autoimmunity dependent on the lipase-like protein ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1), suggesting that CAD1 integrity is monitored by the plant immune system. We further demonstrate that CAD1 localizes to both the cytosol and plasma membrane using confocal microscopy and subcellular fractionation. Our results offer new insights into immune homeostasis and provide tools to further decipher the intriguing role of MACPF proteins in plants.

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.

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.

Cellulose defects in the Arabidopsis secondary cell wall promote early chloroplast development

Lincomycin (LIN)-mediated inhibition of protein synthesis in chloroplasts prevents the greening of seedlings, represses the activity of photosynthesis-related genes in the nucleus, including LHCB1.2, and induces the phenylpropanoid pathway, resulting in the production of anthocyanins. In genomes uncoupled (gun) mutants, LHCB1.2 expression is maintained in the presence of LIN or other inhibitors of early chloroplast development. In a screen using concentrations of LIN lower than those employed to isolate gun mutants, we have identified happy on lincomycin (holi) mutants. Several holi mutants show an increased tolerance to LIN, exhibiting de-repressed LHCB1.2 expression and chlorophyll synthesis in seedlings. The mutations responsible were identified by whole-genome single-nucleotide polymorphism (SNP) mapping, and most were found to affect the phenylpropanoid pathway; however, LHCB1.2 expression does not appear to be directly regulated by phenylpropanoids, as indicated by the metabolic profiling of mutants. The most potent holi mutant is defective in a subunit of cellulose synthase encoded by IRREGULAR XYLEM 3, and comparative analysis of this and other cell-wall mutants establishes a link between secondary cell-wall integrity and early chloroplast development, possibly involving altered ABA metabolism or sensing.

A SNP Mutation of SiCRC Regulates Seed Number Per Capsule and Capsule Length of cs1 Mutant in Sesame

Seed number per capsule (SNC) is a major factor influencing seed yield and is an important trait with complex gene interaction effects. We first performed genetic analysis, gene cloning, and molecular mechanism study for an EMS-induced sesame mutant cs1 with fewer SNC and shorter capsule length (CL). The mutant traits were due to the pleiotropism of a regressive gene (Sics1). Capsule hormone determination showed that five out of 12 hormones, including auxin indole-3-acetic acid (IAA), had significantly different levels between wild type (WT) and mutant type (MT). KEGG pathway analysis showed that plant hormone signal transduction, especially the auxin signal transduction pathway, was the most abundant differentially expressed signaling pathway. After the cross-population association and regional genome screening, we found that three homozygous loci were retained in cs1. Further analysis of these three loci resulted in the identification of SiCRC as the candidate gene for cs1. SiCRC consists of seven exons and six introns encoding 163 amino acids. The SiCRC in cs1 showed a point mutation at intron 5 and exon 6 junction, resulting in the splice site being frame-shifted eight nucleotides further downstream, causing incorrect splicing. Taken together, we assumed the SNP mutation in SiCRC disrupted the function of the transcription factor, which might act downstream of the CRC-auxin signal transduction pathway, resulting in a shorter CL and less SNC mutation of cs1 in sesame. Our results highlight the molecular framework underlying the transcription factor CRC-mediated role of auxin transduction in SNC and CL development.

Sequencing of bulks of segregants allows dissection of genetic control of amylose content in rice

Amylose content (AC) is a key quality trait in rice. A cross between Oryza glaberrima (African rice) and Oryza sativa (Asian rice) segregating for AC was analysed by sequencing bulks of individuals with high and low AC. SNP associated with the granule bound starch synthase (GBSS1) locus on chromosome 6 were polymorphic between the bulks. In particular, a G/A SNP that would result in an Asp to Asn mutation was identified. This amino acid substitution may be responsible for differences in GBSS activity as it is adjacent to a disulphide linkage conserved in all grass GBSS proteins. Other polymorphisms in genomic regions closely surrounding this variation may be the result of linkage drag. In addition to the variant in the starch biosynthesis gene, SNP on chromosomes 1 and 11 linked to AC was also identified. SNP was found in the genes encoding the NAC and CCAAT-HAP5 transcription factors that have previously been linked to starch biosynthesis. This study has demonstrated that the approach of sequencing bulks was able to identify genes on different chromosomes associated with this complex trait.