Allele-specific genome targeting in the development of precision medicine

Allele-specific CRISPR/Cas9 editing inactivates a single nucleotide variant associated with collagen VI muscular dystrophy

The application of allele-specific gene editing tools can expand the therapeutic options for dominant genetic conditions, either via gene correction or via allelic gene inactivation in situations where haploinsufficiency is tolerated. Here, we used allele-targeted CRISPR/Cas9 guide RNAs (gRNAs) to introduce inactivating frameshifting indels at a single nucleotide variant in the COL6A1 gene (c.868G>A; G290R), a variant that acts as dominant negative and that is associated with a severe form of congenital muscular dystrophy. We expressed spCas9 along with allele-targeted gRNAs, without providing a repair template, in primary fibroblasts derived from four patients and one control subject. Amplicon deep-sequencing for two gRNAs tested showed that single nucleotide deletions accounted for the majority of indels introduced. While activity of the two gRNAs was greater at the G290R allele, both gRNAs were also active at the wild-type allele. To enhance allele-selectivity, we introduced deliberate additional mismatches to one gRNA. One of these optimized gRNAs showed minimal activity at the WT allele, while generating productive edits and improving collagen VI matrix in cultured patient fibroblasts. This study strengthens the potential of gene editing to treat dominant-negative disorders, but also underscores the challenges in achieving allele selectivity with gRNAs.

Development of super-specific epigenome editing by targeted allele-specific DNA methylation

BACKGROUND

Epigenome editing refers to the targeted reprogramming of genomic loci using an EpiEditor which may consist of an sgRNA/dCas9 complex that recruits DNMT3A/3L to the target locus. Methylation of the locus can lead to a modulation of gene expression. Allele-specific DNA methylation (ASM) refers to the targeted methylation delivery only to one allele of a locus. In the context of diseases caused by a dominant mutation, the selective DNA methylation of the mutant allele could be used to repress its expression but retain the functionality of the normal gene.

RESULTS

To set up allele-specific targeted DNA methylation, target regions were selected from hypomethylated CGIs bearing a heterozygous SNP in their promoters in the HEK293 cell line. We aimed at delivering maximum DNA methylation with highest allelic specificity in the targeted regions. Placing SNPs in the PAM or seed regions of the sgRNA, we designed 24 different sgRNAs targeting single alleles in 14 different gene loci. We achieved efficient ASM in multiple cases, such as ISG15, MSH6, GPD1L, MRPL52, PDE8A, NARF, DAP3, and GSPT1, which in best cases led to five to tenfold stronger average DNA methylation at the on-target allele and absolute differences in the DNA methylation gain at on- and off-target alleles of > 50%. In general, loci with the allele discriminatory SNP positioned in the PAM region showed higher success rate of ASM and better specificity. Highest DNA methylation was observed on day 3 after transfection followed by a gradual decline. In selected cases, ASM was stable up to 11 days in HEK293 cells and it led up to a 3.6-fold change in allelic expression ratios.

CONCLUSIONS

We successfully delivered ASM at multiple genomic loci with high specificity, efficiency and stability. This form of super-specific epigenome editing could find applications in the treatment of diseases caused by dominant mutations, because it allows silencing of the mutant allele without repression of the expression of the normal allele thereby minimizing potential side-effects of the treatment.

Full-Length Transcript Phasing with Third-Generation Sequencing

Haplotyping individual full-length transcripts can be important in diagnosis and treatment of certain genetic diseases. One set of diseases, repeat expansions of simple tandem repeat sequences are the cause of over 40 neurological disorders. In many of these conditions, expanding a polymorphic repeat beyond a given threshold has been strongly associated with disease onset and severity. Given that most repeat expansions are inherited in an autosomal dominant pattern, repeat expansion disorders are typically characterized by a heterozygous expansion locus associated with a single haplotype. Precision genetic medicines can be used to selectively target expansion-containing sequences in a haplotype-specific manner.However, repeat expansion lengths often exceed the capacity of next-generation sequencing (NGS) reads. Therefore, the accurate length and haplotype determination of repeat expansions requires special considerations and requires the development of custom methods. Here we highlight a method for targeted haplotype phasing of the HTT gene, which can be adopted for use with other full-length transcripts and in other repeat expansion disorders.

CRISPR-Cas9-mediated gene therapy in lung cancer

As the largest cause of cancer-related deaths worldwide, pulmonary cancer is the most common form of the disease. Several genetic, epigenetic, and environmental factors come into play during the multi-step mechanism of tumorigenesis. The heterogeneity that makes discovering successful therapeutics for pulmonary cancer problematic is significantly influenced by the epigenetic landscape, including DNA methylation, chromatin architecture, histone modifications, and noncoding RNA control. Clinical activity of epigenetic-targeted medicines has been reported in hematological tumors, and these compounds may also have therapeutic effects in solid tumors. Over the course of the past few years, some researchers have successfully modified the expression of genes in cells using the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas (CRISPR-associated proteins) technique. The utilization of this technology allows for the induction of site-specific mutagenesis, epigenetic alterations, and the regulation of gene expression. This study will present an overview of the primary epigenetic alterations seen in pulmonary cancer, as well as a summary of therapeutic implications for targeting epigenetics in the management of pulmonary cancer, with a particular emphasis on the technique known as CRISPR/Cas9.

Identification of a Novel Functional Non-synonymous Single Nucleotide Polymorphism in Frizzled Class Receptor 6 Gene for Involvement in Depressive Symptoms

Although numerous susceptibility loci for depression have been identified in recent years, their biological function and molecular mechanism remain largely unknown. By using an exome-wide association study for depressive symptoms assessed by the Center for Epidemiological Studies Depression (CES-D) score, we discovered a novel missense single nucleotide polymorphism (SNP), rs61753730 (Q152E), located in the fourth exon of the frizzled class receptor 6 gene (FZD6), which is a potential causal variant and is significantly associated with the CES-D score. Computer-based in silico analysis revealed that the protein configuration and stability, as well as the secondary structure of FZD6 differed greatly between the wild-type (WT) and Q152E mutant. We further found that rs61753730 significantly affected the luciferase activity and expression of FZD6 in an allele-specific way. Finally, we generated Fzd6-knockin (Fzd6-KI) mice with rs61753730 mutation using the CRISPR/Cas9 genome editing system and found that these mice presented greater immobility in the forced swimming test, less preference for sucrose in the sucrose preference test, as well as decreased center entries, center time, and distance traveled in the open filed test compared with WT mice after exposed to chronic social defeat stress. These results indicate the involvement of rs61753730 in depression. Taken together, our findings demonstrate that SNP rs61753730 is a novel functional variant and plays an important role in depressive symptoms.

Identification and characterization of two novel noncoding tyrosinase (TYR) gene variants leading to oculocutaneous albinism type 1

Oculocutaneous albinism type 1 (OCA1), resulting from pathogenic variants in the tyrosinase (TYR) gene, refers to a group of phenotypically heterogeneous autosomal recessive disorders characterized by a partial or a complete absence of pigment in the skin/hair and is also associated with common developmental eye defects. In this study, we identified two novel compound heterozygous TYR variants from a Chinese hypopigmentary patient by whole-exome sequencing. Specifically, the two variants were c.-89T>G, located at the core of the initiator E-box (Inr E-box) of the TYR promoter, and p.S16Y (c.47C>A), located within the signal sequence. We performed both in silico analysis and experimental validation and verified these mutations as OCA1 variants that caused either impaired or complete loss of function of TYR. Mechanistically, the Inr E-box variant dampened TYR binding to microphthalmia-associated transcription factor, a master transcriptional regulator of the melanocyte development, whereas the S16Y variant contributed to endoplasmic reticulum retention, a common and principal cause of impaired TYR activity. Interestingly, we found that the Inr E-box variant creates novel protospacer adjacent motif sites, recognized by nucleases SpCas9 and SaCas9-KKH, respectively, without compromising the functional TYR coding sequence. We further used allele-specific genomic editing by CRISPR activation to specifically target the variant promoter and successfully activated its downstream gene expression, which could lead to potential therapeutic benefits. In conclusion, this study expands the spectrum of TYR variants, especially those within the promoter and noncoding regions, which can facilitate genetic counseling and clinical diagnosis of OCA1.

Allele-aware chromosome-scale assembly of the allopolyploid genome of hexaploid Ma bamboo (Dendrocalamus latiflorus Munro)

Dendrocalamus latiflorus Munro is a woody clumping bamboo with rapid shoot growth. Both genetic transformation and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) gene editing techniques are available for D. latiflorus, enabling reverse genetic approaches. Thus, D. latiflorus has the potential to be a model bamboo species. However, the genome sequence of D. latiflorus has remained unreported due to its polyploidy and large genome size. Here, we sequenced the D. latiflorus genome and assembled it into three allele-aware subgenomes (AABBCC), representing the largest genome of a major bamboo species. We assembled 70 allelic chromosomes (2, 737 Mb) for hexaploid D. latiflorus using both single-molecule sequencing from the Pacific Biosciences (PacBio) Sequel platform and chromosome conformation capture sequencing (Hi-C). Repetitive sequences comprised 52.65% of the D. latiflorus genome. We annotated 135 231 protein-coding genes in the genome based on transcriptomes from eight different tissues. Transcriptome sequencing using RNA-Seq and PacBio single-molecule real-time long-read isoform sequencing revealed highly differential alternative splicing (AS) between non-abortive and abortive shoots, suggesting that AS regulates the abortion rate of bamboo shoots. This high-quality hexaploid genome and comprehensive strand-specific transcriptome datasets for this Poaceae family member will pave the way for bamboo research using D. latiflorus as a model species.

Systematic decomposition of sequence determinants governing CRISPR/Cas9 specificity

The specificity of CRISPR/Cas9 genome editing is largely determined by the sequences of guide RNA (gRNA) and the targeted DNA, yet the sequence-dependent rules underlying off-target effects are not fully understood. To systematically explore the sequence determinants governing CRISPR/Cas9 specificity, here we describe a dual-target system to measure the relative cleavage rate between off- and on-target sequences (off-on ratios) of 1902 gRNAs on 13,314 synthetic target sequences, and reveal a set of sequence rules involving 2 factors in off-targeting: 1) a guide-intrinsic mismatch tolerance (GMT) independent of the mismatch context; 2) an "epistasis-like" combinatorial effect of multiple mismatches, which are associated with the free-energy landscape in R-loop formation and are explainable by a multi-state kinetic model. These sequence rules lead to the development of MOFF, a model-based predictor of Cas9-mediated off-target effects. Moreover, the "epistasis-like" combinatorial effect suggests a strategy of allele-specific genome editing using mismatched guides. With the aid of MOFF prediction, this strategy significantly improves the selectivity and expands the application domain of Cas9-based allele-specific editing, as tested in a high-throughput allele-editing screen on 18 cancer hotspot mutations.

The dopamine transporter gene SLC6A3: multidisease risks

The human dopamine transporter gene SLC6A3 has been consistently implicated in several neuropsychiatric diseases but the disease mechanism remains elusive. In this risk synthesis, we have concluded that SLC6A3 represents an increasingly recognized risk with a growing number of familial mutants associated with neuropsychiatric and neurological disorders. At least five loci were related to common and severe diseases including alcohol use disorder (high activity variant), attention-deficit/hyperactivity disorder (low activity variant), autism (familial proteins with mutated networking) and movement disorders (both regulatory variants and familial mutations). Association signals depended on genetic markers used as well as ethnicity examined. Strong haplotype selection and gene-wide epistases support multimarker assessment of functional variations and phenotype associations. Inclusion of its promoter region's functional markers such as DNPi (rs67175440) and 5'VNTR (rs70957367) may help delineate condensate-based risk action, testing a locus-pathway-phenotype hypothesis for one gene-multidisease etiology.

Johne's Disease in Dairy Cattle: An Immunogenetic Perspective

Johne's disease (JD), also known as paratuberculosis, is a severe production-limiting disease with significant economic and welfare implications for the global cattle industry. Caused by infection with Mycobacterium avium subspecies paratuberculosis (MAP), JD manifests as chronic enteritis in infected cattle. In addition to the economic losses and animal welfare issues associated with JD, MAP has attracted public health concerns with potential association with Crohn's disease, a human inflammatory bowel disease. The lack of effective treatment options, such as a vaccine, has hampered JD control resulting in its increasing global prevalence. The disease was first reported in 1895, but in recognition of its growing economic impact, extensive recent research facilitated by a revolution in technological approaches has led to significantly enhanced understanding of the immunological, genetic, and pathogen factors influencing disease pathogenesis. This knowledge has been derived from a variety of diverse models to elucidate host-pathogen interactions including in vivo and in vitro experimental infection models, studies measuring immune parameters in naturally-infected animals, and by studies conducted at the population level to enable the estimation of genetic parameters, and the identification of genetic markers and quantitative trait loci (QTL) putatively associated with susceptibility or resistance to JD. The main objectives of this review are to summarize these recent developments from an immunogenetics perspective and attempt to extract the principal and common findings emerging from this wealth of recent information. Based on these analyses, and in light of emerging technologies such as gene-editing, we conclude by discussing potential future avenues for effectively mitigating JD in cattle.

Conformational control of Cas9 by CRISPR hybrid RNA-DNA guides mitigates off-target activity in T cells

The off-target activity of the CRISPR-associated nuclease Cas9 is a potential concern for therapeutic genome editing applications. Although high-fidelity Cas9 variants have been engineered, they exhibit varying efficiencies and have residual off-target effects, limiting their applicability. Here, we show that CRISPR hybrid RNA-DNA (chRDNA) guides provide an effective approach to increase Cas9 specificity while preserving on-target editing activity. Across multiple genomic targets in primary human T cells, we show that 2'-deoxynucleotide (dnt) positioning affects guide activity and specificity in a target-dependent manner and that this can be used to engineer chRDNA guides with substantially reduced off-target effects. Crystal structures of DNA-bound Cas9-chRDNA complexes reveal distorted guide-target duplex geometry and allosteric modulation of Cas9 conformation. These structural effects increase specificity by perturbing DNA hybridization and modulating Cas9 activation kinetics to disfavor binding and cleavage of off-target substrates. Overall, these results pave the way for utilizing customized chRDNAs in clinical applications.

Chromothripsis as an on-target consequence of CRISPR-Cas9 genome editing

Genome editing has therapeutic potential for treating genetic diseases and cancer. However, the currently most practicable approaches rely on the generation of DNA double-strand breaks (DSBs), which can give rise to a poorly characterized spectrum of chromosome structural abnormalities. Here, using model cells and single-cell whole-genome sequencing, as well as by editing at a clinically relevant locus in clinically relevant cells, we show that CRISPR-Cas9 editing generates structural defects of the nucleus, micronuclei and chromosome bridges, which initiate a mutational process called chromothripsis. Chromothripsis is extensive chromosome rearrangement restricted to one or a few chromosomes that can cause human congenital disease and cancer. These results demonstrate that chromothripsis is a previously unappreciated on-target consequence of CRISPR-Cas9-generated DSBs. As genome editing is implemented in the clinic, the potential for extensive chromosomal rearrangements should be considered and monitored.

AsCRISPR: A Web Server for Allele-Specific Single Guide RNA Design in Precision Medicine

Allele-specific genomic targeting by CRISPR is a versatile strategy that has been increasingly exploited not only in treating inherited dominant diseases and mutation-driven cancers, but also in other important fields such as genome imprinting, haploinsufficiency, and genome loci imaging. Despite its tremendous utilities, few bioinformatic tools have been implemented for the allele-specific purpose of CRISPR. We thus developed AsCRISPR (Allele-specific CRISPR), a comprehensive web tool to aid the design of short-guide RNA (sgRNA) sequences that can discriminate between alleles. AsCRISPR allows users to analyze both their own identified variants and heterozygous single nucleotide polymorphisms and, importantly, output the candidate sgRNAs and their quality control information. To facilitate targeting dominant diseases, AsCRISPR analyzed dominant single nucleotide variants (SNVs) retrieved from ClinVar and OMIM databases, and generated a dominant database of candidate-discriminating sgRNAs that may specifically target the alternative allele for each dominant SNV site. Moreover, a validated database was established, which manually curated the discriminating sgRNAs that were experimentally validated in the mounting literature for multiple allele-specific purposes.

An ultrasensitive hybridization chain reaction-amplified CRISPR-Cas12a aptasensor for extracellular vesicle surface protein quantification

Tumor-derived extracellular vesicle (TEV) protein biomarkers facilitate cancer diagnosis and prognostic evaluations. However, the lack of reliable and convenient quantitative methods for evaluating TEV proteins prevents their clinical application.

Methods: Here, based on dual amplification of hybridization chain reaction (HCR) and CRISPR-Cas12a, we developed the apta-HCR-CRISPR assay for direct high-sensitivity detection of TEV proteins. The TEV protein-targeted aptamer was amplified by HCR to produce a long-repeated sequence comprising multiple CRISPR RNA (crRNA) targetable barcodes, and the signals were further amplified by CRISPR-Cas12a collateral cleavage activities, resulting in a fluorescence signal.

Results: The established strategy was verified by detecting the TEV protein markers nucleolin and programmed death ligand 1 (PD-L1). Both achieved limit of detection (LOD) values as low as 10² particles/µL, which is at least 10⁴-fold more sensitive than aptamer-ELISA and 10²-fold more sensitive than apta-HCR-ELISA. We directly applied our assay to a clinical analysis of circulating TEVs from 50 µL of serum, revealing potential applications of nucleolin⁺ TEVs for nasopharyngeal carcinoma cancer (NPC) diagnosis and PD-L1⁺ TEVs for therapeutic monitoring.

Conclusion: The platform was simple and easy to operate, and this approach should be useful for the highly sensitive and versatile quantification of TEV proteins in clinical samples.

Array-Based Epigenetic Aging Indices May Be Racially Biased

Epigenetic aging (EA) indices are frequently used as predictors of mortality and other important health outcomes. However, each of the commonly used array-based indices has significant heritable components which could tag ethnicity and potentially confound comparisons across racial and ethnic groups. To determine if this was possible, we examined the relationship of DNA methylation in cord blood from 203 newborns (112 African American (AA) and 91 White) at the 513 probes from the Levine PhenoAge Epigenetic Aging index to ethnicity. Then, we examined all sites significantly associated with race in the newborn sample to determine if they were also associated with an index of ethnic genetic heritage in a cohort of 505 AA adults. After Bonferroni correction, methylation at 50 CpG sites was significantly associated with ethnicity in the newborn cohort. The five most significant sites predicted ancestry with a receiver operator characteristic area under the curve of 0.97. Examination of the top 50 sites in the AA adult cohort showed that methylation status at 11 of those sites was also associated with percentage European ancestry. We conclude that the Levine PhenoAge Index is influenced by cryptic ethnic-specific genetic influences. This influence may extend to similarly constructed EA indices and bias cross-race comparisons.