Getting started in mapping-by-sequencing

Cloning of an Albino Mutation of Arabidopsis thaliana Using Mapping-by-Sequencing

We report the molecular characterization of an ethyl methanesulfonate (EMS)-induced mutation that causes albinism and lethality at the seedling stage in Arabidopsis thaliana. We identified the mutation using a mapping-by-sequencing approach that uses Fisher's exact tests to detect changes in allele frequencies among the seedlings of an F₂ mapping population, which had been pooled according to their phenotypes (wild-type or mutant). After purifying genomic DNA from the plants of both pools, the two samples were sequenced using the Illumina HiSeq 2500 next-generation sequencing platform. The bioinformatic analysis allowed us to identify a point mutation that damages a conserved residue at the acceptor site of an intron of the At2g04030 gene, which encodes the chloroplast-localized AtHsp90.5 protein, a member of the HSP90 family of heat shock proteins. Our RNA-seq analysis demonstrates that the new allele alters the splicing of At2g04030 transcripts in multiple ways, leading to massive deregulation of genes encoding plastid-localized proteins. A search for protein-protein interactions using the yeast two-hybrid method allowed us to identify two members of the GrpE superfamily as potential interactors of AtHsp90.5, as has previously been reported for green algae.

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.

Cloning the barley nec3 disease lesion mimic mutant using complementation by sequencing

Disease lesion mimic (DLM) or necrotic mutants display necrotic lesions in the absence of pathogen infections. They can show improved resistance to some pathogens and their molecular dissection can contribute to revealing components of plant defense pathways. Although forward-genetics strategies to find genes causal to mutant phenotypes are available in crops, these strategies require the production of experimental cross populations, mutagenesis, or gene editing and are time- and resource-consuming or may have to deal with regulated plant materials. In this study, we described a collection of 34 DLM mutants in barley (Hordeum vulgare L.) and applied a novel method called complementation by sequencing (CBS), which enables the identification of the gene responsible for a mutant phenotype given the availability of two or more chemically mutagenized individuals showing the same phenotype. Complementation by sequencing relies on the feasibility to obtain all induced mutations present in chemical mutants and on the low probability that different individuals share the same mutated genes. By CBS, we identified a cytochrome P450 CYP71P1 gene as responsible for orange blotch DLM mutants, including the historical barley nec3 locus. By comparative phylogenetic analysis we showed that CYP71P1 gene family emerged early in angiosperm evolution but has been recurrently lost in some lineages including Arabidopsis thaliana (L.) Heynh. Complementation by sequencing is a straightforward cost-effective approach to clone genes controlling phenotypes in a chemically mutagenized collection. The TILLMore (TM) collection will be instrumental for understanding the molecular basis of DLM phenotypes and to contribute knowledge about mechanisms of host-pathogen interaction.

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.

ATP-Dependent Clp Protease Subunit C1, HvClpC1, Is a Strong Candidate Gene for Barley Variegation Mutant luteostrians as Revealed by Genetic Mapping and Genomic Re-sequencing

Implementation of next-generation sequencing in forward genetic screens greatly accelerated gene discovery in species with larger genomes, including many crop plants. In barley, extensive mutant collections are available, however, the causative mutations for many of the genes remains largely unknown. Here we demonstrate how a combination of low-resolution genetic mapping, whole-genome resequencing and comparative functional analyses provides a promising path toward candidate identification of genes involved in plastid biology and/or photosynthesis, even if genes are located in recombination poor regions of the genome. As a proof of concept, we simulated the prediction of a candidate gene for the recently cloned variegation mutant albostrians (HvAST/HvCMF7) and adopted the approach for suggesting HvClpC1 as candidate gene for the yellow-green variegation mutant luteostrians.

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.

Genomic diversity analysis and identification of novel SSR markers in four tobacco varieties by high-throughput resequencing

Genome resequencing was carried out on two varieties of flue-cured tobacco (LY1306 and Qinyan 96), one variety of sun-cured tobacco (Wanmao 3), and one variety of air-cured Maryland tobacco (Wufeng 1), for a comparative analysis of genomic variation across the four varieties. Single nucleotide polymorphisms (SNPs), insertions and deletions (InDels), structural variations (SVs), and copy-number variations (CNVs) were then identified in each tobacco variety. Furthermore, a functional analysis of mutated genes was carried out. Through in-depth comparative analysis of genomes of different tobacco varieties, we identified genome variations in a number of SNPs, InDels, SVs, and CNVs, respectively. Computational analysis to predict the function of mutated genes containing these differential SNPs, InDels, SVs, and CNVs showed that they were mainly involved in different functions, such as carbohydrate metabolism and secondary metabolites biosynthesis. We mainly focused on genes that were involved in the biosynthesis of secondary metabolites and nicotine metabolism. In addition, we identified five simple sequence repeat (SSR)-based markers and verified them by PCR amplification in 10 tobacco varieties. Taken together, our study increases the understanding of genetic differences between tobacco types or varieties and identifies five SSR markers to classify tobacco varieties or types.

LMM24 Encodes Receptor-Like Cytoplasmic Kinase 109, Which Regulates Cell Death and Defense Responses in Rice

Lesion mimic mutants are excellent models for research on molecular mechanisms of cell death and defense responses in rice. We identified a new rice lesion mimic mutant lmm24 from a mutant pool of indica rice cultivar “ZhongHui8015”. The LMM24 gene was identified by MutMap, and LMM24 was confirmed as a receptor-like cytoplasmic kinase 109 by amino acid sequence analysis. The lmm24 mutant displayed dark brown lesions in leaves and growth retardation that were not observed in wild-type ZH8015. The results of histochemical staining and TUNEL assays showed enhanced ROS accumulation and cell death in lmm24. Chloroplast degradation was observed in lmm24 leaves, with decreased expression of photosynthesis-related genes and increased expression of the senescence-induced STAYGREEN (SGR) gene and other senescence-associated genes. Furthermore, lmm24 exhibited enhanced resistance to rice blast fungus Magnaporthe oryzae (M. oryzae) and up-regulation of defense response genes. Our data demonstrate that LMM24 regulates cell death and defense responses in rice.

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.

Genome-Wide Analysis of Genetic Variations and the Detection of Rich Variants of NBS-LRR Encoding Genes in Common Wild Rice Lines

Common wild rice (Oryza rufipogon Griff.) is invaluable genetic resource for rice resistance breeding. Whole-genome re-sequencing was conducted to systematically analyze the variations in two new inbred lines (Huaye 3 and Huaye 4) developed from a common wild rice. A total of 4,841,127 SNPs, 1,170,479 InDels, 24,080 structural variations (SVs), and 298 copy number variations (CNVs) were identified in three materials. Approximately 16.24 and 5.64% of the total SNPs and InDels of Huaye 3 and Huaye 4 were located in genic regions, respectively. Together, 12,486 and 15,925 large-effect SNPs, and 12,417 and 14,513 large-effect InDels, which affect the integrity of the encoded protein, were identified in Huaye 3 and Huaye 4, respectively. The distribution map of 194 and 245 NBS-LRR encoding homologs was constructed across 12 rice chromosomes. Further, GO enrichment analysis of the homologs with identical genotype variations in Huaye 3 and Huaye 4 revealed 67, 82, and 58 homologs involved in cell death, response to stress, and both terms, respectively. Comparative analysis displayed that 550 out of 652 SNPs and 129 out of 147 InDels were present in a widely used blast-susceptible rice variety (LTH). Protein-protein interaction analysis revealed a strong interaction between NBS-LRR candidates and several known R genes. One homolog of disease resistance protein (RPM1) was involved in the plant-pathogen interaction pathway. Artificial inoculation of disease/insect displayed resistance phenotypes against rice blast and brown planthopper in two lines. The results will provide allele-specific markers for rice molecular breeding.

A re-sequencing-based ultra-dense genetic map reveals a gummy stem blight resistance-associated gene in Cucumis melo

The melon (Cucumis melo) genome and genetic maps with hundreds to thousands of single nucleotide polymorphism markers were recently released. However, a high-resolution genetic map was lacking. Gummy stem blight (Gsb) is a destructive disease responsible for considerable economic losses during melon production. We herein describe the development of an ultra-dense genetic map consisting of 12,932 recombination bin markers covering 1,818 cM, with an average distance of 0.17 cM between adjacent tags. A comparison of the genetic maps for melon, watermelon, and cucumber revealed chromosome-level syntenic relationships and recombination events among the three Cucurbitaceae species. Our genetic map was useful for re-anchoring the genome scaffolds of melon. More than 92% assembly was anchored to 12 pseudo-chromosomes and 90% of them were oriented. Furthermore, 1,135 recombination hotspots revealed an unbalanced recombination rate across the melon genome. Genetic analyses of the Gsb-resistant and -susceptible lines indicated the resistance phenotype is mediated by a single dominant gene. We identified Gsb-resistance gene candidates in a 108-kb region on pseudo-chromosome 4. Our findings verify the utility of an ultra-dense genetic map for mapping a gene of interest, and for identifying new disease resistant genes.

Rapid identification of lettuce seed germination mutants by bulked segregant analysis and whole genome sequencing

Lettuce (Lactuca sativa) seeds exhibit thermoinhibition, or failure to complete germination when imbibed at warm temperatures. Chemical mutagenesis was employed to develop lettuce lines that exhibit germination thermotolerance. Two independent thermotolerant lettuce seed mutant lines, TG01 and TG10, were generated through ethyl methanesulfonate mutagenesis. Genetic and physiological analyses indicated that these two mutations were allelic and recessive. To identify the causal gene(s), we applied bulked segregant analysis by whole genome sequencing. For each mutant, bulked DNA samples of segregating thermotolerant (mutant) seeds were sequenced and analyzed for homozygous single-nucleotide polymorphisms. Two independent candidate mutations were identified at different physical positions in the zeaxanthin epoxidase gene (ABSCISIC ACID DEFICIENT 1/ZEAXANTHIN EPOXIDASE, or ABA1/ZEP) in TG01 and TG10. The mutation in TG01 caused an amino acid replacement, whereas the mutation in TG10 resulted in alternative mRNA splicing. Endogenous abscisic acid contents were reduced in both mutants, and expression of the ABA1 gene from wild-type lettuce under its own promoter fully complemented the TG01 mutant. Conventional genetic mapping confirmed that the causal mutations were located near the ZEP/ABA1 gene, but the bulked segregant whole genome sequencing approach more efficiently identified the specific gene responsible for the phenotype.

Genomics-based strategies for the use of natural variation in the improvement of crop metabolism

Next-generation genomics holds great potential in the study of plant phenotypic variation. With several crop reference genomes now available, the affordable costs of de novo genome assembly or target resequencing offer the opportunity to mine the enormous amount of genetic diversity hidden in crop wild relatives. Wide introgressions from these wild ancestors species or land races represent a possible strategy to improve cultivated varieties. In this review, we discuss the mechanisms underlying metabolic diversity within plant species and the possible strategies (and barriers) to introgress novel metabolic traits into cultivated varieties. We show how deep genomic surveys uncover various types of structural variants from extended gene pools of major crops and highlight how this variation may be used for the improvement of crop metabolism.

Bioinformatics in Africa: The Rise of Ghana?

Until recently, bioinformatics, an important discipline in the biological sciences, was largely limited to countries with advanced scientific resources. Nonetheless, several developing countries have lately been making progress in bioinformatics training and applications. In Africa, leading countries in the discipline include South Africa, Nigeria, and Kenya. However, one country that is less known when it comes to bioinformatics is Ghana. Here, I provide a first description of the development of bioinformatics activities in Ghana and how these activities contribute to the overall development of the discipline in Africa. Over the past decade, scientists in Ghana have been involved in publications incorporating bioinformatics analyses, aimed at addressing research questions in biomedical science and agriculture. Scarce research funding and inadequate training opportunities are some of the challenges that need to be addressed for Ghanaian scientists to continue developing their expertise in bioinformatics.