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.     

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.     

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.     

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.     

A CRISPR/Cas9 toolkit for efficient targeted base editing to induce genetic variations in rice

正Dear Editor,Gene-targeting technologies using sequence specific nucleases have been widely adopted to induce genetic modifications in both plant molecular biology research and crop improvement,of which generating targeted point mutation for gain-of-function phenotype is of great value.Although tremendous efforts have been made in developing geneinsertion or replacement approaches by harnessing     

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.     

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.

     

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.     

High efficiency of targeted mutagenesis in arabidopsis via meiotic promoter-driven expression of Cas9 endonuclease

Key message

The use of a meiosis I-specific promoter increased the efficiency of targeted mutagenesis and will facilitate the manipulation of homologous recombination.

Abstract

The CRISPR/Cas9 system has been harnessed for targeted engineering of eukaryotic genomes, including plants; however, CRISPR/Cas9 efficiency varies considerably in different plant tissues and species. In Arabidopsis, the generation of homozygous or bi-allelic mutants in the first (T1) generation is inefficient. Here, we used specific promoters to drive the expression of Cas9 during meiosis to maximize the efficiency of recovering heritable mutants in T1 plants. Our data reveal that the use of a promoter active in meiosis I resulted in high-efficiency (28 %) recovery of targeted mutants in the T1 generation. Moreover, this method enabled efficient simultaneous targeting of three genes for mutagenesis. Taken together, our results show that the use of meiosis-specific promoters will improve methods for functional genomic analysis and studying the molecular underpinnings of homologous recombination.

     

下一代测序中ChIP-seq数据的处理与分析

将染色质免疫共沉淀技术(ChIP)与下一代高通量测序技术相结合的染色质免疫共沉淀测序(ChIP-seq),已成为功能基因组学、特别是基因表达调控领域研究的关键技术。ChIP-seq实验带来的海量数据向生物信息学研究人员提出了新的挑战。由于此领域数据处理技术的发展大大滞后于实验技术进步,有必要系统地介绍和回顾ChIP-seq数据处理的各个方面,以便更多研究人员进入此领域设计或改进相应的算法。文章结合实例详细介绍了ChIP-seq数据整个流程,并重点讨论了其中的主要问题和关键环节,为这一研究领域的科研人员提供一个快速而深入的认识。     

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.     

Comparative genomics of first available bovine Anaplasma phagocytophilum genome obtained with targeted sequence capture

Background

Anaplasma phagocytophilum is a zoonotic and obligate intracellular bacterium transmitted by ticks. In domestic ruminants, it is the causative agent of tick-borne fever, which causes significant economic losses in Europe. As A. phagocytophilum is difficult to isolate and cultivate, only nine genome sequences have been published to date, none of which originate from a bovine strain.Our goals were to; 1/ develop a sequencing methodology which efficiently circumvents the difficulties associated with A. phagocytophilum isolation and culture; 2/ describe the first genome of a bovine strain; and 3/ compare it with available genomes, in order to both explore key genomic features at the species level, and to identify candidate genes that could be specific to bovine strains.

Results

DNA was extracted from a bovine blood sample infected by A. phagocytophilum. Following a whole genome capture approach, A. phagocytophilum DNA was enriched 197-fold in the sample and then sequenced using Illumina technology. In total, 58.9% of obtained reads corresponded to the A. phagocytophilum genome, covering 85.3% of the HZ genome. Then by performing comparisons with nine previously-sequenced A. phagocytophilum genomes, we determined the core genome of these ten strains. Following analysis, 1281 coding DNA sequences, including 1001 complete sequences, were detected in the A. phagocytophilum bovine genome, of which four appeared to be unique to the bovine isolate. These four coding DNA sequences coded for "hypothetical proteins of unknown function” and require further analysis. We also identified nine proteins common to both European domestic ruminants tested.

Conclusion

Using a whole genome capture approach, we have sequenced the first A. phagocytophilum genome isolated from a cow. To the best of our knowledge, this is the first time that this method has been used to selectively enrich pathogenic bacterial DNA from samples also containing host DNA. The four proteins unique to the A. phagocytophilum bovine genome could be involved in host tropism, therefore their functions need to be explored.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-973) contains supplementary material, which is available to authorized users.     

A genome-wide gain-of-function analysis of rice genes using the FOX-hunting system

The latest report has estimated the number of rice genes to be approximately 32,000. To elucidate the functions of a large population of rice genes and to search efficiently for agriculturally useful genes, we have been taking advantage of the Full-length cDNA Over-eXpresser (FOX) gene-hunting system. This system is very useful for analyzing various gain-of-function phenotypes from large populations of transgenic plants overexpressing cDNAs of interest and others with unknown or important functions. We collected the plasmid DNAs of 13,980 independent full-length cDNA (FL-cDNA) clones to produce a FOX library by placing individual cDNAs under the control of the maize Ubiquitin-1 promoter. The FOX library was transformed into rice by Agrobacterium-mediated high-speed transformation. So far, we have generated approximately 12,000 FOX-rice lines. Genomic PCR analysis indicated that the average number of FL-cDNAs introduced into individual lines was 1.04. Sequencing analysis of the PCR fragments carrying FL-cDNAs from 8615 FOX-rice lines identified FL-cDNAs in 8225 lines, and a database search classified the cDNAs into 5462 independent ones. Approximately 16.6% of FOX-rice lines examined showed altered growth or morphological characteristics. Three super-dwarf mutants overexpressed a novel gibberellin 2-oxidase gene,confirming the importance of this system. We also show here the other morphological alterations caused by individual FL-cDNA expression. These dominant phenotypes should be valuable indicators for gene discovery and functional analysis.     

Efficient generation of pink-fruited tomatoes using CRISPR/Cas9 system

正Tomato(Solanum lycopersicum)is the leading vegetable crop worldwide and an essential component of a healthy diet(Lin et al.,2014;Du et al.,2017).Fruit color is regarded as one of the most important commercial traits in tomato(The Tomato Genome Consortium,2012).Consumers in different regions have different color preferences.For example,European and American