Phosphonoacetate Modifications Enhance the Stability and Editing Yields of Guide RNAs for Cas9 Editors.

CRISPR gene editing and control systems continue to emerge and inspire novel research and clinical applications. Advances in CRISPR performance such as optimizing the duration of activity in cells, tissues, and organisms, as well as limiting off-target activities, have been extremely important for expanding the utility of CRISPR-based systems. By investigating the effects of various chemical modifications in guide RNAs (gRNAs) at defined positions and combinations, we find that 2'-O-methyl-3'-phosphonoacetate (MP) modifications can be substantially more effective than 2'-O-methyl-3'-phosphorothioate (MS) modifications at the 3' ends of single-guide RNAs (sgRNAs) to promote high editing yields, in some instances showing an order of magnitude higher editing yield in human cells. MP-modified 3' ends are especially effective at promoting the activity of guide RNAs cotransfected with Cas messenger RNA (mRNA), as the gRNA must persist in cells until the Cas protein is expressed. We demonstrate such an MP enhancement for sgRNAs cotransfected with a BE4 mRNA for cytidine base editing and also demonstrate that MP at the 3' ends of prime editing guide RNAs (pegRNAs) cotransfected with PE2 mRNA can promote maximal prime editing yields. In the presence of serum, sgRNAs with MP-modified 3' ends showed marked improvements in editing efficiency over sgRNAs with MS-modified 3' ends codelivered with Cas9 mRNA and showed more modest improvements at enhancing the activity of transfected ribonucleoprotein (RNP) complexes. Our results suggest that MP should be considered as a performance-enhancing modification for the 3' ends of synthetic gRNAs, especially in situations where the guide RNAs may be susceptible to exonuclease-mediated degradation.

Global Transcriptional Response to CRISPR/Cas9-AAV6-Based Genome Editing in CD34+ Hematopoietic Stem and Progenitor Cells.

Genome-editing technologies are currently being translated to the clinic. However, cellular effects of the editing machinery have yet to be fully elucidated. Here, we performed global microarray-based gene expression measurements on human CD34+ hematopoietic stem and progenitor cells that underwent editing. We probed effects of the entire editing process as well as each component individually, including electroporation, Cas9 (mRNA or protein) with chemically modified sgRNA, and AAV6 transduction. We identified differentially expressed genes relative to control treatments, which displayed enrichment for particular biological processes. All editing machinery components elicited immune, stress, and apoptotic responses. Cas9 mRNA invoked the greatest amount of transcriptional change, eliciting a distinct viral response and global transcriptional downregulation, particularly of metabolic and cell cycle processes. Electroporation also induced significant transcriptional change, with notable downregulation of metabolic processes. Surprisingly, AAV6 evoked no detectable viral response. We also found Cas9/sgRNA ribonucleoprotein treatment to be well tolerated, in spite of eliciting a DNA damage signature. Overall, this data establishes a benchmark for cellular tolerance of CRISPR/Cas9-AAV6-based genome editing, ensuring that the clinical protocol is as safe and efficient as possible.

Improving CRISPR-Cas specificity with chemical modifications in single-guide RNAs.

CRISPR systems have emerged as transformative tools for altering genomes in living cells with unprecedented ease, inspiring keen interest in increasing their specificity for perfectly matched targets. We have developed a novel approach for improving specificity by incorporating chemical modifications in guide RNAs (gRNAs) at specific sites in their DNA recognition sequence ('guide sequence') and systematically evaluating their on-target and off-target activities in biochemical DNA cleavage assays and cell-based assays. Our results show that a chemical modification (2'-O-methyl-3'-phosphonoacetate, or 'MP') incorporated at select sites in the ribose-phosphate backbone of gRNAs can dramatically reduce off-target cleavage activities while maintaining high on-target performance, as demonstrated in clinically relevant genes. These findings reveal a unique method for enhancing specificity by chemically modifying the guide sequence in gRNAs. Our approach introduces a versatile tool for augmenting the performance of CRISPR systems for research, industrial and therapeutic applications.

Chemically modified guide RNAs enhance CRISPR-Cas genome editing in human primary cells.

CRISPR-Cas-mediated genome editing relies on guide RNAs that direct site-specific DNA cleavage facilitated by the Cas endonuclease. Here we report that chemical alterations to synthesized single guide RNAs (sgRNAs) enhance genome editing efficiency in human primary T cells and CD34(+) hematopoietic stem and progenitor cells. Co-delivering chemically modified sgRNAs with Cas9 mRNA or protein is an efficient RNA- or ribonucleoprotein (RNP)-based delivery method for the CRISPR-Cas system, without the toxicity associated with DNA delivery. This approach is a simple and effective way to streamline the development of genome editing with the potential to accelerate a wide array of biotechnological and therapeutic applications of the CRISPR-Cas technology.

RNA-protein analysis using a conditional CRISPR nuclease.

RNA-binding proteins control the fate and function of the transcriptome in all cells. Here we present technology for isolating RNA-protein partners efficiently and accurately using an engineered clustered regularly interspaced short palindromic repeats (CRISPR) endoribonuclease. An inactive version of the Csy4 nuclease binds irreversibly to transcripts engineered with a 16-nt hairpin sequence at their 5' ends. Once immobilized by Csy4 on a solid support, contaminating proteins and other molecules can be removed by extensive washing. Upon addition of imidazole, Csy4 is activated to cleave the RNA, removing the hairpin tag and releasing the native transcript along with its specifically bound protein partners. This conditional Csy4 enzyme enables recovery of specific RNA-binding partners with minimal false-positive contamination. We use this method, coupled with quantitative MS, to identify cell type-specific human pre-microRNA-binding proteins. We also show that this technology is suitable for analyzing diverse size transcripts, and that it is suitable for adaptation to a high-throughput discovery format.

lincRNAs act in the circuitry controlling pluripotency and differentiation.

Although thousands of large intergenic non-coding RNAs (lincRNAs) have been identified in mammals, few have been functionally characterized, leading to debate about their biological role. To address this, we performed loss-of-function studies on most lincRNAs expressed in mouse embryonic stem (ES) cells and characterized the effects on gene expression. Here we show that knockdown of lincRNAs has major consequences on gene expression patterns, comparable to knockdown of well-known ES cell regulators. Notably, lincRNAs primarily affect gene expression in trans. Knockdown of dozens of lincRNAs causes either exit from the pluripotent state or upregulation of lineage commitment programs. We integrate lincRNAs into the molecular circuitry of ES cells and show that lincRNA genes are regulated by key transcription factors and that lincRNA transcripts bind to multiple chromatin regulatory proteins to affect shared gene expression programs. Together, the results demonstrate that lincRNAs have key roles in the circuitry controlling ES cell state.

Population-genetic properties of differentiated human copy-number polymorphisms.

Copy-number variants (CNVs) can reach appreciable frequencies in the human population, and recent discoveries have shown that several of these copy-number polymorphisms (CNPs) are associated with human diseases, including lupus, psoriasis, Crohn disease, and obesity. Despite new advances, significant biases remain in terms of CNP discovery and genotyping. We developed a method based on single-channel intensity data and benchmarked against copy numbers determined from sequencing read depth to successfully obtain CNP genotypes for 1495 CNPs from 487 human DNA samples of diverse ethnic backgrounds. This microarray contained CNPs in segmental duplication-rich regions and insertions of sequences not represented in the reference genome assembly or on standard SNP microarray platforms. We observe that CNPs in segmental duplications are more likely to be population differentiated than CNPs in unique regions (p = 0.015) and that biallelic CNPs show greater stratification when compared to frequency-matched SNPs (p = 0.0026). Although biallelic CNPs show a strong correlation of copy number with flanking SNP genotypes, the majority of multicopy CNPs do not (40% with r > 0.8). We selected a subset of CNPs for further characterization in 1876 additional samples from 62 populations; this revealed striking population-differentiated structural variants in genes of clinical significance such as OCLN, a tight junction protein involved in hepatitis C viral entry. Our microarray design allows these variants to be rapidly tested for disease association and our results suggest that CNPs (especially those that cannot be imputed from SNP genotypes) might have contributed disproportionately to human diversity and selection.

Diversity of human copy number variation and multicopy genes.

Copy number variants affect both disease and normal phenotypic variation, but those lying within heavily duplicated, highly identical sequence have been difficult to assay. By analyzing short-read mapping depth for 159 human genomes, we demonstrated accurate estimation of absolute copy number for duplications as small as 1.9 kilobase pairs, ranging from 0 to 48 copies. We identified 4.1 million "singly unique nucleotide" positions informative in distinguishing specific copies and used them to genotype the copy and content of specific paralogs within highly duplicated gene families. These data identify human-specific expansions in genes associated with brain development, reveal extensive population genetic diversity, and detect signatures consistent with gene conversion in the human species. Our approach makes ~1000 genes accessible to genetic studies of disease association.

Characterization of missing human genome sequences and copy-number polymorphic insertions.

The extent of human genomic structural variation suggests that there must be portions of the genome yet to be discovered, annotated and characterized at the sequence level. We present a resource and analysis of 2,363 new insertion sequences corresponding to 720 genomic loci. We found that a substantial fraction of these sequences are either missing, fragmented or misassigned when compared to recent de novo sequence assemblies from short-read next-generation sequence data. We determined that 18-37% of these new insertions are copy-number polymorphic, including loci that show extensive population stratification among Europeans, Asians and Africans. Complete sequencing of 156 of these insertions identified new exons and conserved noncoding sequences not yet represented in the reference genome. We developed a method to accurately genotype these new insertions by mapping next-generation sequencing datasets to the breakpoint, thereby providing a means to characterize copy-number status for regions previously inaccessible to single-nucleotide polymorphism microarrays.

Mapping and sequencing of structural variation from eight human genomes.

Genetic variation among individual humans occurs on many different scales, ranging from gross alterations in the human karyotype to single nucleotide changes. Here we explore variation on an intermediate scale--particularly insertions, deletions and inversions affecting from a few thousand to a few million base pairs. We employed a clone-based method to interrogate this intermediate structural variation in eight individuals of diverse geographic ancestry. Our analysis provides a comprehensive overview of the normal pattern of structural variation present in these genomes, refining the location of 1,695 structural variants. We find that 50% were seen in more than one individual and that nearly half lay outside regions of the genome previously described as structurally variant. We discover 525 new insertion sequences that are not present in the human reference genome and show that many of these are variable in copy number between individuals. Complete sequencing of 261 structural variants reveals considerable locus complexity and provides insights into the different mutational processes that have shaped the human genome. These data provide the first high-resolution sequence map of human structural variation--a standard for genotyping platforms and a prelude to future individual genome sequencing projects.

The fine-scale and complex architecture of human copy-number variation.

Despite considerable excitement over the potential functional significance of copy-number variants (CNVs), we still lack knowledge of the fine-scale architecture of the large majority of CNV regions in the human genome. In this study, we used a high-resolution array-based comparative genomic hybridization (aCGH) platform that targeted known CNV regions of the human genome at approximately 1 kb resolution to interrogate the genomic DNAs of 30 individuals from four HapMap populations. Our results revealed that 1020 of 1153 CNV loci (88%) were actually smaller in size than what is recorded in the Database of Genomic Variants based on previously published studies. A reduction in size of more than 50% was observed for 876 CNV regions (76%). We conclude that the total genomic content of currently known common human CNVs is likely smaller than previously thought. In addition, approximately 8% of the CNV regions observed in multiple individuals exhibited genomic architectural complexity in the form of smaller CNVs within larger ones and CNVs with interindividual variation in breakpoints. Future association studies that aim to capture the potential influences of CNVs on disease phenotypes will need to consider how to best ascertain this previously uncharacterized complexity.

Promiscuous mutations activate the noncanonical NF-kappaB pathway in multiple myeloma.

Activation of NF-kappaB has been noted in many tumor types, however only rarely has this been linked to an underlying genetic mutation. An integrated analysis of high-density oligonucleotide array CGH and gene expression profiling data from 155 multiple myeloma samples identified a promiscuous array of abnormalities contributing to the dysregulation of NF-kappaB in approximately 20% of patients. We report mutations in ten genes causing the inactivation of TRAF2, TRAF3, CYLD, cIAP1/cIAP2 and activation of NFKB1, NFKB2, CD40, LTBR, TACI, and NIK that result primarily in constitutive activation of the noncanonical NF-kappaB pathway, with the single most common abnormality being inactivation of TRAF3. These results highlight the critical importance of the NF-kappaB pathway in the pathogenesis of multiple myeloma.

Array CGH analysis of copy number variation identifies 1284 new genes variant in healthy white males: implications for association studies of complex diseases.

The discovery of copy number variation in healthy individuals is far from complete, and owing to the resolution of detection systems used, the majority of loci reported so far are relatively large ( approximately 65%>10 kb). Applying a two-stage high-resolution array comparative genomic hybridization approach to analyse 50 healthy Caucasian males from northern France, we discovered 2208 copy number variants (CNVs) detected by more than one consecutive probe. These clustered into 1469 CNV regions (CNVRs), of which 721 are thought to be novel. The majority of these are small (median size 4.4 kb) and most have common boundaries, with a coefficient of variation less than 0.1 for 83% of endpoints in those observed in multiple samples. Only 6% of the CNVRs analysed showed evidence of both copy number losses and gains at the same site. A further 6089 variants were detected by single probes: 48% of these were observed in more than one individual. In total, 2570 genes were seen to intersect variants: 1284 in novel loci. Genes involved in differentiation and development were significantly over-represented and approximately half of the genes identified feature in the Online Mendelian Inheritance in Man database. The biological importance of many genes affected, along with the well-conserved nature of the majority of the CNVs, suggests that they could have important implications for phenotype and, thus, be useful for association studies of complex diseases.

High resolution oligonucleotide CGH using DNA from archived prostate tissue.

BACKGROUND: The current focus on biomarker discovery is a result of an improved understanding of the biological basis for carcinogenesis and advances in technology. Biomarkers can aid in diagnosis, prognosis, treatment selection, and drug development. There is an urgent need for high-resolution tools that perform well using archived tissue for biomarker discovery and tools that can translate into the clinic. METHODS: Oligonucleotide array comparative genomic hybridization (oCGH) was compared to BAC-based aCGH using unamplified total genomic DNA from formalin fixed paraffin-embedded (FFPE) prostate tissue. RESULTS: The copy number aberrations detected with the BAC and oligonucleotide arrays were highly correlated in cases where the arrays contained probes in similar genomic locations. The oligonucleotide array platform provided more precise mapping due to the higher density of oligonucleotide probes. CONCLUSIONS: These results demonstrate the utility of high-resolution oligonucleotide arrays designed to use genomic DNA for CGH measurements using archived tissue samples for discovery and clinic based assays.

Robust interlaboratory reproducibility of a gene expression signature measurement consistent with the needs of a new generation of diagnostic tools.

BACKGROUND: The increasing use of DNA microarrays in biomedical research, toxicogenomics, pharmaceutical development, and diagnostics has focused attention on the reproducibility and reliability of microarray measurements. While the reproducibility of microarray gene expression measurements has been the subject of several recent reports, there is still a need for systematic investigation into what factors most contribute to variability of measured expression levels observed among different laboratories and different experimenters. RESULTS: We report the results of an interlaboratory comparison of gene expression array measurements on the same microarray platform, in which the RNA amplification and labeling, hybridization and wash, and slide scanning were each individually varied. Identical input RNA was used for all experiments. While some sources of variation have measurable influence on individual microarray signals, they showed very low influence on sample-to-reference ratios based on averaged triplicate measurements in the two-color experiments. RNA labeling was the largest contributor to interlaboratory variation. CONCLUSION: Despite this variation, measurement of one particular breast cancer gene expression signature in three different laboratories was found to be highly robust, showing a high intralaboratory and interlaboratory reproducibility when using strictly controlled standard operating procedures.

Molecular signatures determining coronary artery and saphenous vein smooth muscle cell phenotypes: distinct responses to stimuli.

OBJECTIVE: Phenotypic differences between vascular smooth muscle cell (VSMC) subtypes lead to diverse pathological processes including atherosclerosis, postangioplasty restenosis and vein graft disease. To better understand the molecular mechanisms underlying functional differences among distinct SMC subtypes, we compared gene expression profiles and functional responses to oxidized low-density lipoprotein (OxLDL) and platelet-derived growth factor (PDGF) between cultured SMCs from human coronary artery (CASM) and saphenous vein (SVSM). METHODS AND RESULTS: OxLDL and PDGF elicited markedly different functional responses and expression profiles between the 2 SMC subtypes. In CASM, OxLDL inhibited cell proliferation and migration and modified gene expression of chemokines (CXCL10, CXCL11 and CXCL12), proinflammatory cytokines (IL-1, IL-6, and IL-18), insulin-like growth factor binding proteins (IGFBPs), and both endothelial and smooth muscle marker genes. In SVSM, OxLDL promoted proliferation partially via IGF1 signaling, activated NF-kappaB and phosphatidylinositol signaling pathways, and upregulated prostaglandin (PG) receptors and synthases. In untreated cells, alpha-chemokines, proinflammatory cytokines, and genes associated with apoptosis, inflammation, and lipid biosynthesis were higher in CASM, whereas some beta-chemokines, metalloproteinase inhibitors, and IGFBPs were higher in SVSM. Interestingly, the basal expression levels of these genes seemed closely related to their responses to OxLDL and PDGF. In summary, our results suggest dramatic differences in gene expression patterns and functional responses to OxLDL and PDGF between venous and arterial SMCs, with venous SMCs having stronger proliferative/migratory responses to stimuli but also higher expression of atheroprotective genes at baseline. CONCLUSIONS: These results reveal molecular signatures that define the distinct phenotypes characteristics of coronary artery and saphenous vein SMC subtypes.

Differences in vascular bed disease susceptibility reflect differences in gene expression response to atherogenic stimuli.

Atherosclerosis occurs predominantly in arteries and only rarely in veins. The goal of this study was to test whether differences in the molecular responses of venous and arterial endothelial cells (ECs) to atherosclerotic stimuli might contribute to vascular bed differences in susceptibility to atherosclerosis. We compared gene expression profiles of primary cultured ECs from human saphenous vein (SVEC) and coronary artery (CAEC) exposed to atherogenic stimuli. In addition to identifying differentially expressed genes, we applied statistical analysis of gene ontology and pathway annotation terms to identify signaling differences related to cell type and stimulus. Differential gene expression of untreated venous and arterial endothelial cells yielded 285 genes more highly expressed in untreated SVEC (P<0.005 and fold change >1.5). These genes represented various atherosclerosis-related pathways including responses to proliferation, oxidoreductase activity, antiinflammatory responses, cell growth, and hemostasis functions. Moreover, stimulation with oxidized LDL induced dramatically greater gene expression responses in CAEC compared with SVEC, relating to adhesion, proliferation, and apoptosis pathways. In contrast, interleukin 1beta and tumor necrosis factor alpha activated similar gene expression responses in both CAEC and SVEC. The differences in functional response and gene expression were further validated by an in vitro proliferation assay and in vivo immunostaining of alphabeta-crystallin protein. Our results strongly suggest that different inherent gene expression programs in arterial versus venous endothelial cells contribute to differences in atherosclerotic disease susceptibility.

Pathway analysis of coronary atherosclerosis.

Large-scale gene expression studies provide significant insight into genes differentially regulated in disease processes such as cancer. However, these investigations offer limited understanding of multisystem, multicellular diseases such as atherosclerosis. A systems biology approach that accounts for gene interactions, incorporates nontranscriptionally regulated genes, and integrates prior knowledge offers many advantages. We performed a comprehensive gene level assessment of coronary atherosclerosis using 51 coronary artery segments isolated from the explanted hearts of 22 cardiac transplant patients. After histological grading of vascular segments according to American Heart Association guidelines, isolated RNA was hybridized onto a customized 22-K oligonucleotide microarray, and significance analysis of microarrays and gene ontology analyses were performed to identify significant gene expression profiles. Our studies revealed that loss of differentiated smooth muscle cell gene expression is the primary expression signature of disease progression in atherosclerosis. Furthermore, we provide insight into the severe form of coronary artery disease associated with diabetes, reporting an overabundance of immune and inflammatory signals in diabetics. We present a novel approach to pathway development based on connectivity, determined by language parsing of the published literature, and ranking, determined by the significance of differentially regulated genes in the network. In doing this, we identify highly connected "nexus" genes that are attractive candidates for therapeutic targeting and followup studies. Our use of pathway techniques to study atherosclerosis as an integrated network of gene interactions expands on traditional microarray analysis methods and emphasizes the significant advantages of a systems-based approach to analyzing complex disease.

Comparative genomic hybridization using oligonucleotide microarrays and total genomic DNA.

Array-based comparative genomic hybridization (CGH) measures copy-number variations at multiple loci simultaneously, providing an important tool for studying cancer and developmental disorders and for developing diagnostic and therapeutic targets. Arrays for CGH based on PCR products representing assemblies of BAC or cDNA clones typically require maintenance, propagation, replication, and verification of large clone sets. Furthermore, it is difficult to control the specificity of the hybridization to the complex sequences that are present in each feature of such arrays. To develop a more robust and flexible platform, we created probe-design methods and assay protocols that make oligonucleotide microarrays synthesized in situ by inkjet technology compatible with array-based comparative genomic hybridization applications employing samples of total genomic DNA. Hybridization of a series of cell lines with variable numbers of X chromosomes to arrays designed for CGH measurements gave median ratios for X-chromosome probes within 6% of the theoretical values (0.5 for XY/XX, 1.0 for XX/XX, 1.4 for XXX/XX, 2.1 for XXXX/XX, and 2.6 for XXXXX/XX). Furthermore, these arrays detected and mapped regions of single-copy losses, homozygous deletions, and amplicons of various sizes in different model systems, including diploid cells with a chromosomal breakpoint that has been mapped and sequenced to a precise nucleotide and tumor cell lines with highly variable regions of gains and losses. Our results demonstrate that oligonucleotide arrays designed for CGH provide a robust and precise platform for detecting chromosomal alterations throughout a genome with high sensitivity even when using full-complexity genomic samples.

Novel role for the potent endogenous inotrope apelin in human cardiac dysfunction.

BACKGROUND: Apelin is among the most potent stimulators of cardiac contractility known. However, no physiological or pathological role for apelin-angiotensin receptor-like 1 (APJ) signaling has ever been described. METHODS AND RESULTS: We performed transcriptional profiling using a spotted cDNA microarray with 12 814 unique clones on paired samples of left ventricle obtained before and after placement of a left ventricular assist device in 11 patients. The significance analysis of microarrays and a novel rank consistency score designed to exploit the paired structure of the data confirmed that natriuretic peptides were among the most significantly downregulated genes after offloading. The most significantly upregulated gene was the G-protein-coupled receptor APJ, the specific receptor for apelin. We demonstrate here using immunoassay and immunohistochemical techniques that apelin is localized primarily in the endothelium of the coronary arteries and is found at a higher concentration in cardiac tissue after mechanical offloading. These findings imply an important paracrine signaling pathway in the heart. We additionally extend the clinical significance of this work by reporting for the first time circulating human apelin levels and demonstrating increases in the plasma level of apelin in patients with left ventricular dysfunction. CONCLUSIONS: The apelin-APJ signaling pathway emerges as an important novel mediator of cardiovascular control.