Haplotype determination of the Bombyx mori nucleopolyhedrovirus by Nanopore sequencing and linkage of single nucleotide variants.

Naturally occurring isolates of baculoviruses, such as the Bombyx mori nucleopolyhedrovirus (BmNPV), usually consist of numerous genetically different haplotypes. Deciphering the different haplotypes of such isolates is hampered by the large size of the dsDNA genome, as well as the short read length of next generation sequencing (NGS) techniques that are widely applied for baculovirus isolate characterization. In this study, we addressed this challenge by combining the accuracy of NGS to determine single nucleotide variants (SNVs) as genetic markers with the long read length of Nanopore sequencing technique. This hybrid approach allowed the comprehensive analysis of genetically homogeneous and heterogeneous isolates of BmNPV. Specifically, this allowed the identification of two putative major haplotypes in the heterogeneous isolate BmNPV-Ja by SNV position linkage. SNV positions, which were determined based on NGS data, were linked by the long Nanopore reads in a Position Weight Matrix. Using a modified Expectation-Maximization algorithm, the Nanopore reads were assigned according to the occurrence of variable SNV positions by machine learning. The cohorts of reads were de novo assembled, which led to the identification of BmNPV haplotypes. The method demonstrated the strength of the combined approach of short- and long-read sequencing techniques to decipher the genetic diversity of baculovirus isolates.

The crane fly glycosylated triketide δ-lactone cornicinine elicits akinete differentiation of the cyanobiont in aquatic Azolla fern symbioses.

The restriction of plant-symbiont dinitrogen fixation by an insect semiochemical had not been previously described. Here we report on a glycosylated triketide δ-lactone from Nephrotoma cornicina crane flies, cornicinine, that causes chlorosis in the floating-fern symbioses from the genus Azolla. Only the glycosylated trans-A form of chemically synthesized cornicinine was active: 500 nM cornicinine in the growth medium turned all cyanobacterial filaments from Nostoc azollae inside the host leaf-cavities into akinetes typically secreting CTB-bacteriocins. Cornicinine further inhibited akinete germination in Azolla sporelings, precluding re-establishment of the symbiosis during sexual reproduction. It did not impact development of the plant Arabidopsis thaliana or several free-living cyanobacteria from the genera Anabaena or Nostoc but affected the fern host without cyanobiont. Fern-host mRNA sequencing from isolated leaf cavities confirmed high NH4-assimilation and proanthocyanidin biosynthesis in this trichome-rich tissue. After cornicinine treatment, it revealed activation of Cullin-RING ubiquitin-ligase-pathways, known to mediate metabolite signaling and plant elicitation consistent with the chlorosis phenotype, and increased JA-oxidase, sulfate transport and exosome formation. The work begins to uncover molecular mechanisms of cyanobiont differentiation in a seed-free plant symbiosis important for wetland ecology or circular crop-production today, that once caused massive CO2 draw-down during the Eocene geological past.

Advancing pathogen surveillance by nanopore sequencing and genotype characterization of Acheta domesticus densovirus in mass-reared house crickets.

Rapid and reliable detection of pathogens is crucial to complement the growing industry of mass-reared insects, in order to safeguard the insect colonies from outbreak of diseases, which may cause significant economic loss. Current diagnostic methods are mainly based on conventional PCR and microscopic examination, requiring prior knowledge of disease symptoms and are limited to identifying known pathogens. Here, we present a rapid nanopore-based metagenomics approach for detecting entomopathogens from the European house cricket (Acheta domesticus). In this study, the Acheta domesticus densovirus (AdDV) was detected from diseased individuals using solely Nanopore sequencing. Virus reads and genome assemblies were obtained within twenty-four hours after sequencing. Subsequently, due to the length of the Nanopore reads, it was possible to reconstruct significantly large parts or even the entire AdDV genome to conduct studies for genotype identification. Variant analysis indicated the presence of three AdDV genotypes within the same house cricket population, with association to the vital status of the diseased crickets. This contrast provided compelling evidence for the existence of non-lethal AdDV genotypes. These findings demonstrated nanopore-based metagenomics sequencing as a powerful addition to the diagnostic tool kit for routine pathogen surveillance and diagnosis in the insect rearing industry.

The structure of the tetraploid sour cherry 'Schattenmorelle' (Prunus cerasus L.) genome reveals insights into its segmental allopolyploid nature.

Sour cherry (Prunus cerasus L.) is an important allotetraploid cherry species that evolved in the Caspian Sea and Black Sea regions from a hybridization of the tetraploid ground cherry (Prunus fruticosa Pall.) and an unreduced pollen of the diploid sweet cherry (P. avium L.) ancestor. Details of when and where the evolution of this species occurred are unclear, as well as the effect of hybridization on the genome structure. To gain insight, the genome of the sour cherry cultivar 'Schattenmorelle' was sequenced using Illumina NovaSeqTM and Oxford Nanopore long-read technologies, resulting in a ~629-Mbp pseudomolecule reference genome. The genome could be separated into two subgenomes, with subgenome PceS_a originating from P. avium and subgenome PceS_f originating from P. fruticosa. The genome also showed size reduction compared to ancestral species and traces of homoeologous sequence exchanges throughout. Comparative analysis confirmed that the genome of sour cherry is segmental allotetraploid and evolved very recently in the past.

Finding needles in a haystack: identification of inter-specific introgressions in wheat genebank collections using low-coverage sequencing data.

Recently, entire genebank collections of wheat have been extensively characterized with sequencing data. We have identified introgressions using these genotyping-by-sequencing and whole-genome sequencing data. On the basis of our results, we provide information about predicted introgressions at 1-Mb resolution for 9,172 wheat samples as a resource for breeders and scientists. We recommend that all plant genetic resources, including genebank collections, be characterized using a combination of variant calling and introgression prediction. This is necessary to identify potential duplicates in collections efficiently and reliably, and to select promising germplasms with potentially beneficial introgressions for further characterization and prospective breeding application.

The Australasian dingo archetype: de novo chromosome-length genome assembly, DNA methylome, and cranial morphology.

BACKGROUND: One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long-read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed. FINDINGS: We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on chromosomes 11, 16, 25, and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and 9 previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mitochondrial DNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified 2 differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphologic data, comprising geometric morphometric assessment of cranial morphology, place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue shows she had a larger cranial capacity than a similar-sized domestic dog. CONCLUSIONS: These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphologic characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

The Australasian dingo archetype: De novo chromosome-length genome assembly, DNA methylome, and cranial morphology.

Background: One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed. Findings: We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on Chromosomes 11, 16, 25 and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and nine previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mtDNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified two differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphological data, comprising geometric morphometric assessment of cranial morphology place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue show she had a larger cranial capacity than a similar-sized domestic dog. Conclusions: These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphological characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

Insights into breeding history, hotspot regions of selection, and untapped allelic diversity for bread wheat breeding.

Breeding has increasingly altered the genetics of crop plants since the domestication of their wild progenitors. It is postulated that the genetic diversity of elite wheat breeding pools is too narrow to cope with future challenges. In contrast, plant genetic resources (PGRs) of wheat stored in genebanks are valuable sources of unexploited genetic diversity. Therefore, to ensure breeding progress in the future, it is of prime importance to identify the useful allelic diversity available in PGRs and to transfer it into elite breeding pools. Here, a diverse collection consisting of modern winter wheat cultivars and genebank accessions was investigated based on reduced-representation genomic sequencing and an iSelect single nucleotide polymorphism (SNP) chip array. Analyses of these datasets provided detailed insights into population structure, levels of genetic diversity, sources of new allelic diversity, and genomic regions affected by breeding activities. We identified 57 regions representing genomic signatures of selection and 827 regions representing private alleles associated exclusively with genebank accessions. The presence of known functional wheat genes, quantitative trait loci, and large chromosomal modifications, i.e., introgressions from wheat wild relatives, provided initial evidence for putative traits associated within these identified regions. These findings were supported by the results of ontology enrichment analyses. The results reported here will stimulate further research and promote breeding in the future by allowing for the targeted introduction of novel allelic diversity into elite wheat breeding pools.

Comparative genomics reveals low levels of inter- and intraspecies diversity in the causal agents of dwarf and common bunt of wheat and hint at conspecificity of Tilletia caries and T. laevis.

Tilletia caries and T. laevis, which are the causal agents of common bunt, as well as T. controversa, which causes dwarf bunt of wheat, threaten especially organic wheat farming. The three closely related fungal species differ in their teliospore morphology and partially in their physiology and infection biology. The gene content as well as intraspecies variation in these species and the genetic basis of their separation is unknown. We sequenced the genome of four T. caries, five T. controversa, and two T. laevis and extended this dataset with five publicly available ones. The genomes of the three species displayed microsynteny with up to 94.3% pairwise aligned regions excluding repetitive regions. The majority of functionally characterized genes involved in pathogenicity, life cycle, and infection of corn smut, Ustilago maydis, were found to be absent or poorly conserved in the draft genomes and the biosynthetic pathway for trimethylamine in Tilletia spp. could be different from bacteria. Overall, 75% of the identified protein-coding genes comprising 84% of the total predicted carbohydrate utilizing enzymes, 72.5% putatively secreted proteins, and 47.4% of effector-like proteins were conserved and shared across all 16 isolates. We predicted nine highly identical secondary metabolite biosynthesis gene clusters comprising in total 62 genes in all species and none were species-specific. Less than 0.1% of the protein-coding genes were species-specific and their function remained mostly unknown. Tilletia controversa had the highest intraspecies genetic variation, followed by T. caries and the lowest in T. laevis. Although the genomes of the three species are very similar, employing 241 single copy genes T. controversa was phylogenetically distinct from T. caries and T. laevis, however these two could not be resolved as individual monophyletic groups. This was in line with the genome-wide number of single nucleotide polymorphisms and small insertions and deletions. Despite the conspicuously different teliospore ornamentation of T. caries and T. laevis, a high degree of genomic identity and scarcity of species-specific genes indicate that the two species could be conspecific.

The Australian dingo is an early offshoot of modern breed dogs.

Dogs are uniquely associated with human dispersal and bring transformational insight into the domestication process. Dingoes represent an intriguing case within canine evolution being geographically isolated for thousands of years. Here, we present a high-quality de novo assembly of a pure dingo (CanFam_DDS). We identified large chromosomal differences relative to the current dog reference (CanFam3.1) and confirmed no expanded pancreatic amylase gene as found in breed dogs. Phylogenetic analyses using variant pairwise matrices show that the dingo is distinct from five breed dogs with 100% bootstrap support when using Greenland wolf as the outgroup. Functionally, we observe differences in methylation patterns between the dingo and German shepherd dog genomes and differences in serum biochemistry and microbiome makeup. Our results suggest that distinct demographic and environmental conditions have shaped the dingo genome. In contrast, artificial human selection has likely shaped the genomes of domestic breed dogs after divergence from the dingo.

Genome wide association study of frost tolerance in wheat.

Winter wheat growing areas in the Northern hemisphere are regularly exposed to heavy frost. Due to the negative impact on yield, the identification of genetic factors controlling frost tolerance (FroT) and development of tools for breeding is of prime importance. Here, we detected QTL associated with FroT by genome wide association studies (GWAS) using a diverse panel of 276 winter wheat genotypes that was phenotyped at five locations in Germany and Russia in three years. The panel was genotyped using the 90 K iSelect array and SNPs in FroT candidate genes. In total, 17,566 SNPs were used for GWAS resulting in the identification of 53 markers significantly associated (LOD ≥ 4) to FroT, corresponding to 23 QTL regions located on 11 chromosomes (1A, 1B, 2A, 2B, 2D, 3A, 3D, 4A, 5A, 5B and 7D). The strongest QTL effect confirmed the importance of chromosome 5A for FroT. In addition, to our best knowledge, eight FroT QTLs were discovered for the first time in this study comprising one QTL on chromosomes 3A, 3D, 4A, 7D and two on chromosomes 1B and 2D. Identification of novel FroT candidate genes will help to better understand the FroT mechanism in wheat and to develop more effective combating strategies.

Detecting major introgressions in wheat and their putative origins using coverage analysis.

Introgressions from crop wild relatives (CWRs) have been used to introduce beneficial traits into cultivated plants. Introgressions have traditionally been detected using cytological methods. Recently, single nucleotide polymorphism (SNP)-based methods have been proposed to detect introgressions in crosses for which both parents are known. However, for unknown material, no method was available to detect introgressions and predict the putative donor species. Here, we present a method to detect introgressions and the putative donor species. We demonstrate the utility of this method using 10 publicly available wheat genome sequences and identify nine major introgressions. We show that the method can distinguish different introgressions at the same locus. We trace introgressions to early wheat cultivars and show that natural introgressions were utilised in early breeding history and still influence elite lines today. Finally, we provide evidence that these introgressions harbour resistance genes.

The draft chromosome-level genome assembly of tetraploid ground cherry (Prunus fruticosa Pall.) from long reads.

Cherries are stone fruits and belong to the economically important plant family of Rosaceae with worldwide cultivation of different species. The ground cherry, Prunus fruticosa Pall., is an ancestor of cultivated sour cherry, an important tetraploid cherry species. Here, we present a long read chromosome-level draft genome assembly and related plastid sequences using the Oxford Nanopore Technology PromethION platform and R10.3 pore type. We generated a final consensus genome sequence of 366 Mb comprising eight chromosomes. The N50 scaffold was ~44 Mb with the longest chromosome being 66.5 Mb. The chloroplast and mitochondrial genomes were 158,217 bp and 383,281 bp long, which is in accordance with previously published plastid sequences. This is the first report of the genome of ground cherry (P. fruticosa) sequenced by long read technology only. The datasets obtained from this study provide a foundation for future breeding, molecular and evolutionary analysis in Prunus studies.

No Evidence of Unexpected Transgenic Insertions in T1190 - A Transgenic Apple Used in Rapid Cycle Breeding - Following Whole Genome Sequencing.

Rapid cycle breeding uses transgenic early flowering plants as crossbreed parents to facilitate the shortening of breeding programs for perennial crops with long-lasting juvenility. Rapid cycle breeding in apple was established using the transgenic genotype T1190 expressing the BpMADS4 gene of silver birch. In this study, the genomes of T1190 and its non-transgenic wild-type PinS (F1-offspring of 'Pinova' and 'Idared') were sequenced by Illumina short-read sequencing in two separate experiments resulting in a mean sequencing depth of 182× for T1190 and 167× for PinS. The sequencing revealed 8,450 reads, which contain sequences of ≥20 bp identical to the plant transformation vector. These reads were assembled into 125 contigs, which were examined to see whether they contained transgenic insertions or if they are not using a five-step procedure. The sequence of one contig represents the known T-DNA insertion on chromosome 4 of T1190. The sequences of the remaining contigs were either equally present in T1190 and PinS, their part with sequence identity to the vector was equally present in apple reference genomes, or they seem to result from endophytic contaminations rather than from additional transgenic insertions. Therefore, we conclude that the transgenic apple plant T1190 contains only one transgenic insertion, located on chromosome 4, and shows no further partial insertions of the transformation vector. Accession Numbers: JQ974028.1.

Transposable elements and introgression introduce genetic variation in the invasive ant Cardiocondyla obscurior.

Introduced populations of invasive organisms have to cope with novel environmental challenges, while having reduced genetic variation caused by founder effects. The mechanisms associated with this "genetic paradox of invasive species" has received considerable attention, yet few studies have examined the genomic architecture of invasive species. Populations of the heart node ant Cardiocondyla obscurior belong to two distinct lineages, a New World lineage so far only found in Latin America and a more globally distributed Old World lineage. In the present study, we use population genomic approaches to compare populations of the two lineages with apparent divergent invasive potential. We find that the strong genetic differentiation of the two lineages began at least 40,000 generations ago and that activity of transposable elements (TEs) has contributed significantly to the divergence of both lineages, possibly linked to the very unusual genomic distribution of TEs in this species. Furthermore, we show that introgression from the Old World lineage is a dominant source of genetic diversity in the New World lineage, despite the lineages' strong genetic differentiation. Our study uncovers mechanisms underlying novel genetic variation in introduced populations of C. obscurior that could contribute to the species' adaptive potential.

Chromosome-length genome assembly and structural variations of the primal Basenji dog (Canis lupus familiaris) genome.

BACKGROUND: Basenjis are considered an ancient dog breed of central African origins that still live and hunt with tribesmen in the African Congo. Nicknamed the barkless dog, Basenjis possess unique phylogeny, geographical origins and traits, making their genome structure of great interest. The increasing number of available canid reference genomes allows us to examine the impact the choice of reference genome makes with regard to reference genome quality and breed relatedness. RESULTS: Here, we report two high quality de novo Basenji genome assemblies: a female, China (CanFam_Bas), and a male, Wags. We conduct pairwise comparisons and report structural variations between assembled genomes of three dog breeds: Basenji (CanFam_Bas), Boxer (CanFam3.1) and German Shepherd Dog (GSD) (CanFam_GSD). CanFam_Bas is superior to CanFam3.1 in terms of genome contiguity and comparable overall to the high quality CanFam_GSD assembly. By aligning short read data from 58 representative dog breeds to three reference genomes, we demonstrate how the choice of reference genome significantly impacts both read mapping and variant detection. CONCLUSIONS: The growing number of high-quality canid reference genomes means the choice of reference genome is an increasingly critical decision in subsequent canid variant analyses. The basal position of the Basenji makes it suitable for variant analysis for targeted applications of specific dog breeds. However, we believe more comprehensive analyses across the entire family of canids is more suited to a pangenome approach. Collectively this work highlights the importance the choice of reference genome makes in all variation studies.

Chromosome-scale genome assembly provides insights into rye biology, evolution and agronomic potential.

Rye (Secale cereale L.) is an exceptionally climate-resilient cereal crop, used extensively to produce improved wheat varieties via introgressive hybridization and possessing the entire repertoire of genes necessary to enable hybrid breeding. Rye is allogamous and only recently domesticated, thus giving cultivated ryes access to a diverse and exploitable wild gene pool. To further enhance the agronomic potential of rye, we produced a chromosome-scale annotated assembly of the 7.9-gigabase rye genome and extensively validated its quality by using a suite of molecular genetic resources. We demonstrate applications of this resource with a broad range of investigations. We present findings on cultivated rye's incomplete genetic isolation from wild relatives, mechanisms of genome structural evolution, pathogen resistance, low-temperature tolerance, fertility control systems for hybrid breeding and the yield benefits of rye-wheat introgressions.

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

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

Transfer of the Dominant Virus Resistance Gene AV-1 pro From Asparagus prostratus to Chromosome 2 of Garden Asparagus A. officinalis L.

An introgression breeding programme was carried out to transfer the virus resistance gene AV-1 pro from the wild species Asparagus prostratus to the garden asparagus Asparagus officinalis. Serious crossing barriers caused by genetic distance and different ploidy levels of the crossing parents have been overcome using embryo rescue for the F1, BC1, and BC2 generations. The male and female fertility was widely restored in BC2 and was shown to be comparable to the cultivated asparagus. Five AV-1 resistant diploid (2n = 2x = 20) BC2 plants were selected and reciprocally backcrossed with asparagus cultivars. Segregation analyses of fourteen seedborne BC3 progenies suggested a monogenic dominant inheritance of the AV-1 resistance. Genotyping by sequencing analysis gave a strong hint for location of the resistance gene on asparagus Chromosome 2. Using an Axiom single nucleotide polymorphism (SNP) genotyping array for the analysis of three BC3 families with 10 AV-1 resistant and 10 AV-1 susceptible plants each, as well as 25 asparagus cultivars, the AV-1 pro locus on Chromosome 2 was further narrowed down. The SNP with the highest LOD score was converted to a kompetitive allele specific PCR (KASP) marker, shown to be useful for the further backcross programme and serving as the starting point for the development of a diagnostic marker.

Which Factors Affect the Occurrence of Off-Target Effects Caused by the Use of CRISPR/Cas: A Systematic Review in Plants.

CRISPR/Cas enables a targeted modification of DNA sequences. Despite their ease and efficient use, one limitation is the potential occurrence of associated off-target effects. This systematic review aims to answer the following research question: Which factors affect the occurrence of off-target effects caused by the use of CRISPR/Cas in plants? Literature published until March 2019 was considered for this review. Articles were screened for relevance based on pre-defined inclusion criteria. Relevant studies were subject to critical appraisal. All studies included in the systematic review were synthesized in a narrative report, but studies rated as high and medium/high validity were reported separately from studies rated as low and medium/low or unclear validity. In addition, we ran a binary logistic regression analysis to verify five factors that may affect the occurrence of off-target effects: (1) Number of mismatches (2) Position of mismatches (3) GC-content of the targeting sequence (4) Altered nuclease variants (5) Delivery methods. In total, 180 relevant articles were included in this review containing 468 studies therein. Seventy nine percentage of these studies were rated as having high or medium/high validity. Within these studies, 6,416 potential off-target sequences were assessed for the occurrence of off-target effects. Results clearly indicate that an increased number of mismatches between the on-target and potential off-target sequence steeply decreases the likelihood of off-target effects. The observed rate of off-target effects decreased from 59% when there is one mismatch between the on-target and off-target sequences toward 0% when four or more mismatches exist. In addition, mismatch/es located within the first eight nucleotides proximal to the PAM significantly decreased the occurrence of off-target effects. There is no evidence that the GC-content significantly affects off-target effects. The database regarding the impact of the nuclease variant and the delivery method is very poor as the majority of studies applied the standard nuclease SpCas9 and the CRISPR/Cas system was stably delivered in the genome. Hence, a general significant impact of these two factors on the occurrence of off-target effects cannot be proved. This identified evidence gap needs to be filled by systematic studies exploring these individual factors in sufficient numbers.