• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    CRISPR/Cas9: Development and Application in Rice Breeding

    2019-10-08 06:45:30FernandoMatasRomeroAndrGaticaArias
    Rice Science 2019年5期

    Fernando Matías Romero, Andrés Gatica-Arias

    CRISPR/Cas9: Development and Application in Rice Breeding

    Fernando Matías Romero, Andrés Gatica-Arias

    (,,,)

    Rice (L.) is an important staple food crop worldwide due to its adaptability to different environmental conditions. Because of its great economic and social importance, there is a constant requirement for new varieties with improved agronomic characteristics, such as tolerance to different biotic (such as bacterium, fungus, insect and virus) and abiotic stresses (such as salinity, drought and temperature), higher yield and better organoleptic and nutritional value. Among the new genome editing technologies, the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas) (CRISPR/Cas) system allows precise and specific edition in a targeted genome region. It is one of the most frequently used techniques for the study of the function of new genes and for the development of mutant lines with enhanced tolerance to biotic and abiotic stresses, herbicide resistance or improved yield. The wide varieties of applications for this technology include simple non-homologous end joining, homologous recombination, gene replacement, and base editing. In this review, we analyzed how some of these applications have been used in rice cultivars to obtain rice varieties better adapted to current environmental conditions and market requirements.

    rice; CRISPR; Cas; genome editing; biotic stress; abiotic stress; yield; off-target; mutation

    Rice (L.) is one of the most important cereal crops in the world. The major area of production is Asia, with more than 6.6 billion tonnes produced in 2016, followed by America (36 million tonnes) and Africa (32 million tonnes)(FAOSTAT, http://www.fao.org/faostat/en/#home). One of the advantages of rice is that it can be grown in a wide range of environments, even in areas not suitable for other crops (Breviario and Genga, 2013). However, because of climate change, biotic stress and stressful environmental conditions such as drought, salinity and extreme temperatures have increased all over the world, affecting rice production and causing enormous losses. In addition, the world’s population is projected to reach 9.1 billion by the year 2050, increasing by 34%. According to the Food and Agriculture Organization of the United Nations (FAO), the global demand for agricultural products will increase by about 70% by 2050. Annual cereal production need to increase from 2.1 billion today to about 3.0 billion tonnes, in order to feed the projected world population.

    Based on these, the main concern in current agricultural research is to improve plant tolerance to biotic and abiotic stresses and increase productivity (Delorge et al, 2014). In recent years, new precise genome editing techniques have been developed and replaced previous techniques such as random mutagenesis, naturally occurring mutations and classical breeding techniques because these methodologies are time-consuming and take too long to obtain individuals with the desired phenotype. New techniques have become important tools for plant science and plant molecular breeding (Endo et al, 2018). Genome editing allows the introduction of deletions, insertions or base substitutions by causing damage, double-strand breaks (DSBs), in targeted deoxyribonucleic acid (DNA) (Fig. 1-A). In order to repair this damage, plant cells follow two strategies. The preferred repair pathway in higher plants is non-homologous end-joining (NHEJ) (Fig. 1-B), which mainly causes insertions or deletions (InDels) and can result in frameshift mutations. A second repair pathway, homologous recombination (HR) (Fig. 1-C), occurs when a template with homologous sequence surrounding the DSB is available and is used for DSB repair, resulting in gene replacement (Endo et al, 2018). Regardless of the repair mechanism used, DSBs are induced by sequence-specific nucleases (SSNs), including zinc finger nucleases (ZFNs) (Pabo et al, 2001) and transcription activator-like effector nucleases (TALENs) (Moscou and Bogdanove, 2009). Both SSNs are artificial bipartite enzymes that consist of a DNA-binding domain and thenuclease domain. In order to target the nuclease to the desired DNA sequence, a modular DNA-binding domain must be designed and assembled. Although the design and construction of large modular proteins are both laborious and expensive (Belhaj et al, 2015), they have been used for molecular breeding with relative success (Sprink et al, 2015).

    The newest and most widely used genome editing technique is based on clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated (Cas) (CRISPR/Cas). This technique is different from ZFN and TALEN in terms of the DNA-binding system. Cas9 can be targeted to a specific genomic sequence by an easily engineered 20 base pair (bp) RNA guide sequence that binds to its DNA target by base-pairing. The use of different genome editing techniques in rice from 2013 until early 2018 was recently reviewed (Mishra et al, 2018). However, in this review, we will focus on the use of CRISPR/Cas technology to improve desirable traits in rice, such as tolerance to biotic and abiotic stresses as well as different parameters involved in determining yield (Table 1). In addition, we will discuss the pre-requisites, limitations and advantages of this technique in rice.

    Fig. 1. CRISPR/Cas9 system components and pathways.

    A, Components required for genome editing using the CRISPR/Cas9 system: a DNA endonuclease (the most commonly used is the Cas9 protein from) and a customizable single guide RNA (sgRNA). B, Non-homologous end-joining (NHEJ) repair pathway. C, Homologous recombination (HR) repair pathway.

    PAM, Protospacer adjacent motif.

    NHEJ, Non-homologous end-joining repair pathway; HR, Homologous recombination repair pathway; PEG, Polyethylene glycol.

    CRISPR/Cas system

    CRISPR/Cas is a type of adaptive immune system, which degrades exogenous DNA, discovered by accident in 1987 in(Ishino et al, 1987). Later it was reported in other bacteria including,and(Nakata et al, 1989; Hermans et al, 1991). It was firstly reported in archaea by Mojica et al (1993). But it was not popular until 2012 when the potential to exploit this system for genome editing was suggested (Jinek et al, 2012).

    The immune memory of this system is stored in the form of spacer sequences from foreign genomes inserted into CRISPR arrays (Koonin et al, 2017). These spacers, along with Cas proteins, act as a surveillance system to recognize and degrade foreign nucleic acids. This process has three stages. The first stage, called adaptation, immunization or spacer acquisition, consists of the recognition and integration of foreign DNA (spacers) into the CRISPR locus. The sequence of viral or plasmid DNA in the spacer is called the protospacer. Usually there is a short conserved sequence in the immediate vicinity of the protospacer, referred to as the protospacer adjacent motif (PAM). The second stage consists of the expression of the system. A primary transcript (pre-crRNA) is transcribed from the CRISPR locus and processed into small CRISPR RNAs (crRNA). In the final stage, called interference or immunity, crRNAs, along with trans-activating crRNA (tracrRNA), form a ribonucleoprotein complex with Cas proteins. This complex recognizes foreign DNA by base pairing and then degrades it (Bhaya et al, 2011).

    CRISPR/Cas system has been engineered and used for genome editing in wide varieties of organisms such as mammalian cells, bacteria, fungi and plants (Jiang W Y et al, 2013; Krappmann, 2017). Two components are required for genome editing using this system: a DNA endonuclease (the most commonly used is the Cas9 protein from) and a customizable single guide RNA (sgRNA) (Fig. 1-A). Cas9 has a bi-lobed architecture, with a large globular recognition lobe (REC) and a small nuclease lobe (NUC) with two nuclease domains, RuvC and HNH, which each cut a specific DNA strand. The sgRNA is a small non-coding RNA that is the fusion of crRNA and tracrRNA (Jinek et al, 2012). CRISPR/Cas system can potentially be used to edit any sequence in the genome of any organism. The only requisite is that the target sequence must be adjacent to a PAM sequence. Cas9 fromrecognizes the sequence 5′-NGG-3′ as PAM. However, there are other Cas proteins with different PAM sequences, such as Cpf1 (or Cas12) that recognizes the sequence 5′-TTTN-3′ or 5′-TTN-3′ as PAM. Moreover, new Cas9 variants (VQR, EQR and VRER) have been developed to recognize alternative PAMs, thereby increasing the possibilities for modifying any target sequence in a genome (Anders et al, 2016). Use of CRISPR/Cas9 system in plants (rice and wheat) was firstly reported by Shan et al (2013). Since then, there have been several reports of successful plant genome editing using this system. One important feature of the system is the ability to act in-trans on the respective target, allowing the separation of the mutagenesis agent (CRISPR protein/gene and sgRNA) from the sequence with the desired modification (Ricroch et al, 2017). In order to do this, it is necessary to obtain stable production of CRISPR/Cas-mutated lines without CRISPR/Cas expression cassettes in the final plants. This can be achieved by several methods such as genetic segregation, transient expression of the CRISPR components, or transformation by ribonucleoprotein complexes (Pyott et al, 2016; Zhang et al, 2016; Liang et al, 2017). The easiest way to achieve this in rice is the genetic segregation of transgene in the next generation, which is the most common method to obtain transgene-free plants. However, transformation with ribonucleoprotein complexes has been also achieved in rice (Woo et al, 2015).

    Screening of CRISPR/Cas-generated mutant plants

    The increasing interest in this technology has led to the development and/or adaptation of new techniques for screening CRISPR/Cas-generated mutant plants. Chen et al (2018) developed a rapid, cost-effective, high-throughput mutant screening protocol based on Illumina sequencing followed by high-resolution melting (HRM) analysis. First, pooled DNA samples are amplified using primers targeting the sgRNA region, and the PCR products are sequenced. Since nucleotide variant frequency (NVF) is a measure of the frequency of abnormal nucleotides relative to the wild type, samples with high NVF values are considered likely to contain mutants. In this way, it is possible to detect one mutant plant in a 42-plant pooled sample. HRM analysis was then used to identify individual mutant line within the pooled samples. Liu et al (2019) developed a high throughput method to genotype mutant rice plants, based on protocols described earlier (Bell et al, 2014). Briefly, DNA samples are arranged in 96-well microplates, and two PCR steps are performed: the first one with specific primers designed to amplify the expected mutation region, and the second one with general primers coupled to barcodes in order to identify each sample. These PCR products are pooled together and sequenced by any next generation technique. Thus, the sequences generated are analyzed by the pipeline high-throughput tracking of mutation (Hi-Tom), which allows the extraction of primers and barcodes, filters low quality sequences, compares the sequences with the reference genome (wild type) and reports any mutation detected in each sample. Along the same lines, Peng et al (2018) developed a simple high-throughput qPCR-based method that can distinguish wild types from mutant samples and allow the identification of heterozygous and homozygous mutations in several plant species, including,,and. This technique uses two differently labelled probes for the detection of the same PCR product with one probe located outside the mutation position and the other within the mutation position. In this way, any mutation at the target site will block the binding of the second probe. Using qPCR, wild type, heterozygous and homozygous samples are easily differentiated.

    In rice, two new techniques for the identification of CRISPR/Cas-induced mutants have been reported. One uses non-invasive near-infrared (NIR) hyperspectral imaging. Radiation from light sources in the NIR region provides a unique spectral signature of organic molecules. Because modifications in DNA affect this unique spectral signature, wild type and mutant seeds can be differentiated (Feng et al, 2017). The other new technique combines CRISPR/Cas with the RNA interference (RNAi) technique, and takes advantage of the natural resistance of rice plants to the herbicide benzaton (Lu H P et al, 2017). In this technique, an RNAi expression element targeting the gene responsible for benzaton resistance is incorporated into a CRISPR/Cas9 vector. After-mediated transformation, tillers of T0plants are separated and grown as two subplants, one of which is sprayed with benzaton. Then, genome editing is analyzed in tolerant and susceptible plants. Plants carrying the RNAi expression element along with the CRISPR/Cas9 cassette show the mutated genotype (sensitivity to the herbicide). Thus, bentazon susceptibility is 100% correlated with the targeted mutations in T0plants (Lu H P et al, 2017).

    Pioneer reports of rice genome editing by CRISPR/Cas

    CRISPR/Cas technology for genome editing of rice and other plants (, tobacco, sorghum and wheat) was first reported in 2013 (Jiang W Z et al, 2013; Miao et al, 2013; Xie and Yang, 2013). Shan et al (2013) reported that two sgRNAs targeting the rice phytoene desaturase geneare designed and an efficiency of mutagenesis (15%)is obtained in transformed rice (Nipponbare) protoplasts. Nipponbare calliis bombarded with Cas9 plasmid and sgRNA expression plasmids designed to cleave eitheror(). Mutations inare identified in 9 out of 96 independent transgenic plants (9.4%), and mutations inare detected in 7 out of 98 transgenic plants (7.1%). HR repair system was evaluated using a single-stranded oligo introduced intoduring protoplast transformation, and2 out of 29 colonies incorporate the designed oligo in the targeted site. Miao et al (2013)used the reporter genewhich restores β-glucuronidase (GUS) activity through single strand annealing upon DNA cleavage between the repeat regions in thereporter. A sgRNA and a dual-crRNA:tracrRNA targeting this region are designed. Restoration of GUS activity was evaluated in rice calli from the Nipponbare. Both the sgRNA and dual-crRNA:tracrRNA are effective in this assay. Furthermore, the CRISPR/Cas technology was used to disrupt endogenous genes that allow easy identification of mutant plants. The() mutant, which shows a pale green leaf phenotype, and thegene mutant were chosen because bothmutants exhibit a pronounced tiller-spreading phenotype after the tillering stage. In this analysis, a high mutation efficiency ranging from 83% to 91% is observed in transformed lines of rice variety Kitaake (Miao et al, 2013). Other works demonstrated the ability to edit the rice genome using CRISPR/Cas technology by targeting genes involved in defence responses. Jiang W Z et al (2013) designed sgRNAs to edit the genesand, which are involved in resistance to bacterial blight caused by. After protoplast transformation, edition of the target sites was confirmed by PCR and sequencing. Xie and Yang (2013) used three sgRNAs targeting the rice() gene, a negative regulator of the plant defence response. Protoplasts from the rice variety Nipponbare are transformed and efficiencies ranging from 2% to 10% for different sgRNAs are observed.

    Since the publications of these reports validate the use of CRISPR/Cas genome editing in rice, there have been many studies implementing it for different purposes (Table 1).

    Pre-requisites and limitations of CRISPR/ Cas9 technology for rice improvement

    Genome editing in plants requires the introduction of the editing system inside the plant cells. In recent years, there has been a remarkable breakthrough in genetic transformation, and the introduction of foreign DNA is now possible for most major crops, including rice (Jung and Seo, 2017). However, since these techniques involve tissue culture, which is difficult in some rice varieties, especially, the application of genome editing could be limited in important commercial varieties. Moreover, the development and optimization of tissue culture protocols are labor-intensive and time-consuming procedures. Tissue culture can lead to somaclonal variations that can compromise the overall of regenerated plants (Kaeppler et al, 2000; Sarmast, 2016). Minimizing or eliminating the tissue culture requirement would improve the efficiency of crop transformation and genome editing (Jung and Seo, 2017).

    Delivery of CRISPR/Cas elements (Cas9 and sgRNAs) is another step during genome editing in plants. In order to express CRISPR/Cas elements and select genome-edited plants, vectors containing all of the elements and a selectable marker are delivered into the cells. These genes are usually integrated into the cell genome to allow the expression of Cas9 and the sgRNAs. Once the genome-edited plants are identified, these elements can be removed by segregation in the next generation or by site-specific recombination. An alternative is the transient expression of these elements without incorporation into the host genome. Transient expression could be sufficient for genome editing by CRISPR/Cas elements (Hamada et al, 2018). Another alternative is the use of ribonucleoprotein complexes instead of the expression cassette of the elements. In this way, transgene-free mutant plants can be obtained from the beginning. This strategy has been tested successfully in some crops including rice (Woo et al, 2015; Svitashev et al, 2016; Liang et al, 2017, 2018). Regarding delivery systems,- and particle bombardment-mediated transformation are the most widely utilized. However, both systems have their pros and cons. Efficient-mediated transformation is limited to a narrow range of genotypes within a species. Particle bombardment-mediated gene transfer can be applied to a wider range of genotypes than-mediated gene transfer, but plant regeneration after bombardment can be limiting (Altpeter et al, 2016). Therefore, the improvement of these techniques or the development of new technical advances can help in the development of genome-edited lines.

    The main concern about genome editing using CRISPR/Cas technology is the occurrence of undesired off-target events. This has been an important barrier to the application of genome editing systems in plant breeding. The use of new nucleases with higher specificity than Cas9, such as Cpf1, has improved this aspect of genome editing. Yin et al (2017) comparedthe two nucleases for genome editing in thevariety IR64 using CRISPR/Cas9 and CRISPR/Cpf1 and indicated that both systems are equally efficient for editing the() gene. T0plants generated by-mediated transformation of immature embryos show mutation frequencies (heterozygous) of 4% and 10% for Cas9 and Cpf1, respectively. However, T1plants show a higher number of homozygous mutant plants with Cas9 (25%–28%) than with Cpf1 (16%–20%). Off-target effects are analyzed in T2homozygous transgene-free plants using the Cas9 system and no off-target effects are observed (Yin et al, 2017).

    Recently, whole genome sequences of 34 edited rice (Nipponbare) plants using Cas9 (12 sgRNAs) and 15 edited plants using Cpf1 (3 crRNAs) were analyzed for off-targets. Both T0and T1plants are analyzed and only one Cas9 sgRNA (about 8%) results in off-target mutations in T0lines at predicted sites. None of the T1plants shows off-target mutations with Cas9 or Cpf1. Moreover, several mutations due to the tissue culture process are observed (approximately 102–148 single nucleotide variations and 32–83 inserts/deletions per plant). These results indicate that both enzymes are very specific in generating target mutations (Tang et al, 2018).

    Advantages of CRISPR/Cas9 system

    One of the most appealing characteristics of CRISPR/ Cas system is its multiplex editing capability, which leads to the mutation of several genes in the same transformation event, thereby reducing the time and cost of generating new varieties. This trait has been applied for the simultaneous mutation of an entire gene family in rice (Nipponbare and Kitaake). Eight sgRNAs targeting four members of the mitogen-activated protein kinase gene family have been designed. Each pair targets two genomic sites within a gene locus, 350–750 bp apart, so that excision could be observed in edited cells. Numerous sgRNAs are produced from a single polycistronic gene using the endogenous tRNA-processing system, which precisely cleaves both ends of the tRNA precursor. After protoplast transformation, analysis of the target genes reveals efficiently excised chromosomal fragments at all the four targeted loci with 4%–20% frequency. When the same construct is used to transform embryogenic calli by, transformed plants show mutations at all the eight sites. However, fragment deletion is detected only at two loci (Xie et al, 2015). With this system, simple, double and quadruple mutants are obtained. Mutation efficiency is 86%–100% for simple mutants, 67%–100% for double mutants and 86% for quadruple mutants. Phenotypic analysis of the mutants reveals severely dwarfed and sterilemutants and defective embryo development in homozygousseeds from heterozygous parents. In contrast, heterozygousmutant plants completely fail to produce homozygousseeds (Minkenberg et al, 2017).

    Another important advantage of CRISPR/Cas is that transgene-free mutant plants can be obtained in the first generation. This can be achieved by segregation of the transgene. In rice, there are several examples of T1transgene-free mutant plants (Li J et al, 2016; Wang et al, 2016). Another way to obtain transgene-free mutant plants is by transformation with ribonucleoprotein complexes (Cas9 and sgRNAs) for DNA-free genome editing. In this way, rice (Dongjin) protoplasts are transformed with Cas9 coupled to a sgRNA targeting the cytochromegene, and a mutation frequency of about 20% is achieved (Woo et al, 2015).

    Recently, a new rice transformation system was developed in order to obtain T1transgene-free plants. This system takes advantage of a pair of suicide transgenes that effectively kills all CRISPR/Cas-containing pollen and embryos produced by T0plants. One of these genes is the bacterialgene, which encodes a toxic protein that kills plant cells. This gene is under the control of thepromoter, which is expressed during early embryo development. This system contains the gene, encoding a rice male gametophyte specific lethal protein under thepromoter. In this way, T0plants containing these expression cassettes along with the CRISPR/Cas9 system will produce toxic proteins that kill male gametophytes and embryos, and only the transgene-free seeds will be produced. This construct is named Transgene Killer CRISPR (TKC) (He et al, 2018). In order to test this system, a sgRNA is designed to target the genebecausemutants show a large tiller angle and allow an easy way to identifymutant individuals. Aftertransformation of rice calli, 65 T0plants are obtained and 29 showphenotype. Seeds from 10 T0plants are evaluated for transgenes and none of them show the presence of TKC or CRISPR/Cas systems. The mutationsegregation pattern does not correspond to the expected Mendelian ratio. However, every single T1plant contains mutations at the target site and these plants were either homozygous or bi-allelic, demonstrating the 100% efficiency in eliminating both the LAZY1 function and the transgenes (He et al, 2018).

    Applications of CRISPR/Cas in rice

    A highly interesting application of this technology is the generation of a genome-wide mutant library that can be used for identifying gene functions and for genetic improvement. Meng et al (2017) searched the rice genome for every potential sgRNA, taking into consideration that the seed region (12 bp adjacent to PAM) should match only once in the genome in order to avoid off-targets. They found 1535852 target sites located in the exon regions of 52916 rice genes. In order to enrich loss-of-function mutations, sgRNAs located in exons near the beginning of ORFs are selected and two sgRNAs in each candidate gene are chosen. In this way, 12802 genes and 25604 sgRNAs are generated. The plasmid library is checked for quality, gene coverage and sgRNA accuracy by sequencing. Then, it is transformed intofor rice (Zhonghua11) transformation and more than 14000 independent T0lines are obtained. A total of 182 T0plants are randomly chosen and sequenced to evaluate the quality and coverage of the mutants. Of these, 139 plants harbouring single correct sgRNAs are found, of which 136 are different, indicating a good representation of every sgRNA. Finally, predicted off-targets were checked for some sgRNAs and no off-target mutations are found (Meng et al, 2017). In a similar way, Lu Y M et al (2017) designed a sgRNA library with 88541 members, targeting 34234 genes. There is an average of 2.59 sgRNA per gene in the rice cultivar MSU7. The library construction is checked by sequencing. Thus, more than 90% of the plasmid is confirmed and the coverage of 99% is observed. It is transformed intofor rice (MSU7) transformation. A total of 84384 transgenic plants are obtained with a mutation frequency of 83.9%. A high-throughput genotyping is performed as described by Bell et al (2014). In this way, 9216 plants are analyzed and 7004 are successfully identified (PCR positive), from which 86.5% contains a sgRNA. After all, 5541 plants are identified and 2326 loci are covered (Lu Y M et al, 2017). These results showed that the CRISPR/Cas9 system is a powerful tool for constructing mutant libraries in rice and can significantly accelerate the breeding process for the first time. Moreover, it will help in the functional characterization of unknown genes, not only in rice, but in other crops in general.

    Another interesting use of CRISPR/Cas is the precise elimination of marker genes in transgenic plants. Because of biosafety concerns, this is especially important when the selection marker gene encodes antibiotic and/or herbicide resistance. Moreover,the presence of selection markers makes the introduction of new genes into transgenic lines more difficult. Therefore, removing and recycling selection markers are highly desirable. For this purpose, two sgRNAs are designed targeting both ends of thegene in the B1transgenic rice (Nipponbare) line (Nandy and Srivastava, 2012). With this strategy, 5 out of 113 transgenic events show the expected excision of the full gene. Two of these are biallelic excisions. Sequence analysis of these events indicated that the excision is extremely precise since no mutations are observed in the two cutting sites (Srivastava et al, 2017). This new application of CRISPR/Cas could be used for selectable marker gene excision in rice or other plant species to improve breeding techniques.

    Genome editing to improve resistance to biotic and abiotic stresses

    Global climate change has increased the areas where adverse biotic and abiotic conditions seriously threat rice production and cause huge loss around the world. The International Food Policy Research Institute (IFPRI) forecasts that by 2050 rice yield loss could be about 13% worldwide, with the most affected areas being middle east-north Africa, Latin America and the Caribbean. The loss will cause rice price increased by32%to 37%. These adverse conditions include stresses caused by drought, salinity, extreme temperature, flooding, low nutrient availability, virus, bacterium, fungus, nematode, insect and herbivore. Drought and salinity are the main abiotic causes of rice yield loss worldwide. Blast disease, caused by, and bacterial blight, caused by, are the most destructive rice diseases.

    Resistance to

    utilizes a group of type III TAL (transcription activator-like) effectors to induce host gene expression and condition host susceptibility. Some of the targets of these effectors are SWEET genes, which have been shown to be sugar transporters. Zhou et al (2015) tested the role of thegene in rice resistance to. A sgRNA is designed to target the first exon of the gene. After-mediated transformation of rice calli (Kitaake), two mutant lines with deletions of 4 and 11 bases in the coding region ofare selected. These lines show the same phenotypic characteristics as control plants in normal growth conditions.However, they show increased resistance toandlesion length is reduced by about 90% compared to wild type plants in disease occurrence (Zhou et al, 2015).

    Resistance to

    Rice blast is the most devastating disease in all rice-growing countries, with typical loss of 10%–30%, although regional epidemics can be devastating (Dean et al, 2012). Therefore, augmenting resistance tois one of the most effective approaches for managing this disease. There is evidence that plant ethylene responsive factors (ERF) are involved in the modulation of multiple stress tolerance (Müller and Munné-Bosch, 2015). Moreover, knockdown expression of ricefrom accession Zhonghua 17 by RNAi enhances resistance to, indicating that this gene may act as the negative regulator of resistance (Liu et al, 2012). CRISPR/Cas9 is used to mutate thegene. Throughtransformation of rice calli (Kuiku131) using one sgRNA targeting the first exon of the gene, 21 mutant lines are obtained with different types of mutations, including deletions (64.3%), insertions (23.8%) and both (11.9%). Due to segregation, mutant homozygous transgene-free rice plants are obtained in the T1generation. The analysis of different agronomic traits (plant height, flag leaf length and width, number of productive panicles, panicle length, number of grains per panicle, seed-setting rate and 1000-grain weight) of these mutants shows no differences with the wild type. However, mutant lines show enhanced resistance tocompared to the wild type. Lesion length in the mutant lines is about 66% smaller than that in the wild type (Wang et al, 2016). The possibility of using several sgRNAs targeting different regions of the same gene is also analyzed. Thus, using three sgRNAs, 90% of the plants carry mutations in the three different sites at the same time, which suggests that this system allows highly efficient multilocus mutation in rice plants.

    In a different line of evidence, Ma et al (2018) explored the role of thegene in the defense response to.is an important subunit of the exocyst complex in rice. In this case, the authors designed two sgRNAs targeting the third and tenth exons of this gene. Five rice (Kitaake) mutant lines are obtained showing deletions, insertions or both in the targeted sites. These lines show altered growth and agronomic traits compared to the wild type, with smaller seedlings, shorter main roots, and decreased or impaired agricultural traits (plant height, panicle length, tiller number, 1000-grain weight and spikelet fertility). One line is chosen for further analysis; however, no transgene-free plants are obtained. This mutant line shows enhanced activation of the defence response, as evidenced by up-regulation of pathogenesis-related proteins and salicylic acid synthesis-related genes and increased levels of salicylic acid. Due to this augmented defence state,mutants are more tolerant toinfection (Ma et al, 2018).

    There are some reports about the use of CRISPR/Cas to study the role of genes involved in the rice response to different abiotic stresses.

    Salinity tolerance

    Duan et al (2016) studied the promoter of thegene, a transcription factor involved in the response to saline stress. A specific region of this promoter,, which is essential for the salt induction of this gene, is identified. In order to elucidate the role of this element, a sgRNA is designed to target theregion of the promoter. Two out of twelve rice (Nipponbare) lines are selected because they have lost theelement. The inability of the mutant lines to overexpress thegene under conditions of high salinity confirms the importance of this region.

    Drought tolerance

    Lou et al (2017) elucidated the role of(osmotic stress/ABA-activated protein kinase 2) using loss-of-function mutants by CRISPR/Cas. A sgRNA is designed to target the third exon of the gene. After-mediated transformation, 20 transgenic T0lines are obtained. Two homozygous T1mutant plants are obtained from these lines and characterised by sequencing the target region. Themutants exhibit an ABA-insensitive phenotype and are more sensitive to drought stress than the wildtype, indicating thatis important for the response to drought conditions in rice.

    Cold tolerance

    The role of the rice annexinduring cold stress was studied by Shen et al (2017). In this case, the second exon of the gene is targeted by the sgRNA. After-mediated transformation, the target regions are sequenced in 19 transgenic T0rice (Taipei309) lines to confirm mutations. In this way, six mutants are detected with four types of NHEJ mutations, 1-bp insertion, 1-bp deletion, 3-bp deletion and 4-bp deletion. Of these, 3 lines are homozygous mutants and 1 line is a biallelic mutant. Three T1lines from the T0homozygous and biallelic mutants are evaluated for cold tolerance, and all the mutant lines tested are more susceptible to cold stress. Interestingly, theauthors also analyzed off-target effects in other members of the annexin family and no off-target mutations are detected, confirming the high specificity of the CRISPR/Cas system in rice.

    Genome editing to improve herbicide tolerance

    Another important desirable trait that researchers have attempted to introduce in rice lines is resistance to herbicides. Several attempts have been made using transgenic lines (Inui et al, 2001; Zhao et al, 2011; Dong et al, 2017; Fartyal et al, 2018), however, due to the legal and commercial limitations regarding genetically modified organisms, they are moderately successful. The use of genome editing techniques allows the generation of transgene-free edited plants. There are several reports of adaptations of CRISPR/Cas technology to obtain new non-transgenic lines resistant to herbicides. In many cases, this agricultural trait is conferred by point mutations and is best achieved by precise genome editing.

    In this sense, a P178S mutation in() lead to herbicide resistance in soybean (Li Z S et al, 2015). ALS1is a key enzyme for the biosynthesis of branched chain amino acids and is a major target for important herbicides including chlorsulfuron and bispyribac sodium. Sun et al (2016) used the CRISPR/Cas technique with HR repair to replace the wild typegene with one carrying two point mutations (W548 L and S627 I) that can lead to herbicide resistance. Two sgRNAs are designed to target the region of interest (about 1625–1888 bp of the gene). A single plasmid carrying the Cas9 expressioncassette, both sgRNAs and the template for homologous recombination (with the point mutations of interest),is constructed. The template also contains the target site for both sgRNAs, and therefore Cas9 could cut both the genomic DNA and the plasmid with the HR template. The template is flanked by a left arm (100 bp) and a right arm (46 bp) homologous to the target region to improve recombination. In order to provide enough donor fragments for HR, the free DNA donor fragment and the plasmid are co-introduced into rice (Nipponbare) calli through particle bombardment. With this strategy, HR occurs in all the plants analyzed, although in different HR events. Every plant is sequenced to confirm the correct edition in the target region (Sun et al, 2016). It is important to point out that a similar strategy is used previously in maize (Svitashev et al, 2015), but when it is tested in rice, this strategy works much less efficiently, suggesting that a successful strategy in one plant species needs to be modified for other plant species.

    Gene replacement can be used to introduce point mutations. The enzyme 5-enolpyruvylshikimate-3-phosphate synthase () has a conserved motif in all plants that is the binding site for its substrate pyruvate and the herbicide glyphosate. A double amino acid substitution (T102I + P106S) in this domain leads to resistance to herbicide (Yu et al, 2015). In rice (Nipponbare), the second exon ofis replaced using CRISPR/Cas9 by designing two sgRNAs to target introns 1 and 2. The second exon with the required mutations flanked by the recognition sites for both sgRNAs is used as donor DNA. With this strategy, exon replacement is obtained in 2% of the analyzed plants. Since sgRNAs are directed to introns, the splicing pattern is also analyzed to make sure that it is not modified by mutations in introns. No alterations in splicing are found (Li J et al, 2016).

    Another strategy to induce point mutations through CRISPR/Cas is the fusion of Cas9 to a base editor, such as Petromyzon marinus cytidine deaminase (PmCDA1). This is called target-activation-induced cytidine deaminase (Target-AID). Basically, two mutated forms of Cas9 (dCas9 with no nuclease activity or nCas9 that has nickase activity) are fused to the base editor, so that base edition can be directed to a specific site through sgRNAs. The point mutation C287T of thegene endows rice plants with resistance to the herbicide imazamox. Shimatani et al (2017) designed a sgRNA targeting the C287 base, transformed rice (Nipponbare) calli using this system with both dCas9 and nCas9, and selected transformed calli on medium containing imazamox. As a result, 3 and 14 resistant lines from dCas9 and nCas9 transformants are obtained, respectively, and7out of the 14 lines obtained with nCas9 contain the expected C287T mutation and no off-target mutations are reported. Moreover, they demonstrated that using this system, point mutations could be transmitted from calli to regenerates and their progenies to obtain selectable marker-free herbicide tolerant plants through segregation (Shimatani et al, 2018).

    Another base editor used to create herbicide-tolerant rice plants using CRISPR/Cas is the adenine base editor (ABE), which consists of tRNA adenosine deaminase fused with nCas9. It generates A-T to C-G conversions when directed by sgRNAs. In order to create a suitable ABE system for rice, seven ABE fusion proteins are created that vary in the position of adenosine deaminase and the number and location of the nuclear localization sequences (Li J et al, 2018). These authors chose the most efficient configuration and designed a sgRNA to achieve the point mutation C2186R in the gene coding for acetyl-coenzyme A carboxylase (ACC) in rice (Zhonghua11) with resistance to herbicides across the aryloxyphenoxypropionate, cyclohexanedione and phenylpyrazoline chemical groups. Analysis of 160 transformed lines shows a mutation efficiency of 20.6% in the target region. Importantly, no off-target mutations are found among the T1mutant lines analysed. As expected, rice plants carrying the C2186R mutation are tolerant to the herbicide haloxyfop-R-methyl (Li J et al, 2018).

    Yield improvement by genome editing

    Grain yield is a complex trait governed by many different factors including agronomical practices and many genes known as quantitative trait loci (QTLs). The most popular method for plant breeding in essence is to generate various QTL combinations and select the elite ones (Shen et al, 2018). Some of these QTLs have been selected as target genes in studies using CRISPR/Cas to obtain new rice lines with improved yield. Li M R et al (2016) mutated the rice (Zhonghua11) genes,,and, which function as regulators of grain number, panicle architecture, grain size and plant architecture, respectively. The first exon ofandand the third exon ofandare chosen as targets for the sgRNAs because mutations in those regions have been proven to originate the desired phenotype (higher yield). And the mutation frequencies in T0transgenic lines are 42.5%, 67.5%, 57.5% and 27.5% for,,and, respectively. Themutation increases plant height, panicle size and number of flowers per panicle (about 90% more than the control). Themutants show decreased plant height and short panicles (about 20% less than the control), but an increased number of flowers per panicle (about 50% more than the control). Mutant lines ofshow larger grain size and longer awns on the husks. Finally, mutations ofresult in three phenotypes depending on the nature of the mutation. The sgRNA is designed to mutate the miR156 target site. If the mutation did not originate a frame shift or modify the miR156 target site, the phenotype was the same as the wild type. If the mutation changed the miR156 target site but did not produce a frame shift of the protein, mutant plants had 2–4 tillers, which is fewer than the 7–10 tillers typical of the wildtype, but plant height, flower number and panicle length were all increased compared with those of the wild type. Finally, if the mutation provoked a frame shift in the protein, the mutant plants had a dwarf phenotype with an increased number of tillers (Li M R et al, 2016). Three out of four sgRNAs produce off-target events in the predicted genomic regions. Off-target mutation frequenciesare 67% for the sgRNA of, 47.5% forand 2.5% for. This reinforces the idea that sgRNAs should be carefully designed in order to avoid unwanted genomic modifications. In this way, several bioinformatics tools have been created to design highly efficient and specific sgRNA in order to increase mutation efficiency and avoid unwanted off-target effects (Bae et al, 2014; Doench et al, 2016; Chari et al, 2017; Zhao et al, 2017).

    andmutations also result in higher grain number or increased grain size (Shen et al, 2018). However, rice yield in those mutants vary depending on the genetic background. Fivevarieties,Nanjing 9108, Wuyunjing 27, Yangjing 4227, Zhejing 22 and Zhejing 88 are chosen for mutagenesis. A single plasmid containing two sgRNAs targetingandis used for-mediated transformation of rice calli. Simpleand doublemutants are obtained for each variety and 10 new genotypes were obtained. Seven out of the ten novel genotypes show decreased grain yields per plant, ranging from 1% to 30% less than wild type. Only three genotypes show grain yields higher than the wild type,with increases of 3%–7%. The explanation is that despite the higher grain number and size of the mutant, some backgrounds show fewer effective tillers and lower grain number per tiller (Shen et al, 2018). Similarly, Zhou et al (2018) observed that mutations in the yield-related genes,andresult in increased yields per panicle in the varieties Jijing809 and Liaojing237, but not in Chuannongxiang jing. Multiplex sgRNA is used to generate double and triple mutants of these genes, and an additive effect of these genes is observed in Jijing 809 and Liaojing 237 (Zhou et al, 2018). Importantly, these authors analyzed the two most likely off-targets of each sgRNA and no mutations at those sites are found. It is important to highlight that analysing mutations in the same genes in different cultivars gives contrasting results. This shows the need to test those mutations in different agronomic cultivars around the globe.

    The CRISPR/Cas9-mediated multiplex genome editing system is one of the most attractive characteristics of this technology, especially for modifying QTLs to improve rice yield. In this way, Xu et al (2016) edited three QTLs related to grain weight,,and. The functions of,andare well characterized, and null mutations in any of them lead to significant increases in grain weight. Moreover,andmutants have larger grains than the wild type. The double mutant exhibits increases of about 12%, 8% and 13% in grain length, grain width and 1000-grain weight, respectively. The triple mutants display increases of about 25%, 20% and 28% in grain length, grain width and 1000-grain weight, respectively. The results show that this system is well suited for the rapid generation and pyramiding of beneficial alleles in rice.

    The possibility of deleting a precise region of a gene with CRISPR/Cas has also been used to increase yield. A 625 bp-deletion ingene confers dense and erect panicles with a higher grain number and lower plant height than the wild type (Huang et al, 2009). Taking this into consideration, Wang et al (2017) designed four sgRNAs (S1–S4) to achieve the same deletion inrice (inbred line IR58025B). Thefour sgRNAs are distributed in four different constructs, S1/S4, S1/S3, S2/S4 and S1/S2/S3/S4. Deletion frequencies of upto 21% for a 430 bp target and 9% for a 10 kb target, among 96 T0events per construct, are achieved. Themutants produced using CRISPR/Cas technology show the desired phenotype (dense and erect panicles and reduced plant height).

    Rice heading date, one of the most important agronomic traits determining rice distribution and production, is controlled by both genetic and environmental factors (Matsubara et al, 2014).,andare major genes for negative control of the heading date of rice varieties (Li X F et al, 2015), and mutation of those genes could lead to early-maturing varieties. Different sgRNAs (two forandand one for) are designed by Li et al (2017). Seven rice cultivars, Longdao16 (),Longdao18 (), Daohuaxiang2 (),Songjing19 (), Dongnong430 (), Dongnong429 () and Longqingdao2 (), were transformed by. A total of 18 independent transgenic lines are generated from the five varieties with genotype, among which 14 seedlings show mutations (single-base insertions or deletions) within all three genes, indicating a high efficiency (77.8%) of simultaneous editing of multiple genes. T1transgenic lines from the seven varieties are analyzed and nine homozygous T2lines are chosen for evaluation. The heading date in all of the lines is 5–30 d earlier than that in the wild type (Li et al, 2017).

    Hybrid rice plays an advantage of 10%–20% over conventional rice in yield, and therefore it has a key role in rice production worldwide. Hybrid lines can be developed using the three- and two-line hybrid breeding systems (Cheng et al, 2007). The two-line breeding system uses either photoperiod-sensitive genic male-sterile (PGMS) or thermo-sensitive genic male-sterile (TGMS) lines as sterility or maintainer lines under restrictive or permissive conditions, respectively. Thus, development of new P/TGMS lines would improve and accelerate the breeding process. Zhou et al (2016) modified the geneby CRISPR/Cas to obtain a new TGMS line. Ten sgRNAsare designed targeting this gene and used to transform rice calli (Zhonghua 11) through. T0plants are obtained with mutations in nine out of ten sites, with a homozygous mutation rate of about 32%. Among T0plants, 30%–85% were pollen-sterile individuals under restrictive conditions. To obtain TGMS plants, T0sterile plants are treated at a low temperature and transgene-free TGMS lines are obtained in T1plants. Furthermore, the most efficient sgRNA is chosen to transform 11 cultivars to obtain TGMS plants, including tencultivars (Zhenshan97B, ZhongzheB, TiangfengB, YixiangB, ReB, HuahuiB, HuanongB, Yuejingsimiao, Yuenongsimiao and Wushansimiao) and onecultivar GAZ. In this case, the mutation frequency ranges from 72% to 100% with a homozygous mutation rate of 11%–54%. Mutant transgene-free plants are obtained in T1generation and the mutations are stable in T2and T3generations. To obtain TGMS lines, these lines are treated with differenttemperature gradients and found that fertility is gradually reduced with increasing temperature. The use of CRISPR/Cas to produce new TGMS lines can facilitate the exploitation of heterosis in breeding process. Moreover, this technology can be applied potentially for breeding other hybrid crops not only rice.

    Genome editing application in rice quality

    Another goal of rice improvement is to create new varieties with increased nutritional and commercial value. Of particular interest is the modification of starch and amylose content. A cereal grain higher in amylose is a good source of resistant starch, which is a special kind of starch that is not digested in the stomach or small intestine, but passes directly to the large intestine. Resistant starch helps to improve human health and reduce the risk of non-infectious diseases (Regina et al, 2006). Amylose and resistant starch contents can be increased by modifying the starch branching enzymes (SBEs).is mainly expressed in rice grains, which is why it is chosen as a target for CRISPR/Cas mutation. After-mediatedtransformation of rice (Kitaake) calli, 22 of 30 transgenic plants show mutations in the expected region, with frequencies of 6.7%, 36.6% and 26.7% for heterozygous, bi-allelic and homozygous lines, respectively. No off-target effects are detected in the predicted sites. After the generation of transgenic-free mutant lines through segregation,mutants show significantly increased amylose and resistant starch contents (upto 25.0% and 9.8%, respectively) (Sun et al, 2017).

    Amylose content is also the most important quality indicator in rice, especially for cooking and eating quality. The quality of some rice varieties, especiallyhybrids, is considered poor in some markets owing to their high amylose content. In this case, reducing the amylose content to obtain glutinous rice is a highly desirable characteristic. A single dominant gene calledcontrols the amylose content in rice. Therefore, in order to obtain a loss-of-function mutant through CRISPR/Cas, a sgRNA targeting the first exon of thegene is designed by Zhang et al (2018). Following-mediated calli transformation of two cultivars (Xiushui134 and Wuyunjing 7), a high mutagenesis efficiency(82%– 87%) is observed in T0transformants. About 4%–15% of the mutants were homozygous. Interestingly, no off-target mutations are observed in any of the mutants analyzed. Analysis of agronomic traits in T1generation plants shows no differences between wild type andmutants regarding plant height, grain number per panicle, panicle number per plant, yield per plot, grain width, grain length and 1000-grain weight. The amylose content is significantly lower (85%) inmutants than in the wild type plants (Zhang et al, 2018).

    Rice bran oil (RBO) is a commercial derivative of rice that is produced and consumed mainly in Asian countries. RBO contains components with valuable health-promoting properties, one of which is oleic acid, which helps prevent lifestyle diseases. Modifications in fatty acid synthesis pathways could lead to higher levels of oleic acid and more valuable RBO. The enzyme fatty acid desaturase 2 (FAD2) is responsible for the conversion of oleic acid to linoleic acid. The rice genome has fourgenes, of whichis the most highly expressed in rice seeds. Thus, the knockout ofby CRISPR/Cas could lead to a variety with high oleic acid content. Through-mediated transformation of rice (Nipponbare) calli, six independent mutant lines are regenerated. After evaluating the progeny to confirm the mutation, fatty acid profiles of the T2lines are analyzed. Inhomozygous mutants, the oleic acid content increases to more than twice of that in wild type, and linoleic acid is not detected (Abe et al, 2018).

    Cadmium (Cd) is highly toxic heavy metal for living organisms. Rice grains with excesive Cd content are a serious threat for people who consume rice as a staple food. Persistent intake of Cd could lead to different health problems (Bertin and Averbeck, 2006; Clemens et al, 2013). The development of elite rice cultivars with low Cd accumulation is very important, especially forvarieties, which can accumulate more Cd thanones. Tang et al (2017) created a newrice line with low Cd accumulation by mutating the metal transporter geneusing CRISPR/Cas9 system. This transporter mediates the root uptake of Cd, and therefore mutating this gene causes a great reduction in Cd content. Two sgRNAs are designed to target the ninth exon of the gene.-transformation of embriogenic calli leads to a mutation efficiency of 82.4% for Huazhan and 80% for Longke6385. No off-targets effects are observed. Five (three from Huazhan and two from Longke6385) transgene-free homozygous knockout lines are obtained in T1generation by segregation.lines show less Cd accumulation (0.05 mg/kg) when grown on Cd-contaminated paddy fields compared to the wild type (0.33–2.90 mg/kg). Moreover, mutant plants show no differences in agronomic traits such as grain yield, straw weight or grain quality compared with wild type plants.

    Genome editing application in nutrient use

    As mentioned in previous sections, the CRISPR/Cas system can accelerate breeding processes that would otherwise take years using traditional crossing/backcrossing techniques. An example of this is the introduction of an elite allele into commercial cultivars. For instance, a single nucleotide polymorphism in the nitrate transportergene is responsible for improved nitrogen use efficiency inrice. Therefore, the replacement of this gene could lead to commercial cultivars with improved nitrogen use efficiency. Li J Y et al (2018) employed the strategy described earlier (Sun et al, 2016) with some modifications to replace theallele with theallele. First, the two sgRNAs target different strands of DNA to minimize the non-homologous sequences when released from the vector. Second, the PAM site is mutated in the template sequence to prevent further cuts after HR is successful. Finally, a longer right arm is used to enhance HR frequency. With this improved strategy, the precise replacement of the allele, either homozygous or heterozygous, is obtained in 6.72% of the transformed plants. Partial HR and reverse complementation are also observed, but no off-target mutations are found in the plants analyzed (Li J Y et al, 2018),which demonstrated that CRISPR/Cas technology can be used for efficient allelic replacement by HR in rice in one generation, thereby accelerating the crop improvement process.

    PERSPECTS

    Since the CRISPR/Cas genome editing system was first implemented in 2013, its use and efficiency for site-directed mutation have been tested in plants. This system has become the most common tool for functional and/or agronomic studies in different plant species, especially in rice, which is not only very important at the agronomic level, but is also a model species for physiological and molecular studies in monocots. The growing interest and constant search for new applications has led to the generation of new knowledge and improvements that facilitate the application of this technology in rice and other crops of agronomic interest. The research demonstrates the wide variety of applications of this technology: the generation of small InDels that cause mutation by frame shift changes, large gene deletions, the precise edition of certain bases to generate point mutations using the base-editing systems coupled to CRISPR/Cas and gene substitution by homologous recombination at specific sites in the genome. Even the generation of a collection of mutants in more than 12000 genes is possible thanks to this technology. This technology and the growing advances in massive sequencing techniques have made possible the creation of new rice lines with specific site directed mutations without off-target mutations. More importantly,because of the strict regulation of transgenic organisms, this technique allows the generation of transgene-free mutant plants. Altogether, this technique facilitates not only the analysis of the function of different genes, but also the process of genetic improvement of this important species. It is because of all these advantages that this technology is preferred by researchers dedicated to genetic improvement of rice worldwide.

    ACKNOWLEDGEMENTS

    The authors thank Dr. Isabel Mora Ramírez (Leibniz Institute of Plant Genetics and Crop Plant Research, Gatersleben, Germany) and Dr. Marco Herde (Leibniz University of Hanover, Germany) for their helpful comments on this review.

    Abe K, Araki E, Suzuki Y, Toki S, Saika H. 2018. Production of high oleic/low linoleic rice by genome editing., 131: 58–62.

    Altpeter F, Springer N M, Bartley L E, Blechl A E, Brutnell T P, Citovsky V, Conrad L J, Gelvin S B, Jackson D P, Kausch A P, Lemaux P G, Medford J I, Orozco-Cárdenas M L, Tricoli D M, Van Eck J, Voytas D F, Walbot V, Wang K, Zhang Z J, Stewart C N Jr. 2016. Advancing crop transformation in the era of genome editing., 28(7): 1510–1520.

    Anders C, Bargsten K, Jinek M. 2016. Structural plasticity of PAM recognition by engineered variants of the RNA-guided endonuclease Cas9., 61(6): 895–902.

    Bae S, Park J, Kim JS. 2014. Cas-OFFinder: A fast and versatile algorithm that searches for potential off-target sites of Cas9 RNA-guided endonucleases., 30(10): 1473–1475.

    Belhaj K, Chaparro-Garcia A, Kamoun S, Patron N J, Nekrasov V. 2015. Editing plant genomes with CRISPR/Cas9., 32: 76–84.

    Bell C C, Magor G W, Gillinder K R, Perkins A C. 2014. A high-throughput screening strategy for detecting CRISPR-Cas9 induced mutations using next-generation sequencing., 15(1): 1002.

    Bertin G, Averbeck D. 2006. Cadmium: Cellular effects, modifications ofbiomolecules, modulation ofDNA repair andgenotoxic consequences., 88(11): 1549–1559.

    Bhaya D, Davison M, Barrangou R. 2011. CRISPR-Cas systems in bacteria and archaea: Versatile small RNAs for adaptive defense and regulation., 45(1): 273–297.

    Breviario D, Genga A. 2013. Stress response in rice., 2: e104.

    Clemens S, Aarts M G M, Thomine S, Verbruggen N. 2013. Plant science: The key to preventing slow cadmium poisoning., 18(2): 92–99.

    Chari R, Yeo N C, Chavez A, Church G M. 2017. sgRNA scorer 2.0: A species-independent model to predict CRISPR/Cas9 activity., 6(5): 902–904.

    Chen L Z, Li W, Katin-Grazzini L, Ding J, Gu X B, Li Y J, Gu T T, Wang R, Lin X C, Deng Z N, McAvoy R J, Gmitter F G, Deng Z N, Zhao Y D, Li Y. 2018. A method for the production and expedient screening of CRISPR/Cas9-mediated non-transgenic mutant plants., 5: 13.

    Cheng SH, Zhuang JY, Fan YY, Du JH, Cao LY. 2007. Progress in research and development on hybrid rice: A super-domesticate in China., 100(5): 959–966.

    Dean R, Van Kan J A L, Pretorius Z A, Hammond-Kosack K E, Di Pietro A, Spanu P D, Rudd J J, Dickman M, Kahmann R, Ellis J, Foster G D. 2012. The top 10 fungal pathogens in molecular plant pathology., 13(4): 414–430.

    Delorge I, Janiak M, Carpentier S, van Dijck P. 2014. Fine tuning of trehalose biosynthesis and hydrolysis as novel tools for the generation of abiotic stress tolerant plants., 5: 147.

    Doench J G, Fusi N, Sullender M, Hegde M, Vaimberg E W, Donovan K F, Smith I, Tothova Z, Wilen C, Orchard R, Virgin H W, Listgarten J, Root D E. 2016. Optimized sgRNA design to maximize activity and minimize off-target effects of CRISPR-Cas9., 34: 184–191.

    Dong Y F, Jin X, Tang Q L, Zhang X, Yang J T, Liu X J, Cai J F, Zhang X B, Wang X J, Wang Z X. 2017. Development and event-specific detection of transgenic glyphosate-resistant rice expressing thegene., 8: 885.

    Duan YB, Li J, Qin RY, Xu RF, Li H, Yang YC, Ma H, Li L, Wei PC, Yang JB. 2016. Identification of a regulatory element responsible for salt induction of ricethroughandpromoter analysis., 90(1): 49–62.

    Endo M, Nishizawa-Yokoi A, Toki S. 2018. Rice genome editing.: Sasaki T, Ashikari M. Rice Genomics, Genetics and Breeding.Singapore: Springer:523–539.

    Fartyal D, Agarwal A, James D, Borphukan B, Ram B, Sheri V, Agrawal P K, Achary V M M, Reddy M K. 2018. Developing dual herbicide tolerant transgenic rice plants for sustainable weed management., 8(1): 11598.

    Feng X P, Peng C, Chen Y, Liu X D, Feng X J, He Y. 2017. Discrimination of CRISPR/Cas9-induced mutants of rice seeds using near-infrared hyperspectral imaging., 7: 15934.

    Hamada H, Liu Y L, Nagira Y, Miki R, Taoka N, Imai R. 2018. Biolistic-delivery-based transient CRISPR/Cas9 expression enables in planta genome editing in wheat., 8: 14422.

    He Y B, Zhu M, Wang L H, Wu J H, Wang Q Y, Wang R C, Zhao Y D. 2018. Programmed self-elimination of the CRISPR/Cas9 construct greatly accelerates the isolation of edited and transgene-free rice plants., 11(9): 1210–1213.

    Hermans P W, van Soolingen D, Bik E M, de Haas P E, Dale J W, van Embden J D. 1991. Insertion element IS987 fromBCG is located in a hot-spot integration region for insertion elements incomplex strains., 59(8): 2695–2705.

    Huang X Z, Qian Q, Liu Z B, Sun H Y, He S Y, Luo D, Xia G M, Chu C C, Li J Y, Fu X D. 2009. Natural variation at thelocus enhances grain yield in rice., 41: 494–497.

    Inui H, Shiota N, Ido Y, Inoue T, Hirose S, Kawahigashi H, Ohkawa Y, Ohkawa H. 2001. Herbicide metabolism and tolerance in the transgenic rice plants expressing human CYP2C9 and CYP2C19., 71(3): 156–169.

    Ishino Y, Shinagawa H, Makino K, Amemura M, Nakata A. 1987. Nucleotide sequence of thegene, responsible for alkaline phosphatase isozyme conversion in, and identification of the gene product., 169(12): 5429–5433.

    Jiang W Y, Bikard D, Cox D, Zhang F, Marraffini L A. 2013. RNA-guided editing of bacterial genomes using CRISPR-Cas systems., 31(3): 233–239.

    Jiang W Z, Zhou H B, Bi H H, Fromm M, Yang B, Weeks D P. 2013. Demonstration of CRISPR/Cas9/sgRNA-mediated targeted gene modification in, tobacco, sorghum and rice., 41(20): e188.

    Jinek M, Chylinski K, Fonfara I, Hauer M, Doudna J A, Charpentier E. 2012. A programmable dual-RNA: Guided DNA endonuclease in adaptive bacterial immunity., 337: 816–821.

    Jung J H, Seo Y W. 2017. Challenges in wide implementation of genome editing for crop improvement., 20(2): 129–135.

    Kaeppler S M, Kaeppler H F, Rhee Y. 2000. Epigenetic aspects of somaclonal variation in plants., 43: 179–188.

    Koonin E V, Makarova K S, Zhang F. 2017. Diversity, classification and evolution of CRISPR-Cas systems., 37: 67–78.

    Krappmann S. 2017. CRISPR-Cas9, the new kid on the block of fungal molecular biology., 55(1): 16–23.

    Li C, Zong Y, Wang Y P, Jin S, Zhang D B, Song Q N, Zhang R, Gao C X. 2018. Expanded base editing in rice and wheat using a Cas9-adenosine deaminase fusion., 19: 59.

    Li J, Meng X B, Zong Y, Chen K L, Zhang H W, Liu J X, Li J Y, Gao C X. 2016. Gene replacements and insertions in rice by intron targeting using CRISPR-Cas9., 2: 16139.

    Li J Y, Zhang X, Sun Y W, Zhang J H, Du W M, Guo X P, Li S Y, Zhao Y D, Xia L Q. 2018. Efficient allelic replacement in rice by gene editing: A case study of thegene., 60(7): 536–540.

    Li M R, Li X X, Zhou Z J, Wu P Z, Fang M C, Pan X P, Lin Q P, Luo W B, Wu G J, Li H Q. 2016. Reassessment of the four yield-related genes,,,andin rice using a CRISPR/Cas9 system., 7: 377.

    Li X F, Liu H Z, Wang M Q, Liu H L, Tian X J, Zhou W J, Lü T X, Wang Z Y, Chu C C, Fang J, Bu Q Y. 2015. Combinations ofandgenes determine rice adaptability to Heilongjiang Province, northern limit of China., 57(8): 698–707.

    Li X F, Zhou W J, Ren Y K, Tian X J, Lv T X, Wang Z Y, Fang J, Chu C C, Yang J, Bu Q Y. 2017. High-efficiency breeding of early-maturing rice cultivars via CRISPR/Cas9-mediated genome editing., 44(3): 175–178.

    Li Z S, Liu ZB, Xing A Q, Moon B P, Koellhoffer J P, Huang L X, Ward R T, Clifton E, Falco S C, Cigan A M. 2015. Cas9-guide RNA directed genome editing in soybean., 169(2): 960–970.

    Liang Z, Chen K L, Li T D, Zhang Y, Wang Y P, Zhao Q, Liu J X, Zhang H W, Liu C M, Ran Y D, Gao C X. 2017. Efficient DNA-free genome editing of bread wheat using CRISPR/Cas9 ribonucleoprotein complexes., 8: 14261.

    Liang Z, Chen K L, Zhang Y, Liu J X, Yin K Q, Qiu JL, Gao C X. 2018. Genome editing of bread wheat using biolistic delivery of CRISPR/Cas9transcripts or ribonucleoproteins., 13: 413–430.

    Liu D F, Chen X J, Liu J Q, Ye J C, Guo Z J. 2012. The rice ERF transcription factornegatively regulates resistance toand salt tolerance., 63(10): 3899–3911.

    Liu Q, Wang C, Jiao X Z, Zhang H W, Song L L, Li Y X, Gao C X, Wang K J. 2019. Hi-TOM: A platform for high-throughput tracking of mutations induced by CRISPR/Cas systems., 62(1): 1–7.

    Lou D J, Wang H P, Liang G, Yu D Q. 2017.confers abscisic acid sensitivity and tolerance to drought stress in rice., 8: 993.

    Lu HP, Liu SM, Xu SL, Chen WY, Zhou X, Tan YY, Huang JZ, Shu QY. 2017. CRISPR-S: An active interference element for a rapid and inexpensive selection of genome-edited, transgene-free rice plants., 15(11): 1371–1373.

    Lu Y M, Ye X, Guo R M, Huang J, Wang W, Tang J Y, Tan L T, Zhu J K, Chu C C, Qian Y W. 2017. Genome-wide targeted mutagenesis in rice using the CRISPR/Cas9 system., 10(9): 1242–1245.

    Ma J, Chen J, Wang M, Ren Y L, Wang S, Lei C L, Cheng Z J, Sodmergen. 2018. Disruption ofincreases the content of salicylic acid and induces plant defense responses in rice., 69(5): 1051–1064.

    Matsubara K, Hori K, Ogiso-Tanaka E, Yano M. 2014. Cloning of quantitative trait genes from rice reveals conservation and divergence of photoperiod flowering pathways inand rice., 5: 193.

    Meng X B, Yu H, Zhang Y, Zhuang F F, Song X G, Gao S S, Gao C X, Li J Y. 2017. Construction of a genome-wide mutant library inrice using CRISPR/Cas9., 10(9): 1238–1241.

    Miao J, Guo D S, Zhang J Z, Huang Q P, Qin G J, Zhang X, Wan J M, Gu H Y, Qu LJ. 2013. Targeted mutagenesis in rice using CRISPR-Cas system., 23(10): 1233–1236.

    Minkenberg B, Xie K B, Yang Y N. 2017. Discovery of rice essential genes by characterizing a CRISPR-edited mutation of closely related rice MAP kinase genes., 89(3): 636–648.

    Mishra R, Joshi R K, Zhao K J. 2018. Genome editing in rice: Recent advances, challenges, and future implications., 9: 1361.

    Mojica F J M, Juez G, Rodriguez-Valera F. 1993. Transcription at different salinities ofsequences adjacent to partially modifiedsites., 9(3): 613–621.

    Moscou M J, Bogdanove A J. 2009. A simple cipher governs DNA recognition by TAL effectors., 326: 1501.

    Müller M, Munné-Bosch S. 2015. Ethylene response factors: A key regulatory hub in hormone and stress signaling., 169(1): 32–41.

    Nakata A, Amemura M, Makino K. 1989. Unusual nucleotide arrangement with repeated sequences in theK-12 chromosome., 171(6): 3553–3556.

    Nandy S, Srivastava V. 2012. Marker-free site-specific gene integration in rice based on the use of two recombination systems., 10(8): 904–912.

    Pabo C O, Peisach E, Grant R A. 2001. Design and selection of novel Cys2His2zinc finger proteins., 70(1): 313–340.

    Peng C, Wang H, Xu X L, Wang X F, Chen X Y, Wei W, Lai Y M, Liu G Q, Godwin I D, Li J Q, Zhang L, Xu J F. 2018. High-throughput detection and screening of plants modified by gene editing using quantitative real-time polymerase chain reaction., 95(3): 557–567.

    Pyott D E, Sheehan E, Molnar A. 2016. Engineering of CRISPR/Cas9-mediated potyvirus resistance in transgene-freeplants., 17(8): 1276–1288.

    Regina A, Bird A, Topping D, Bowden S, Freeman J, Barsby T, Kosar-Hashemi B, Li Z Y, Rahman S, Morell M. 2006. High-amylose wheat generated by RNA interference improves indices of large-bowel health in rats., 103(10): 3546–3551.

    Ricroch A, Clairand P, Harwood W. 2017. Use of CRISPR systems in plant genome editing: Toward new opportunities in agriculture., 1(2): 169–182.

    Sarmast M K. 2016. Genetic transformation and somaclonal variation in conifers., 10(6): 309–325.

    Shan Q W, Wang Y P, Li J, Zhang Y, Chen K L, Liang Z, Zhang K, Liu J X, Xi J J, Qiu JL, Gao C X. 2013. Targeted genome modification of crop plants using a CRISPR-Cas system., 31(8): 686–688.

    Shen C X, Que Q Z, Xia Y M, Tang N, Li D, He R H, Cao M L. 2017. Knock out of the annexin genevia CRISPR/Cas9-mediated genome editing decreased cold tolerance in rice., 60(6): 539–547.

    Shen L, Wang C, Fu Y P, Wang J J, Liu Q, Zhang X M, Yan C J, Qian Q, Wang K J. 2018. QTL editing confers opposing yield performance in different rice varieties., 60(2): 89–93.

    Shimatani Z, Kashojiya S, Takayama M, Terada R, Arazoe T, Ishii H, Teramura H, Yamamoto T, Komatsu H, Miura K, Ezura H, Nishida K, Ariizumi T, Kondo A. 2017. Targeted base editing in rice and tomato using a CRISPR-Cas9 cytidine deaminase fusion., 35: 441–443.

    Shimatani Z, Fujikura U, Ishii H, Matsui Y, Suzuki M, Ueke Y, Taoka K, Terada R, Nishida K, Kondo A. 2018. Inheritance of co-edited genes by CRISPR-based targeted nucleotide substitutions in rice., 131: 78–83.

    Sprink T, Metje J, Hartung F. 2015. Plant genome editing by novel tools: TALEN and other sequence specific nucleases., 32: 47–53.

    Srivastava V, Underwood J L, Zhao S. 2017. Dual-targeting by CRISPR/Cas9 for precise excision of transgenes from rice genome., 129(1): 153–160.

    Sun Y W, Zhang X, Wu C Y, He Y B, Ma Y Z, Hou H, Guo X P, Du W M, Zhao Y D, Xia L Q. 2016. Engineering herbicide-resistant rice plants through CRISPR/Cas9-mediated homologous recombination of acetolactate synthase., 9(4): 628–631.

    Sun Y W, Jiao G A, Liu Z P, Zhang X, Li J Y, Guo X P, Du W M, Du J L, Francis F, Zhao Y D, Xia L Q. 2017. Generation of high-amylose rice through CRISPR/Cas9-mediated targeted mutagenesis of starch branching enzymes., 8: 298.

    Svitashev S, Schwartz C, Lenderts B, Young J K, Mark Cigan A. 2016. Genome editing in maize directed by CRISPR-Cas9 ribonucleoprotein complexes., 7: 13274.

    Tang L, Mao B G, Li Y K, Lv Q M, Zhang L P, Chen C Y, He H J, Wang W P, Zeng X F, Shao Y, Pan Y L, Hu Y Y, Peng Y, Fu X Q, Li H Q, Xia S T, Zhao B R. 2017. Knockout ofusing the CRISPR/Cas9 system produces low Cd-accumulatingrice without compromising yield., 7: 14438.

    Tang X, Liu G Q, Zhou J P, Ren Q R, You Q, Tian L, Xin X H, Zhong Z H, Liu B L, Zheng X L, Zhang D W, Malzahn A, Gong Z Y, Qi Y P, Zhang T, Zhang Y. 2018. A large-scale whole-genome sequencing analysis reveals highly specific genome editing by both Cas9 and Cpf1 (Cas12a) nucleases in rice., 19: 84.

    Wang F J, Wang C L, Liu P Q, Lei C L, Hao W, Gao Y, Liu YG, Zhao K J. 2016. Enhanced rice blast resistance by CRISPR/Cas9-targeted mutagenesis of the ERF transcription factor gene., 11(4): e0154027.

    Wang Y, Geng L Z, Yuan M L, Wei J, Jin C, Li M, Yu K, Zhang Y, Jin H B, Wang E, Chai Z J, Fu X D, Li X G. 2017. Deletion of a target gene inrice via CRISPR/Cas9., 36(8): 1333–1343.

    Woo J W, Kim J, Kwon S I, Corvalán C, Cho S W, Kim H, Kim SG, Kim ST, Choe S, Kim JS. 2015. DNA-free genome editing in plants with preassembled CRISPR-Cas9 ribonucleoproteins., 33: 1162–1164.

    Xie K B, Yang Y N. 2013. RNA-guided genome editing in plants using a CRISPR-Cas system., 6(6): 1975–1983.

    Xie K B, Minkenberg B, Yang Y N. 2015. Boosting CRISPR/Cas9 multiplex editing capability with the endogenous tRNA-processing system., 112(11): 3570–3575.

    Xu R F, Yang Y C, Qin R Y, Li H, Qiu C H, Li L, Wei P C, Yang J B. 2016. Rapid improvement of grain weight via highly efficient CRISPR/Cas9-mediated multiplex genome editing in rice., 43(8): 529–532.

    Yin X J, Biswal A K, Dionora J, Perdigon K M, Balahadia C P, Mazumdar S, Chater C, Lin HC, Coe R A, Kretzschmar T, Gray J E, Quick P W, Bandyopadhyay A. 2017. CRISPR-Cas9 and CRISPR-Cpf1 mediated targeting of a stomatal developmental genein rice., 36(5): 745–757.

    Yu Q, Jalaludin A, Han H P, Chen M, Sammons R D, Powles S B. 2015. Evolution of a double amino acid substitution in the 5-enolpyruvylshikimate-3-phosphate synthase inconferring high-level glyphosate resistance., 167(4): 1440–1447.

    Zhang J S, Zhang H, Botella J R, Zhu JK. 2018. Generation of new glutinous rice by CRISPR/Cas9-targeted mutagenesis of thegene in elite rice varieties., 60(5): 369–375.

    Zhang Y, Liang Z, Zong Y, Wang Y P, Liu J X, Chen K L, Qiu JL, Gao C X. 2016. Efficient and transgene-free genome editing in wheat through transient expression of CRISPR/Cas9 DNA or RNA., 7: 12617.

    Zhao C Z, Zheng X G, Qu W B, Li G L, Li X Y, Miao YL, Han X S, Liu X D, Li Z H, Ma Y L, Shao Q Z, Li H W, Sun F, Xie S S, Zhao S H. 2017. CRISPR-offinder: A CRISPR guide RNA design and off-target searching tool for user-defined protospacer adjacent motif., 13(12): 1470–1478.

    Zhao T, Lin C Y, Shen Z C. 2011. Development of transgenic glyphosate-resistant rice withgene encoding 5-enolpyruvylshikimate-3-phosphate synthase., 10(9): 1307–1312.

    Zhou H, He M, Li J, Chen L, Huang Z F, Zheng S Y, Zhu L Y, Ni E D, Jiang D G, Zhao B R, Zhuang C X. 2016. Development of commercial thermo-sensitive genic male sterile rice accelerates hybrid rice breeding using the CRISPR/Cas9-mediatedediting system., 6: 37395.

    Zhou J H, Peng Z, Long J Y, Sosso D, Liu B, Eom JS, Huang S, Liu S Z, Vera Cruz C, Frommer W B, White F F, Yang B. 2015. Gene targeting by the TAL effector PthXo2 reveals cryptic resistance gene for bacterial blight of rice., 82(4): 632–643.

    Zhou J P, Xin X H, He Y, Chen H Q, Li Q, Tang X, Zhong Z H, Deng K J, Zheng X L, Akher S A, Cai G Z, Qi Y P, Zhang Y. 2018. Multiplex QTL editing of grain-related genes improves yield in elite rice varieties., 38(4): 475–485.

    15 December 2018;

    7 March 2019

    Andrés Gatica-Arias(andres.gatica@ucr.ac.cr)

    Copyright ? 2019, China National Rice Research Institute. Hosting by Elsevier B V

    This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/)

    Peer review under responsibility of China National Rice Research Institute

    http://dx.doi.org/10.1016/j.rsci.2019.08.001

    (Managing Editor: Wang Caihong)

    黄色片一级片一级黄色片| 三级毛片av免费| 18禁裸乳无遮挡免费网站照片| 色播亚洲综合网| 岛国在线免费视频观看| 男女视频在线观看网站免费| 亚洲国产精品成人综合色| 亚洲av免费高清在线观看| 99精品久久久久人妻精品| 国产中年淑女户外野战色| 色av中文字幕| 欧美激情久久久久久爽电影| 搡老岳熟女国产| 黄色视频,在线免费观看| 国产精品美女特级片免费视频播放器| 欧美一级a爱片免费观看看| 成人国产综合亚洲| 国产蜜桃级精品一区二区三区| 亚洲av第一区精品v没综合| 日本撒尿小便嘘嘘汇集6| 国产成人av激情在线播放| 在线观看免费视频日本深夜| 草草在线视频免费看| 91在线精品国自产拍蜜月 | 亚洲人成网站在线播放欧美日韩| 在线观看一区二区三区| 国产精品 国内视频| 欧美一级毛片孕妇| 成年版毛片免费区| 国产欧美日韩一区二区三| 免费搜索国产男女视频| 丰满人妻熟妇乱又伦精品不卡| 少妇高潮的动态图| 久久99热这里只有精品18| 女人十人毛片免费观看3o分钟| 日本黄色片子视频| 欧美精品啪啪一区二区三区| 18禁黄网站禁片免费观看直播| 国内毛片毛片毛片毛片毛片| 狂野欧美激情性xxxx| 欧美乱码精品一区二区三区| 久久国产乱子伦精品免费另类| 露出奶头的视频| 熟女人妻精品中文字幕| 在线观看舔阴道视频| 男人舔奶头视频| 女生性感内裤真人,穿戴方法视频| 一级作爱视频免费观看| 脱女人内裤的视频| 人人妻,人人澡人人爽秒播| 色综合站精品国产| 性欧美人与动物交配| 搡老熟女国产l中国老女人| 好男人电影高清在线观看| 九色成人免费人妻av| 国产精品嫩草影院av在线观看 | 高清毛片免费观看视频网站| 脱女人内裤的视频| 男女午夜视频在线观看| 99国产精品一区二区蜜桃av| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品一卡2卡三卡4卡5卡| 亚洲片人在线观看| 久久九九热精品免费| 亚洲人成网站高清观看| 欧美一区二区精品小视频在线| 有码 亚洲区| 日韩有码中文字幕| 搡老熟女国产l中国老女人| 2021天堂中文幕一二区在线观| 国产一区二区三区在线臀色熟女| 特大巨黑吊av在线直播| 午夜两性在线视频| 亚洲av日韩精品久久久久久密| 午夜激情福利司机影院| 欧美绝顶高潮抽搐喷水| 国产精品av视频在线免费观看| 日本 欧美在线| 丰满人妻熟妇乱又伦精品不卡| 日本成人三级电影网站| 久久香蕉国产精品| 国产成人欧美在线观看| 99久久九九国产精品国产免费| 久久亚洲精品不卡| 国模一区二区三区四区视频| 欧美大码av| 亚洲国产高清在线一区二区三| 男女床上黄色一级片免费看| 非洲黑人性xxxx精品又粗又长| 91在线观看av| 日本 av在线| 久久精品91无色码中文字幕| 怎么达到女性高潮| 亚洲av成人不卡在线观看播放网| 91在线精品国自产拍蜜月 | АⅤ资源中文在线天堂| aaaaa片日本免费| 亚洲av成人不卡在线观看播放网| 超碰av人人做人人爽久久 | 精品电影一区二区在线| 露出奶头的视频| 欧美3d第一页| 欧美最新免费一区二区三区 | 久久久国产精品麻豆| 老司机午夜十八禁免费视频| 欧美日本视频| 亚洲自拍偷在线| 午夜亚洲福利在线播放| 午夜久久久久精精品| 亚洲一区二区三区色噜噜| 亚洲精品粉嫩美女一区| 国产精品久久久久久久久免 | 亚洲专区国产一区二区| 欧美一区二区国产精品久久精品| 国产v大片淫在线免费观看| 老司机在亚洲福利影院| 国产黄a三级三级三级人| 国产成人av教育| 欧美乱色亚洲激情| 五月伊人婷婷丁香| 成年人黄色毛片网站| 啦啦啦韩国在线观看视频| 97碰自拍视频| 色播亚洲综合网| 日韩欧美国产在线观看| 97碰自拍视频| 18美女黄网站色大片免费观看| 一本精品99久久精品77| 日韩欧美国产一区二区入口| 免费观看人在逋| 精品久久久久久成人av| 国产高潮美女av| 俄罗斯特黄特色一大片| 日本在线视频免费播放| 大型黄色视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲第一电影网av| 日韩欧美国产一区二区入口| 色精品久久人妻99蜜桃| 午夜日韩欧美国产| 国产精品久久久久久精品电影| 国产亚洲av嫩草精品影院| 亚洲av中文字字幕乱码综合| 91av网一区二区| 久久精品91无色码中文字幕| 午夜日韩欧美国产| 精品国产三级普通话版| 国产精品美女特级片免费视频播放器| 久久精品国产亚洲av涩爱 | 亚洲五月天丁香| 99热这里只有是精品50| 国产成人系列免费观看| 日韩亚洲欧美综合| 在线十欧美十亚洲十日本专区| 久久久成人免费电影| 久久久色成人| 亚洲国产日韩欧美精品在线观看 | 夜夜躁狠狠躁天天躁| 91字幕亚洲| 搡老妇女老女人老熟妇| 国产探花在线观看一区二区| 国产一区在线观看成人免费| 欧美日韩福利视频一区二区| 一个人看的www免费观看视频| 亚洲av美国av| 欧美一级毛片孕妇| 国产高潮美女av| 少妇裸体淫交视频免费看高清| 成人18禁在线播放| 久久人人精品亚洲av| 中亚洲国语对白在线视频| 日韩大尺度精品在线看网址| 国产男靠女视频免费网站| 99国产精品一区二区蜜桃av| 一级毛片高清免费大全| 久久久成人免费电影| 亚洲人成电影免费在线| 又紧又爽又黄一区二区| 噜噜噜噜噜久久久久久91| 欧美极品一区二区三区四区| 国产爱豆传媒在线观看| 最近最新免费中文字幕在线| 天堂av国产一区二区熟女人妻| 一级a爱片免费观看的视频| 高清毛片免费观看视频网站| 99国产精品一区二区三区| 精品国产美女av久久久久小说| 国产真实乱freesex| 亚洲国产欧美网| 久久亚洲精品不卡| 熟女人妻精品中文字幕| 老熟妇乱子伦视频在线观看| 淫秽高清视频在线观看| 亚洲av成人精品一区久久| 国产野战对白在线观看| 日本一二三区视频观看| 日本熟妇午夜| 深夜精品福利| 国产精华一区二区三区| 免费观看精品视频网站| 国产免费男女视频| 日本一本二区三区精品| 三级男女做爰猛烈吃奶摸视频| 91麻豆精品激情在线观看国产| 麻豆国产av国片精品| 啦啦啦观看免费观看视频高清| 国产免费一级a男人的天堂| 亚洲国产欧洲综合997久久,| 久9热在线精品视频| 久久久久久久久中文| 18禁黄网站禁片免费观看直播| 男女视频在线观看网站免费| 黄色丝袜av网址大全| 成人永久免费在线观看视频| 欧美av亚洲av综合av国产av| 欧美成人性av电影在线观看| 色综合欧美亚洲国产小说| 黄色日韩在线| 丰满人妻一区二区三区视频av | 欧美xxxx黑人xx丫x性爽| 免费在线观看亚洲国产| 看片在线看免费视频| 日本免费一区二区三区高清不卡| 国产精品一及| 母亲3免费完整高清在线观看| 国产精品亚洲av一区麻豆| 久久精品亚洲精品国产色婷小说| 日本五十路高清| 偷拍熟女少妇极品色| 亚洲成a人片在线一区二区| 9191精品国产免费久久| 人人妻,人人澡人人爽秒播| 国产精品自产拍在线观看55亚洲| 日本撒尿小便嘘嘘汇集6| 亚洲午夜理论影院| 99国产综合亚洲精品| 亚洲性夜色夜夜综合| 欧美中文日本在线观看视频| 变态另类成人亚洲欧美熟女| 国产精品免费一区二区三区在线| 亚洲av免费高清在线观看| 国产激情偷乱视频一区二区| 久久久久久久久大av| 亚洲成a人片在线一区二区| 欧美三级亚洲精品| 免费人成视频x8x8入口观看| 成人永久免费在线观看视频| 久久久国产精品麻豆| 一本精品99久久精品77| 国产真实伦视频高清在线观看 | 亚洲av一区综合| 亚洲在线自拍视频| 欧美午夜高清在线| 校园春色视频在线观看| 亚洲国产精品999在线| 国产精品嫩草影院av在线观看 | 亚洲国产精品久久男人天堂| 日韩欧美精品v在线| 1024手机看黄色片| 久久香蕉精品热| 69人妻影院| 看片在线看免费视频| 久久精品国产清高在天天线| 精品久久久久久成人av| 精华霜和精华液先用哪个| 18禁黄网站禁片午夜丰满| 九九在线视频观看精品| 国产精品久久久久久人妻精品电影| 天堂网av新在线| 最新中文字幕久久久久| 少妇的丰满在线观看| 欧美一区二区精品小视频在线| 精品久久久久久久末码| 国产精品一区二区三区四区免费观看 | x7x7x7水蜜桃| 一区二区三区高清视频在线| 亚洲五月婷婷丁香| 波野结衣二区三区在线 | 欧美最黄视频在线播放免费| 丰满人妻一区二区三区视频av | 人人妻人人澡欧美一区二区| 日韩有码中文字幕| 不卡一级毛片| 男插女下体视频免费在线播放| 少妇高潮的动态图| 国产高清有码在线观看视频| 欧美日韩国产亚洲二区| 亚洲av电影在线进入| 国产色婷婷99| 亚洲av美国av| 九色国产91popny在线| 最新中文字幕久久久久| 色综合欧美亚洲国产小说| 国产精品久久久久久亚洲av鲁大| bbb黄色大片| 夜夜看夜夜爽夜夜摸| 亚洲黑人精品在线| 香蕉久久夜色| 深爱激情五月婷婷| 国产熟女xx| 麻豆国产97在线/欧美| 欧美大码av| 国产精品香港三级国产av潘金莲| 精品日产1卡2卡| 日本五十路高清| 在线观看免费午夜福利视频| 每晚都被弄得嗷嗷叫到高潮| 午夜精品在线福利| 国产极品精品免费视频能看的| 天天躁日日操中文字幕| 国产精品99久久99久久久不卡| 精品免费久久久久久久清纯| 亚洲va日本ⅴa欧美va伊人久久| 久久久国产精品麻豆| 免费观看精品视频网站| 国产国拍精品亚洲av在线观看 | 中文字幕人成人乱码亚洲影| 久久草成人影院| 嫩草影视91久久| 99在线视频只有这里精品首页| 在线观看免费视频日本深夜| a在线观看视频网站| 国产精品1区2区在线观看.| 国产黄色小视频在线观看| 中文字幕av在线有码专区| 亚洲熟妇熟女久久| 99久久精品一区二区三区| 天堂动漫精品| 欧美一级a爱片免费观看看| 午夜免费观看网址| 欧洲精品卡2卡3卡4卡5卡区| 男女做爰动态图高潮gif福利片| 99久国产av精品| 免费看a级黄色片| 中文字幕精品亚洲无线码一区| 两个人的视频大全免费| 亚洲人成电影免费在线| 91字幕亚洲| 欧美日韩福利视频一区二区| 中文字幕av成人在线电影| 精品一区二区三区视频在线 | 很黄的视频免费| 大型黄色视频在线免费观看| 美女 人体艺术 gogo| 老熟妇仑乱视频hdxx| 毛片女人毛片| 身体一侧抽搐| 欧美日韩福利视频一区二区| 亚洲精品亚洲一区二区| 日本撒尿小便嘘嘘汇集6| 男人舔奶头视频| 热99re8久久精品国产| 悠悠久久av| 亚洲欧美日韩高清专用| 午夜福利高清视频| 一级黄片播放器| 悠悠久久av| 久9热在线精品视频| 色哟哟哟哟哟哟| 给我免费播放毛片高清在线观看| 国产免费一级a男人的天堂| 男人和女人高潮做爰伦理| 熟女电影av网| 久久性视频一级片| 不卡一级毛片| 熟女人妻精品中文字幕| 夜夜看夜夜爽夜夜摸| 国产黄色小视频在线观看| 国产av麻豆久久久久久久| 国产精品久久久久久精品电影| 欧美丝袜亚洲另类 | 99视频精品全部免费 在线| 老鸭窝网址在线观看| 老司机在亚洲福利影院| 亚洲精华国产精华精| 免费av不卡在线播放| 免费看十八禁软件| 日本黄色视频三级网站网址| 精品一区二区三区视频在线 | 国产精品久久久久久久久免 | 欧美黑人欧美精品刺激| 午夜精品一区二区三区免费看| www.999成人在线观看| 国内精品久久久久精免费| 亚洲av成人av| 成人精品一区二区免费| 噜噜噜噜噜久久久久久91| 国产亚洲精品av在线| 国产精品女同一区二区软件 | 日日摸夜夜添夜夜添小说| 亚洲内射少妇av| 看免费av毛片| 一个人免费在线观看的高清视频| 久久久久免费精品人妻一区二区| 深爱激情五月婷婷| 少妇高潮的动态图| 一本久久中文字幕| 日韩欧美 国产精品| 一进一出抽搐动态| 男女午夜视频在线观看| 午夜福利成人在线免费观看| 国产精品av视频在线免费观看| 久久久久精品国产欧美久久久| 国产精品久久久久久精品电影| 亚洲第一电影网av| 亚洲精品色激情综合| 怎么达到女性高潮| 露出奶头的视频| 国产亚洲精品一区二区www| 亚洲av一区综合| 成人高潮视频无遮挡免费网站| 最近最新中文字幕大全免费视频| 婷婷精品国产亚洲av| 成人无遮挡网站| 成年女人看的毛片在线观看| 一级作爱视频免费观看| 免费看a级黄色片| 亚洲第一欧美日韩一区二区三区| 国产蜜桃级精品一区二区三区| 亚洲天堂国产精品一区在线| 十八禁网站免费在线| 天堂av国产一区二区熟女人妻| 噜噜噜噜噜久久久久久91| 熟妇人妻久久中文字幕3abv| 亚洲精华国产精华精| 一区二区三区激情视频| 亚洲人与动物交配视频| 女人高潮潮喷娇喘18禁视频| 精品午夜福利视频在线观看一区| 国内少妇人妻偷人精品xxx网站| 欧美zozozo另类| 亚洲中文日韩欧美视频| 88av欧美| 欧美国产日韩亚洲一区| 精品久久久久久,| 亚洲人与动物交配视频| 精品一区二区三区视频在线观看免费| 亚洲黑人精品在线| 国产精品亚洲美女久久久| 琪琪午夜伦伦电影理论片6080| 欧美大码av| 久久精品夜夜夜夜夜久久蜜豆| 日本撒尿小便嘘嘘汇集6| 国产高潮美女av| 午夜免费激情av| 国产精品1区2区在线观看.| 性色av乱码一区二区三区2| 国产欧美日韩精品一区二区| 亚洲国产欧洲综合997久久,| 中文字幕人妻熟人妻熟丝袜美 | 欧美成人免费av一区二区三区| 国产男靠女视频免费网站| 日韩有码中文字幕| 一进一出抽搐gif免费好疼| 麻豆成人午夜福利视频| 人妻夜夜爽99麻豆av| 好男人在线观看高清免费视频| 欧美日韩黄片免| 亚洲最大成人手机在线| 成年女人看的毛片在线观看| 免费观看的影片在线观看| 动漫黄色视频在线观看| 性色avwww在线观看| 精品不卡国产一区二区三区| 国产精品 国内视频| 激情在线观看视频在线高清| 日本a在线网址| 欧美日韩瑟瑟在线播放| 亚洲av熟女| 色老头精品视频在线观看| 69人妻影院| 757午夜福利合集在线观看| 狠狠狠狠99中文字幕| 免费看光身美女| 国产高清有码在线观看视频| 婷婷亚洲欧美| 日韩av在线大香蕉| 免费无遮挡裸体视频| 一区二区三区免费毛片| 亚洲中文字幕日韩| 午夜老司机福利剧场| 日本黄色视频三级网站网址| 久久久色成人| 久久久久久九九精品二区国产| 欧美色视频一区免费| 全区人妻精品视频| 18美女黄网站色大片免费观看| 黄色成人免费大全| 欧美xxxx黑人xx丫x性爽| 一个人看视频在线观看www免费 | a级一级毛片免费在线观看| 色综合站精品国产| 99精品在免费线老司机午夜| 亚洲,欧美精品.| xxx96com| 欧美午夜高清在线| 男女那种视频在线观看| 极品教师在线免费播放| 在线观看免费视频日本深夜| 在线免费观看不下载黄p国产 | 婷婷精品国产亚洲av在线| 欧美日韩亚洲国产一区二区在线观看| 亚洲熟妇中文字幕五十中出| 久久精品国产亚洲av涩爱 | 欧美性猛交黑人性爽| 免费在线观看日本一区| 中文字幕人妻丝袜一区二区| 亚洲av成人精品一区久久| 国产91精品成人一区二区三区| 国产久久久一区二区三区| 一级a爱片免费观看的视频| 欧美中文日本在线观看视频| 欧美日本亚洲视频在线播放| 国产精品久久久久久久久免 | 亚洲av日韩精品久久久久久密| 757午夜福利合集在线观看| 欧美成狂野欧美在线观看| 午夜福利18| 久久人妻av系列| or卡值多少钱| 精品不卡国产一区二区三区| aaaaa片日本免费| 国产精品久久久久久精品电影| 成人精品一区二区免费| 最近最新免费中文字幕在线| 97碰自拍视频| 高清在线国产一区| 嫩草影院精品99| 欧美日韩综合久久久久久 | av在线天堂中文字幕| 精品不卡国产一区二区三区| 老司机午夜十八禁免费视频| 国产免费av片在线观看野外av| 一区二区三区国产精品乱码| 精品久久久久久久末码| 国产精品影院久久| 国产高清videossex| a级毛片a级免费在线| 亚洲专区国产一区二区| 久久久久性生活片| 精品国产美女av久久久久小说| 18+在线观看网站| 欧美成人一区二区免费高清观看| 人人妻,人人澡人人爽秒播| 尤物成人国产欧美一区二区三区| 精品久久久久久久毛片微露脸| 一级黄色大片毛片| 桃色一区二区三区在线观看| 久久久久久久久久黄片| 两性午夜刺激爽爽歪歪视频在线观看| 日韩亚洲欧美综合| 国产不卡一卡二| 国产精品野战在线观看| 国产精品久久久久久人妻精品电影| 日本免费一区二区三区高清不卡| 国产精品 欧美亚洲| or卡值多少钱| 日韩中文字幕欧美一区二区| 又爽又黄无遮挡网站| 亚洲精品美女久久久久99蜜臀| 欧美日韩瑟瑟在线播放| 操出白浆在线播放| 国产蜜桃级精品一区二区三区| 9191精品国产免费久久| 久久性视频一级片| 欧美在线黄色| 亚洲av中文字字幕乱码综合| 国模一区二区三区四区视频| 婷婷精品国产亚洲av在线| 真实男女啪啪啪动态图| 69人妻影院| 日韩免费av在线播放| 国产精品久久电影中文字幕| 女生性感内裤真人,穿戴方法视频| 国产精品1区2区在线观看.| 在线看三级毛片| 88av欧美| 99久久九九国产精品国产免费| 少妇人妻精品综合一区二区 | 久久久久久国产a免费观看| av福利片在线观看| 欧美丝袜亚洲另类 | 丰满的人妻完整版| 中文资源天堂在线| 久久久色成人| 国产高清激情床上av| 久久精品国产综合久久久| 嫁个100分男人电影在线观看| 搡老岳熟女国产| 亚洲av中文字字幕乱码综合| 18禁裸乳无遮挡免费网站照片| 日韩欧美精品v在线| 国产激情欧美一区二区| 老鸭窝网址在线观看| 最好的美女福利视频网| 日韩大尺度精品在线看网址| 欧美在线一区亚洲| 亚洲av成人不卡在线观看播放网| 国产色爽女视频免费观看| 极品教师在线免费播放| 色av中文字幕| 91麻豆精品激情在线观看国产| 99久久成人亚洲精品观看| 免费电影在线观看免费观看| 天堂网av新在线| 麻豆国产97在线/欧美| 日韩欧美在线乱码| 亚洲无线在线观看| 欧美日本视频| 色播亚洲综合网| 亚洲电影在线观看av| 在线观看日韩欧美| 亚洲真实伦在线观看| 美女高潮的动态| avwww免费|