GABI-DUPLO: a collection of double mutants to overcome genetic redundancy in Arabidopsis thaliana

Establishment of single-cell transcriptional states during seed germination

Germination involves highly dynamic transcriptional programs as the cells of seeds reactivate and express the functions necessary for establishment in the environment. Individual cell types have distinct roles within the embryo, so must therefore have cell type-specific gene expression and gene regulatory networks. We can better understand how the functions of different cell types are established and contribute to the embryo by determining how cell type-specific transcription begins and changes through germination. Here we describe a temporal analysis of the germinating Arabidopsis thaliana embryo at single-cell resolution. We define the highly dynamic cell type-specific patterns of gene expression and how these relate to changing cellular function as germination progresses. Underlying these are unique gene regulatory networks and transcription factor activity. We unexpectedly discover that most embryo cells transition through the same initial transcriptional state early in germination, even though cell identity has already been established during embryogenesis. Cells later transition to cell type-specific gene expression patterns. Furthermore, our analyses support previous findings that the earliest events leading to the induction of seed germination take place in the vasculature. Overall, our study constitutes a general framework with which to characterize Arabidopsis cell transcriptional states through seed germination, allowing investigation of different genotypes and other plant species whose seed strategies may differ.

AGC kinases OXI1 and AGC2-2 regulate camalexin secretion and disease resistance by phosphorylating transporter PDR6

Plant transporters regulating the distribution of secondary metabolites play critical roles in defending against pathogens, insects, and interacting with beneficial microbes. The phosphorylation of these transporters can alter their activity, stability, and intracellular protein trafficking. However, the regulatory mechanism underlying this modification remains elusive. In this study, we discovered two orthologs of mammalian PKA, PKG, and PKC (AGC) kinases, oxidative signal-inducible 1 (OXI1) and its closest homologue, AGC subclass 2 member 2 (AGC2-2; 75% amino acid sequence identity with OXI1), associated with the extracellular secretion of camalexin and Arabidopsis (Arabidopsis thaliana) resistance to Pseudomonas syringae, and Botrytis cinerea. These kinases can undergo in vitro kinase reactions with three pleiotropic drug resistance (PDR) transporters: PDR6, PDR8, and PDR12. Moreover, our investigation confirmed PDR6 interaction with OXI1 and AGC2-2. By performing LC-MS/MS and parallel reaction monitoring, we identified the phosphorylation sites on PDR6 targeted by these kinases. Notably, chitin-induced PDR6 phosphorylation at specific residues, namely S31, S33, S827, and T832. Additional insights emerged by expressing dephosphorylated PDR6 variants in a pdr6 mutant background, revealing that the target residues S31, S33, and S827 promote PDR6 efflux activity, while T832 potentially contributes to PDR6 stability within the plasma membrane. The findings of this study elucidate partial mechanisms involved in the activity regulation of PDR-type transporters, providing valuable insights for their potential application in future plant breeding endeavors.

The prodomain of Arabidopsis metacaspase 2 positively regulates immune signaling mediated by pattern-recognition receptors

Metacaspases (MCs) are structural homologs of mammalian caspases found in plants, fungi, and protozoa. Type-I MCs carry an N-terminal prodomain, the function of which is unclear. Through genetic analysis of Arabidopsis mc2-1, a T-DNA insertion mutant of MC2, we demonstrated that the prodomain of metacaspase 2 (MC2) promotes immune signaling mediated by pattern-recognition receptors (PRRs). In mc2-1, immune responses are constitutively activated. The receptor-like kinases (RLKs) BAK1/BKK1 and SOBIR1 are required for the autoimmune phenotype of mc2-1, suggesting that immune signaling mediated by the receptor-like protein (RLP)-type PRRs is activated in mc2-1. A suppressor screen identified multiple mutations in the first exon of MC2, which suppress the autoimmunity in mc2-1. Further analysis revealed that the T-DNA insertion at the end of exon 1 of MC2 causes elevated expression of the MC2 prodomain, and overexpression of the MC2 prodomain in wild-type (WT) plants results in the activation of immune responses. The MC2 prodomain interacts with BIR1, which inhibits RLP-mediated immune signaling by interacting with BAK1, suggesting that the MC2 prodomain promotes plant defense responses by interfering with the function of BIR1. Our study uncovers an unexpected function of the prodomain of a MC in plant immunity.

Exploiting plant transcriptomic databases: Resources, tools, and approaches

There are now more than 300 000 RNA sequencing samples available, stemming from thousands of experiments capturing gene expression in organs, tissues, developmental stages, and experimental treatments for hundreds of plant species. The expression data have great value, as they can be re-analyzed by others to ask and answer questions that go beyond the aims of the study that generated the data. Because gene expression provides essential clues to where and when a gene is active, the data provide powerful tools for predicting gene function, and comparative analyses allow us to study plant evolution from a new perspective. This review describes how we can gain new knowledge from gene expression profiles, expression specificities, co-expression networks, differential gene expression, and experiment correlation. We also introduce and demonstrate databases that provide user-friendly access to these tools.

TOR coordinates nucleotide availability with ribosome biogenesis in plants

TARGET OF RAPAMYCIN (TOR) is a conserved eukaryotic Ser/Thr protein kinase that coordinates growth and metabolism with nutrient availability. We conducted a medium-throughput functional genetic screen to discover essential genes that promote TOR activity in plants, and identified a critical regulatory enzyme, cytosolic phosphoribosyl pyrophosphate (PRPP) synthetase (PRS4). PRS4 synthesizes cytosolic PRPP, a key upstream metabolite in nucleotide synthesis and salvage pathways. We found that prs4 knockouts are embryo-lethal in Arabidopsis thaliana, and that silencing PRS4 expression in Nicotiana benthamiana causes pleiotropic developmental phenotypes, including dwarfism, aberrant leaf shape, and delayed flowering. Transcriptomic analysis revealed that ribosome biogenesis is among the most strongly repressed processes in prs4 knockdowns. Building on these results, we discovered that TOR activity is inhibited by chemical or genetic disruption of nucleotide biosynthesis, but that this effect can be reversed by supplying plants with nucleobases. Finally, we show that TOR transcriptionally promotes nucleotide biosynthesis to support the demands of ribosomal RNA synthesis. We propose that TOR coordinates ribosome biogenesis with nucleotide availability in plants to maintain metabolic homeostasis and support growth.

Predictive Models of Genetic Redundancy in Arabidopsis thaliana

Genetic redundancy refers to a situation where an individual with a loss-of-function mutation in one gene (single mutant) does not show an apparent phenotype until one or more paralogs are also knocked out (double/higher-order mutant). Previous studies have identified some characteristics common among redundant gene pairs, but a predictive model of genetic redundancy incorporating a wide variety of features derived from accumulating omics and mutant phenotype data is yet to be established. In addition, the relative importance of these features for genetic redundancy remains largely unclear. Here, we establish machine learning models for predicting whether a gene pair is likely redundant or not in the model plant Arabidopsis thaliana based on six feature categories: functional annotations, evolutionary conservation including duplication patterns and mechanisms, epigenetic marks, protein properties including posttranslational modifications, gene expression, and gene network properties. The definition of redundancy, data transformations, feature subsets, and machine learning algorithms used significantly affected model performance based on holdout, testing phenotype data. Among the most important features in predicting gene pairs as redundant were having a paralog(s) from recent duplication events, annotation as a transcription factor, downregulation during stress conditions, and having similar expression patterns under stress conditions. We also explored the potential reasons underlying mispredictions and limitations of our studies. This genetic redundancy model sheds light on characteristics that may contribute to long-term maintenance of paralogs, and will ultimately allow for more targeted generation of functionally informative double mutants, advancing functional genomic studies.

Biotin attachment domain-containing proteins mediate hydroxy fatty acid-dependent inhibition of acetyl CoA carboxylase

Hundreds of naturally occurring specialized fatty acids (FAs) have potential as desirable chemical feedstocks if they could be produced at large scale by crop plants; however, transgenic expression of their biosynthetic genes has generally been accompanied by dramatic reductions in oil yield. For example, expression of castor (Ricinus communis) FA hydroxylase (FAH) in the Arabidopsis thaliana FA elongation mutant fae1 resulted in a 50% reduction of FA synthesis rate that was attributed to inhibition of acetyl-CoA carboxylase (ACCase) by an undefined mechanism. Here, we tested the hypothesis that the ricinoleic acid-dependent decrease in ACCase activity is mediated by biotin attachment domain-containing (BADC) proteins. BADCs are inactive homologs of biotin carboxy carrier protein that lack a biotin cofactor and can inhibit ACCase. Arabidopsis contains three BADC genes. To reduce expression levels of BADC1 and BADC3 in fae1/FAH plants, a homozygous badc1,3/fae1/FAH line was created. The rate of FA synthesis in badc1,3/fae1/FAH seeds doubled relative to fae1/FAH, restoring it to fae1 levels, increasing both native FA and HFA accumulation. Total FA per seed, seed oil content, and seed yield per plant all increased in badc1,3/fae1/FAH, to 5.8 µg, 37%, and 162 mg, respectively, relative to 4.9 µg, 33%, and 126 mg, respectively, for fae1/FAH. Transcript levels of FA synthesis-related genes, including those encoding ACCase subunits, did not significantly differ between badc1,3/fae1/FAH and fae1/FAH. These results demonstrate that BADC1 and BADC3 mediate ricinoleic acid-dependent inhibition of FA synthesis. We propose that BADC-mediated FAS inhibition as a general mechanism that limits FA accumulation in specialized FA-accumulating seeds.

Functional and evolutionary analysis of the Arabidopsis 4R-MYB protein SNAPc4 as part of the SNAP complex

Transcription initiation of the genes coding for small nuclear RNA (snRNA) has been extensively analyzed in humans and fruit fly, but only a single ortholog of a snRNA-activating protein complex (SNAPc) subunit has so far been characterized in plants. The genome of the model plant Arabidopsis thaliana encodes orthologs of all three core SNAPc subunits, including A. thaliana SNAP complex 4 (AtSNAPc4)-a 4R-MYB-type protein with four-and-a-half adjacent MYB repeat units. We report the conserved role of AtSNAPc4 as subunit of a protein complex involved in snRNA gene transcription and present genetic evidence that AtSNAPc4 is an essential gene in gametophyte and zygote development. We present experimental evidence that the three A. thaliana SNAPc subunits assemble into a SNAP complex and demonstrate the binding of AtSNAPc4 to snRNA promoters. In addition, co-localization studies show a link between AtSNAPc4 accumulation and Cajal bodies, known to aggregate at snRNA gene loci in humans. Moreover, we show the strong evolutionary conservation of single-copy 4R-MYB/SNAPc4 genes in a broad range of eukaryotes and present additional shared protein features besides the MYB domain, suggesting a conservation of the snRNA transcription initiation machinery along the course of the eukaryotic evolution.

Highly efficient multiplex editing: one-shot generation of 8× Nicotiana benthamiana and 12× Arabidopsis mutants

Genome editing by RNA-guided nucleases, such as SpCas9, has been used in numerous different plant species. However, to what extent multiple independent loci can be targeted simultaneously by multiplexing has not been well documented. Here, we developed a toolkit, based on a highly intron-optimized zCas9i gene, which allows assembly of nuclease constructs expressing up to 32 single guide RNAs (sgRNAs). We used this toolkit to explore the limits of multiplexing in two major model species, and report on the isolation of transgene-free octuple (8×) Nicotiana benthamiana and duodecuple (12×) Arabidopsis thaliana mutant lines in a single generation (T1 and T2 , respectively). We developed novel counter-selection markers for N. benthamiana, most importantly Sl-FAST2, comparable to the well-established Arabidopsis seed fluorescence marker, and FCY-UPP, based on the production of toxic 5-fluorouracil in the presence of a precursor. Targeting eight genes with an array of nine different sgRNAs and relying on FCY-UPP for selection of non-transgenic T1 , we identified N. benthamiana mutant lines with astonishingly high efficiencies: All analyzed plants carried mutations in all genes (approximately 112/116 target sites edited). Furthermore, we targeted 12 genes by an array of 24 sgRNAs in A. thaliana. Efficiency was significantly lower in A. thaliana, and our results indicate Cas9 availability is the limiting factor in such higher-order multiplexing applications. We identified a duodecuple mutant line by a combination of phenotypic screening and amplicon sequencing. The resources and results presented provide new perspectives for how multiplexing can be used to generate complex genotypes or to functionally interrogate groups of candidate genes.

CAFRI-Rice: CRISPR applicable functional redundancy inspector to accelerate functional genomics in rice

Rice (Oryza sativa L.) is a staple crop with agricultural traits that have been intensively investigated. However, despite the variety of mutant population and multi-omics data that have been generated, rice functional genomic research has been bottlenecked due to the functional redundancy in the genome. This phenomenon has masked the phenotypes of knockout mutants by functional compensation and redundancy. Here, we present an intuitive tool, CRISPR applicable functional redundancy inspector to accelerate functional genomics in rice (CAFRI-Rice; cafri-rice.khu.ac.kr). To create this tool, we generated a phylogenetic heatmap that can estimate the similarity between protein sequences and expression patterns, based on 2,617 phylogenetic trees and eight tissue RNA-sequencing datasets. In this study, 33,483 genes were sorted into 2,617 families, and about 24,980 genes were tested for functional redundancy using a phylogenetic heatmap approach. It was predicted that 7,075 genes would have functional redundancy, according to the threshold value validated by an analysis of 111 known genes functionally characterized using knockout mutants and 5,170 duplicated genes. In addition, our analysis demonstrated that an anther/pollen-preferred gene cluster has more functional redundancy than other clusters. Finally, we showed the usefulness of the CAFRI-Rice-based approach by overcoming the functional redundancy between two root-preferred genes via loss-of-function analyses as well as confirming the functional dominancy of three genes through a literature search. This CAFRI-Rice-based target selection for CRISPR/Cas9-mediated mutagenesis will not only accelerate functional genomic studies in rice but can also be straightforwardly expanded to other plant species.

Genome-wide identification of EMBRYO-DEFECTIVE (EMB) genes required for growth and development in Arabidopsis

With the emergence of high-throughput methods in plant biology, the importance of long-term projects characterized by incremental advances involving multiple laboratories can sometimes be overlooked. Here, I highlight my 40-year effort to isolate and characterize the most common class of mutants encountered in Arabidopsis (Arabidopsis thaliana): those defective in embryo development. I present an updated dataset of 510 EMBRYO-DEFECTIVE (EMB) genes identified throughout the Arabidopsis community; include important details on 2200 emb mutants and 241 pigment-defective embryo (pde) mutants analyzed in my laboratory; provide curated datasets with key features and publication links for each EMB gene identified; revisit past estimates of 500-1000 total EMB genes in Arabidopsis; document 83 double mutant combinations reported to disrupt embryo development; emphasize the importance of following established nomenclature guidelines and acknowledging allele history in research publications; and consider how best to extend community-based curation and screening efforts to approach saturation for this diverse class of mutants in the future. Continued advances in identifying EMB genes and characterizing their loss-of-function mutant alleles are needed to understand genotype-to-phenotype relationships in Arabidopsis on a broad scale, and to document the contributions of large numbers of essential genes to plant growth and development.

Distributed phenomics with the unPAK project reveals the effects of mutations

Determining how genes are associated with traits in plants and other organisms is a major challenge in modern biology. The unPAK project - undergraduates phenotyping Arabidopsis knockouts - has generated phenotype data for thousands of non-lethal insertion mutation lines within a single Arabidopsis thaliana genomic background. The focal phenotypes examined by unPAK are complex macroscopic fitness-related traits, which have ecological, evolutionary and agricultural importance. These phenotypes are placed in the context of the wild-type and also natural accessions (phytometers), and standardized for environmental differences between assays. Data from the unPAK project are used to describe broad patterns in the phenotypic consequences of insertion mutation, and to identify individual mutant lines with distinct phenotypes as candidates for further study. Inclusion of undergraduate researchers is at the core of unPAK activities, and an important broader impact of the project is providing students an opportunity to obtain research experience.

A set of Arabidopsis genes involved in the accommodation of the downy mildew pathogen Hyaloperonospora arabidopsidis

The intracellular accommodation structures formed by plant cells to host arbuscular mycorrhiza fungi and biotrophic hyphal pathogens are cytologically similar. Therefore we investigated whether these interactions build on an overlapping genetic framework. In legumes, the malectin-like domain leucine-rich repeat receptor kinase SYMRK, the cation channel POLLUX and members of the nuclear pore NUP107-160 subcomplex are essential for symbiotic signal transduction and arbuscular mycorrhiza development. We identified members of these three groups in Arabidopsis thaliana and explored their impact on the interaction with the oomycete downy mildew pathogen Hyaloperonospora arabidopsidis (Hpa). We report that mutations in the corresponding genes reduced the reproductive success of Hpa as determined by sporangiophore and spore counts. We discovered that a developmental transition of haustorial shape occurred significantly earlier and at higher frequency in the mutants. Analysis of the multiplication of extracellular bacterial pathogens, Hpa-induced cell death or callose accumulation, as well as Hpa- or flg22-induced defence marker gene expression, did not reveal any traces of constitutive or exacerbated defence responses. These findings point towards an overlap between the plant genetic toolboxes involved in the interaction with biotrophic intracellular hyphal symbionts and pathogens in terms of the gene families involved.

Our work reveals genetic commonalities between biotrophic intracellular interactions with symbiotic and pathogenic hyphal microbes. The majority of land plants engages in arbuscular mycorrhiza (AM) symbiosis with phosphate-acquiring arbuscular mycorrhizal fungi to avoid phosphate starvation. Nutrient exchange in this interaction occurs via arbuscules, tree-shaped fungal structures, hosted within plant root cells. A series of plant genes including the Symbiosis Receptor-like kinase (SYMRK), members of the NUP107-160 subcomplex and nuclear envelope localised cation channels are required for a signalling process leading to the development of AM. The model plant Arabidopsis thaliana lost the ability to form AM. Although the ortholog of SYMRK was deleted during evolution, members of the malectin-like domain leucine-rich repeat receptor kinase (MLD-LRR-RK) gene family, components of the NUP107-160 subcomplex, and an ortholog of the nuclear envelope-localized cation channel POLLUX, are still present in the Arabidopsis genome, and Arabidopsis leaf cells retained the ability to accommodate haustoria, presumed feeding structures of the obligate biotrophic downy mildew pathogen Hyaloperonospora arabidopsidis. We discovered that both of these plant-microbe interactions utilize a corresponding set of genes including the ortholog of POLLUX, members of the NUP107-160 subcomplex and members of the MLD-LRR-RK gene family, thus revealing similarities in the plant program for the intracellular accommodation of biotrophic organisms in symbiosis and disease.

Transcript Profiling Identifies NAC-Domain Genes Involved in Regulating Wall Ingrowth Deposition in Phloem Parenchyma Transfer Cells of Arabidopsis thaliana

Transfer cells (TCs) play important roles in facilitating enhanced rates of nutrient transport at key apoplasmic/symplasmic junctions along the nutrient acquisition and transport pathways in plants. TCs achieve this capacity by developing elaborate wall ingrowth networks which serve to increase plasma membrane surface area thus increasing the cell's surface area-to-volume ratio to achieve increased flux of nutrients across the plasma membrane. Phloem parenchyma (PP) cells of Arabidopsis leaf veins trans-differentiate to become PP TCs which likely function in a two-step phloem loading mechanism by facilitating unloading of photoassimilates into the apoplasm for subsequent energy-dependent uptake into the sieve element/companion cell (SE/CC) complex. We are using PP TCs in Arabidopsis as a genetic model to identify transcription factors involved in coordinating deposition of the wall ingrowth network. Confocal imaging of pseudo-Schiff propidium iodide-stained tissue revealed different profiles of temporal development of wall ingrowth deposition across maturing cotyledons and juvenile leaves, and a basipetal gradient of deposition across mature adult leaves. RNA-Seq analysis was undertaken to identify differentially expressed genes common to these three different profiles of wall ingrowth deposition. This analysis identified 68 transcription factors up-regulated two-fold or more in at least two of the three experimental comparisons, with six of these transcription factors belonging to Clade III of the NAC-domain family. Phenotypic analysis of these NAC genes using insertional mutants revealed significant reductions in levels of wall ingrowth deposition, particularly in a double mutant of NAC056 and NAC018, as well as compromised sucrose-dependent root growth, indicating impaired capacity for phloem loading. Collectively, these results support the proposition that Clade III members of the NAC-domain family in Arabidopsis play important roles in regulating wall ingrowth deposition in PP TCs.

Chlorophyll Fluorescence Video Imaging: A Versatile Tool for Identifying Factors Related to Photosynthesis

Measurements of chlorophyll fluorescence provide an elegant and non-invasive means of probing the dynamics of photosynthesis. Advances in video imaging of chlorophyll fluorescence have now made it possible to study photosynthesis at all levels from individual cells to entire crop populations. Since the technology delivers quantitative data, is easily scaled up and can be readily combined with other approaches, it has become a powerful phenotyping tool for the identification of factors relevant to photosynthesis. Here, we review genetic chlorophyll fluorescence-based screens of libraries of Arabidopsis and Chlamydomonas mutants, discuss its application to high-throughput phenotyping in quantitative genetics and highlight potential future developments.

Upstream kinases of plant SnRKs are involved in salt stress tolerance

Sucrose non-fermenting 1-related protein kinases (SnRKs) are important for plant growth and stress responses. This family has three clades: SnRK1, SnRK2 and SnRK3. Although plant SnRKs are thought to be activated by upstream kinases, the overall mechanism remains obscure. Geminivirus Rep-Interacting Kinase (GRIK)1 and GRIK2 phosphorylate SnRK1s, which are involved in sugar/energy sensing, and the grik1-1 grik2-1 double mutant shows growth retardation under regular growth conditions. In this study, we established another Arabidopsis mutant line harbouring a different allele of gene GRIK1 (grik1-2 grik2-1) that grows similarly to the wild-type, enabling us to evaluate the function of GRIKs under stress conditions. In the grik1-2 grik2-1 double mutant, phosphorylation of SnRK1.1 was reduced, but not eliminated, suggesting that the grik1-2 mutation is a weak allele. In addition to high sensitivity to glucose, the grik1-2 grik2-1 mutant was sensitive to high salt, indicating that GRIKs are also involved in salinity signalling pathways. Salt Overly Sensitive (SOS)2, a member of the SnRK3 subfamily, is a critical mediator of the response to salinity. GRIK1 phosphorylated SOS2 in vitro, resulting in elevated kinase activity of SOS2. The salt tolerance of sos2 was restored to normal levels by wild-type SOS2, but not by a mutated form of SOS2 lacking the T168 residue phosphorylated by GRIK1. Activation of SOS2 by GRIK1 was also demonstrated in a reconstituted system in yeast. Our results indicate that GRIKs phosphorylate and activate SnRK1 and other members of the SnRK3 family, and that they play important roles in multiple signalling pathways in vivo.

Fitness effects of mutation: testing genetic redundancy in Arabidopsis thaliana

Screens of organisms with disruptive mutations in a single gene often fail to detect phenotypic consequences for the majority of mutants. One explanation for this phenomenon is that the presence of paralogous loci provides genetic redundancy. However, it is also possible that the assayed traits are affected by few loci, that effects could be subtle or that phenotypic effects are restricted to certain environments. We assayed a set of T-DNA insertion mutant lines of Arabidopsis thaliana to determine the frequency with which mutation affected fitness-related phenotypes. We found that between 8% and 42% of the assayed lines had altered fitness from the wild type. Furthermore, many of these lines exhibited fitness greater than the wild type. In a second experiment, we grew a subset of the lines in multiple environments and found whether a T-DNA insert increased or decreased fitness traits depended on the assay environment. Overall, our evidence contradicts the hypothesis that genetic redundancy is a common phenomenon in A. thaliana for fitness traits. Evidence for redundancy from prior screens of knockout mutants may often be an artefact of the design of the phenotypic assays which have focused on less complex phenotypes than fitness and have used single environments. Finally, our study adds to evidence that beneficial mutations may represent a significant component of the mutational spectrum of A. thaliana.

The impact of Arabidopsis thaliana SNF1-related-kinase 1 (SnRK1)-activating kinase 1 (SnAK1) and SnAK2 on SnRK1 phosphorylation status: characterization of a SnAK double mutant

Arabidopsis thaliana SNF1-related-kinase 1 (SnRK1)-activating kinase 1 (AtSnAK1) and AtSnAK2 have been shown to phosphorylate in vitro and activate the energy signalling integrator, SnRK1. To clarify this signalling cascade in planta, a genetic- and molecular-based approach was developed. Homozygous single AtSnAK1 and AtSnAK2 T-DNA insertional mutants did not display an apparent phenotype. Crossing of the single mutants did not allow the isolation of double-mutant plants, whereas self-pollinating the S1-/- S2+/- sesquimutant specifically gave approximatively 22% individuals in their offspring that, when rescued on sugar-supplemented media in vitro, were shown to be AtSnAK1 AtSnAK2 double mutants. Interestingly, this was not obtained in the case of the other sesquimutant, S1+/- S2-/-. Although reduced in size, the double mutant had the capacity to produce flowers, but not seeds. Immunological characterization established the T-loop of the SnRK1 catalytic subunit to be non-phosphorylated in the absence of both SnAKs. When the double mutant was complemented with a DNA construct containing an AtSnAK2 open reading frame driven by its own promoter, a normal phenotype was restored. Therefore, wild-type plant growth and development is dependent on the presence of SnAK in vivo, and this is correlated with SnRK1 phosphorylation. These data show that both SnAKs are kinases phosphorylating SnRK1, and thereby they contribute to energy signalling in planta.

Enhancing the GABI-Kat Arabidopsis thaliana T-DNA Insertion Mutant Database by Incorporating Araport11 Annotation

SimpleSearch provides access to a database containing information about T-DNA insertion lines of the GABI-Kat collection of Arabidopsis thaliana mutants. These mutants are an important tool for reverse genetics, and GABI-Kat is the second largest collection of such T-DNA insertion mutants. Insertion sites were deduced from flanking sequence tags (FSTs), and the database contains information about mutant plant lines as well as insertion alleles. Here, we describe improvements within the interface (available at http://www.gabi-kat.de/db/genehits.php) and with regard to the database content that have been realized in the last five years. These improvements include the integration of the Araport11 genome sequence annotation data containing the recently updated A. thaliana structural gene descriptions, an updated visualization component that displays groups of insertions with very similar insertion positions, mapped confirmation sequences, and primers. The visualization component provides a quick way to identify insertions of interest, and access to improved data about the exact structure of confirmed insertion alleles. In addition, the database content has been extended by incorporating additional insertion alleles that were detected during the confirmation process, as well as by adding new FSTs that have been produced during continued efforts to complement gaps in FST availability. Finally, the current database content regarding predicted and confirmed insertion alleles as well as primer sequences has been made available as downloadable flat files.

A double-mutant collection targeting MAP kinase related genes in Arabidopsis for studying genetic interactions

Mitogen-activated protein kinase cascades are conserved in all eukaryotes. In Arabidopsis thaliana there are approximately 80 genes encoding MAP kinase kinase kinases (MAP3K), 10 genes encoding MAP kinase kinases (MAP2K), and 20 genes encoding MAP kinases (MAPK). Reverse genetic analysis has failed to reveal abnormal phenotypes for a majority of these genes. One strategy for uncovering gene function when single-mutant lines do not produce an informative phenotype is to perform a systematic genetic interaction screen whereby double-mutants are created from a large library of single-mutant lines. Here we describe a new collection of 275 double-mutant lines derived from a library of single-mutants targeting genes related to MAP kinase signaling. To facilitate this study, we developed a high-throughput double-mutant generating pipeline using a system for growing Arabidopsis seedlings in 96-well plates. A quantitative root growth assay was used to screen for evidence of genetic interactions in this double-mutant collection. Our screen revealed four genetic interactions, all of which caused synthetic enhancement of the root growth defects observed in a MAP kinase 4 (MPK4) single-mutant line. Seeds for this double-mutant collection are publicly available through the Arabidopsis Biological Resource Center. Scientists interested in diverse biological processes can now screen this double-mutant collection under a wide range of growth conditions in order to search for additional genetic interactions that may provide new insights into MAP kinase signaling.

The Structural Features of Thousands of T-DNA Insertion Sites Are Consistent with a Double-Strand Break Repair-Based Insertion Mechanism

Transformation by Agrobacterium tumefaciens, an important tool in modern plant research, involves the integration of T-DNA initially present on a plasmid in agrobacteria into the genome of plant cells. The process of attachment of the agrobacteria to plant cells and the transport of T-DNA into the cell and further to the nucleus has been well described. However, the exact mechanism of integration into the host's DNA is still unclear, although several models have been proposed. During confirmation of T-DNA insertion alleles from the GABI-Kat collection of Arabidopsis thaliana mutants, we have generated about 34,000 sequences from the junctions between inserted T-DNA and adjacent genome regions. Here, we describe the evaluation of this dataset with regard to existing models for T-DNA integration. The results suggest that integration into the plant genome is mainly mediated by the endogenous plant DNA repair machinery. The observed integration events showed characteristics highly similar to those of repair sites of double-strand breaks with respect to microhomology and deletion sizes. In addition, we describe unexpected integration events, such as large deletions and inversions at the integration site that are relevant for correct interpretation of results from T-DNA insertion mutants in reverse genetics experiments.

Genes and networks regulating root anatomy and architecture

The root is an excellent model for studying developmental processes that underlie plant anatomy and architecture. Its modular structure, the lack of cell movement and relative accessibility to microscopic visualization facilitate research in a number of areas of plant biology. In this review, we describe several examples that demonstrate how cell type-specific developmental mechanisms determine cell fate and the formation of defined tissues with unique characteristics. In the last 10 yr, advances in genome-wide technologies have led to the sequencing of thousands of plant genomes, transcriptomes and proteomes. In parallel with the development of these high-throughput technologies, biologists have had to establish computational, statistical and bioinformatic tools that can deal with the wealth of data generated by them. These resources provide a foundation for posing more complex questions about molecular interactions, and have led to the discovery of new mechanisms that control phenotypic differences. Here we review several recent studies that shed new light on developmental processes, which are involved in establishing root anatomy and architecture. We highlight the power of combining large-scale experiments with classical techniques to uncover new pathways in root development.

A Genome-Wide Association Analysis Reveals Epistatic Cancellation of Additive Genetic Variance for Root Length in Arabidopsis thaliana

Efforts to identify loci underlying complex traits generally assume that most genetic variance is additive. Here, we examined the genetics of Arabidopsis thaliana root length and found that the genomic narrow-sense heritability for this trait in the examined population was statistically zero. The low amount of additive genetic variance that could be captured by the genome-wide genotypes likely explains why no associations to root length could be found using standard additive-model-based genome-wide association (GWA) approaches. However, as the broad-sense heritability for root length was significantly larger, and primarily due to epistasis, we also performed an epistatic GWA analysis to map loci contributing to the epistatic genetic variance. Four interacting pairs of loci were revealed, involving seven chromosomal loci that passed a standard multiple-testing corrected significance threshold. The genotype-phenotype maps for these pairs revealed epistasis that cancelled out the additive genetic variance, explaining why these loci were not detected in the additive GWA analysis. Small population sizes, such as in our experiment, increase the risk of identifying false epistatic interactions due to testing for associations with very large numbers of multi-marker genotypes in few phenotyped individuals. Therefore, we estimated the false-positive risk using a new statistical approach that suggested half of the associated pairs to be true positive associations. Our experimental evaluation of candidate genes within the seven associated loci suggests that this estimate is conservative; we identified functional candidate genes that affected root development in four loci that were part of three of the pairs. The statistical epistatic analyses were thus indispensable for confirming known, and identifying new, candidate genes for root length in this population of wild-collected A. thaliana accessions. We also illustrate how epistatic cancellation of the additive genetic variance explains the insignificant narrow-sense and significant broad-sense heritability by using a combination of careful statistical epistatic analyses and functional genetic experiments.

Complex traits, such as many human diseases or climate adaptation and production traits in crops, arise through the action and interaction of many genes and environmental factors. Classic approaches to identify contributing genes generally assume that these factors contribute mainly additive genetic variance. Recent methods, such as genome-wide association studies, often adhere to this additive genetics paradigm. However, additive models of complex traits do not reflect that genes can also contribute with non-additive genetic variance. In this study, we use Arabidopsis thaliana to determine the additive and non-additive genetic contributions to the phenotypic variation in root length. Surprisingly, much of the observed phenotypic variation in root length across genetically divergent strains was explained by epistasis. We mapped seven loci contributing to the epistatic genetic variance and validated four genes in these loci with mutant analysis. For three of these genes, this is their first implication in root development. Together, our results emphasize the importance of considering both non-additive and additive genetic variance when dissecting complex trait variation, in order not to lose sensitivity in genetic analyses.

Unmasking alternative splicing inside protein-coding exons defines exitrons and their role in proteome plasticity

Alternative splicing (AS) diversifies transcriptomes and proteomes and is widely recognized as a key mechanism for regulating gene expression. Previously, in an analysis of intron retention events in Arabidopsis, we found unusual AS events inside annotated protein-coding exons. Here, we also identify such AS events in human and use these two sets to analyse their features, regulation, functional impact, and evolutionary origin. As these events involve introns with features of both introns and protein-coding exons, we name them exitrons (exonic introns). Though exitrons were detected as a subset of retained introns, they are clearly distinguishable, and their splicing results in transcripts with different fates. About half of the 1002 Arabidopsis and 923 human exitrons have sizes of multiples of 3 nucleotides (nt). Splicing of these exitrons results in internally deleted proteins and affects protein domains, disordered regions, and various post-translational modification sites, thus broadly impacting protein function. Exitron splicing is regulated across tissues, in response to stress and in carcinogenesis. Intriguingly, annotated intronless genes can be also alternatively spliced via exitron usage. We demonstrate that at least some exitrons originate from ancestral coding exons. Based on our findings, we propose a “splicing memory” hypothesis whereby upon intron loss imprints of former exon borders defined by vestigial splicing regulatory elements could drive the evolution of exitron splicing. Altogether, our studies show that exitron splicing is a conserved strategy for increasing proteome plasticity in plants and animals, complementing the repertoire of AS events.

Comprehensive identification of mutations induced by heavy-ion beam irradiation in Arabidopsis thaliana

Heavy-ion beams are widely used for mutation breeding and molecular biology. Although the mutagenic effects of heavy-ion beam irradiation have been characterized by sequence analysis of some restricted chromosomal regions or loci, there have been no evaluations at the whole-genome level or of the detailed genomic rearrangements in the mutant genomes. In this study, using array comparative genomic hybridization (array-CGH) and resequencing, we comprehensively characterized the mutations in Arabidopsis thaliana genomes irradiated with Ar or Fe ions. We subsequently used this information to investigate the mutagenic effects of the heavy-ion beams. Array-CGH demonstrated that the average number of deleted areas per genome were 1.9 and 3.7 following Ar-ion and Fe-ion irradiation, respectively, with deletion sizes ranging from 149 to 602,180 bp; 81% of the deletions were accompanied by genomic rearrangements. To provide a further detailed analysis, the genomes of the mutants induced by Ar-ion beam irradiation were resequenced, and total mutations, including base substitutions, duplications, in/dels, inversions, and translocations, were detected using three algorithms. All three resequenced mutants had genomic rearrangements. Of the 22 DNA fragments that contributed to the rearrangements, 19 fragments were responsible for the intrachromosomal rearrangements, and multiple rearrangements were formed in the localized regions of the chromosomes. The interchromosomal rearrangements were detected in the multiply rearranged regions. These results indicate that the heavy-ion beams led to clustered DNA damage in the chromosome, and that they have great potential to induce complicated intrachromosomal rearrangements. Heavy-ion beams will prove useful as unique mutagens for plant breeding and the establishment of mutant lines.

AtNIGT1/HRS1 integrates nitrate and phosphate signals at the Arabidopsis root tip

Nitrogen and phosphorus are among the most widely used fertilizers worldwide. Nitrate (NO₃⁻) and phosphate (PO₄³⁻) are also signaling molecules whose respective transduction pathways are being intensively studied. However, plants are continuously challenged with combined nutritional deficiencies, yet very little is known about how these signaling pathways are integrated. Here we report the identification of a highly NO₃⁻-inducible NRT1.1-controlled GARP transcription factor, HRS1, document its genome-wide transcriptional targets, and validate its cis-regulatory-elements. We demonstrate that this transcription factor and a close homolog repress primary root growth in response to P deficiency conditions, but only when NO₃⁻ is present. This system defines a molecular logic gate integrating P and N signals. We propose that NO₃⁻ and P signaling converge via double transcriptional and post-transcriptional control of the same protein, HRS1

Molecular systems governing leaf growth: from genes to networks

Arabidopsis leaf growth consists of a complex sequence of interconnected events involving cell division and cell expansion, and requiring multiple levels of genetic regulation. With classical genetics, numerous leaf growth regulators have been identified, but the picture is far from complete. With the recent advances made in quantitative phenotyping, the study of the quantitative, dynamic, and multifactorial features of leaf growth is now facilitated. The use of high-throughput phenotyping technologies to study large numbers of natural accessions or mutants, or to screen for the effects of large sets of chemicals will allow for further identification of the additional players that constitute the leaf growth regulatory networks. Only a tight co-ordination between these numerous molecular players can support the formation of a functional organ. The connections between the components of the network and their dynamics can be further disentangled through gene-stacking approaches and ultimately through mathematical modelling. In this review, we describe these different approaches that should help to obtain a holistic image of the molecular regulation of organ growth which is of high interest in view of the increasing needs for plant-derived products.

Update on transparent testa mutants from Arabidopsis thaliana: characterisation of new alleles from an isogenic collection

We present a comprehensive overview on flavonoid-related phenotypes of A. thaliana tt and tds mutants, provide tools for their characterisation, increase the number of available alleles and demonstrate that tds3 is allelic to tt12 and tds5 to aha10. Flavonoid biosynthesis is one of the best-studied secondary metabolite pathways in plants. In the model system Arabidopsis thaliana it leads to the synthesis of three phenolic compound classes: flavonol glycosides, anthocyanins and proanthocyanidins (PAs). PAs appear brown in their oxidised polymeric forms, and most A. thaliana mutants impaired in flavonoid accumulation were identified through screens for lack of this seed coat pigmentation. These mutants are referred to as transparent testa (tt) or tannin-deficient seed (tds). More than 20 mutants of these types have been published, probably representing most of the genes relevant for PA accumulation in A. thaliana. However, data about the genes involved in PA deposition or oxidation are still rather scarce. Also, for some of the known mutants it is unclear if they represent additional loci or if they are allelic to known genes. For the present study, we have performed a systematic phenotypic characterisation of almost all available tt and tds mutants and built a collection of mutants in the genetic background of the accession Columbia to minimise effects arising from ecotype variation. We have identified a novel tt6 allele from a forward genetic screen and demonstrated that tds3 is allelic to tt12 and tds5 to aha10.

Multi-gene silencing in Arabidopsis: a collection of artificial microRNAs targeting groups of paralogs encoding transcription factors

Functional redundancy often hampers the analysis of gene families. To overcome this difficulty, we constructed Arabidopsis thaliana lines that expressed artificial microRNAs designed to simultaneously target two to six paralogous genes encoding members of transcription factor families. Of the 576 genes that we chose as targets, only 122 had already been functionally studied at some level. As a simple indicator of the inhibitory effects of our amiRNAs on their targets, we examined the amiRNA-expressing transgenic lines for morphological phenotypes at the rosette stage. Of 338 transgenes tested, 21 caused a visible morphological phenotype in leaves, a proportion that is much higher than that expected as a result of insertional mutagenesis. Also, our collection probably represents many other mutant phenotypes, not just those in leaves. This robust, versatile method enables functional examination of redundant transcription factor paralogs, and is particularly useful for genes that occur in tandem.

Leaf phenomics: a systematic reverse genetic screen for Arabidopsis leaf mutants

The study and eventual manipulation of leaf development in plants requires a thorough understanding of the genetic basis of leaf organogenesis. Forward genetic screens have identified hundreds of Arabidopsis mutants with altered leaf development, but the genome has not yet been saturated. To identify genes required for leaf development we are screening the Arabidopsis Salk Unimutant collection. We have identified 608 lines that exhibit a leaf phenotype with full penetrance and almost constant expressivity and 98 additional lines with segregating mutant phenotypes. To allow indexing and integration with other mutants, the mutant phenotypes were described using a custom leaf phenotype ontology. We found that the indexed mutation is present in the annotated locus for 78% of the 553 mutants genotyped, and that in half of these the annotated T-DNA is responsible for the phenotype. To quickly map non-annotated T-DNA insertions, we developed a reliable, cost-effective and easy method based on whole-genome sequencing. To enable comprehensive access to our data, we implemented a public web application named PhenoLeaf (http://genetics.umh.es/phenoleaf) that allows researchers to query the results of our screen, including text and visual phenotype information. We demonstrated how this new resource can facilitate gene function discovery by identifying and characterizing At1g77600, which we found to be required for proximal-distal cell cycle-driven leaf growth, and At3g62870, which encodes a ribosomal protein needed for cell proliferation and chloroplast function. This collection provides a valuable tool for the study of leaf development, characterization of biomass feedstocks and examination of other traits in this fundamental photosynthetic organ.

Defined single-gene and multi-gene deletion mutant collections in Salmonella enterica sv Typhimurium

We constructed two collections of targeted single gene deletion (SGD) mutants and two collections of targeted multi-gene deletion (MGD) mutants in Salmonella enterica sv Typhimurium 14028s. The SGD mutant collections contain (1), 3517 mutants in which a single gene is replaced by a cassette containing a kanamycin resistance (KanR) gene oriented in the sense direction (SGD-K), and (2), 3376 mutants with a chloramphenicol resistance gene (CamR) oriented in the antisense direction (SGD-C). A combined total of 3773 individual genes were deleted across these SGD collections. The MGD collections contain mutants bearing deletions of contiguous regions of three or more genes and include (3), 198 mutants spanning 2543 genes replaced by a KanR cassette (MGD-K), and (4), 251 mutants spanning 2799 genes replaced by a CamR cassette (MGD-C). Overall, 3476 genes were deleted in at least one MGD collection. The collections with different antibiotic markers permit construction of all viable combinations of mutants in the same background. Together, the libraries allow hierarchical screening of MGDs for different phenotypic followed by screening of SGDs within the target MGD regions. The mutants of these collections are stored at BEI Resources (www.beiresources.org) and publicly available.

Genome-enabled plant metabolomics

The grand challenge currently facing metabolomics is that of comprehensitivity whilst next generation sequencing and advanced proteomics methods now allow almost complete and at least 50% coverage of their respective target molecules, metabolomics platforms at best offer coverage of just 10% of the small molecule complement of the cell. Here we discuss the use of genome sequence information as an enabling tool for peak identity and for translational metabolomics. Whilst we argue that genome information is not sufficient to compute the size of a species metabolome it is highly useful in predicting the occurrence of a wide range of common metabolites. Furthermore, we describe how via gene functional analysis in model species the identity of unknown metabolite peaks can be resolved. Taken together these examples suggest that genome sequence information is current (and likely will remain), a highly effective tool in peak elucidation in mass spectral metabolomics strategies.

An easy-to-use primer design tool to address paralogous loci and T-DNA insertion sites in the genome of Arabidopsis thaliana

BACKGROUND

More than 90% of the Arabidopsis thaliana genes are members of multigene families. DNA sequence similarities present in such related genes can cause trouble, e.g. when molecularly analysing mutant alleles of these genes. Also, flanking-sequence-tag (FST) based predictions of T-DNA insertion positions are often located within paralogous regions of the genome. In such cases, the prediction of the correct insertion site must include careful sequence analyses on the one hand and a paralog specific primer design for experimental confirmation of the prediction on the other hand.

RESULTS

GABI-Kat is a large A. thaliana insertion line resource, which uses in-house confirmation to provide highly reliable access to T-DNA insertion alleles. To offer trustworthy mutant alleles of paralogous loci, we considered multiple insertion site predictions for single FSTs and implemented this 1-to-N relation in our database. The resulting paralogous predictions were addressed experimentally and the correct insertion locus was identified in most cases, including cases in which there were multiple predictions with identical prediction scores. A newly developed primer design tool that takes paralogous regions into account was developed to streamline the confirmation process for paralogs. The tool is suitable for all parts of the genome and is freely available at the GABI-Kat website. Although the tool was initially designed for the analysis of T-DNA insertion mutants, it can be used for any experiment that requires locus-specific primers for the A. thaliana genome. It is easy to use and also able to design amplimers with two genome-specific primers as required for genotyping segregating families of insertion mutants when looking for homozygous offspring.

CONCLUSIONS

The paralog-aware confirmation process significantly improved the reliability of the insertion site assignment when paralogous regions of the genome were affected. An automatic online primer design tool that incorporates experience from the in-house confirmation of T-DNA insertion lines has been made available. It provides easy access to primers for the analysis of T-DNA insertion alleles, but it is also beneficial for other applications as well.