Fine mapping of a male sterility gene <Emphasis Type="Italic">ms</Emphasis>-<Emphasis Type="Italic">3</Emphasis> in a novel cucumber (<Emphasis Type="Italic">Cucumis sativus</Emphasis> L.) mutant

Key message

The cucumber male sterility gene ms - 3 was fine mapped in a 76 kb region harboring an MMD1 -like gene Csa3M006660 that may be responsible for the male sterile in cucumber.

Abstract

A cucumber (Cucumis sativus L.) male sterile mutant (ms-3) in an advanced-generation inbred line was identified, and genetic analysis revealed that the male sterility trait was controlled by a recessive nuclear gene, ms-3, which was stably inherited. Histological studies suggested that the main cause of the male sterility was defective microsporogenesis, resulting in no tetrad or microspores being formed. Bulked segregant analysis (BSA) and genotyping of an F₂ population of 2553 individuals were employed used to fine map ms-3, which was delimited to a 76 Kb region. In this region, a single non-synonymous SNP was found in the Csa3M006660 gene locus, which was predicted to result in an amino acid change. Quantitative RT-PCR analysis of Csa3M006660 was consistent with the fact that it plays a role in the early development of cucumber pollen. The protein encoded by Csa3M006660 is predicted to be homeodomain (PHD) finger protein, and the high degree of sequence conservation with homologs from a range of plant species further suggested the importance of the ms-3 non-synonymous mutation. The data presented here provide support for Csa3M006660 as the most likely candidate gene for Ms-3.

     

Genetic composition of interspecific potato somatic hybrids and autofused 4<Emphasis Type="Italic">x</Emphasis> plants evaluated by DArT and cytoplasmic DNA markers

Key message

Using DArT analysis, we demonstrated that all Solanum × michoacanum (+) S. tuberosum somatic hybrids contained all parental chromosomes. However, from 13.9 to 29.6 % of the markers from both parents were lost in the hybrids.

Abstract

Somatic hybrids are an interesting material for research of nucleus-cytoplasm interaction and sources of new nuclear and cytoplasmic combinations. Analyses of genomes of somatic hybrids are essential for studies on genome compatibility between species, its evolution and are important for their efficient exploitation. Diversity array technology (DArT) permits analysis of the composition of nuclear DNA of somatic hybrids. The nuclear genome compositions of 97 Solanum × michoacanum (+) S. tuberosum [mch (+) tbr] somatic hybrids from five fusion combinations and 11 autofused 4x mch were analyzed for the first time based on DArT markers. Out of 5358 DArT markers generated in a single assay, greater than 2000 markers were polymorphic between parents, of which more than 1500 have a known chromosomal location on potato genetic or physical map. DArT markers were distributed along the entire length of 12 chromosomes. We noticed elimination of markers of wild and tbr fusion components. The nuclear genome of individual somatic hybrids was diversified. Mch is a source of resistance to Phytophthora infestans. From 97 mch (+) tbr somatic hybrids, two hybrids and all 11 autofused 4x mch were resistant to P. infestans. The analysis of the structure of particular hybrids’ chromosomes indicated the presence of markers from both parental genomes as well as missing markers spread along the full length of the chromosome. Markers specific to chloroplast DNA and mitochondrial DNA were used for analysis of changes within the organellar genomes of somatic hybrids. Random and non-random segregations of organellar DNA were noted.

     

Development of <Emphasis Type="Italic">Brassica oleracea</Emphasis>-<Emphasis Type="Italic">nigra</Emphasis> monosomic alien addition lines: genotypic,cytological and morphological analyses

Key message

We report the development and characterization of Brassica oleracea - nigra monosomic alien addition lines (MAALs) to dissect the Brassica B genome.

Abstract

Brassica nigra (2n = 16, BB) represents the diploid Brassica B genome which carries many useful genes and traits for breeding but received limited studies. To dissect the B genome from B. nigra, the triploid F₁ hybrid (2n = 26, CCB) obtained previously from the cross B. oleracea var. alboglabra (2n = 18, CC) × B. nigra was used as the maternal parent and backcrossed successively to parental B. oleracea. The progenies in BC₁ to BC₃ generations were analyzed by the methods of FISH and SSR markers to screen the monosomic alien addition lines (MAALs) with each of eight different B-genome chromosomes added to C genome (2n = 19, CC + 1B_1?8), and seven different MAALs were established, except for the one with chromosome B2 which existed in one triple addition. Most of these MAALs were distinguishable morphologically from each other, as they expressed the characters from B. nigra differently and at variable extents. The alien chromosome remained unpaired as a univalent in 86.24% pollen mother cells at diakinesis or metaphase I, and formed a trivalent with two C-genome chromosomes in 13.76% cells. Transmission frequency of all the added chromosomes was far higher through the ovules (averagely 14.40%) than the pollen (2.64%). The B1, B4 and B5 chromosomes were transmitted by female at much higher rates (22.38–30.00%) than the other four (B3, B6, B7, B8) (5.04–8.42%). The MAALs should be valuable for exploiting the genome structure and evolution of B. nigra.

     

The seed dormancy allele <Emphasis Type="Italic">TaSdr</Emphasis>-<Emphasis Type="Italic">A1a</Emphasis> associated with pre-harvest sprouting tolerance is mainly present in Chinese wheat landraces

Key message

We cloned TaSdr - A1 gene, and developed a gene-specific marker for TaSdr - A1 . A QTL for germination index at the TaSdr - A1 locus was identified in the Yangxiaomai/Zhongyou 9507 RIL population.

Abstract

Pre-harvest sprouting (PHS) affects yield and end-use quality in bread wheat (Triticum aestivum L.). In the present study we found an association between the TaSdr-A1 gene and PHS tolerance in bread wheat. TaSdr-A1 on chromosome 2A was cloned using a homologous cloning approach. Sequence analysis of TaSdr-A1 revealed an SNP at position 643, with the G allele being present in genotypes with lower germination index (GI) values and A in those with higher GI. These alleles were designated as TaSdr-A1a and TaSdr-A1b, respectively. A cleaved amplified polymorphism sequence (CAPS) marker Sdr2A based on the SNP was developed, and linkage mapping and QTL analysis were conducted to confirm the association between TaSdr-A1 and seed dormancy. Sdr2A was located in a 2.9 cM interval between SSR markers Xgwm95 and Xgwm372. A QTL for GI at the TaSdr-A1 locus explained 6.6, 7.3, and 8.2 % of the phenotypic variances in a Yangxiaomai/Zhongyou 9507 RIL population grown at Beijing, Shijiazhuang, and the averaged data from the two environments, respectively. Two sets of Chinese wheat cultivars used for validating the TaSdr-A1 polymorphism and the corresponding gene-specific marker Sdr2A showed that TaSdr-A1 was significantly associated with GI. Among 29 accessions with TaSdr-A1a, 24 (82.8 %) were landraces, indicating the importance of Chinese wheat landraces as sources of PHS tolerance. This study identified a novel PHS resistance allele TaSdr-A1a mainly presented in Chinese landraces and it is likely to be the causal gene for QPhs.ccsu-2A.3, providing new information for an understanding of seed dormancy.

     

Growth and nodulation in <Emphasis Type="Italic">Alnus sieboldiana</Emphasis> in response to <Emphasis Type="Italic">Frankia</Emphasis> inoculation and nitrogen treatments

Key message

We investigated a Frankia – Alnus sieboldiana symbiosis, including the minimum inoculum dose for constant nodulation, the period of time to nodulation after inoculation, and the effects of N on nodulation.

Abstract

Frankia is a nitrogen-fixing actinomycete that forms root nodules in some dicotyledonous plants, which are called actinorhizal. We studied nodule formation in Alnus sieboldiana, an actinorhizal plant, after inoculation with a Frankia isolate to establish techniques for Frankia inoculation and the cultivation of inoculated plants. Root nodules formed on seedlings of A. sieboldiana by 2 weeks after inoculation, and N₂ fixation measured by acetylene reduction activity started 3 weeks after inoculation. Nodulation was observed after inoculation with a Frankia isolate at 0.001 μL packed cell volume (pcv). The number of nodules formed on the seedlings inoculated with Frankia at more than 0.05 μL pcv was not significantly different. Nodule development and N₂ fixation were reduced when inoculated seedlings were treated weekly with 15 mM NH₄NO₃-N. In contrast, treatment with 3.75 or 0.9375 mM NH₄NO₃-N did not inhibit nodule development or N₂ fixation of inoculated seedlings by 15 weeks of N treatment.

     

Genetic diversity and population structure of <Emphasis Type="Italic">Melia azedarach</Emphasis> in North-Western Plains of India

Key message

Genetic structure among M. azedarach populations was detected and two subpopulations were present among them. A significant ‘isolation by distance’ was found in M. azedarach population in North-Western Plains of India.

Abstract

Melia azedarach is an important forest tree with pharmaceutical, insecticidal, pesticidal, and commercial significance. It is a good reforestation tree because of its fast growth and drought hardy nature. Genetic variation in a species allows itself to adapt, evolve and respond to environmental stress. It provides the basis for survival of a species and critically influences its evolutionary potential. Assessment of genetic diversity is necessary for improvement and conservation of a species. For this, microsatellite markers are of particular interest given the attributes like co-dominance, reproducibility, hyper variability and abundance throughout the genome. In the present study, we analyzed the genetic diversity and population structure of M. azedarach, an ecologically imperative species growing in the North-Western Plains of India. We developed 43 microsatellite markers, of which 20 were subsequently employed for analysis of diversity and population structure among 33 populations encompassing 318 genotypes representing North-Western Plains of India. A moderate level of diversity (Na = 5.1, Ho = 0.506, He = 0.712, I = 1.386) was assessed. The highest value of ΔK estimated using STRUCTURE indicated 2 subpopulations (K = 2). AMOVA exhibited 73 % variation within populations and 12 % variation was found among regions. Significant positive correlation between geographical and genetic distance was found (Rxy = 0.365, P = 0.010). The present study lays a foundation on a better understanding of genetic dynamics of the species and reveals its diversity and population structure in North-Western Plains of India.

     

Rapid and efficient <Emphasis Type="Italic">Agrobacterium</Emphasis>-mediated transformation of sorghum (<Emphasis Type="Italic">Sorghum bicolor</Emphasis>) employing standard binary vectors and <Emphasis Type="Italic">bar</Emphasis> gene as a selectable marker

Key message

A rapid and efficient Agrobacterium -mediated transformation system in sorghum has been developed employing standard binary vectors and bar gene as a selectable marker.

Abstract

Sorghum (Sorghum bicolor) is an important food and biofuel crop worldwide, for which improvements in genetic transformation are needed to study its biology and facilitate agronomic and commercial improvement. Here, we report optimization of regeneration and transformation of public sorghum genotype P898012 using standard binary vectors and bar gene as a selectable marker. The tissue culture regeneration time frame has been reduced to 7–12 weeks with a yield of over 18 plants per callus, and the optimized transformation system employing Agrobacterium tumefaciens strain AGL1 and the bar with a MAS promoter achieved an average frequency over 14 %. Of randomly analyzed independent transgenic events, 40–50 % carry single copy of integrated T-DNA. Some independent transgenic events were derived from the same embryogenic callus lines, but a 3:1 Mendelian segregation ratio was found in all transgenic events with single copy as estimated by Southern blots. The system described here should facilitate studies of sorghum biology and agronomic improvement.

     

Allometric biomass,nutrient and carbon stock models for <Emphasis Type="Italic">Kandelia candel</Emphasis> of the Sundarbans,Bangladesh

Key message

Present study recommends DBH as independent variable of the derived allometric models and Biomass = a + b DBH ² has been selected for total above-ground biomass, nutrients and carbon stock.

Abstract

Kandelia candel (L.) Druce is a shrub to small tree of the Sundarbans mangrove forest of Bangladesh. The aim of the study was to derive the allometric models for estimating above-ground biomass, nutrient and carbon stock in K. candel. A total of eight linear models with 64 regression equations were tested to derive the allometric models for biomass of each part of plant; and nutrients and carbon stock in total above-ground biomass. The best fitted allometric models were selected by considering the values of R ², CV, R _mse, MS_error, S _a, S _b, F value, AICc and Furnival Index. The selected allometric models were Biomass = 0.014 DBH² + 0.03; √Biomass = 0.29 DBH ? 0.21; √Biomass = 0.66 √DBH ? 0.57; √Biomass = 1.19 √DBH ? 1.02; Biomass = 0.21 DBH² + 0.12 for leaves, branches, bark, stem without bark and total above-ground biomass, respectively. The selected allometric models for Nitrogen, Phosphorous, Potassium and Carbon stock in total above-ground biomass were N = 0.39 DBH² + 0.49, P = 0.77 DBH² + 0.14, K = 0.87 DBH² + 0.07 and C = 0.09 DBH² + 0.05, respectively. The derived allometric models have included DBH as a single independent variable, which may give quick and accurate estimation of the above-ground biomass, nutrient and carbon stock in this species. This information may also contribute to a broader study of nutrient cycling, nutrient budgeting and carbon sequestration of the studied forest.

     

Revision of the <Emphasis Type="Italic">Rhynchospora glomerata</Emphasis> species complex,focusing on the taxonomic status of <Emphasis Type="Italic">R. leptocarpa</Emphasis> (Cyperaceae)

Rhynchospora glomerata and its closest relatives comprise a group of beakesedges widespread and frequent in much of North America. The classification of the R. glomerata complex remains unresolved and controversial. The goals of this study are to determine the number of taxa in the complex and their ranks, and identify their best diagnostic characters. Measurements of eight characters from each of 101 specimens from throughout the geographic range of the complex furnished data for morphometric analyses. These analyses reveal the R. glomerata complex contains three species and no infraspecific taxa: R. capitellata, R. glomerata, and R. leptocarpa. We detected 10 validly published basionyms in the complex, five of which required lectotypification. Accordingly, we designated lectotypes for R. glomerata var. discutiens, R. glomerata var. minor, R. glomerata var. paniculata, and R. glomerata var. robustior, and the second-step lectotype for R. capitellata var. controversa.     

Characterisation of a novel quantitative trait locus, <Emphasis Type="Italic">GN4</Emphasis>-<Emphasis Type="Italic">1</Emphasis>, for grain number and yield in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Key message

A novel QTL for grain number, GN4-1, was identified and fine-mapped to an ~ 190-kb region on the long arm of rice chromosome 4.

Abstract

Rice grain yield is primarily determined by three components: number of panicles per plant, grain number per panicle and grain weight. Among these traits, grain number per panicle is the major contributor to grain yield formation and is a crucial trait for yield improvement. In this study, we identified a major quantitative trait locus (QTL) responsible for rice grain number on chromosome 4, designated GN4-1 (a QTL for Grain Number on chromosome 4), using advanced segregating populations derived from the crosses between an elite indica cultivar ‘Zhonghui 8006’ (ZH8006) and a japonica rice ‘Wuyunjing 8’ (WYJ8). GN4-1 was delimited to an ~ 190-kb region on chromosome 4. The genetic effect of GN4-1 was estimated using a pair of near-isogenic lines. The GN4-1 gene from WYJ8 promoted accumulation of cytokinins in the inflorescence and increased grain number per panicle by ~ 17%. More importantly, introduction of the WYJ8 GN4-1 gene into ZH8006 increased grain yield by ~ 14.3 and ~ 11.5% in the experimental plots in 2014 and 2015, respectively. In addition, GN4-1 promoted thickening of the culm and may enhance resistance to lodging. These results demonstrate that GN4-1 is a potentially valuable gene for improvement of yield and lodging resistance in rice breeding.

     

Update on <Emphasis Type="Italic">transparent testa</Emphasis> mutants from <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>: characterisation of new alleles from an isogenic collection

Main conclusion

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

     

Fine mapping of the chromosome 5B region carrying closely linked rust resistance genes <Emphasis Type="Italic">Yr47</Emphasis> and <Emphasis Type="Italic">Lr52</Emphasis> in wheat

Key message

Fine mapping of Yr47 and Lr52 in chromosome arm 5BS of wheat identified close linkage of the marker sun180 to both genes and its robustness for marker-assisted selection was demonstrated.

Abstract

The widely effective and genetically linked rust resistance genes Yr47 and Lr52 have previously been mapped in the short arm of chromosome 5B in two F₃ populations (Aus28183/Aus27229 and Aus28187/Aus27229). The Aus28183/Aus27229 F₃ population was advanced to generate an F₆ recombinant inbred line (RIL) population to identify markers closely linked with Yr47 and Lr52. Diverse genomic resources including flow-sorted chromosome survey sequence contigs representing the orthologous region in Brachypodium distachyon, the physical map of chromosome arm 5BS, expressed sequence tags (ESTs) located in the 5BS6-0.81-1.00 deletion bin and resistance gene analog contigs of chromosome arm 5BS were used to develop markers to saturate the target region. Selective genotyping was also performed using the iSelect 90 K Infinium wheat SNP assay. A set of SSR, STS, gene-based and SNP markers were developed and genotyped on the Aus28183/Aus27229 RIL population. Yr47 and Lr52 are genetically distinct genes that mapped 0.4 cM apart in the RIL population. The SSR marker sun180 co-segregated with Lr52 and mapped 0.4 cM distal to Yr47. In a high resolution mapping population of 600 F₂ genotypes Yr47 and Lr52 mapped 0.2 cM apart and marker sun180 was placed 0.4 cM distal to Lr52. The amplification of a different sun180 amplicon (195 bp) than that linked with Yr47 and Lr52 (200 bp) in 204 diverse wheat genotypes demonstrated its robustness for marker-assisted selection of these genes.

     

Glandular trichomes and essential oil characteristics of in vitro propagated <Emphasis Type="Italic">Micromeria pulegium</Emphasis> (Rochel) Benth. (Lamiaceae)

Main conclusion

In vitro conditions and benzyladenine influenced both content and composition of micropropagated Micromeria pulegium essential oils, with pulegone and menthone being the main essential oil components. The content and chemical composition of Micromeria pulegium (Rochel) Benth. essential oils were studied in native plant material at vegetative stage and in micropropagated plants, obtained from nodal segments cultured on solid MS medium supplemented with N⁶–benzyladenine (BA) or kinetin at different concentrations, alone or in combination with indole-3-acetic acid. Shoot proliferation was achieved in all treatments, but the highest biomass production was obtained after treatment with 10 μM BA. Phytochemical analysis identified up to 21 compounds in the essential oils of wild-growing and in vitro cultivated plants, both showing very high percentages of total monoterpenoids dominated by oxygenated monoterpenes of the menthane type. Pulegone and menthone were the main essential oil components detected in both wild-growing plants (60.07 and 26.85 %, respectively) and micropropagated plants grown on either plant growth regulator-free medium (44.57 and 29.14 %, respectively) or BA-supplemented medium (50.77 and 14.45 %, respectively). The percentage of total sesquiterpenoids increased in vitro, particularly owing to sesquiterpene hydrocarbons that were not found in wild-growing plants. Differences in both content and the composition of the essential oils obtained from different samples indicated that in vitro culture conditions and plant growth regulators significantly influence the essential oils properties. In addition, the morphology and structure of M. pulegium glandular trichomes in relation to the secretory process were characterized for the first time using SEM and light microscopy, and their secretion was histochemically analyzed.

     

Complete chloroplast genome sequences of <Emphasis Type="Italic">Solanum commersonii</Emphasis> and its application to chloroplast genotype in somatic hybrids with <Emphasis Type="Italic">Solanum tuberosum</Emphasis>

Key message

Chloroplast genome of Solanum commersonii and S olanum tuberosum were completely sequenced, and Indel markers were successfully applied to distinguish chlorotypes demonstrating the chloroplast genome was randomly distributed during protoplast fusion.

Abstract

Somatic hybridization has been widely employed for the introgression of resistance to several diseases from wild Solanum species to overcome sexual barriers in potato breeding. Solanum commersonii is a major resource used as a parent line in somatic hybridization to improve bacterial wilt resistance in interspecies transfer to cultivated potato (S. tuberosum). Here, we sequenced the complete chloroplast genomes of Lz3.2 (S. commersonii) and S. tuberosum (PT56), which were used to develop fusion products, then compared them with those of five members of the Solanaceae family, S. tuberosum, Capsicum annum, S. lycopersicum, S. bulbocastanum and S. nigrum and Coffea arabica as an out-group. We then developed Indel markers for application in chloroplast genotyping. The complete chloroplast genome of Lz3.2 is composed of 155,525 bp, which is larger than the PT56 genome with 155,296 bp. Gene content, order and orientation of the S. commersonii chloroplast genome were highly conserved with those of other Solanaceae species, and the phylogenetic tree revealed that S. commersonii is located within the same node of S. tuberosum. However, sequence alignment revealed nine Indels between S. commersonii and S. tuberosum in their chloroplast genomes, allowing two Indel markers to be developed. The markers could distinguish the two species and were successfully applied to chloroplast genotyping (chlorotype) in somatic hybrids and their progenies. The results obtained in this study confirmed the random distribution of the chloroplast genome during protoplast fusion and its maternal inheritance and can be applied to select proper plastid genotypes in potato breeding program.

     

<Emphasis Type="Italic">qRfg3</Emphasis>, a novel quantitative resistance locus against <Emphasis Type="Italic">Gibberella</Emphasis> stalk rot in maize

Key message

A quantitative trait locus  qRfg3 imparts recessive resistance to maize Gibberella stalk rot. qRfg3 has been mapped into a 350-kb interval and could reduce the disease severity index by ~26.6%.

Abstract

Gibberella stalk rot, caused by the fungal pathogen Fusarium graminearum, severely affects maize yield and grain quality worldwide. To identify more resistance quantitative trait loci (QTLs) against this disease, we analyzed a recombinant inbred line (RIL) population derived from a cross between resistant H127R and susceptible C7-2 inbred lines. Within this population, maize resistance to Gibberella stalk rot had high broad-sense heritability. A major QTL, qRfg3, on chromosome 3 was consistently detected across three field trials, accounting for 10.7–19.4% of the total phenotypic variation. Using a progeny-based sequential fine-mapping strategy, we narrowed qRfg3 down to an interval of ~350 kb. We further demonstrated that qRfg3 is a recessive resistance locus to Gibberella stalk rot that reduced the disease severity index by ~26.6%. Both the gene location and recessive genetic mode distinguish qRfg3 from other stalk rot resistance loci. Hence, qRfg3 is valuable as a complement to existing resistance QTLs to improve maize resistance to Gibberella stalk rot.

     

Site-directed mutagenesis in <Emphasis Type="Italic">Petunia</Emphasis> × <Emphasis Type="Italic">hybrida</Emphasis> protoplast system using direct delivery of purified recombinant Cas9 ribonucleoproteins

Key message

Site-directed mutagenesis of nitrate reductase genes using direct delivery of purified Cas9 protein preassembled with guide RNA produces mutations efficiently in Petunia × hybrida protoplast system.

Abstract

The clustered, regularly interspaced, short palindromic repeat (CRISPR)-CRISPR associated endonuclease 9 (CRISPR/Cas9) system has been recently announced as a powerful molecular breeding tool for site-directed mutagenesis in higher plants. Here, we report a site-directed mutagenesis method targeting Petunia nitrate reductase (NR) gene locus. This method could create mutations efficiently using direct delivery of purified Cas9 protein and single guide RNA (sgRNA) into protoplast cells. After transient introduction of RNA-guided endonuclease (RGEN) ribonucleoproteins (RNPs) with different sgRNAs targeting NR genes, mutagenesis at the targeted loci was detected by T7E1 assay and confirmed by targeted deep sequencing. T7E1 assay showed that RGEN RNPs induced site-specific mutations at frequencies ranging from 2.4 to 21 % at four different sites (NR1, 2, 4 and 6) in the PhNR gene locus with average mutation efficiency of 14.9 ± 2.2 %. Targeted deep DNA sequencing revealed mutation rates of 5.3–17.8 % with average mutation rate of 11.5 ± 2 % at the same NR gene target sites in DNA fragments of analyzed protoplast transfectants. Further analysis from targeted deep sequencing showed that the average ratio of deletion to insertion produced collectively by the four NR-RGEN target sites (NR1, 2, 4, and 6) was about 63:37. Our results demonstrated that direct delivery of RGEN RNPs into protoplast cells of Petunia can be exploited as an efficient tool for site-directed mutagenesis of genes or genome editing in plant systems.