The art and design of genetic screens: Arabidopsis thaliana

The intron binding protein EMB-4 is an opposite regulator of cold and high temperature tolerance in Caenorhabditis elegans

Adaptation and tolerance to changes in heat and cold temperature are essential for survival and proliferation in plants and animals. However, there is no clear information regarding the common molecules between animals and plants. In this study, we found that heat, and cold tolerance of the nematode Caenorhabditis elegans is oppositely regulated by the RNA-binding protein EMB-4, whose plant homolog contains polymorphism causing heat tolerance diversity. Caenorhabditis elegans alters its cold and heat tolerance depending on the previous cultivation temperature, wherein EMB-4 respectively acts as a positive and negative controller of heat and cold tolerance by altering gene expression. Among the genes whose expression is regulated by EMB-4, a phospholipid scramblase, and an acid sphingomyelinase, which are involved in membrane lipid metabolism, were found to play essential roles in the negative regulation of heat tolerance.

Transposase-assisted tagmentation: an economical and scalable strategy for single-worm whole-genome sequencing

AlphaMissense identifies 23 million human missense variants as likely pathogenic, but only 0.1% have been clinically classified. To experimentally validate these predictions, chemical mutagenesis presents a rapid, cost-effective method to produce billions of mutations in model organisms. However, the prohibitive costs and limitations in the throughput of whole-genome sequencing (WGS) technologies, crucial for variant identification, constrain its widespread application. Here, we introduce a Tn5 transposase-assisted tagmentation technique for conducting WGS in Caenorhabditis elegans, Escherichia coli, Saccharomyces cerevisiae, and Chlamydomonas reinhardtii. This method, demands merely 20 min of hands-on time for a single-worm or single-cell clones and incurs a cost below 10 US dollars. It effectively pinpoints causal mutations in mutants defective in cilia or neurotransmitter secretion and in mutants synthetically sterile with a variant analogous to the B-Raf Proto-oncogene, Serine/Threonine Kinase (BRAF) V600E mutation. Integrated with chemical mutagenesis, our approach can generate and identify missense variants economically and efficiently, facilitating experimental investigations of missense variants in diverse species.

Dynamics of ER stress-induced gene regulation in plants

Endoplasmic reticulum (ER) stress is a potentially lethal condition that is induced by the abnormal accumulation of unfolded or misfolded secretory proteins in the ER. In eukaryotes, ER stress is managed by the unfolded protein response (UPR) through a tightly regulated, yet highly dynamic, reprogramming of gene transcription. Although the core principles of the UPR are similar across eukaryotes, unique features of the plant UPR reflect the adaptability of plants to their ever-changing environments and the need to balance the demands of growth and development with the response to environmental stressors. The past decades have seen notable progress in understanding the mechanisms underlying ER stress sensing and signalling transduction pathways, implicating the UPR in the effects of physiological and induced ER stress on plant growth and crop yield. Facilitated by sequencing technologies and advances in genetic and genomic resources, recent efforts have driven the discovery of transcriptional regulators and elucidated the mechanisms that mediate the dynamic and precise gene regulation in response to ER stress at the systems level.

ZmPILS6 is an auxin efflux carrier required for maize root morphogenesis

Plant root systems play a pivotal role in plant physiology and exhibit diverse phenotypic traits. Understanding the genetic mechanisms governing root growth and development in model plants like maize is crucial for enhancing crop resilience to drought and nutrient limitations. This study focused on identifying and characterizing ZmPILS6, an annotated auxin efflux carrier, as a key regulator of various crown root traits in maize. ZmPILS6-modified roots displayed reduced network area and suppressed lateral root formation, which are desirable traits for the "steep, cheap, and deep" ideotype. The research revealed that ZmPILS6 localizes to the endoplasmic reticulum and plays a vital role in controlling the spatial distribution of indole-3-acetic acid (IAA or "auxin") in primary roots. The study also demonstrated that ZmPILS6 can actively efflux IAA when expressed in yeast. Furthermore, the loss of ZmPILS6 resulted in significant proteome remodeling in maize roots, particularly affecting hormone signaling pathways. To identify potential interacting partners of ZmPILS6, a weighted gene coexpression analysis was performed. Altogether, this research contributes to the growing knowledge of essential genetic determinants governing maize root morphogenesis, which is crucial for guiding agricultural improvement strategies.

A K homology (KH) domain protein identified by a forward genetic screen affects bundle sheath anatomy in Arabidopsis thaliana

Because of their photosynthetic capacity, leaves function as solar panels providing the basis for the growth of the entire plant. Although the molecular mechanisms of leaf development have been well studied in model dicot and monocot species, a lot of information is still needed about the interplay of the genes that regulate cell division and differentiation and thereby affect the photosynthetic performance of the leaf. We were specifically interested in understanding the differentiation of mesophyll and bundle sheath cells in Arabidopsis thaliana and aimed to identify genes that are involved in determining bundle sheath anatomy. To this end, we established a forward genetic screen by using ethyl methanesulfonate (EMS) for mutagenizing a reporter line expressing a chloroplast-targeted green fluorescent protein (sGFP) under the control of a bundle sheath-specific promoter. Based on the GFP fluorescence phenotype, numerous mutants were produced, and by pursuing a mapping-by-sequencing approach, the genomic segments containing mutated candidate genes were identified. One of the lines with an enhanced GFP fluorescence phenotype (named ELEVATED BUNDLE SHEATH CELLS SIGNAL 1 [ebss1]) was selected for further study, and the responsible gene was verified by CRISPR/Cas9-based mutagenesis of candidate genes located in the mapped genomic segment. The verified gene, At2g25970, encodes a K homology (KH) domain-containing protein.

The unequal functional redundancy of Arabidopsis INCURVATA11 and CUPULIFORMIS2 is not dependent on genetic background

The paralogous genes INCURVATA11 (ICU11) and CUPULIFORMIS2 (CP2) encode components of the epigenetic machinery in Arabidopsis and belong to the 2-oxoglutarate and Fe (II)-dependent dioxygenase superfamily. We previously inferred unequal functional redundancy between ICU11 and CP2 from a study of the synergistic phenotypes of the double mutant and sesquimutant combinations of icu11 and cp2 mutations, although they represented mixed genetic backgrounds. To avoid potential confounding effects arising from different genetic backgrounds, we generated the icu11-5 and icu11-6 mutants via CRISPR/Cas genome editing in the Col-0 background and crossed them to cp2 mutants in Col-0. The resulting mutants exhibited a postembryonic-lethal phenotype reminiscent of strong embryonic flower (emf) mutants. Double mutants involving icu11-5 and mutations affecting epigenetic machinery components displayed synergistic phenotypes, whereas cp2-3 did not besides icu11-5. Our results confirmed the unequal functional redundancy between ICU11 and CP2 and demonstrated that it is not allele or genetic background specific. An increase in sucrose content in the culture medium partially rescued the post-germinative lethality of icu11 cp2 double mutants and sesquimutants, facilitating the study of their morphological phenotypes throughout their life cycle, which include floral organ homeotic transformations. We thus established that the ICU11-CP2 module is required for proper flower organ identity.

The Ka /Ks and πa /πs Ratios under Different Models of Gametophytic and Sporophytic Selection

Alternation of generations in plant life cycle provides a biological basis for natural selection occurring in either the gametophyte or the sporophyte phase or in both. Divergent biphasic selection could yield distinct evolutionary rates for phase-specific or pleiotropic genes. Here, we analyze models that deal with antagonistic and synergistic selection between alternative generations in terms of the ratio of nonsynonymous to synonymous divergence (Ka/Ks). Effects of biphasic selection are opposite under antagonistic selection but cumulative under synergistic selection for pleiotropic genes. Under the additive and comparable strengths of biphasic allelic selection, the absolute Ka/Ks for the gametophyte gene is equal to in outcrossing but smaller than, in a mixed mating system, that for the sporophyte gene under antagonistic selection. The same pattern is predicted for Ka/Ks under synergistic selection. Selfing reduces efficacy of gametophytic selection. Other processes, including pollen and seed flow and genetic drift, reduce selection efficacy. The polymorphism (πa) at a nonsynonymous site is affected by the joint effects of selfing with gametophytic or sporophytic selection. Likewise, the ratio of nonsynonymous to synonymous polymorphism (πa/πs) is also affected by the same joint effects. Gene flow and genetic drift have opposite effects on πa or πa/πs in interacting with gametophytic and sporophytic selection. We discuss implications of this theory for detecting natural selection in terms of Ka/Ks and for interpreting the evolutionary divergence among gametophyte-specific, sporophyte-specific, and pleiotropic genes.

Current insights and advances into plant male sterility: new precision breeding technology based on genome editing applications

Plant male sterility (MS) represents the inability of the plant to generate functional anthers, pollen, or male gametes. Developing MS lines represents one of the most important challenges in plant breeding programs, since the establishment of MS lines is a major goal in F1 hybrid production. For these reasons, MS lines have been developed in several species of economic interest, particularly in horticultural crops and ornamental plants. Over the years, MS has been accomplished through many different techniques ranging from approaches based on cross-mediated conventional breeding methods, to advanced devices based on knowledge of genetics and genomics to the most advanced molecular technologies based on genome editing (GE). GE methods, in particular gene knockout mediated by CRISPR/Cas-related tools, have resulted in flexible and successful strategic ideas used to alter the function of key genes, regulating numerous biological processes including MS. These precision breeding technologies are less time-consuming and can accelerate the creation of new genetic variability with the accumulation of favorable alleles, able to dramatically change the biological process and resulting in a potential efficiency of cultivar development bypassing sexual crosses. The main goal of this manuscript is to provide a general overview of insights and advances into plant male sterility, focusing the attention on the recent new breeding GE-based applications capable of inducing MS by targeting specific nuclear genic loci. A summary of the mechanisms underlying the recent CRISPR technology and relative success applications are described for the main crop and ornamental species. The future challenges and new potential applications of CRISPR/Cas systems in MS mutant production and other potential opportunities will be discussed, as generating CRISPR-edited DNA-free by transient transformation system and transgenerational gene editing for introducing desirable alleles and for precision breeding strategies.

Spatial Lipidomics of EPSPS and PAT Transgenic and Non-Transgenic Soybean Seeds Using Matrix-Assisted Laser Desorption/Ionization Mass Spectrometry Imaging

Herbicide-resistant soybeans are among the most widely planted transgenic crops. The in situ evaluation of spatial lipidomics in transgenic and non-transgenic soybeans is important for directly assessing the unintended effects of exogenous gene introduction. In this study, matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI)-based non-targeted analytical strategies were used for the first time for in situ detection and imaging of endogenous lipid distributions in transgenic (EPSPS and PAT genes) herbicide-resistant soybean (Glycine max Merrill) (S4003.14) and non-transgenic soybean (JACK) seeds. Statistical analysis revealed significant differences in lipids between S4003.14 and JACK seeds. The variable importance of projection analysis further revealed that 18 identified lipids, including six phosphatidylcholines (PCs), four phosphatidylethanolamines (PEs), five triacylglycerols (TAGs), and three cytidine diphosphate-diacylglycerols (CDP-DAGs), had the strongest differential expression between S4003.14 and JACK seeds. Among those, the upregulated expressions of PC(P-36:1), PC(36:2), PC(P-36:0), PC(37:5), PE(40:2), TAG(52:1), TAG(55:5), and CDP-DAG(37:2) and the downregulated expressions of PC (36:1), TAG(43:0), and three PEs (i.e., PE(P-38:1), PE(P-38:0), and PE(P-40:3)) were successfully found in the S4003.14 seeds, compared to these lipids detected in the JACK seeds. Meanwhile, the lipids of PC (44:8), CDP-DAG(38:0), and CDP-DAG(42:0) were uniquely detected in the S4003.14 soybean seeds, and TAG(45:2) and TAG(57:10) were detected as the unique lipids in the JACK seeds. The heterogeneous distribution of these lipids in the soybean seeds was also clearly visualized using MALDI-MSI. MSI results showed that lipid expression was significantly up/downregulated in S4003.14 seeds, compared to that in JACK seeds. This study improves our understanding of the unintended effects of herbicide-resistant EPSPS and PAT gene transfers on spatial lipidomes in soybean seeds and enables the continued progression of MALDI-MSI as an emerging, reliable, and rapid molecular imaging tool for evaluating unintended effects in transgenic plants.

Genetic Screens for Floral Mutants in Arabidopsis thaliana: Enhancers and Suppressors

The flower is a hallmark feature that has contributed to the evolutionary success of land plants. Diverse mutagenic agents have been employed as a tool to genetically perturb flower development and identify genes involved in floral patterning and morphogenesis. Since the initial studies to identify genes governing processes such as floral organ specification, mutagenesis in sensitized backgrounds has been used to isolate enhancers and suppressors to further probe the molecular basis of floral development. Here, we first describe two commonly employed methods for mutagenesis (using ethyl methanesulfonate (EMS) or T-DNAs as mutagens), and then describe three methods for identifying a mutation that leads to phenotypic alterations: traditional map-based cloning, modified high-efficiency thermal asymmetric interlaced PCR (mhiTAIL-PCR), and deep sequencing in the plant model Arabidopsis thaliana.

FIND-IT: Accelerated trait development for a green evolution

Improved agricultural and industrial production organisms are required to meet the future global food demands and minimize the effects of climate change. A new resource for crop and microbe improvement, designated FIND-IT (Fast Identification of Nucleotide variants by droplet DigITal PCR), provides ultrafast identification and isolation of predetermined, targeted genetic variants in a screening cycle of less than 10 days. Using large-scale sample pooling in combination with droplet digital PCR (ddPCR) greatly increases the size of low-mutation density and screenable variant libraries and the probability of identifying the variant of interest. The method is validated by screening variant libraries totaling 500,000 barley (Hordeum vulgare) individuals and isolating more than 125 targeted barley gene knockout lines and miRNA or promoter variants enabling functional gene analysis. FIND-IT variants are directly applicable to elite breeding pipelines and minimize time-consuming technical steps to accelerate the evolution of germplasm.

The need for unbiased genetic screens to dissect aggression in Drosophila melanogaster

Aggression is an evolutionarily conserved behavior present in most animals and is necessary for survival when competing for limited resources and mating partners. Studies have shown that aggression is modulated both genetically and epigenetically, but details of how the molecular and cellular mechanisms interact to determine aggressive behavior remain to be elucidated. In recent decades, Drosophila melanogaster has emerged as a powerful model system to understand the mechanisms that regulate aggression. Surprisingly most of the findings discovered to date have not come from genetic screens despite the fly's long and successful history of using screens to unravel its biology. Here, we highlight the tools and techniques used to successfully screen for aggression-linked behavioral elements in Drosophila and discuss the potential impact future screens have in advancing our knowledge of the underlying genetic and neural circuits governing aggression.

Induced variations of ethyl methane sulfonate mutagenized cowpea (Vigna unguiculata L. walp) plants

Unique variants are desired in the development of genetically improved crops to meet farmer and market needs hence ethyl methane sulfonate (EMS) was used to induce genetic variability in cowpea (Vigna unguiculata cv. Asontem). The main objective of this research was to characterize induced variations in EMS chemically mutagenized population of cowpea (Vigna unguiculata L. Walp Var. Asontem) in the M₁ and M₂ generations. The optimum concentration (LD50) of EMS for generating the mutagenized population was determined by treating seeds with different concentrations of EMS (0.0, 0.2, 0.4, 0.6, and 0.8% v/v) and observing the germination count after 5 days of planting the seeds in Petri dishes. Three thousand cowpea seeds were treated with the 0.4% EMS to generate the M₁ and M₂ populations that were evaluated for agronomic and morphological traits with untreated seeds serving as control. Data analysis involved distribution of qualitative and quantitative traits. Germination was significantly reduced in the mutagenized population (17.8%) and compared with that of the wild type (61.6%). Percentage survival was significantly higher in wild type (98.38%) as compared with the M₁ population (78.46%). Percentage germination in the M₂ population (74.03%) was lower than the wild type (80%). A wide spectrum of agro-morphological abnormalities was observed in the M₂ population. Wide variations and uniquely different phenotypic classes were observed in leaf color, leaf shape, growth habit, plant pigmentation, twining tendency, pod curvature, seed shape, and seed coat color. M₂ individuals were widely distributed for days to flowering, number of pods per plant, number of seeds per pod, number of locules per pods, percentage seed set, pod length and number of seeds per plant. In conclusion, the EMS mutagenesis was effective in inducing the unique variations that will be useful for breeding and development of new farmer preferred varieties.

Developmental genetics with model organisms

In Darwin's and Mendel's times, researchers investigated a wealth of organisms, chosen to solve particular problems for which they seemed especially well suited. Later, a focus on a few organisms, which are accessible to systematic genetic investigations, resulted in larger repertoires of methods and applications in these few species. Genetic animal model organisms with large research communities are the nematode Caenorhabditis elegans, the fly Drosophila melanogaster, the zebrafish Danio rerio, and the mouse Mus musculus. Due to their specific strengths, these model organisms have their strongest impacts in rather different areas of biology. C. elegans is unbeatable in the analysis of cell-to-cell contacts by saturation mutagenesis, as worms can be grown very fast in very high numbers. In Drosophila, a rich pattern is generated in the embryo as well as in adults that is used to unravel the underlying mechanisms of morphogenesis. The transparent larvae of zebrafish are uniquely suited to study organ development in a vertebrate, and the superb versatility of reverse genetics in the mouse made it the model organism to study human physiology and diseases. The combination of these models allows the in-depth genetic analysis of many fundamental biological processes using a plethora of different methods, finally providing many specific approaches to combat human diseases. The plant model Arabidopsis thaliana provides an understanding of many aspects of plant biology that might ultimately be useful for breeding crops.

A Forward Genetic Approach to Identify Plasmodesmal Trafficking Regulators Based on Trichome Rescue

Plasmodesmata (PD) are channels in the walls of plant cells which enable cell-to-cell information transfer. This includes the selective transport of specific transcription factors that control cell fate during plant development. KNOTTED1 (KN1) homeobox (KNOX) family transcription factors that are essential for the maintenance and function of stem cells in shoot meristems use this trafficking pathway, but its mechanism is largely unknown. Here we describe a forward genetic approach to the identification of regulators of selective KN1 trafficking through PD, using a trichome rescue system that permits simple visual analysis in Arabidopsis leaves. A KN1 trafficking regulator identified in this approach had the capacity to regulate the transport not only of KN1 but also of another mobile regulatory protein, TRANSPARENT TESTA GLABRA1 (TTG1). Our system could be easily adapted to reveal the mechanism underlying the selective transport of additional mobile signals through PD.

SciApps: An Automated Platform for Processing and Distribution of Plant Genomics Data

SciApps is an open-source, web-based platform for processing, storing, visualizing, and distributing genomic data and analysis results. Built upon the Tapis (formerly Agave) platform, SciApps brings users TB-scale of data storage via CyVerse Data Store and over one million CPUs via the Extreme Science and Engineering Discovery Environment (XSEDE) resources at Texas Advanced Computing Center (TACC). SciApps provides users ways to chain individual jobs into automated and reproducible workflows in a distributed cloud and provides a management system for data, associated metadata, individual analysis jobs, and multi-step workflows. This chapter provides examples of how to (1) submitting, managing, constructing workflows, (2) using public workflows for Bulked Segregant Analysis (BSA), (3) constructing a Data Analysis Center (DAC), and Data Coordination Center (DCC) for the plant ENCODE project.

The genome-wide rate and spectrum of EMS-induced heritable mutations in the microcrustacean Daphnia: on the prospect of forward genetics

Forward genetic screening using the alkylating mutagen ethyl methanesulfonate (EMS) is an effective method for identifying phenotypic mutants of interest, which can be further genetically dissected to pinpoint the causal genetic mutations. An accurate estimate of the rate of EMS-induced heritable mutations is fundamental for determining the mutant sample size of a screening experiment that aims to saturate all the genes in a genome with mutations. This study examines the genome-wide EMS-induced heritable base-substitutions in three species of the freshwater microcrustacean Daphnia to help guide screening experiments. Our results show that the 10 mM EMS treatment induces base substitutions at an average rate of 1.17 × 10^-6/site/generation across the three species, whereas a significantly higher average mutation rate of 1.75 × 10^-6 occurs at 25 mM. The mutation spectrum of EMS-induced base substitutions at both concentration is dominated by G:C to A:T transitions. Furthermore, we find that female Daphnia exposed to EMS (F₀ individuals) can asexually produce unique mutant offspring (F₁) for at least 3 consecutive broods, suggestive of multiple broods as F₁ mutants. Lastly, we estimate that about 750 F₁s are needed for all genes in the Daphnia genome to be mutated at least once with a 95% probability. We also recommend 4-5 F₂s should be collected from each F₁ mutant through sibling crossing so that all induced mutations could appear in the homozygous state in the F₂ population at 70-80% probability.

mRNA surveillance complex PELOTA-HBS1 regulates phosphoinositide-dependent protein kinase1 and plant growth

The quality control system for messenger RNA (mRNA) is fundamental for cellular activities in eukaryotes. To elucidate the molecular mechanism of 3'-phosphoinositide-dependent protein kinase1 (PDK1), a master regulator that is essential throughout eukaryotic growth and development, we employed a forward genetic approach to screen for suppressors of the loss-of-function T-DNA insertion double mutant pdk1.1 pdk1.2 in Arabidopsis thaliana. Notably, the severe growth attenuation of pdk1.1 pdk1.2 was rescued by sop21 (suppressor of pdk1.1 pdk1.2), which harbors a loss-of-function mutation in PELOTA1 (PEL1). PEL1 is a homolog of mammalian PELOTA and yeast (Saccharomyces cerevisiae) DOM34p, which each form a heterodimeric complex with the GTPase HBS1 (HSP70 SUBFAMILY B SUPPRESSOR1, also called SUPERKILLER PROTEIN7, SKI7), a protein that is responsible for ribosomal rescue and thereby assures the quality and fidelity of mRNA molecules during translation. Genetic analysis further revealed that a dysfunctional PEL1-HBS1 complex failed to degrade the T-DNA-disrupted PDK1 transcripts, which were truncated but functional, and thus rescued the growth and developmental defects of pdk1.1 pdk1.2. Our studies demonstrated the functionality of a homologous PELOTA-HBS1 complex and identified its essential regulatory role in plants, providing insights into the mechanism of mRNA quality control.

slim shady is a novel allele of PHYTOCHROME B present in the T-DNA line SALK_015201

Auxin is a hormone that is required for hypocotyl elongation during seedling development. In response to auxin, rapid changes in transcript and protein abundance occur in hypocotyls, and some auxin responsive gene expression is linked to hypocotyl growth. To functionally validate proteomic studies, a reverse genetics screen was performed on mutants in auxin-regulated proteins to identify novel regulators of plant growth. This uncovered a long hypocotyl mutant, which we called slim shady, in an annotated insertion line in IMMUNOREGULATORY RNA-BINDING PROTEIN (IRR). Overexpression of the IRR gene failed to rescue the slim shady phenotype and characterization of a second T-DNA allele of IRR found that it had a wild-type (WT) hypocotyl length. The slim shady mutant has an elevated expression of numerous genes associated with the brassinosteroid-auxin-phytochrome (BAP) regulatory module compared to WT, including transcription factors that regulate brassinosteroid, auxin, and phytochrome pathways. Additionally, slim shady seedlings fail to exhibit a strong transcriptional response to auxin. Using whole genome sequence data and genetic complementation analysis with SALK_015201C, we determined that a novel single nucleotide polymorphism in PHYTOCHROME B was responsible for the slim shady phenotype. This is predicted to induce a frameshift and premature stop codon at leucine 1125, within the histidine kinase-related domain of the carboxy terminus of PHYB, which is required for phytochrome signaling and function. Genetic complementation analyses with phyb-9 confirmed that slim shady is a mutant allele of PHYB. This study advances our understanding of the molecular mechanisms in seedling development, by furthering our understanding of how light signaling is linked to auxin-dependent cell elongation. Furthermore, this study highlights the importance of confirming the genetic identity of research material before attributing phenotypes to known mutations sourced from T-DNA stocks.

CRISPR screens in plants: approaches, guidelines, and future prospects

Clustered regularly interspaced short palindromic repeat (CRISPR)-associated systems have revolutionized genome engineering by facilitating a wide range of targeted DNA perturbations. These systems have resulted in the development of powerful new screens to test gene functions at the genomic scale. While there is tremendous potential to map and interrogate gene regulatory networks at unprecedented speed and scale using CRISPR screens, their implementation in plants remains in its infancy. Here we discuss the general concepts, tools, and workflows for establishing CRISPR screens in plants and analyze the handful of recent reports describing the use of this strategy to generate mutant knockout collections or to diversify DNA sequences. In addition, we provide insight into how to design CRISPR knockout screens in plants given the current challenges and limitations and examine multiple design options. Finally, we discuss the unique multiplexing capabilities of CRISPR screens to investigate redundant gene functions in highly duplicated plant genomes. Combinatorial mutant screens have the potential to routinely generate higher-order mutant collections and facilitate the characterization of gene networks. By integrating this approach with the numerous genomic profiles that have been generated over the past two decades, the implementation of CRISPR screens offers new opportunities to analyze plant genomes at deeper resolution and will lead to great advances in functional and synthetic biology.

Enhancer-mediated reporter gene expression in Arabidopsis thaliana: a forward genetic screen

Gene expression is controlled and regulated by interactions between cis-regulatory DNA elements (CREs) and regulatory proteins. Enhancers are one of the most important classes of CREs in eukaryotes. Eukaryotic genes, especially those related to development or responses to environmental cues, are often regulated by multiple enhancers in different tissues and/or at different developmental stages. Remarkably, little is known about the molecular mechanisms by which enhancers regulate gene expression in plants. We identified a distal enhancer, CREβ, which regulates the expression of AtDGK7, which encodes a diacylglycerol kinase in Arabidopsis. We developed a transgenic line containing the luciferase reporter gene (LUC) driven by CREβ fused with a minimal cauliflower mosaic virus (CaMV) 35S promoter. The CREβ enhancer was shown to play a role in the response to osmotic pressure of the LUC reporter gene. A forward genetic screen pipeline based on the transgenic line was established to generate mutations associated with altered expression of the LUC reporter gene. We identified a suite of mutants with variable LUC expression levels as well as different segregation patterns of the mutations in populations. We demonstrate that this pipeline will allow us to identify trans-regulatory factors associated with CREβ function as well as those acting in the regulation of the endogenous AtDGK7 gene.

BSAseq: an interactive and integrated web-based workflow for identification of causal mutations in bulked F2 populations

SUMMARY

With the advance of next-generation sequencing technologies and reductions in the costs of these techniques, bulked segregant analysis (BSA) has become not only a powerful tool for mapping quantitative trait loci but also a useful way to identify causal gene mutations underlying phenotypes of interest. However, due to the presence of background mutations and errors in sequencing, genotyping, and reference assembly, it is often difficult to distinguish true causal mutations from background mutations. In this study, we developed the BSAseq workflow, which includes an automated bioinformatics analysis pipeline with a probabilistic model for estimating the linked region (the region linked to the causal mutation) and an interactive Shiny web application for visualizing the results. We deeply sequenced a sorghum male-sterile parental line (ms8) to capture the majority of background mutations in our bulked F2 data. We applied the workflow to 11 bulked sorghum F2 populations and 1 rice F2 population and identified the true causal mutation in each population. The workflow is intuitive and straightforward, facilitating its adoption by users without bioinformatics analysis skills. We anticipate that the BSAseq workflow will be broadly applicable to the identification of causal mutations for many phenotypes of interest.

AVAILABILITY AND IMPLEMENTATION

BSAseq is freely available on https://www.sciapps.org/page/bsa.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Can genetic engineering-based methods for gene function identification be eclipsed by genome editing in plants? A comparison of methodologies

Finding and explaining the functions of genes in plants have promising applications in crop improvement and bioprospecting and hence, it is important to compare various techniques available for gene function identification in plants. Today, the most popular technology among researchers to identify the functions of genes is the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 (Cas9)-based genome editing method. But by no means can we say that CRISPR/Cas9 is the go-to method for all purposes. It comes with its own baggage. Researchers will agree and have lived through at least seven more technologies deployed to find the functions of genes, which come under three umbrellas: 1. genetic engineering, 2. transient expression, and 3. chemical/physical mutagenesis. Each of the methods evolved when the previous one ran into an insurmountable problem. In this review, we compare the eight technologies against one another on 14 parameters. This review lays bare the pros and cons, and similarities and dissimilarities of various methods. Every method comes with its advantages and disadvantages. For example, the CRISPR/Cas9-based genome editing is an excellent method for modifying gene sequences, creating allelic versions of genes, thereby aiding the understanding of gene function. But it comes with the baggage of unwanted or off-target mutations. Then, we have methods based on random or targeted knockout of the gene, knockdown, and overexpression of the gene. Targeted disruption of genes is required for complete knockout of gene function, which may not be accomplished by editing. We have also discussed the strategies to overcome the shortcomings of the targeted gene-knockout and the CRISPR/Cas9-based methods. This review serves as a comprehensive guide towards the understanding and comparison of various technologies available for gene function identification in plants and hence, it will find application for crop improvement and bioprospecting related research.

Distinguishing HapMap Accessions Through Recursive Set Partitioning in Hierarchical Decision Trees

The HapMap (haplotype map) projects have produced valuable genetic resources in life science research communities, allowing researchers to investigate sequence variations and conduct genome-wide association study (GWAS) analyses. A typical HapMap project may require sequencing hundreds, even thousands, of individual lines or accessions within a species. Due to limitations in current sequencing technology, the genotype values for some accessions cannot be clearly called. Additionally, allelic heterozygosity can be very high in some lines, causing genetic and sometimes phenotypic segregation in their descendants. Genetic and phenotypic segregation degrades the original accession's specificity and makes it difficult to distinguish one accession from another. Therefore, it is vitally important to determine and validate HapMap accessions before one conducts a GWAS analysis. However, to the best of our knowledge, there are no prior methodologies or tools that can readily distinguish or validate multiple accessions in a HapMap population. We devised a bioinformatics approach to distinguish multiple HapMap accessions using only a minimum number of genetic markers. First, we assign each candidate marker with a distinguishing score (DS), which measures its capability in distinguishing accessions. The DS score prioritizes those markers with higher percentages of homozygous genotypes (allele combinations), as they can be stably passed on to offspring. Next, we apply the "set-partitioning" concept to select optimal markers by recursively partitioning accession sets. Subsequently, we build a hierarchical decision tree in which a specific path represents the selected markers and the homogenous genotypes that can be used to distinguish one accession from others in the HapMap population. Based on these algorithms, we developed a web tool named MAD-HiDTree (Multiple Accession Distinguishment-Hierarchical Decision Tree), designed to analyze a user-input genotype matrix and construct a hierarchical decision tree. Using genetic marker data extracted from the Medicago truncatula HapMap population, we successfully constructed hierarchical decision trees by which the original 262 M. truncatula accessions could be efficiently distinguished. PCR experiments verified our proposed method, confirming that MAD-HiDTree can be used for the identification of a specific accession. MAD-HiDTree was developed in C/C⁺⁺ in Linux. Both the source code and test data are publicly available at https://bioinfo.noble.org/MAD-HiDTree/.

Overexpression of the Auxin Receptor AFB3 in Arabidopsis Results in Salt Stress Resistance and the Modulation of NAC4 and SZF1

Soil salinity is a key problem for crop production worldwide. High salt concentration in soil negatively modulates plant growth and development. In roots, salinity affects the growth and development of both primary and lateral roots. The phytohormone auxin regulates various developmental processes during the plant’s life cycle, including several aspects of root architecture. Auxin signaling involves the perception by specialized receptors which module several regulatory pathways. Despite their redundancy, previous studies have shown that their functions can also be context-specific depending on tissue, developmental or environmental cues. Here we show that the over-expression of Auxin Signaling F-Box 3 receptor results in an increased resistance to salinity in terms of root architecture and germination. We also studied possible downstream signaling components to further characterize the role of auxin in response to salt stress. We identify the transcription factor SZF1 as a key component in auxin-dependent salt stress response through the regulation of NAC4. These results give lights of an auxin-dependent mechanism that leads to the modulation of root system architecture in response to salt identifying a hormonal cascade important for stress response.

Whole-Seedling-Based Chemical Genetic Screens in Arabidopsis

Forward genetics has been extremely powerful for dissecting biological pathways in various model organisms. However, it is limited by the fact that redundant gene families and essential genes cannot be readily uncovered through such methods. Chemical genetics, on the other hand, provides a valuable complementary approach to probe biological processes and is suitable for not only genetic model organisms but also genetically less tractable species. We describe here a high-throughput chemical genetic screening method simply based on plant growth and developmental phenotypes in Arabidopsis. It was successfully utilized to study plant immunity and can be easily adapted for dissecting other plant signal transduction pathways.

A Genetic Screen to Identify New Molecular Players Involved in Photoprotection qH in Arabidopsis thaliana

Photosynthesis is a biological process which converts light energy into chemical energy that is used in the Calvin-Benson cycle to produce organic compounds. An excess of light can induce damage to the photosynthetic machinery. Therefore, plants have evolved photoprotective mechanisms such as non-photochemical quenching (NPQ). To focus molecular insights on slowly relaxing NPQ processes in Arabidopsis thaliana, previously, a qE-deficient line-the PsbS mutant-was mutagenized and a mutant with high and slowly relaxing NPQ was isolated. The mutated gene was named suppressor of quenching 1, or SOQ1, to describe its function. Indeed, when present, SOQ1 negatively regulates or suppresses a form of antenna NPQ that is slow to relax and is photoprotective. We have now termed this component qH and identified the plastid lipocalin, LCNP, as the effector for this energy dissipation mode to occur. Recently, we found that the relaxation of qH1, ROQH1, protein is required to turn off qH. The aim of this study is to identify new molecular players involved in photoprotection qH by a whole genome sequencing approach of chemically mutagenized Arabidopsis thaliana. We conducted an EMS-mutagenesis on the soq1 npq4 double mutant and used chlorophyll fluorescence imaging to screen for suppressors and enhancers of qH. Out of 22,000 mutagenized plants screened, the molecular players cited above were found using a mapping-by-sequencing approach. Here, we describe the phenotypic characterization of the other mutants isolated from this genetic screen and an additional 8000 plants screened. We have classified them in several classes based on their fluorescence parameters, NPQ kinetics, and pigment content. A high-throughput whole genome sequencing approach on 65 mutants will identify the causal mutations thanks to allelic mutations from having reached saturation of the genetic screen. The candidate genes could be involved in the formation or maintenance of quenching sites for qH, in the regulation of qH at the transcriptional level, or be part of the quenching site itself.

Identification and characterization of unconventional membrane protein trafficking regulators in Arabidopsis: A genetic approach

Proper trafficking and subcellular localization of membrane proteins are essential for plant growth and development. The plant endomembrane system contains several membrane-bound organelles with distinct functions including the endoplasmic reticulum (ER), Golgi apparatus, trans-Golgi network (TGN) or early endosome, prevacuolar compartment (PVC) or multivesicular body (MVB) and vacuole. Multiple approaches have been successfully used to identify and study the regulators and components important for signal transduction, growth and development, as well as membrane trafficking in the endomembrane system in plants. These include the homologous characterization of the counterparts in mammals or yeast employing both reverse genetic as well as the forward genetic screen approaches. However, the deletion or mutation of membrane trafficking related proteins usually leads to seedling lethality due to their essential roles in plant development and organelle biogenesis. To overcome the limitation of lethal phenotype of the target proteins, we used DEX-inducible RNAi knock-down lines to study their function in plants. More recently, we developed and used both RNAi knock-down and T-DNA insertional lines as starting materials to screen for mutations that could suppress and rescue the lethal phenotype, or a suppressor screening. Further characterization of the newly identified suppressor mutants has resulted in the identification of novel negative regulators in mediating membrane trafficking and organelle biogenesis in plants. In this review, we summarize the current approaches in studying protein trafficking in the endomembrane system. We then describe three examples of suppressor screening with distinct starting materials (i.e. FREE1, MON1, and SH3P2 that are regulators of MVB, vacuole, and autophagosomes, respectively) to discuss the rationale, procedures, advantages and disadvantages, and possible outcomes of such a suppressor screening. We finally propose that these novel screening approaches will lead to the identification of new unconventional players in regulating protein trafficking and organelle biogenesis in plants and discuss their impact on plant cell biology research.

Phenotyping in Arabidopsis and Crops-Are We Addressing the Same Traits? A Case Study in Tomato

The convenient model Arabidopsis thaliana has allowed tremendous advances in plant genetics and physiology, in spite of only being a weed. It has also unveiled the main molecular networks governing, among others, abiotic stress responses. Through the use of the latest genomic tools, Arabidopsis research is nowadays being translated to agronomically interesting crop models such as tomato, but at a lagging pace. Knowledge transfer has been hindered by invariable differences in plant architecture and behaviour, as well as the divergent direct objectives of research in Arabidopsis versus crops compromise transferability. In this sense, phenotype translation is still a very complex matter. Here, we point out the challenges of "translational phenotyping" in the case study of drought stress phenotyping in Arabidopsis and tomato. After briefly defining and describing drought stress and survival strategies, we compare drought stress protocols and phenotyping techniques most commonly used in the two species, and discuss their potential to gain insights, which are truly transferable between species. This review is intended to be a starting point for discussion about translational phenotyping approaches among plant scientists, and provides a useful compendium of methods and techniques used in modern phenotyping for this specific plant pair as a case study.

Photosynthesis in a Changing Global Climate: Scaling Up and Scaling Down in Crops

Photosynthesis is the major process leading to primary production in the Biosphere. There is a total of 7000bn tons of CO2 in the atmosphere and photosynthesis fixes more than 100bn tons annually. The CO2 assimilated by the photosynthetic apparatus is the basis of crop production and, therefore, of animal and human food. This has led to a renewed interest in photosynthesis as a target to increase plant production and there is now increasing evidence showing that the strategy of improving photosynthetic traits can increase plant yield. However, photosynthesis and the photosynthetic apparatus are both conditioned by environmental variables such as water availability, temperature, [CO2], salinity, and ozone. The "omics" revolution has allowed a better understanding of the genetic mechanisms regulating stress responses including the identification of genes and proteins involved in the regulation, acclimation, and adaptation of processes that impact photosynthesis. The development of novel non-destructive high-throughput phenotyping techniques has been important to monitor crop photosynthetic responses to changing environmental conditions. This wealth of data is being incorporated into new modeling algorithms to predict plant growth and development under specific environmental constraints. This review gives a multi-perspective description of the impact of changing environmental conditions on photosynthetic performance and consequently plant growth by briefly highlighting how major technological advances including omics, high-throughput photosynthetic measurements, metabolic engineering, and whole plant photosynthetic modeling have helped to improve our understanding of how the photosynthetic machinery can be modified by different abiotic stresses and thus impact crop production.

Development and Characterization of an Ethyl Methane Sulfonate (EMS) Induced Mutant Population in Capsicum annuum L

Plant breeding explores genetic diversity in useful traits to develop new, high-yielding, and improved cultivars. Ethyl methane sulfonate (EMS) is a chemical widely used to induce mutations at loci that regulate economically essential traits. Additionally, it can knock out genes, facilitating efforts to elucidate gene functions through the analysis of mutant phenotypes. Here, we developed a mutant population using the small and pungent ornamental Capsicum annuum pepper “Micro-Pep”. This accession is particularly suitable for mutation studies and molecular research due to its compact growth habit and small size. We treated 9500 seeds with 1.3% EMS and harvested 3996 M₂ lines. We then selected 1300 (32.5%) independent M₂ families and evaluated their phenotypes over four years. The mutants displayed phenotypic variations in plant growth, habit, leaf color and shape, and flower and fruit morphology. An experiment to optimize Targeting Induced Local Lesions IN Genomes (TILLING) in pepper detected nine EMS-induced mutations in the eIF4E gene. The M₂ families developed here exhibited broad phenotypic variation and should be valuable genetic resources for functional gene analysis in pepper molecular breeding programs using reverse genetics tools, including TILLING.

Can-Seq: a PCR and DNA sequencing strategy for identifying new alleles of known and candidate genes

BACKGROUND

Forward genetic screens are a powerful approach for identifying the genes contributing to a trait of interest. However, mutants arising in genes already known can obscure the identification of new genes contributing to the trait. Here, we describe a strategy called Candidate gene-Sequencing (Can-Seq) for rapidly identifying and filtering out mutants carrying new alleles of known and candidate genes.

RESULTS

We carried out a forward genetic screen and identified 40 independent Arabidopsis mutants with defects in systemic spreading of RNA interference (RNAi), or more specifically in root-to-shoot transmission of post-transcriptional gene silencing (rtp). To classify the mutants as either representing a new allele of a known or candidate gene versus carrying a mutation in an undiscovered gene, bulk genomic DNA from up to 23 independent mutants was used as template to amplify a collection of 47 known or candidate genes. These amplified sequences were combined into Can-Seq libraries and deep sequenced. Subsequently, mutations in the known and candidate genes were identified using a custom Snakemake script (https://github.com/Carroll-Lab/can_seq), and PCR zygosity tests were then designed and used to identify the individual mutants carrying each mutation. Using this approach, we showed that 28 of the 40 rtp mutants carried homozygous nonsense, missense or splice site mutations in one or more of the 47 known or candidate genes. We conducted complementation tests to demonstrate that several of the candidate mutations were responsible for the rtp defect. Importantly, by exclusion, the Can-Seq pipeline also identified rtp mutants that did not carry a causative mutation in any of the 47 known and candidate genes, and these mutants represent an undiscovered gene(s) required for systemic RNAi.

CONCLUSIONS

Can-Seq offers an accurate, cost-effective method for classifying new mutants into known versus unknown genes. It has several advantages over existing genetic and DNA sequencing approaches that are currently being used in forward genetic screens for gene discovery. Using Can-Seq in conjunction with map-based gene cloning is a cost-effective approach towards identifying the full complement of genes contributing to a trait of interest.

Pentatricopeptide repeat protein MID1 modulates nad2 intron 1 splicing and Arabidopsis development

As one of the best-studied RNA binding proteins in plant, pentatricopeptide repeats (PPRs) protein are mainly targeted to mitochondria and/or chloroplasts for RNA processing to regulate the biogenesis and function of the organelles, but its molecular mechanism and role in development remain to be further revealed. Here, we identified a mitochondria-localized P-type small PPR protein, MITOCHONDRION-MEDIATED GROWTH DEFECT 1 (MID1) that is crucial for Arabidopsis development. Mutation in MID1 causes retarded embryo development and stunted plant growth with defects in cell expansion and proliferation. Molecular experiments showed that MID1 is required for the splicing of the nad2 intron 1 in mitochondria. Consistently, mid1 plants display significant reduction in the abundance and activity of mitochondrial respiration complex I, accompanied by abnormal mitochondrial morphology and energy metabolism. Furthermore, MID1 is associated with other trans-factors involved in NICOTINAMIDE ADENINE DINUCLEOTIDE HYDROGEN (NADH) DEHYDROGENASE SUBUNIT 2 (nad2) intron 1 splicing, and interacts directly with itself and MITOCHONDRIAL STABILITY FACTOR 1 (MTSF1). This suggests that MID1 most likely functions as a dimer for nad2 intron 1 splicing. Together, we characterized a novel PPR protein MID1 for nad2 intron 1 splicing.

Genetic Screens to Target Embryo and Endosperm Pathways in Arabidopsis and Maize

The major tissue types and stem-cell niches of plants are established during embryogenesis, and thus knowledge of embryo development is essential for a full understanding of plant development. Studies of seed development are also important for human health, because the nutrients stored in both the embryo and endosperm of plant seeds provide an essential part of our diet. Arabidopsis and maize have evolved different types of seeds, opening a range of experimental opportunities. Development of the Arabidopsis embryo follows an almost invariant pattern, while cell division patterns of maize embryos are variable. Embryo-endosperm interactions are also different between the two species: in Arabidopsis, the endosperm is consumed during seed development, while mature maize seeds contain an enormous endosperm. Genetic screens have provided important insights into seed development in both species. In the genomic era, genetic analysis will continue to provide important tools for understanding embryo and endosperm biology in plants, because single gene functional studies can now be integrated with genome-wide information. Here, we lay out important factors to consider when designing genetic screens to identify new genes or to probe known pathways in seed development. We then highlight the technical details of two previous genetic screens that may serve as useful examples for future experiments.

CrossPlan: systematic planning of genetic crosses to validate mathematical models

Motivation

Mathematical models of cellular processes can systematically predict the phenotypes of novel combinations of multi-gene mutations. Searching for informative predictions and prioritizing them for experimental validation is challenging since the number of possible combinations grows exponentially in the number of mutations. Moreover, keeping track of the crosses needed to make new mutants and planning sequences of experiments is unmanageable when the experimenter is deluged by hundreds of potentially informative predictions to test.

Results

We present CrossPlan, a novel methodology for systematically planning genetic crosses to make a set of target mutants from a set of source mutants. We base our approach on a generic experimental workflow used in performing genetic crosses in budding yeast. We prove that the CrossPlan problem is NP-complete. We develop an integer-linear-program (ILP) to maximize the number of target mutants that we can make under certain experimental constraints. We apply our method to a comprehensive mathematical model of the protein regulatory network controlling cell division in budding yeast. We also extend our solution to incorporate other experimental conditions such as a delay factor that decides the availability of a mutant and genetic markers to confirm gene deletions. The experimental flow that underlies our work is quite generic and our ILP-based algorithm is easy to modify. Hence, our framework should be relevant in plant and animal systems as well.

Availability and implementation

CrossPlan code is freely available under GNU General Public Licence v3.0 at https://github.com/Murali-group/crossplan.

Supplementary information

Supplementary data are available at Bioinformatics online.

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

Reporter-based forward genetic screen to identify bundle sheath anatomy mutants in A. thaliana

The evolution of C₄ photosynthesis proceeded stepwise with each small step increasing the fitness of the plant. An important pre-condition for the introduction of a functional C₄ cycle is the photosynthetic activation of the C₃ bundle sheath by increasing its volume and organelle number. Therefore, to engineer C₄ photosynthesis into existing C₃ crops, information about genes that control the bundle sheath cell size and organelle content is needed. However, very little information is known about the genes that could be manipulated to create a more C₄ -like bundle sheath. To this end, an ethylmethanesulfonate (EMS)-based forward genetic screen was established in the Brassicaceae C₃  species Arabidopsis thaliana. To ensure a high-throughput primary screen, the bundle sheath cells of A. thaliana were labeled using a luciferase (LUC68) or by a chloroplast-targeted green fluorescent protein (sGFP) reporter using a bundle sheath specific promoter. The signal strengths of the reporter genes were used as a proxy to search for mutants with altered bundle sheath anatomy. Here, we show that our genetic screen predominantly identified mutants that were primarily affected in the architecture of the vascular bundle, and led to an increase in bundle sheath volume. By using a mapping-by-sequencing approach the genomic segments that contained mutated candidate genes were identified.

Sequencing-Based Bin Map Construction of a Tomato Mapping Population, Facilitating High-Resolution Quantitative Trait Loci Detection

Genotyping-by-sequencing (GBS) was employed to construct a highly saturated genetic linkage map of a tomato ( L.) recombinant inbred line (RIL) population, derived from a cross between cultivar NC EBR-1 and the wild tomato L. accession LA2093. A pipeline was developed to convert single nucleotide polymorphism (SNP) data into genomic bins, which could be used for fine mapping of quantitative trait loci (QTL) and identification of candidate genes. The pipeline, implemented in a python script named SNPbinner, adopts a hidden Markov model approach for calculation of recombination breakpoints followed by genomic bins construction. The total length of the newly developed high-resolution genetic map was 1.2-fold larger than previously estimated based on restriction fragment length polymorphism (RFLP) and polymerase chain reaction (PCR)-based markers. The map was used to verify and refine QTL previously identified for two fruit quality traits in the RIL population, fruit weight (FW) and fruit lycopene content (LYC). Two well-described FW QTL ( and ) were localized precisely at their known underlying causative genes, and the QTL intervals were decreased by two- to tenfold. A major QTL for LYC content () was verified at high resolution and its underlying causative gene was determined to be ζ (). The RIL population, the high resolution genetic map, and the easy-to-use genotyping pipeline, SNPbinner, are made publicly available.

Rapid fine mapping of causative mutations from sets of unordered, contig-sized fragments of genome sequence

BACKGROUND

Traditional Map based Cloning approaches, used for the identification of desirable alleles, are extremely labour intensive and years can elapse between the mutagenesis and the detection of the polymorphism. High throughput sequencing based Mapping-by-sequencing approach requires an ordered genome assembly and cannot be used with fragmented, un-scaffolded draft genomes, limiting its application to model species and precluding many important organisms.

RESULTS

We addressed this gap in resource and presented a computational method and software implementations called CHERIPIC (Computing Homozygosity Enriched Regions In genomes to Prioritise Identification of Candidate variants). We have successfully validated implementation of CHERIPIC using three different types of bulk segregant sequence data from Arabidopsis, maize and barley, respectively.

CONCLUSIONS

CHERIPIC allows users to rapidly analyse bulk segregant sequence data and we have made it available as a pre-packaged binary with all dependencies for Linux and MacOS and as Galaxy tool.

Deciphering genetic factors that determine melon fruit-quality traits using RNA-Seq-based high-resolution QTL and eQTL mapping

Combined quantitative trait loci (QTL) and expression-QTL (eQTL) mapping analysis was performed to identify genetic factors affecting melon (Cucumis melo) fruit quality, by linking genotypic, metabolic and transcriptomic data from a melon recombinant inbred line (RIL) population. RNA sequencing (RNA-Seq) of fruit from 96 RILs yielded a highly saturated collection of > 58 000 single-nucleotide polymorphisms, identifying 6636 recombination events that separated the genome into 3663 genomic bins. Bin-based QTL analysis of 79 RILs and 129 fruit-quality traits affecting taste, aroma and color resulted in the mapping of 241 QTL. Thiol acyltransferase (CmThAT1) gene was identified within the QTL interval of its product, S-methyl-thioacetate, a key component of melon fruit aroma. Metabolic activity of CmThAT1-encoded protein was validated in bacteria and in vitro. QTL analysis of flesh color intensity identified a candidate white-flesh gene (CmPPR1), one of two major loci determining fruit flesh color in melon. CmPPR1 encodes a member of the pentatricopeptide protein family, involved in processing of RNA in plastids, where carotenoid and chlorophyll pigments accumulate. Network analysis of > 12 000 eQTL mapped for > 8000 differentially expressed fruit genes supported the role of CmPPR1 in determining the expression level of plastid targeted genes. We highlight the potential of RNA-Seq-based QTL analysis of small to moderate size, advanced RIL populations for precise marker-assisted breeding and gene discovery. We provide the following resources: a RIL population genotyped with a unique set of SNP markers, confined genomic segments that harbor QTL governing 129 traits and a saturated set of melon eQTLs.

Whole-Genome Sequence Accuracy Is Improved by Replication in a Population of Mutagenized Sorghum

The accurate detection of induced mutations is critical for both forward and reverse genetics studies. Experimental chemical mutagenesis induces relatively few single base changes per individual. In a complex eukaryotic genome, false positive detection of mutations can occur at or above this mutagenesis rate. We demonstrate here, using a population of ethyl methanesulfonate (EMS)-treated Sorghum bicolor BTx623 individuals, that using replication to detect false positive-induced variants in next-generation sequencing (NGS) data permits higher throughput variant detection with greater accuracy. We used a lower sequence coverage depth (average of 7×) from 586 independently mutagenized individuals and detected 5,399,493 homozygous single nucleotide polymorphisms (SNPs). Of these, 76% originated from only 57,872 genomic positions prone to false positive variant calling. These positions are characterized by high copy number paralogs where the error-prone SNP positions are at copies containing a variant at the SNP position. The ability of short stretches of homology to generate these error-prone positions suggests that incompletely assembled or poorly mapped repeated sequences are one driver of these error-prone positions. Removal of these false positives left 1,275,872 homozygous and 477,531 heterozygous EMS-induced SNPs, which, congruent with the mutagenic mechanism of EMS, were >98% G:C to A:T transitions. Through this analysis, we generated a collection of sequence indexed mutants of sorghum. This collection contains 4035 high-impact homozygous mutations in 3637 genes and 56,514 homozygous missense mutations in 23,227 genes. Each line contains, on average, 2177 annotated homozygous SNPs per genome, including seven likely gene knockouts and 96 missense mutations. The number of mutations in a transcript was linearly correlated with the transcript length and also the G+C count, but not with the GC/AT ratio. Analysis of the detected mutagenized positions identified CG-rich patches, and flanking sequences strongly influenced EMS-induced mutation rates. This method for detecting false positive-induced mutations is generally applicable to any organism, is independent of the choice of in silico variant-calling algorithm, and is most valuable when the true mutation rate is likely to be low, such as in laboratory-induced mutations or somatic mutation detection in medicine.

Generation of EMS-Mutagenized Populations of Arabidopsis thaliana for Polyamine Genetics

In the recent years, genetic engineering of polyamine biosynthetic genes has provided evidence for their involvement in plant stress responses and different aspects of plant development. Such approaches are being complemented with the use of reverse genetics, in which mutants affected on a particular trait, tightly associated with polyamines, are isolated and the causal genes mapped. Reverse genetics enables the identification of novel genes in the polyamine pathway, which may be involved in downstream signaling, transport, homeostasis, or perception. Here, we describe a basic protocol for the generation of ethyl methanesulfonate (EMS) mutagenized populations of Arabidopsis thaliana for its use in reverse genetics applied to polyamines.

The Application of a Meiocyte-Specific CRISPR/Cas9 (MSC) System and a Suicide-MSC System in Generating Inheritable and Stable Mutations in Arabidopsis

The CRISPR/Cas9 system has been widely used for generating targeted mutations in various species. In Arabidopsis, it largely relies on the edited cells where the Cas9 protein performs its activity to obtain heritable and stable mutated lines. Here, we designed an improved CRISPR/Cas9 system, named as the MSC (meiocyte-specific CRISPR/Cas9) system, in which the Cas9 expression is driven by an experimentally approved meiocyte-specific promoter (AtDMC1 promoter). Two endogenous genes, including vegetative gene AtDET2 and reproductive gene AtDMC1, were targeted. We obtained heterozygous T1 plants for targeted genes with high efficiency (64%). In the T2 generation, the homozygous plants were abundant with high efficiency (37%). Analysis of Sanger sequencing results of T2 generation revealed that heritable gene mutations were high (52%). Moreover, we showed that the MSC system could sufficiently delete a middle size DNA fragment (∼500 bp) between two cleavage sites with a high rate (64.15%) in the T1 plants, providing direct evidence for making complete knock-out or certain domain-depletion mutations. In addition, we further made a suicide-MSC system, which can edit the targeted endogenous gene and the exogenous Cas9 gene simultaneously, not only successfully avoiding the further destroy of alleles brought in by molecular complementary or genic allelic test, but also maintaining the stable mutated alleles for functional studies. In short, the two systems provide new approaches to generate mutations for gene functional studies.

The HEM Lines: A New Library of Homozygous Arabidopsis thaliana EMS Mutants and its Potential to Detect Meiotic Phenotypes

Genetic screens have been crucial for deciphering many important biological processes, including meiosis. In Arabidopsis thaliana, previous forward screens have likely identified almost all the meiotic genes that when mutated lead to a pronounced decrease in fertility. However, the increasing number of genes identified in reverse genetics studies that play crucial roles in meiosis, but do not exhibit strong phenotypes when mutated, suggests that there are still many genes with meiotic function waiting to be discovered. In this study, we produced 897 A. thaliana homozygous mutant lines using Ethyl Methyl Sulfonate (EMS) mutagenesis followed by either single seed descent or haploid doubling. Whole genome sequencing of a subset of lines showed an average of 696 homozygous mutations per line, 195 of which (28%) modify a protein sequence. To test the power of this library, we carried out a forward screen looking for meiotic defects by observing chromosomes at metaphase I of male meiosis. Among the 649 lines analyzed, we identified 43 lines with meiotic defects. Of these, 21 lines had an obvious candidate causal mutation, namely a STOP or splicing site mutation in a gene previously shown to play a role in meiosis (ATM, MLH3, MLH1, MER3, HEI10, FLIP, ASY4, FLIP, PRD2, REC8, FANCL, and PSS1). Interestingly, this was the first time that six of these genes were identified in a forward screen in Arabidopsis (MLH3, MLH1, SGO1, PSS1, FANCL, and ASY4). These results illustrate the potential of this mutant population for screening for any qualitative or quantitative phenotype. Thus, this new mutant library is a powerful tool for functional genomics in A. thaliana. The HEM (Homozygote EMS Mutants) lines are available at the Versailles Arabidopsis stock center.

Mutation rate analysis via parent-progeny sequencing of the perennial peach. I. A low rate in woody perennials and a higher mutagenicity in hybrids

Mutation rates vary between species, between strains within species and between regions within a genome. What are the determinants of these forms of variation? Here, via parent–offspring sequencing of the peach we ask whether (i) woody perennials tend to have lower per unit time mutation rates compared to annuals, and (ii) hybrid strains have high mutation rates. Between a leaf from a low heterozygosity individual, derived from an intraspecific cross, to a leaf of its selfed progeny, the mutation rate is 7.77 × 10⁻⁹ point mutations per bp per generation, similar to Arabidopsis thaliana (7.0–7.4 × 10⁻⁹ point mutations per bp per generation). This suggests a low per unit time mutation rate as the generation time is much longer in peach. This is supported by our estimate of 9.48 × 10⁻⁹ point mutations per bp per generation from a 200-year-old low heterozygosity peach to its progeny. From a more highly heterozygous individual derived from an interspecific cross to its selfed progeny, the mutation rate is 1.38 × 10⁻⁸ mutations per site per generation, consistent with raised rates in hybrids. Our data thus suggest that (i) peach has an approximately order of magnitude lower mutation rate per unit time than Arabidopsis, consistent with reports of low evolutionary rates in woody perennials, and (ii) hybridization may, indeed, be associated with increased mutation rates as considered over a century ago.

The ancestral levels of transcription and the evolution of sexual phenotypes in filamentous fungi

Changes in gene expression have been hypothesized to play an important role in the evolution of divergent morphologies. To test this hypothesis in a model system, we examined differences in fruiting body morphology of five filamentous fungi in the Sordariomycetes, culturing them in a common garden environment and profiling genome-wide gene expression at five developmental stages. We reconstructed ancestral gene expression phenotypes, identifying genes with the largest evolved increases in gene expression across development. Conducting knockouts and performing phenotypic analysis in two divergent species typically demonstrated altered fruiting body development in the species that had evolved increased expression. Our evolutionary approach to finding relevant genes proved far more efficient than other gene deletion studies targeting whole genomes or gene families. Combining gene expression measurements with knockout phenotypes facilitated the refinement of Bayesian networks of the genes underlying fruiting body development, regulation of which is one of the least understood processes of multicellular development.