The sonication-assisted whisker method enables CRISPR-Cas9 ribonucleoprotein delivery to induce genome editing in rice.

CRISPR/Cas9-based genome editing represents an unprecedented potential for plant breeding. Unlike animal cells, plant cells contain a rigid cell wall, genome editing tool delivery into plant cells is thus challenging. In particular, the delivery of the Cas9-gRNA ribonucleoprotein (RNP) into plant cells is desired since the transgene insertion into the genome should be avoided for industrial applications in plants. In this study, we present a novel RNP delivery approach in rice. We applied the sonication-assisted whisker method, conventionally developed for DNA delivery in plants, for RNP delivery in rice. Combined with marker gene delivery, we successfully isolated OsLCYß genome-edited lines generated by RNPs. The calli and regenerated shoot of the OsLCYß mutant showed abnormal carotenoid accumulation. In addition, we also detected, although at a low frequency, genome editing events in rice calli cells by RNP delivery using the sonication-assisted whisker method without any additional. Therefore, the sonication-assisted whisker method could be an attractive way to create RNP-based genome-edited lines in plants.

LysM-mediated signaling in Marchantia polymorpha highlights the conservation of pattern-triggered immunity in land plants.

Pattern-recognition receptor (PRR)-triggered immunity (PTI) wards off a wide range of pathogenic microbes, playing a pivotal role in angiosperms. The model liverwort Marchantia polymorpha triggers defense-related gene expression upon sensing components of bacterial and fungal extracts, suggesting the existence of PTI in this plant model. However, the molecular components of the putative PTI in M. polymorpha and the significance of PTI in bryophytes have not yet been described. We here show that M. polymorpha has four lysin motif (LysM)-domain-containing receptor homologs, two of which, LysM-receptor-like kinase (LYK) MpLYK1 and LYK-related (LYR) MpLYR, are responsible for sensing chitin and peptidoglycan fragments, triggering a series of characteristic immune responses. Comprehensive phosphoproteomic analysis of M. polymorpha in response to chitin treatment identified regulatory proteins that potentially shape LysM-mediated PTI. The identified proteins included homologs of well-described PTI components in angiosperms as well as proteins whose roles in PTI are not yet determined, including the blue-light receptor phototropin MpPHOT. We revealed that MpPHOT is required for negative feedback of defense-related gene expression during PTI. Taken together, this study outlines the basic framework of LysM-mediated PTI in M. polymorpha and highlights conserved elements and new aspects of pattern-triggered immunity in land plants.

Zinc deficiency-induced defensin-like proteins are involved in the inhibition of root growth in Arabidopsis.

The depletion of cellular zinc (Zn) adversely affects plant growth. Plants have adaptation mechanisms for Zn-deficient conditions, inhibiting growth through the action of transcription factors and metal transporters. We previously identified three defensin-like (DEFL) proteins (DEFL203, DEFL206 and DEFL208) that were induced in Arabidopsis thaliana roots under Zn-depleted conditions. DEFLs are small cysteine-rich peptides involved in defense responses, development and excess metal stress in plants. However, the functions of DEFLs in the Zn-deficiency response are largely unknown. Here, phylogenetic tree analysis revealed that seven DEFLs (DEFL202-DEFL208) were categorized into one subgroup. Among the seven DEFLs, the transcripts of five (not DEFL204 and DEFL205) were upregulated by Zn deficiency, consistent with the presence of cis-elements for basic-region leucine-zipper 19 (bZIP19) or bZIP23 in their promoter regions. Microscopic observation of GFP-tagged DEFL203 showed that DEFL203-sGFP was localized to the apoplast and plasma membrane. Whereas a single mutation of the DEFL202 or DEFL203 genes only slightly affected root growth, defl202 defl203 double mutants showed enhanced root growth under all growth conditions. We also showed that the size of the root meristem was increased in the double mutants compared with the wild type. Our results suggest that DEFL202 and DEFL203 are redundantly involved in the inhibition of root growth under Zn-deficient conditions through a reduction in root meristem length and cell number.

Single-Cell RNA Sequencing of Arabidopsis Leaf Tissues Identifies Multiple Specialized Cell Types: Idioblast Myrosin Cells and Potential Glucosinolate-Producing Cells.

The glucosinolate-myrosinase defense system (GMDS), characteristic of Brassicales, is involved in plant defense. Previous single-cell transcriptomic analyses have reported the expression profiles of multiple GMDS-related cell types (i.e. myrosinase-rich myrosin idioblasts and multiple types of potential glucosinolate synthetic cells as well as a candidate S-cell for glucosinolate accumulation). However, differences in plant stages and cell-type annotation methods have hindered comparisons among studies. Here, we used the single-cell transcriptome profiles of extended Arabidopsis leaves and verified the distribution of previously used markers to refine the expression profiles of GMDS-associated cell types. Moreover, we performed beta-glucuronidase promoter assays to confirm the histological expression patterns of newly obtained markers for GMDS-associated candidates. As a result, we found a set of new specific reporters for myrosin cells and potential glucosinolate-producing cells.

Efficacy of the PRESERFLO MicroShunt and a Meta-Analysis of the Literature.

Background: Recent studies on the PRESERFLO MicroShunt suggest that it may be effective in lowering intraocular pressure (IOP); however, the number of studies on this device remains limited. Therefore, we assessed the efficacy of the PRESERFLO MicroShunt in patients with glaucoma and performed a meta-analysis of published results. Methods: Prospective study including all patients that underwent PRESERFLO MicroShunt surgery from 2018 onwards. Sub-analyses were performed for cataract-combined procedures. To compare our results, we performed a systematic review and meta-analysis. IOP, IOP-lowering medication and surgical complications reported in the retrieved studies were assessed. Results: A total of 72 eyes underwent PRESERFLO-implant surgery (59 as standalone procedure and 13 as cataract-combined procedure). No significant differences were found in IOP and IOP-lowering medication between both groups. The mean ± standard deviation IOP and IOP-lowering medications of both groups taken together declined from 21.72 ± 8.35 to 15.92 ± 8.54 mmHg (p < 0.001, 26.7% reduction) and 3.40 to 0.93 (p < 0.001, 72.6% reduction) at 1 year follow-up, respectively. Secondary surgeries were required in 19.4% of eyes, the majority (71.4%) within 6 months. The meta-analysis including 14 studies (totaling 1213 PRESERFLO MicroShunt surgeries) from the systematic review showed a mean preoperative IOP and IOP-lowering medication of 22.28 ± 5.38 and 2.97 ± 1.07, respectively. The three-years postoperative pooled mean was (weighted mean difference, 95% CI) 11.07 (10.27 [8.23−12.32], p < 0.001) mmHg and 0.91 (1.77 [1.26−2.28], p < 0.001) for IOP and IOP-lowering medication, respectively. The most common reported complication was hypotony (2−39%). Conclusion: The PRESERFLO MicroShunt is effective and safe in lowering IOP and the number of IOP-lowering medications.

Viscosity-aided electromechanical poration of cells for transfecting molecules.

Cell poration technologies offer opportunities not only to understand the activities of biological molecules but also to investigate genetic manipulation possibilities. Unfortunately, transferring large molecules that can carry huge genomic information is challenging. Here, we demonstrate electromechanical poration using a core-shell-structured microbubble generator, consisting of a fine microelectrode covered with a dielectric material. By introducing a microcavity at its tip, we could concentrate the electrical field with the application of electric pulses and generate microbubbles for electromechanical stimulation of cells. Specifically, the technology enables transfection with molecules that are thousands of kDa even into osteoblasts and Chlamydomonas, which are generally considered to be difficult to inject. Notably, we found that the transfection efficiency can be enhanced by adjusting the viscosity of the cell suspension, which was presumably achieved by remodeling of the membrane cytoskeleton. The applicability of the approach to a variety of cell types opens up numerous emerging gene engineering applications.

Atezolizumab plus bevacizumab versus lenvatinib or sorafenib in non-viral unresectable hepatocellular carcinoma: an international propensity score matching analysis.

BACKGROUND: A growing body of evidence suggests that non-viral hepatocellular carcinoma (HCC) might benefit less from immunotherapy. MATERIALS AND METHODS: We carried out a retrospective analysis of prospectively collected data from consecutive patients with non-viral advanced HCC, treated with atezolizumab plus bevacizumab, lenvatinib, or sorafenib, in 36 centers in 4 countries (Italy, Japan, Republic of Korea, and UK). The primary endpoint was overall survival (OS) with atezolizumab plus bevacizumab versus lenvatinib. Secondary endpoints were progression-free survival (PFS) with atezolizumab plus bevacizumab versus lenvatinib, and OS and PFS with atezolizumab plus bevacizumab versus sorafenib. For the primary and secondary endpoints, we carried out the analysis on the whole population first, and then we divided the cohort into two groups: non-alcoholic fatty liver disease (NAFLD)/non-alcoholic steatohepatitis (NASH) population and non-NAFLD/NASH population. RESULTS: One hundred and ninety patients received atezolizumab plus bevacizumab, 569 patients received lenvatinib, and 210 patients received sorafenib. In the whole population, multivariate analysis showed that treatment with lenvatinib was associated with a longer OS [hazard ratio (HR) 0.65; 95% confidence interval (CI) 0.44-0.95; P = 0.0268] and PFS (HR 0.67; 95% CI 0.51-0.86; P = 0.002) compared to atezolizumab plus bevacizumab. In the NAFLD/NASH population, multivariate analysis confirmed that lenvatinib treatment was associated with a longer OS (HR 0.46; 95% CI 0.26-0.84; P = 0.0110) and PFS (HR 0.55; 95% CI 0.38-0.82; P = 0.031) compared to atezolizumab plus bevacizumab. In the subgroup of non-NAFLD/NASH patients, no difference in OS or PFS was observed between patients treated with lenvatinib and those treated with atezolizumab plus bevacizumab. All these results were confirmed following propensity score matching analysis. By comparing patients receiving atezolizumab plus bevacizumab versus sorafenib, no statistically significant difference in survival was observed. CONCLUSIONS: The present analysis conducted on a large number of advanced non-viral HCC patients showed for the first time that treatment with lenvatinib is associated with a significant survival benefit compared to atezolizumab plus bevacizumab, in particular in patients with NAFLD/NASH-related HCC.

DNA-functionalized colloidal crystals for macromolecular encapsulation.

Novel DNA-based structures with the ability to encapsulate nanoscale molecules, such as proteins, can be applied to a wide range of areas, including reaction fields and micro/nano drug carriers. DNA-functionalized nanoparticle (DNA-NP) colloidal crystals have emerged as a new class of programmable DNA-based structures harboring metal nanoparticles with improved mechanical properties. The encapsulation of guest molecules into empty spaces in lattice structures is theoretically possible. However, due to the lack of a strategy for versatile encapsulation of guest molecules, the feasibility of nanoscale encapsulation by DNA-NP crystals is unclear. In this study, we developed DNA-functionalized gold nanoparticle (DNA-AuNP) crystals with tunable interparticle spacing for molecular encapsulation. We demonstrated that the modification of DNA-AuNP crystals with functional moieties, that is, biotin molecules, was effective in retaining molecules in the crystals. The crystallinities before and after encapsulation of the molecules were confirmed using small-angle X-ray scattering. We also succeeded in encapsulating CRISPR/Cas9 ribonucleoproteins into DNA-AuNP crystals by harnessing their affinity for target molecules. These findings demonstrated the potential use of metal-DNA hybrid crystals as carriers for direct protein delivery via biolistic bombardment. Thus, this study provides an attractive strategy for creating a new class of DNA-based structures for macromolecular encapsulation, and an alternative research direction toward colloidal crystal engineering using DNA.

Candidate Tidal Disruption Event AT2019fdr Coincident with a High-Energy Neutrino.

The origins of the high-energy cosmic neutrino flux remain largely unknown. Recently, one high-energy neutrino was associated with a tidal disruption event (TDE). Here we present AT2019fdr, an exceptionally luminous TDE candidate, coincident with another high-energy neutrino. Our observations, including a bright dust echo and soft late-time x-ray emission, further support a TDE origin of this flare. The probability of finding two such bright events by chance is just 0.034%. We evaluate several models for neutrino production and show that AT2019fdr is capable of producing the observed high-energy neutrino, reinforcing the case for TDEs as neutrino sources.

Comprehensive analysis of peptide-coding genes and initial characterization of an LRR-only microprotein in Marchantia polymorpha.

In the past two decades, many plant peptides have been found to play crucial roles in various biological events by mediating cell-to-cell communications. However, a large number of small open reading frames (sORFs) or short genes capable of encoding peptides remain uncharacterized. In this study, we examined several candidate genes for peptides conserved between two model plants: Arabidopsis thaliana and Marchantia polymorpha. We examined their expression pattern in M. polymorpha and subcellular localization using a transient assay with Nicotiana benthamiana. We found that one candidate, MpSGF10B, was expressed in meristems, gemma cups, and male reproductive organs called antheridiophores. MpSGF10B has an N-terminal signal peptide followed by two leucine-rich repeat (LRR) domains and was secreted to the extracellular region in N. benthamiana and M. polymorpha. Compared with the wild type, two independent Mpsgf10b mutants had a slightly increased number of antheridiophores. It was revealed in gene ontology enrichment analysis that MpSGF10B was significantly co-expressed with genes related to cell cycle and development. These results suggest that MpSGF10B may be involved in the reproductive development of M. polymorpha. Our research should shed light on the unknown role of LRR-only proteins in land plants.

Nonalcoholic steatohepatitis in hepatocarcinoma: new insights about its prognostic role in patients treated with lenvatinib.

BACKGROUND: Hepatocellular carcinoma (HCC) treatment remains a big challenge in the field of oncology. The liver disease (viral or not viral) underlying HCC turned out to be crucial in determining the biologic behavior of the tumor, including its response to treatment. The aim of this analysis was to investigate the role of the etiology of the underlying liver disease in survival outcomes. PATIENTS AND METHODS: We conducted a multicenter retrospective study on a large cohort of patients treated with lenvatinib as first-line therapy for advanced HCC from both Eastern and Western institutions. Univariate and multivariate analyses were performed. RESULTS: Among the 1232 lenvatinib-treated HCC patients, 453 (36.8%) were hepatitis C virus positive, 268 hepatitis B virus positive (21.8%), 236 nonalcoholic steatohepatitis (NASH) correlate (19.2%) and 275 had other etiologies (22.3%). The median progression-free survival (mPFS) was 6.2 months [95% confidence interval (CI) 5.9-6.7 months] and the median overall survival (mOS) was 15.8 months (95% CI 14.9-17.2 months). In the univariate analysis for OS NASH-HCC was associated with longer mOS [22.2 versus 15.1 months; hazard ratio (HR) 0.69; 95% CI 0.56-0.85; P = 0.0006]. In the univariate analysis for PFS NASH-HCC was associated with longer mPFS (7.5 versus 6.5 months; HR 0.84; 95% CI 0.71-0.99; P = 0.0436). The multivariate analysis confirmed NASH-HCC (HR 0.64; 95% CI 0.48-0.86; P = 0.0028) as an independent prognostic factor for OS, along with albumin-bilirubin (ALBI) grade, extrahepatic spread, neutrophil-to-lymphocyte ratio, portal vein thrombosis, Eastern Cooperative Oncology Group (ECOG) performance status and alpha-fetoprotein. An interaction test was performed between sorafenib and lenvatinib cohorts and the results highlighted the positive predictive role of NASH in favor of the lenvatinib arm (P = 0.0047). CONCLUSION: NASH has been identified as an independent prognostic factor in a large cohort of patients with advanced HCC treated with lenvatinib, thereby suggesting the role of the etiology in the selection of patients for tyrosine kinase treatment. If validated, this result could provide new insights useful to improve the management of these patients.

Soft gamma rays from low accreting supermassive black holes and connection to energetic neutrinos.

The Universe is filled with a diffuse background of MeV gamma-rays and PeV neutrinos, whose origins are unknown. Here, we propose a scenario that can account for both backgrounds simultaneously. Low-luminosity active galactic nuclei have hot accretion flows where thermal electrons naturally emit soft gamma rays via Comptonization of their synchrotron photons. Protons there can be accelerated via turbulence or reconnection, producing high-energy neutrinos via hadronic interactions. We demonstrate that our model can reproduce the gamma-ray and neutrino data. Combined with a contribution by hot coronae in luminous active galactic nuclei, these accretion flows can explain the keV - MeV photon and TeV - PeV neutrino backgrounds. This scenario can account for the MeV background without non-thermal electrons, suggesting a higher transition energy from the thermal to nonthermal Universe than expected. Our model is consistent with X-ray data of nearby objects, and testable by future MeV gamma-ray and high-energy neutrino detectors.

Higher Stomatal Density Improves Photosynthetic Induction and Biomass Production in Arabidopsis Under Fluctuating Light.

Stomatal density (SD) is closely associated with photosynthetic and growth characteristics in plants. In the field, light intensity can fluctuate drastically within a day. The objective of the present study is to examine how higher SD affects stomatal conductance (g s ) and CO2 assimilation rate (A) dynamics, biomass production and water use under fluctuating light. Here, we compared the photosynthetic and growth characteristics under constant and fluctuating light among three lines of Arabidopsis thaliana (L.): the wild type (WT), STOMAGEN/EPFL9-overexpressing line (ST-OX), and EPIDERMAL PATTERNING FACTOR 1 knockout line (epf1). ST-OX and epf1 showed 268.1 and 46.5% higher SD than WT (p < 0.05). Guard cell length of ST-OX was 10.0% lower than that of WT (p < 0.01). There were no significant variations in gas exchange parameters at steady state between WT and ST-OX or epf1, although these parameters tended to be higher in ST-OX and epf1 than WT. On the other hand, ST-OX and epf1 showed faster A induction than WT after step increase in light owing to the higher g s under initial dark condition. In addition, ST-OX and epf1 showed initially faster g s induction and, at the later phase, slower g s induction. Cumulative CO2 assimilation in ST-OX and epf1 was 57.6 and 78.8% higher than WT attributable to faster A induction with reduction of water use efficiency (WUE). epf1 yielded 25.6% higher biomass than WT under fluctuating light (p < 0.01). In the present study, higher SD resulted in faster photosynthetic induction owing to the higher initial g s . epf1, with a moderate increase in SD, achieved greater biomass production than WT under fluctuating light. These results suggest that higher SD can be beneficial to improve biomass production in plants under fluctuating light conditions.

Author Correction: Rapid breeding of parthenocarpic tomato plants using CRISPR/Cas9.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Hidden Cores of Active Galactic Nuclei as the Origin of Medium-Energy Neutrinos: Critical Tests with the MeV Gamma-Ray Connection.

Mysteries about the origin of high-energy cosmic neutrinos have deepened by the recent IceCube measurement of a large diffuse flux in the 10-100 TeV range. Based on the standard disk-corona picture of active galactic nuclei (AGN), we present a phenomenological model enabling us to systematically calculate the spectral sequence of multimessenger emission from the AGN coronae. We show that protons in the coronal plasma can be stochastically accelerated up to PeV energies by plasma turbulence, and find that the model explains the large diffuse flux of medium-energy neutrinos if the cosmic rays carry only a few percent of the thermal energy. We find that the Bethe-Heitler process plays a crucial role in connecting these neutrinos and cascaded MeV gamma rays, and point out that the gamma-ray flux can even be enhanced by the reacceleration of secondary pairs. Critical tests of the model are given by its prediction that a significant fraction of the MeV gamma-ray background correlates with ∼10 TeV neutrinos, and nearby Seyfert galaxies including NGC 1068 are promising targets for IceCube, KM3Net, IceCube-Gen2, and future MeV gamma-ray telescopes.

Manganese Treatment Alleviates Zinc Deficiency Symptoms in Arabidopsis Seedlings.

Plant phenotypes caused by mineral deficiencies differ depending on growth conditions. We recently reported that the growth of Arabidopsis thaliana was severely inhibited on MGRL-based zinc (Zn)-deficient medium but not on Murashige-Skoog-based Zn-deficient medium. Here, we explored the underlying reason for the phenotypic differences in Arabidopsis grown on the different media. The root growth and chlorophyll contents reduced by Zn deficiency were rescued by the addition of extra manganese (Mn) during short-term growth (10 or 14 d). However, this treatment did not affect the growth recovery after long-term growth (38 d). To investigate the reason for plant recovery from Zn deficiency, we performed the RNA-seq analysis of the roots grown on the Zn-basal medium and the Zn-depleted medium with/without additional Mn. Principal component analysis of the RNA-seq data showed that the gene expression patterns of plants on the Zn-basal medium were similar to those on the Zn-depleted medium with Mn, whereas those on the Zn-depleted medium without Mn were different from the others. The expression of several transcription factors and reactive oxygen species (ROS)-related genes was upregulated in only plants on the Zn-depleted medium without Mn. Consistent with the gene expression data, ROS accumulation in the roots grown on this medium was higher than those grown in other conditions. These results suggest that plants accumulate ROS and reduce their biomass under undesirable growth conditions, such as Zn depletion. Taken together, this study shows that the addition of extra Mn to the Zn-depleted medium induces transcriptional changes in ROS-related genes, thereby alleviating short-term growth inhibition due to Zn deficiency.

Two Y-chromosome-encoded genes determine sex in kiwifruit.

Dioecy, the presence of male and female individuals, has evolved independently in multiple flowering plant lineages1-3. Although theoretical models for the evolution of dioecy, such as the 'two-mutations' model, are well established4,5, little is known about the specific genes determining sex and their evolutionary history3. Kiwifruit, a major tree crop consumed worldwide, is a dioecious species. In kiwifruit we previously identified a Y-encoded sex-determinant candidate gene acting as the suppressor of feminization (SuF), named Shy Girl (SyGI)6. Here, we identify a second Y-encoded sex-determinant that we named Friendly Boy (FrBy), which exhibits strong expression in tapetal cells. Gene-editing and complementation analyses in Arabidopsis thaliana and Nicotiana tabacum indicated that FrBy acts for the maintenance of male (M) functions, independently of SyGI, and that these functions are conserved across angiosperm species. We further characterized the genomic architecture of the small (<1 megabase pairs (Mb)) male-specific region of the Y chromosome (MSY), which harbours only two genes expressed extensively in developing gynoecia and androecia, respectively: SyGI and FrBy. Re-sequencing of the genome of a natural hermaphrodite kiwifruit revealed that this individual is genetically male but carries deletion(s) of parts of the Y chromosome, including SyGI. Additionally, expression of FrBy in female kiwifruit resulted in hermaphrodite plants. These results clearly indicate that Y-encoded SyGI and FrBy act independently as the SuF and M factors in kiwifruit, respectively, and provide insight into not only the evolutionary path leading to a two-factor sex-determination system, but also a new breeding approach for dioecious species.

Chromatin-mediated feed-forward auxin biosynthesis in floral meristem determinacy.

In flowering plants, the switch from floral stem cell maintenance to gynoecium (female structure) formation is a critical developmental transition for reproductive success. In Arabidopsis thaliana, AGAMOUS (AG) terminates floral stem cell activities to trigger this transition. Although CRABS CLAW (CRC) is a direct target of AG, previous research has not identified any common targets. Here, we identify an auxin synthesis gene, YUCCA4 (YUC4) as a common direct target. Ectopic YUC4 expression partially rescues the indeterminate phenotype and cell wall defects that are caused by the crc mutation. The feed-forward YUC4 activation by AG and CRC directs a precise change in chromatin state for the shift from floral stem cell maintenance to gynoecium formation. We also showed that two auxin-related direct CRC targets, YUC4 and TORNADO2, cooperatively contribute to the termination of floral stem cell maintenance. This finding provides new insight into the CRC-mediated auxin homeostasis regulation for proper gynoecium formation.

Efficient CRISPR/Cas9-based genome editing and its application to conditional genetic analysis in Marchantia polymorpha.

Marchantia polymorpha is one of the model species of basal land plants. Although CRISPR/Cas9-based genome editing has already been demonstrated for this plant, the efficiency was too low to apply to functional analysis. In this study, we show the establishment of CRISPR/Cas9 genome editing vectors with high efficiency for both construction and genome editing. Codon optimization of Cas9 to Arabidopsis achieved over 70% genome editing efficiency at two loci tested. Systematic assessment revealed that guide sequences of 17 nt or shorter dramatically decreased this efficiency. We also demonstrated that a combinatorial use of this system and a floxed complementation construct enabled conditional analysis of a nearly essential gene. This study reports that simple, rapid, and efficient genome editing is feasible with the series of developed vectors.

CRISPR/Cas9-Based Genome Editing of Transcription Factor Genes in Marchantia polymorpha.

Plant transcription factors (TFs) belong to a wide variety of gene families. The systematic and rapid establishment of knockout lines of TF genes is critical for functional genetics. Genome engineering techniques for dissecting out the molecular function of TFs have been dramatically improved by CRISPR/Cas9-based genome editing technology. In the CRISPR/Cas9 system, Cas9 functions as a Cas9-gRNA ribonucleoprotein complex and uses its DNA endonuclease activity to induce the cleavage of the genome, which is targeted by gRNA. Double-strand breaks sometimes induce insertions and deletions at the target site, leading to frameshift mutations of TF genes. In this chapter, we describe a detailed protocol for the targeted mutagenesis of TFs using the CRISPR/Cas9 system, specifically in the case of Marchantia polymorpha, an emerging model plant for functional genomics. The CRISPR/Cas9 system is highly versatile for targeting genomic sequences because only an alteration of the gRNA sequence is needed to change target sequences. The labor and cost required to establish genome-edited lines are low enough that multiple mutants of TF genes can be generated in one laboratory. The CRISPR/Cas9-based genome editing technique consists of four steps: (1) gRNA design; (2) vector construction; (3) transformation; and (4) isolation of genome-edited lines. This manuscript focuses mainly on the strategy of gRNA design, the workflow for off-target searches, and the selection and identification of genome-edited lines by genotyping. Although we describe a protocol for M. polymorpha, the basic strategy of generating genome-edited lines of TF genes should be applicable widely to other plants.