CRISPR-Cas9-mediated efficient directed mutagenesis and RAD51-dependent and RAD51-independent gene targeting in the moss Physcomitrella patens

A CRISPR/LbCas12a-based method for highly efficient multiplex gene editing in Physcomitrella patens

Due to their high efficiency, specificity, and flexibility, programmable nucleases, such as those of the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas12a (Cpf1) system, have greatly expanded the applicability of editing the genomes of various organisms. Genes from different gene families or genes with redundant functions in the same gene family can be examined by assembling multiple CRISPR RNAs (crRNAs) in a single vector. However, the activity and efficiency of CRISPR/Cas12a in the non-vascular plant Physcomitrella patens are largely unknown. Here, we demonstrate that LbCas12a together with its mature crRNA can target multiple loci simultaneously in P. patens with high efficiency via co-delivery of LbCas12a and a crRNA expression cassette in vivo. The mutation frequencies induced by CRISPR/LbCas12a at a single locus ranged from 26.5 to 100%, with diverse deletions being the most common type of mutation. Our method expands the repertoire of genome editing tools available for P. patens and facilitates the creation of loss-of-function mutants of multiple genes from different gene families.

Heterogeneous Nuclear Ribonucleoprotein H1 Coordinates with Phytochrome and the U1 snRNP Complex to Regulate Alternative Splicing in Physcomitrella patens

Plant photoreceptors tightly regulate gene expression to control photomorphogenic responses. Although gene expression is modulated by photoreceptors at various levels, the regulatory mechanism at the pre-mRNA splicing step remains unclear. Alternative splicing, a widespread mechanism in eukaryotes that generates two or more mRNAs from the same pre-mRNA, is largely controlled by splicing regulators, which recruit spliceosomal components to initiate pre-mRNA splicing. The red/far-red light photoreceptor phytochrome participates in light-mediated splicing regulation, but the detailed mechanism remains unclear. Here, using protein-protein interaction analysis, we demonstrate that in the moss Physcomitrella patens, phytochrome4 physically interacts with the splicing regulator heterogeneous nuclear ribonucleoprotein H1 (PphnRNP-H1) in the nucleus, a process dependent on red light. We show that PphnRNP-H1 is involved in red light-mediated phototropic responses in P. patens and that it binds with higher affinity to the splicing factor pre-mRNA-processing factor39-1 (PpPRP39-1) in the presence of red light-activated phytochromes. Furthermore, PpPRP39-1 associates with the core component of U1 small nuclear RNP in P. patens Genome-wide analyses demonstrated the involvement of both PphnRNP-H1 and PpPRP39-1 in light-mediated splicing regulation. Our results suggest that phytochromes target the early step of spliceosome assembly via a cascade of protein-protein interactions to control pre-mRNA splicing and photomorphogenic responses.

Efficient and modular CRISPR-Cas9 vector system for Physcomitrella patens

CRISPR-Cas9 has been shown to be a valuable tool in recent years, allowing researchers to precisely edit the genome using an RNA-guided nuclease to initiate double-strand breaks. Until recently, classical RAD51-mediated homologous recombination has been a powerful tool for gene targeting in the moss Physcomitrella patens. However, CRISPR-Cas9-mediated genome editing in P. patens was shown to be more efficient than traditional homologous recombination (Plant Biotechnology Journal, 15, 2017, 122). CRISPR-Cas9 provides the opportunity to efficiently edit the genome at multiple loci as well as integrate sequences at precise locations in the genome using a simple transient transformation. To fully take advantage of CRISPR-Cas9 genome editing in P. patens, here we describe the generation and use of a flexible and modular CRISPR-Cas9 vector system. Without the need for gene synthesis, this vector system enables editing of up to 12 loci simultaneously. Using this system, we generated multiple lines that had null alleles at four distant loci. We also found that targeting multiple sites within a single locus can produce larger deletions, but the success of this depends on individual protospacers. To take advantage of homology-directed repair, we developed modular vectors to rapidly generate DNA donor plasmids to efficiently introduce DNA sequences encoding for fluorescent proteins at the 5' and 3' ends of gene coding regions. With regard to homology-directed repair experiments, we found that if the protospacer sequence remains on the DNA donor plasmid, then Cas9 cleaves the plasmid target as well as the genomic target. This can reduce the efficiency of introducing sequences into the genome. Furthermore, to ensure the generation of a null allele near the Cas9 cleavage site, we generated a homology plasmid harboring a "stop codon cassette" with downstream near-effortless genotyping.

Distinct types of short open reading frames are translated in plant cells

Genomes contain millions of short (<100 codons) open reading frames (sORFs), which are usually dismissed during gene annotation. Nevertheless, peptides encoded by such sORFs can play important biological roles, and their impact on cellular processes has long been underestimated. Here, we analyzed approximately 70,000 transcribed sORFs in the model plant Physcomitrella patens (moss). Several distinct classes of sORFs that differ in terms of their position on transcripts and the level of evolutionary conservation are present in the moss genome. Over 5000 sORFs were conserved in at least one of 10 plant species examined. Mass spectrometry analysis of proteomic and peptidomic data sets suggested that tens of sORFs located on distinct parts of mRNAs and long noncoding RNAs (lncRNAs) are translated, including conserved sORFs. Translational analysis of the sORFs and main ORFs at a single locus suggested the existence of genes that code for multiple proteins and peptides with tissue-specific expression. Functional analysis of four lncRNA-encoded peptides showed that sORFs-encoded peptides are involved in regulation of growth and differentiation in moss. Knocking out lncRNA-encoded peptides resulted in a decrease of moss growth. In contrast, the overexpression of these peptides resulted in a diverse range of phenotypic effects. Our results thus open new avenues for discovering novel, biologically active peptides in the plant kingdom.

The XPF-ERCC1 Complex Is Essential for Genome Stability and Is Involved in the Mechanism of Gene Targeting in Physcomitrella patens

The XPF-ERCC1 complex, a highly conserved structure-specific endonuclease, functions in multiple DNA repair pathways that are pivotal for maintaining genome stability, including nucleotide excision repair, interstrand crosslink repair, and homologous recombination. XPF-ERCC1 incises double-stranded DNA at double-strand/single-strand junctions, making it an ideal enzyme for processing DNA structures that contain partially unwound strands. Here, we have examined the role of the XPF-ERCC1 complex in the model bryophyte Physcomitrella patens which exhibits uniquely high gene targeting frequencies. We undertook targeted knockout of the Physcomitrella ERCC1 and XPF genes. Mutant analysis shows that the endonuclease complex is essential for resistance to UV-B and to the alkylating agent MMS, and contributes to the maintenance of genome integrity but is also involved in gene targeting in this model plant. Using different constructs we determine whether the function of the XPF-ERCC1 endonuclease complex in gene targeting was removal of 3' non-homologous termini, similar to SSA, or processing of looped-out heteroduplex intermediates. Interestingly, our data suggest a role of the endonuclease in both pathways and have implications for the mechanism of targeted gene replacement in plants and its specificities compared to yeast and mammalian cells.

POLQ plays a key role in the repair of CRISPR/Cas9-induced double-stranded breaks in the moss Physcomitrella patens

Double-stranded breaks can be repaired by different mechanisms such as homologous recombination (HR), classical nonhomologous end joining (C-NHEJ) and alternative end joining (Alt-EJ). Polymerase Q (POLQ) has been proposed to be the main factor involved in Alt-EJ-mediated DNA repair. Here we describe the role of POLQ in DNA repair and gene targeting in Physcomitrella patens. The disruption of the POLQ gene does not influence the genetic stability of P. patens nor its development. The polq mutant shows the same sensitivity as wild-type towards most of the genotoxic agents tested (ultraviolet (UV), methyl methanesulfonate (MMS) and cisplatin) with the notable exception of bleomycin for which it shows less sensitivity than the wild-type. Furthermore, we show that POLQ is involved in the repair of CRISPR-Cas9-induced double-stranded breaks in P. patens. We also demonstrate that POLQ is a potential competitor and/or inhibitor of the HR repair pathway. This finding has a consequence in terms of genetic engineering, as in the absence of POLQ the frequency of gene targeting is significantly increased and the number of clean two-sided HR-mediated insertions is enhanced. Therefore, the control of POLQ activity in plants could be a useful strategy to optimize the tools of genome engineering for plant breeding.

Efficient Editing of the Nuclear APT Reporter Gene in Chlamydomonas reinhardtii via Expression of a CRISPR-Cas9 Module

The clustered regularly interspaced short palindromic repeat/CRISPR-associated protein 9 (CRISPR/Cas9) technology is a versatile and useful tool to perform genome editing in different organisms ranging from bacteria and yeast to plants and mammalian cells. For a couple of years, it was believed that the system was inefficient and toxic in the alga Chlamydomonas reinhardtii. However, recently the system has been successfully implemented in this model organism, albeit relying mostly on the electroporation of ribonucleoproteins (RNPs) into cell wall deficient strains. This requires a constant source of RNPs and limits the application of the technology to strains that are not necessarily the most relevant from a biotechnological point of view. Here, we show that transient expression of the Streptococcus pyogenes Cas9 gene and sgRNAs, targeted to the single-copy nuclear apt9 gene, encoding an adenine phosphoribosyl transferase (APT), results in efficient disruption at the expected locus. Introduction of indels to the apt9 locus results in cell insensitivity to the otherwise toxic compound 2-fluoroadenine (2-FA). We have used agitation with glass beads and particle bombardment to introduce the plasmids carrying the coding sequences for Cas9 and the sgRNAs in a cell-walled strain of C. reinhardtii (CC-125). Using sgRNAs targeting exons 1 and 3 of apt9, we obtained disruption efficiencies of 3 and 30% on preselected 2-FA resistant colonies, respectively. Our results show that transient expression of Cas9 and a sgRNA can be used for editing of the nuclear genome inexpensively and at high efficiency. Targeting of the APT gene could potentially be used as a pre-selection marker for multiplexed editing or disruption of genes of interest.

CRISPR/Cas9-Mediated Knockout of Physcomitrella patens Phytochromes

Here we describe procedures for gene disruption and excision in Physcomitrella using CRISPR/Cas9 (clustered regularly interspaced short palindromic repeat/CRISPR-associated 9) methods, exemplarily targeting phytochrome (PHY) gene loci. Thereby double-strand breaks (DSBs) are induced using a single guide RNA (sgRNA) with the Cas9 nuclease, leading to insertions or deletions (indels) due to incorrect repair by the nonhomologous-end joining (NHEJ) mechanism. We also include protocols for excision of smaller genomic fragments or whole genes either with or without homologous recombination-assisted repair. The protocol can be adapted to target several loci simultaneously, thereby allowing the physiological analysis of phenotypes that would be masked by functional redundancy. In our particular case, multiple PHY gene knockouts would likely be valuable in understanding phytochrome functions in mosses and, perhaps, higher plants too. Target sites for site-directed induction of DSBs are predicted with the CRISPOR online-tool and are inserted in silico into sequence matrices for the design of sgRNA expression cassettes. The resulting DNAs are cloned into Gateway DONOR vectors and the respective expression plasmids used for moss cotransformation with a Cas9 expression plasmid and a selectable marker (either on a separate plasmid or on one of the other plasmids). After the selection process, genomic DNA is extracted and transformants are analyzed by PCR fingerprinting.

Cytoskeletal discoveries in the plant lineage using the moss Physcomitrella patens

Advances in cell biology have been largely driven by pioneering work in model systems, the majority of which are from one major eukaryotic lineage, the opisthokonts. However, with the explosion of genomic information in many lineages, it has become clear that eukaryotes have incredible diversity in many cellular systems, including the cytoskeleton. By identifying model systems in diverse lineages, it may be possible to begin to understand the evolutionary origins of the eukaryotic cytoskeleton. Within the plant lineage, cell biological studies in the model moss, Physcomitrella patens, have over the past decade provided key insights into how the cytoskeleton drives cell and tissue morphology. Here, we review P. patens attributes that make it such a rich resource for cytoskeletal cell biological inquiry and highlight recent key findings with regard to intracellular transport, microtubule-actin interactions, and gene discovery that promises for many years to provide new cytoskeletal players.

One-step generation of multiple gene knock-outs in the diatom Phaeodactylum tricornutum by DNA-free genome editing

Recently developed transgenic techniques to explore and exploit the metabolic potential of microalgae present several drawbacks associated with the delivery of exogenous DNA into the cells and its subsequent integration at random sites within the genome. Here, we report a highly efficient multiplex genome-editing method in the diatom Phaeodactylum tricornutum, relying on the biolistic delivery of CRISPR-Cas9 ribonucleoproteins coupled with the identification of two endogenous counter-selectable markers, PtUMPS and PtAPT. First, we demonstrate the functionality of RNP delivery by positively selecting the disruption of each of these genes. Then, we illustrate the potential of the approach for multiplexing by generating double-gene knock-out strains, with 65% to 100% efficiency, using RNPs targeting one of these markers and PtAureo1a, a photoreceptor-encoding gene. Finally, we created triple knock-out strains in one step by delivering six RNP complexes into Phaeodactylum cells. This approach could readily be applied to other hard-to-transfect organisms of biotechnological interest.

The manipulation of diatom genomes is essential for industrial applications based on their metabolic abilities. Here the authors present an efficient multiplex DNA-free gene editing method using CRISPR-Cas9 and counter-selectable markers.

PpSARK Regulates Moss Senescence and Salt Tolerance through ABA Related Pathway

Senescence-associated receptor-like kinase (SARK) family members in Arabidopsis, soybean, and rice are known to be positive regulators of leaf senescence. In the meantime, SARKs are extensively involved in stress response. However, their function and underlying molecular mechanism in stress responses in moss are not well known. Here, we investigated functional roles of SARK isolated from Physcomitrella patens (PpSARK) in salt stress response and senescence. PpSARK transcripts significantly accumulated under NaCl and abscisic acid (ABA) treatments, with higher expression in the moss gametophyte stage. Insertional gain-of-function mutants of PpSARK (PpSARKg) were more tolerant to salt stress and ABA than wild type (WT), whereas senescence of mutants was delayed during the protonema stage. Expression of stress-responsive genes in the ABA related pathway, such as PpABI3, PpABI5, PpPP2C, and PpLEA were significantly higher in PpSARKg and WT under salt stress conditions, suggesting that PpSARK might positively regulate salt tolerance via an ABA-related pathway. Endogenous ABA contents also increased 3-fold under salt stress conditions. These results indicate that PpSARK functions as a positive regulator in salt stress responses, while possibly functioning as a negative regulator in senescence in moss.

CRISPR/Cas9-mediated efficient targeted mutagenesis in grape in the first generation

The clustered regularly interspaced short palindromic repeats-associated protein 9 (CRISPR/Cas9) system is a powerful tool for editing plant genomes. Efficient genome editing of grape (Vitis vinifera) suspension cells using the type II CRISPR/Cas9 system has been demonstrated; however, it has not been established whether this system can be applied to get biallelic mutations in the first generation of grape. In this current study, we designed four guide RNAs for the VvWRKY52 transcription factor gene for using with the CRISPR/Cas9 system, and obtained transgenic plants via Agrobacterium-mediated transformation, using somatic embryos of the Thompson Seedless cultivar. Analysis of the first-generation transgenic plants verified 22 mutant plants of the 72 T-DNA-inserted plants. Of these, 15 lines carried biallelic mutations and seven were heterozygous. A range of RNA-guided editing events, including large deletions, were found in the mutant plants, while smaller deletions comprised the majority of the detected mutations. Sequencing of potential off-target sites for all four targets revealed no off-target events. In addition, knockout of VvWRKY52 in grape increased the resistance to Botrytis cinerea. We conclude that the CRISPR/Cas9 system allows precise genome editing in the first generation of grape and represents a useful tool for gene functional analysis and grape molecular breeding.

Genetic redundancy of senescence-associated transcription factors in Arabidopsis

Leaf senescence is a genetically programmed process that constitutes the last stage of leaf development, and involves massive changes in gene expression. As a result of the intensive efforts that have been made to elucidate the molecular genetic mechanisms underlying leaf senescence, 184 genes that alter leaf senescence phenotypes when mutated or overexpressed have been identified in Arabidopsis thaliana over the past two decades. Concurrently, experimental evidence on functional redundancy within senescence-associated genes (SAGs) has increased. In this review, we focus on transcription factors that play regulatory roles in Arabidopsis leaf senescence, and describe the relationships among gene duplication, gene expression level, and senescence phenotypes. Previous findings and our re-analysis demonstrate the widespread existence of duplicate SAG pairs and a correlation between gene expression levels in duplicate genes and senescence-related phenotypic severity of the corresponding mutants. We also highlight effective and powerful tools that are available for functional analyses of redundant SAGs. We propose that the study of duplicate SAG pairs offers a unique opportunity to understand the regulation of leaf senescence and can guide the investigation of the functions of redundant SAGs via reverse genetic approaches.

Deciphering the Evolution and Development of the Cuticle by Studying Lipid Transfer Proteins in Mosses and Liverworts

When plants conquered land, they developed specialized organs, tissues, and cells in order to survive in this new and harsh terrestrial environment. New cell polymers such as the hydrophobic lipid-based polyesters cutin, suberin, and sporopollenin were also developed for protection against water loss, radiation, and other potentially harmful abiotic factors. Cutin and waxes are the main components of the cuticle, which is the waterproof layer covering the epidermis of many aerial organs of land plants. Although the in vivo functions of the group of lipid binding proteins known as lipid transfer proteins (LTPs) are still rather unclear, there is accumulating evidence suggesting a role for LTPs in the transfer and deposition of monomers required for cuticle assembly. In this review, we first present an overview of the data connecting LTPs with cuticle synthesis. Furthermore, we propose liverworts and mosses as attractive model systems for revealing the specific function and activity of LTPs in the biosynthesis and evolution of the plant cuticle.

Precise Editing at DNA Replication Forks Enables Multiplex Genome Engineering in Eukaryotes

We describe a multiplex genome engineering technology in Saccharomyces cerevisiae based on annealing synthetic oligonucleotides at the lagging strand of DNA replication. The mechanism is independent of Rad51-directed homologous recombination and avoids the creation of double-strand DNA breaks, enabling precise chromosome modifications at single base-pair resolution with an efficiency of >40%, without unintended mutagenic changes at the targeted genetic loci. We observed the simultaneous incorporation of up to 12 oligonucleotides with as many as 60 targeted mutations in one transformation. Iterative transformations of a complex pool of oligonucleotides rapidly produced large combinatorial genomic diversity >10⁵. This method was used to diversify a heterologous β-carotene biosynthetic pathway that produced genetic variants with precise mutations in promoters, genes, and terminators, leading to altered carotenoid levels. Our approach of engineering the conserved processes of DNA replication, repair, and recombination could be automated and establishes a general strategy for multiplex combinatorial genome engineering in eukaryotes.

An intragenic mutagenesis strategy in Physcomitrella patens to preserve intron splicing

Gene targeting is a powerful reverse genetics technique for site-specific genome modification. Intrinsic homologous recombination in the moss Physcomitrella patens permits highly effective gene targeting, a characteristic that makes this organism a valuable model for functional genetics. Functional characterization of domains located within a multi-domain protein depends on the ability to generate mutants harboring genetic modifications at internal gene positions while maintaining the reading-frames of the flanking exons. In this study, we designed and evaluated different gene targeting constructs for targeted gene manipulation of sequences corresponding to internal domains of the DEFECTIVE KERNEL1 protein in Physcomitrella patens. Our results show that gene targeting-associated mutagenesis of introns can have adverse effects on splicing, corrupting the normal reading frame of the transcript. We show that successful genetic modification of internal sequences of multi-exon genes depends on gene-targeting strategies which insert the selection marker cassette into the 5' end of the intron and preserve the nucleotide sequence of the targeted intron.

Non-model model organisms

Model organisms are widely used in research as accessible and convenient systems to study a particular area or question in biology. Traditionally only a handful of organisms have been widely studied, but modern research tools are enabling researchers to extend the set of model organisms to include less-studied and more unusual systems. This Forum highlights a range of 'non-model model organisms' as emerging systems for tackling questions across the whole spectrum of biology (and beyond), the opportunities and challenges, and the outlook for the future.

Selective gene dosage by CRISPR-Cas9 genome editing in hexaploid Camelina sativa

In many plant species, gene dosage is an important cause of phenotype variation. Engineering gene dosage, particularly in polyploid genomes, would provide an efficient tool for plant breeding. The hexaploid oilseed crop Camelina sativa, which has three closely related expressed subgenomes, is an ideal species for investigation of the possibility of creating a large collection of combinatorial mutants. Selective, targeted mutagenesis of the three delta-12-desaturase (FAD2) genes was achieved by CRISPR-Cas9 gene editing, leading to reduced levels of polyunsaturated fatty acids and increased accumulation of oleic acid in the oil. Analysis of mutations over four generations demonstrated the presence of a large variety of heritable mutations in the three isologous CsFAD2 genes. The different combinations of single, double and triple mutants in the T3 generation were isolated, and the complete loss-of-function mutants revealed the importance of delta-12-desaturation for Camelina development. Combinatorial association of different alleles for the three FAD2 loci provided a large diversity of Camelina lines with various lipid profiles, ranging from 10% to 62% oleic acid accumulation in the oil. The different allelic combinations allowed an unbiased analysis of gene dosage and function in this hexaploid species, but also provided a unique source of genetic variability for plant breeding.

Plant genome editing with TALEN and CRISPR

Genome editing promises giant leaps forward in advancing biotechnology, agriculture, and basic research. The process relies on the use of sequence specific nucleases (SSNs) to make DNA double stranded breaks at user defined genomic loci, which are subsequently repaired by two main DNA repair pathways: non-homologous end joining (NHEJ) and homology directed repair (HDR). NHEJ can result in frameshift mutations that often create genetic knockouts. These knockout lines are useful for functional and reverse genetic studies but also have applications in agriculture. HDR has a variety of applications as it can be used for gene replacement, gene stacking, and for creating various fusion proteins. In recent years, transcription activator-like effector nucleases and clustered regularly interspaced palindromic repeats (CRISPR) and CRISPR associated protein 9 or CRISPR from Prevotella and Francisella 1 have emerged as the preferred SSNs for research purposes. Here, we review their applications in plant research, discuss current limitations, and predict future research directions in plant genome editing.

Hornworts: An Overlooked Window into Carbon-Concentrating Mechanisms

Hornworts are the only land plant lineage harboring a biophysical carbon-concentrating mechanism (CCM). Here, we argue that hornworts are a promising, yet currently overlooked, model system to study the evolution and genetic basis of CCMs. The results of such studies could have translational values toward engineering a CCM in crop plants.

Mitotic Spindle Assembly in Land Plants: Molecules and Mechanisms

In textbooks, the mitotic spindles of plants are often described separately from those of animals. How do they differ at the molecular and mechanistic levels? In this chapter, we first outline the process of mitotic spindle assembly in animals and land plants. We next discuss the conservation of spindle assembly factors based on database searches. Searches of >100 animal spindle assembly factors showed that the genes involved in this process are well conserved in plants, with the exception of two major missing elements: centrosomal components and subunits/regulators of the cytoplasmic dynein complex. We then describe the spindle and phragmoplast assembly mechanisms based on the data obtained from robust gene loss-of-function analyses using RNA interference (RNAi) or mutant plants. Finally, we discuss future research prospects of plant spindles.

Efficient and Heritable Targeted Mutagenesis in Mosses Using the CRISPR/Cas9 System

Targeted genome modification by RNA-guided nucleases derived from the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated nuclease 9 (Cas9) system has seen rapid development in many organisms, including several plant species. In the present study, we succeeded in introducing the CRISPR/Cas9 system into the non-model organism Scopelophila cataractae, a moss that exhibits heavy metal tolerance, and the model organism Physcomitrella patens Utilizing the process by which moss plants regenerate from protoplasts, we conducted targeted mutagenesis by expression of single-chain guide RNA (sgRNA) and Cas9 in protoplasts. Using this method, the acquisition rate of strains exhibiting phenotypic changes associated with the target genes was approximately 45-69%, and strains with phenotypic changes exhibited various insertion and deletion mutations. In addition, we report that our method is capable of multiplex targeted mutagenesis (two independent genes) and also permits the efficient introduction of large deletions (∼3 kbp). These results demonstrate that the CRISPR/Cas9 system can be used to accelerate investigations of bryology and land plant evolution.

Simple and Efficient Targeting of Multiple Genes Through CRISPR-Cas9 in Physcomitrella patens

Powerful genome editing technologies are needed for efficient gene function analysis. The CRISPR-Cas9 system has been adapted as an efficient gene-knock-out technology in a variety of species. However, in a number of situations, knocking out or modifying a single gene is not sufficient; this is particularly true for genes belonging to a common family, or for genes showing redundant functions. Like many plants, the model organism Physcomitrella patens has experienced multiple events of polyploidization during evolution that has resulted in a number of families of duplicated genes. Here, we report a robust CRISPR-Cas9 system, based on the codelivery of a CAS9 expressing cassette, multiple sgRNA vectors, and a cassette for transient transformation selection, for gene knock-out in multiple gene families. We demonstrate that CRISPR-Cas9-mediated targeting of five different genes allows the selection of a quintuple mutant, and all possible subcombinations of mutants, in one experiment, with no mutations detected in potential off-target sequences. Furthermore, we confirmed the observation that the presence of repeats in the vicinity of the cutting region favors deletion due to the alternative end joining pathway, for which induced frameshift mutations can be potentially predicted. Because the number of multiple gene families in Physcomitrella is substantial, this tool opens new perspectives to study the role of expanded gene families in the colonization of land by plants.