Structure and specificity of the RNA-guided endonuclease Cas9 during DNA interrogation,target binding and cleavage

CRISPR-associated endonuclease Cas9 cuts DNA at variable target sites designated by a Cas9-bound RNA molecule. Cas9''s ability to be directed by single ‘guide RNA’ molecules to target nearly any sequence has been recently exploited for a number of emerging biological and medical applications. Therefore, understanding the nature of Cas9''s off-target activity is of paramount importance for its practical use. Using atomic force microscopy (AFM), we directly resolve individual Cas9 and nuclease-inactive dCas9 proteins as they bind along engineered DNA substrates. High-resolution imaging allows us to determine their relative propensities to bind with different guide RNA variants to targeted or off-target sequences. Mapping the structural properties of Cas9 and dCas9 to their respective binding sites reveals a progressive conformational transformation at DNA sites with increasing sequence similarity to its target. With kinetic Monte Carlo (KMC) simulations, these results provide evidence of a ‘conformational gating’ mechanism driven by the interactions between the guide RNA and the 14th–17th nucleotide region of the targeted DNA, the stabilities of which we find correlate significantly with reported off-target cleavage rates. KMC simulations also reveal potential methodologies to engineer guide RNA sequences with improved specificity by considering the invasion of guide RNAs into targeted DNA duplex.     

Increased efficiency of targeted mutagenesis by CRISPR/Cas9 in plants using heat stress

The CRISPR/Cas9 system has greatly improved our ability to engineer targeted mutations in eukaryotic genomes. While CRISPR/Cas9 appears to work universally, the efficiency of targeted mutagenesis and the adverse generation of off‐target mutations vary greatly between different organisms. In this study, we report that Arabidopsis plants subjected to heat stress at 37°C show much higher frequencies of CRISPR‐induced mutations compared to plants grown continuously at the standard temperature (22°C). Using quantitative assays relying on green fluorescent protein (GFP) reporter genes, we found that targeted mutagenesis by CRISPR/Cas9 in Arabidopsis is increased by approximately 5‐fold in somatic tissues and up to 100‐fold in the germline upon heat treatment. This effect of temperature on the mutation rate is not limited to Arabidopsis, as we observed a similar increase in targeted mutations by CRISPR/Cas9 in Citrus plants exposed to heat stress at 37°C. In vitro assays demonstrate that Cas9 from Streptococcus pyogenes (SpCas9) is more active in creating double‐stranded DNA breaks at 37°C than at 22°C, thus indicating a potential contributing mechanism for the in vivo effect of temperature on CRISPR/Cas9. This study reveals the importance of temperature in modulating SpCas9 activity in eukaryotes, and provides a simple method to increase on‐target mutagenesis in plants using CRISPR/Cas9.     

Optimizing copy number variation analysis using genome-wide short sequence oligonucleotide arrays

The detection of copy number variants (CNV) by array-based platforms provides valuable insight into understanding human diversity. However, suboptimal study design and data processing negatively affect CNV assessment. We quantitatively evaluate their impact when short-sequence oligonucleotide arrays are applied (Affymetrix Genome-Wide Human SNP Array 6.0) by evaluating 42 HapMap samples for CNV detection. Several processing and segmentation strategies are implemented, and results are compared to CNV assessment obtained using an oligonucleotide array CGH platform designed to query CNVs at high resolution (Agilent). We quantitatively demonstrate that different reference models (e.g. single versus pooled sample reference) used to detect CNVs are a major source of inter-platform discrepancy (up to 30%) and that CNVs residing within segmental duplication regions (higher reference copy number) are significantly harder to detect (P < 0.0001). After adjusting Affymetrix data to mimic the Agilent experimental design (reference sample effect), we applied several common segmentation approaches and evaluated differential sensitivity and specificity for CNV detection, ranging 39–77% and 86–100% for non-segmental duplication regions, respectively, and 18–55% and 39–77% for segmental duplications. Our results are relevant to any array-based CNV study and provide guidelines to optimize performance based on study-specific objectives.     

Activities and specificities of CRISPR/Cas9 and Cas12a nucleases for targeted mutagenesis in maize

CRISPR/Cas9 and Cas12a (Cpf1) nucleases are two of the most powerful genome editing tools in plants. In this work, we compared their activities by targeting maize glossy2 gene coding region that has overlapping sequences recognized by both nucleases. We introduced constructs carrying SpCas9‐guide RNA (gRNA) and LbCas12a‐CRISPR RNA (crRNA) into maize inbred B104 embryos using Agrobacterium‐mediated transformation. On‐target mutation analysis showed that 90%–100% of the Cas9‐edited T0 plants carried indel mutations and 63%–77% of them were homozygous or biallelic mutants. In contrast, 0%–60% of Cas12a‐edited T0 plants had on‐target mutations. We then conducted CIRCLE‐seq analysis to identify genome‐wide potential off‐target sites for Cas9. A total of 18 and 67 potential off‐targets were identified for the two gRNAs, respectively, with an average of five mismatches compared to the target sites. Sequencing analysis of a selected subset of the off‐target sites revealed no detectable level of mutations in the T1 plants, which constitutively express Cas9 nuclease and gRNAs. In conclusion, our results suggest that the CRISPR/Cas9 system used in this study is highly efficient and specific for genome editing in maize, while CRISPR/Cas12a needs further optimization for improved editing efficiency.     

An unbiased genome-wide analysis of zinc-finger nuclease specificity

Zinc-finger nucleases (ZFNs) allow gene editing in live cells by inducing a targeted DNA double-strand break (DSB) at a specific genomic locus. However, strategies for characterizing the genome-wide specificity of ZFNs remain limited. We show that nonhomologous end-joining captures integrase-defective lentiviral vectors at DSBs, tagging these transient events. Genome-wide integration site analysis mapped the actual in vivo cleavage activity of four ZFN pairs targeting CCR5 or IL2RG. Ranking loci with repeatedly detectable nuclease activity by deep-sequencing allowed us to monitor the degree of ZFN specificity in vivo at these positions. Cleavage required binding of ZFNs in specific spatial arrangements on DNA bearing high homology to the intended target site and only tolerated mismatches at individual positions of the ZFN binding sites. Whereas the consensus binding sequence derived in vivo closely matched that obtained in biochemical experiments, the ranking of in vivo cleavage sites could not be predicted in silico. Comprehensive mapping of ZFN activity in vivo will facilitate the broad application of these reagents in translational research.     

A targeted capture approach to generating reference sequence databases for chloroplast gene regions

Metabarcoding has improved the way we understand plants within our environment, from their ecology and conservation to invasive species management. The notion of identifying plant taxa within environmental samples relies on the ability to match unknown sequences to known reference libraries. Without comprehensive reference databases, species can go undetected or be incorrectly assigned, leading to false‐positive and false‐negative detections. To improve our ability to generate reference sequence databases, we developed a targeted capture approach using the OZBaits_CP V1.0 set, designed to capture chloroplast gene regions across the entirety of flowering plant diversity. We focused on generating a reference database for coastal temperate plant species given the lack of reference sequences for these taxa. Our approach was successful across all specimens with a target gene recovery rate of 92%, which was achieved in a single assay (i.e., samples were pooled), thus making this approach much faster and more efficient than standard barcoding. Further testing of this database highlighted 80% of all samples could be discriminated to family level across all gene regions with some genes achieving greater resolution than others—which was also dependent on the taxon of interest. Thus, we demonstrate the importance of generating reference sequences across multiple chloroplast gene regions as no single loci are sufficient to discriminate across all plant groups. The targeted capture approach outlined in this study provides a way forward to achieve this.     

Quantitative assessment of engineered Cas9 variants for target specificity enhancement by single-molecule reaction pathway analysis

There have been many engineered Cas9 variants that were developed to minimize unintended cleavage of off-target DNAs, but detailed mechanism for the way they regulate the target specificity through DNA:RNA heteroduplexation remains poorly understood. We used single-molecule FRET assay to follow the dynamics of DNA:RNA heteroduplexation for various engineered Cas9 variants with respect to on-target and off-target DNAs. Just like wild-type Cas9, these engineered Cas9 variants exhibit a strong correlation between their conformational structure and nuclease activity. Compared with wild-type Cas9, the fraction of the cleavage-competent state dropped more rapidly with increasing base-pair mismatch, which gives rise to their enhanced target specificity. We proposed a reaction model to quantitatively analyze the degree of off-target discrimination during the successive process of R-loop expansion. We found that the critical specificity enhancement step is activated during DNA:RNA heteroduplexation for evoCas9 and HypaCas9, while it occurs in the post-heteroduplexation stage for Cas9-HF1, eCas9, and Sniper-Cas9. This study sheds new light on the conformational dynamics behind the target specificity of Cas9, which will help strengthen its rational designing principles in the future.     

The genome-wide DNA sequence specificity of the anti-tumour drug bleomycin in human cells

The cancer chemotherapeutic agent, bleomycin, cleaves DNA at specific sites. For the first time, the genome-wide DNA sequence specificity of bleomycin breakage was determined in human cells. Utilising Illumina next-generation DNA sequencing techniques, over 200 million bleomycin cleavage sites were examined to elucidate the bleomycin genome-wide DNA selectivity. The genome-wide bleomycin cleavage data were analysed by four different methods to determine the cellular DNA sequence specificity of bleomycin strand breakage. For the most highly cleaved DNA sequences, the preferred site of bleomycin breakage was at 5′-GT* dinucleotide sequences (where the asterisk indicates the bleomycin cleavage site), with lesser cleavage at 5′-GC* dinucleotides. This investigation also determined longer bleomycin cleavage sequences, with preferred cleavage at 5′-GT*A and 5′- TGT* trinucleotide sequences, and 5′-TGT*A tetranucleotides. For cellular DNA, the hexanucleotide DNA sequence 5′-RTGT*AY (where R is a purine and Y is a pyrimidine) was the most highly cleaved DNA sequence. It was striking that alternating purine–pyrimidine sequences were highly cleaved by bleomycin. The highest intensity cleavage sites in cellular and purified DNA were very similar although there were some minor differences. Statistical nucleotide frequency analysis indicated a G nucleotide was present at the ?3 position (relative to the cleavage site) in cellular DNA but was absent in purified DNA.     

The genome-wide sequence specificity of DNA cleavage by bleomycin analogues in human cells

Bleomycin (BLM) is a cancer chemotherapeutic agent that cleaves cellular DNA at specific sequences. Using next-generation Illumina sequencing, the genome-wide sequence specificity of DNA cleavage by two BLM analogues, 6′-deoxy-BLM Z and zorbamycin (ZBM), was determined in human HeLa cells and compared with BLM. Over 200 million double-strand breaks were examined for each sample, and the 50,000 highest intensity cleavage sites were analysed. It was found that the DNA sequence specificity of the BLM analogues in human cells was different to BLM, especially at the cleavage site (position “0”) and the “+1” position. In human cells, the 6′-deoxy-BLM Z had a preference for 5′-GTGY*MC (where * is the cleavage site, Y is C or T, M is A or C); it was 5′-GTGY*MCA for ZBM; and 5′-GTGT*AC for BLM. With cellular DNA, the highest ranked tetranucleotides were 5′-TGC*C and 5′-TGT*A for 6′-deoxy-BLM Z; 5′-TGC*C, 5′-TGT*A and 5′-TGC*A for ZBM; and 5′-TGT*A for BLM. In purified human genomic DNA, the DNA sequence preference was 5′-TGT*A for 6′-deoxy-BLM, 5′-RTGY*AYR (where R is G or A) for ZBM, and 5′-TGT*A for BLM. Thus, the sequence specificity of the BLM analogue, 6′-deoxy-BLM Z, was similar to BLM in purified human DNA, while ZBM was different.     

Single-step generation of rabbits carrying a targeted allele of the tyrosinase gene using CRISPR/Cas9

Targeted genome editing of nonrodent mammalian species has provided the potential for highly accurate interventions into gene function in humans and the generation of useful animal models of human diseases. Here we show successful clustered regularly interspaced short palindromic repeat (CRISPR) and CRISPR-associated (Cas)-mediated gene targeting via circular plasmid injection in rabbits. The rabbit tyrosinase gene (TYR) was effectively disrupted, and we confirmed germline transmission by pronuclear injection of a circular plasmid expressing humanized Cas9 (hCas9) and single-guide RNA. Direct injection into pronuclear stage zygotes was possible following an in vitro validation assay. Neither off-target mutagenesis nor hCas9 transgenesis was detected in any of the genetically targeted pups and embryos examined. Gene targeting with this rapid and simplified strategy will help accelerate the development of translational research using other nonrodent mammalian species.     

Computational analysis of the receptor binding specificity of novel influenza A/H7N9 viruses

Background

Influenza viruses are undergoing continuous and rapid evolution. The fatal influenza A/H7N9 has drawn attention since the first wave of infections in March 2013, and raised more grave concerns with its increased potential to spread among humans. Experimental studies have revealed several host and virulence markers, indicating differential host binding preferences which can help estimate the potential of causing a pandemic. Here we systematically investigate the sequence pattern and structural characteristics of novel influenza A/H7N9 using computational approaches.

Results

The sequence analysis highlighted mutations in protein functional domains of influenza viruses. Molecular docking and molecular dynamics simulation revealed that the hemagglutinin (HA) of A/Taiwan/1/2017(H7N9) strain enhanced the binding with both avian and human receptor analogs, compared with the previous A/Shanghai/02/2013(H7N9) strain. The Molecular Mechanics - Poisson Boltzmann Surface Area (MM-PBSA) calculation revealed the change of residue-ligand interaction energy and detected the residues with conspicuous binding preference.

Conclusion

The results are novel and specific to the emerging influenza A/Taiwan/1/2017(H7N9) strain compared with A/Shanghai/02/2013(H7N9). Its enhanced ability to bind human receptor analogs, which are abundant in the human upper respiratory tract, may be responsible for the recent outbreak. Residues showing binding preference were detected, which could facilitate monitoring the circulating influenza viruses.

     

A theoretical study of the sequence specificity in binding of lexitropsins to B-DNA

A theoretical study is presented on the binding to B-DNA of a series of lexitropsins, these ligands being netropsin derivatives in which one or both of the pyrrole rings have been replaced by imidazoles. The best complexes have been located by energy minimisation taking into account nucleic acid flexibility, ligand flexibility, explicit, mobile counterions and solvent dielectric effects. Calculations have been performed for two homopolymeric DNA receptor sequences, AT base sequence, which only decreases in the imidazole derivatives. These results emphasize the decisive role of the molecular electrostatic potential of the nucleic acid in determining the sequence selectivity of these ligands, as opposed to the postulated role of adenine C2 - pyrrole beta hydrogen contacts.