Diffuse large B-cell lymphoma (DLBCL) is characterized by extensive genetic heterogeneity, and this results in unpredictable responses to the current treatment, R-CHOP, which consists of a cancer drug combination supplemented with the humanized CD20-targeting monoclonal antibody rituximab. Despite improvements in the patient response rate through rituximab addition to the treatment plan, up to 40% of DLBCL patients end in a relapsed or refractory state due to inherent or acquired resistance to the regimen. Here, we employ a lentiviral genome-wide clustered regularly interspaced short palindromic repeats library screening approach to identify genes involved in facilitating the rituximab response in cancerous B cells. Along with the CD20-encoding MS4A1 gene, we identify genes related to B-cell receptor (BCR) signaling as mediators of the intracellular signaling response to rituximab. More specifically, the B-cell linker protein (BLNK) and Bruton's tyrosine kinase (BTK) genes stand out as pivotal genes in facilitating direct rituximab-induced apoptosis through mechanisms that occur alongside complement-dependent cytotoxicity (CDC). Our findings demonstrate that rituximab triggers BCR signaling in a BLNK- and BTK-dependent manner and support the existing notion that intertwined CD20 and BCR signaling pathways in germinal center B-cell-like-subtype DLBCL lead to programmed cell death.

Allele-specific protospacer adjacent motif (asPAM)-positioning SNPs and CRISPRs are valuable resources for gene therapy of dominant disorders. However, one technical hurdle is to identify the haplotype comprising the disease-causing allele and the distal asPAM SNPs. Here, we describe a novel CRISPR-based method (CRISPR-hapC) for haplotyping. Based on the generation (with a pair of CRISPRs) of extrachromosomal circular DNA in cells, the CRISPR-hapC can map haplotypes from a few hundred bases to over 200 Mb. To streamline and demonstrate the applicability of the CRISPR-hapC and asPAM CRISPR for allele-specific gene editing, we reanalyzed the 1000 human pan-genome and generated a high frequency asPAM SNP and CRISPR database (www.crispratlas.com/knockout) for four CRISPR systems (SaCas9, SpCas9, xCas9 and Cas12a). Using the huntingtin (HTT) CAG expansion and transthyretin (TTR) exon 2 mutation as examples, we showed that the asPAM CRISPRs can specifically discriminate active and dead PAMs for all 23 loci tested. Combination of the CRISPR-hapC and asPAM CRISPRs further demonstrated the capability for achieving highly accurate and haplotype-specific deletion of the HTT CAG expansion allele and TTR exon 2 mutation in human cells. Taken together, our study provides a new approach and an important resource for genome research and allele-specific (haplotype-specific) gene therapy.

There is a need for functional genome-wide annotation of the protein-coding genes to get a deeper understanding of mammalian biology. Here, a new annotation strategy is introduced based on dimensionality reduction and density-based clustering of whole-body co-expression patterns. This strategy has been used to explore the gene expression landscape in pig, and we present a whole-body map of all protein-coding genes in all major pig tissues and organs.

The transcription factor cyclic-AMP response element-binding protein 1 (CREB1) responds to cAMP level and controls the expression of target genes, which regulates nutrition partitioning. The promoters of CREB1-targeted genes responsive to cAMP have been extensively investigated and characterized with the presence of both cAMP response element and TATA box. Compelling evidence demonstrates that CREB1 also plays an essential role in promoting tumor development. However, only very few genes required for cell survival, proliferation and migration are known to be constitutively regulated by CREB1 in tumors. Their promoters mostly do not harbor any cAMP response element. Thus, it is very likely that CREB1 regulates the expressions of distinct sets of target genes in normal tissues and tumors. The whole gene network constitutively regulated by CREB1 in tumors has remained unrevealed. Here, we employ a systematical and integrative approach to decipher this gene network in the context of both tissue cultured cancer cells and patient samples. We combine transcriptomic, Rank-Rank Hypergeometric Overlap, and Chipseq analysis, to define and characterize CREB1-regulated genes in a multidimensional fashion. A strong cancer relevance of those top-ranked targets, which meet the most stringent criteria, is eventually verified by overall survival analysis of cancer patients. These findings strongly suggest the importance of genes constitutively regulated by CREB1 for their implicative involvement in promoting tumorigenesis.

CRISPR gene therapy is one promising approach for treatment of Duchenne muscular dystrophy (DMD), which is caused by a large spectrum of mutations in the dystrophin gene. To broaden CRISPR gene editing strategies for DMD treatment, we report the efficient restoration of dystrophin expression in induced myotubes by SpCas9 and dual guide RNAs (gRNAs). We first sequenced 32 deletion junctions generated by this editing method and revealed that non-homologous blunt-end joining represents the major indel type. Based on this predictive repair outcome, efficient in-frame deletion of a part of DMD exon 51 was achieved in HEK293T cells with plasmids expressing SpCas9 and dual gRNAs. More importantly, we further corrected a frameshift mutation in human DMD (exon45del) fibroblasts with SpCas9-dual gRNA ribonucleoproteins. The edited DMD fibroblasts were transdifferentiated into myotubes by lentiviral-mediated overexpression of a human MYOD transcription factor. Restoration of DMD expression at both the mRNA and protein levels was confirmed in the induced myotubes. With further development, the combination of SpCas9-dual gRNA-corrected DMD patient fibroblasts and transdifferentiation may provide a valuable therapeutic strategy for DMD.

Extrachromosomal circular DNA (eccDNA) has emerged as a promising biomarker for disease diagnosis and prognosis prediction. However, its role in type 2 diabetes remains unexplored.

Circumstantial evidence suggests that B cells may instruct T cells to break tolerance. Here, to test this hypothesis, we used a murine model in which a single B cell clone precipitates an autoreactive response resembling systemic lupus erythematosus (SLE). The initiating clone did not need to enter germinal centers to precipitate epitope spreading. Rather, it localized to extrafollicular splenic bridging channels early in the response. Autoantibody produced by the initiating clone was not sufficient to drive the autoreactive response. Subsequent epitope spreading depended on antigen presentation and was compartmentalized by major histocompatibility complex (MHC). B cells carrying two MHC haplotypes could bridge the MHC barrier between B cells that did not share MHC. Thus, B cells directly relay autoreactivity between two separate compartments of MHC-restricted T cells, leading to inclusion of distinct B cell populations in germinal centers. Our findings demonstrate that B cells initiate and propagate the autoimmune response.

FUS (TLS) and EWS (EWSR1) belong to the FET-protein family of RNA and DNA binding proteins. FUS and EWS are structurally and functionally related and participate in transcriptional regulation and RNA processing. FUS and EWS are identified in translocation generated cancer fusion proteins and involved in the human neurological diseases amyotrophic lateral sclerosis and fronto-temporal lobar degeneration.

EGFR-mutated non-small cell lung cancer patients experience relapse within 1-2 years of treatment with EGFR-inhibitors, such as erlotinib. Multiple resistance mechanisms have been identified including secondary EGFR-mutations, MET-amplification, and epithelial-mesenchymal transition (EMT). Previous studies have indicated a role of Insulin-like growth factor 1 receptor (IGF1R) in acquired resistance to EGFR-directed drugs as well as in EMT. In the present study, we have investigated the involvement of IGF1R in acquired high-dose erlotinib resistance in the EGFR-mutated lung adenocarcinoma cell line HCC827. We observed that IGF1R was upregulated in the immediate response to erlotinib and hyperactivated in erlotinib resistant HCC827 cells. Resistant cells additionally acquired features of EMT, whereas MET-amplification and secondary EGFR-mutations were absent. Using CRISPR/Cas9, we generated a HCC827(IGFR1-/-) cell line and subsequently investigated resistance development in response to high-dose erlotinib. Interestingly, HCC827(IGFR1-/-) cells were now observed to specifically amplify the MET gene. Additionally, we observed a reduced level of mesenchymal markers in HCC827(IGFR1-/-) indicating an intrinsic enhanced epithelial signature compared to HCC827 cells. In conclusion, our data show that IGF1R have an important role in defining selected resistance mechanisms in response to high doses of erlotinib.

Fusion of DNA methyltransferase domains to the nuclease-deficient clustered regularly interspaced short palindromic repeat (CRISPR) associated protein 9 (dCas9) has been used for epigenome editing, but the specificities of these dCas9 methyltransferases have not been fully investigated.

The P-type ATPases family consists of ion and lipid transporters. Their unique diversity in function and expression is critical for normal development. In this study we investigated human pluripotent stem cells (hPSC) and different neural progenitor states to characterize the expression of the plasma membrane calcium ATPases (PMCAs) during human neural development and in mature mesencephalic dopaminergic (mesDA) neurons. Our RNA sequencing data identified a dynamic change in ATPase expression correlating with the differentiation time of the neural progenitors, which was independent of the neuronal progenitor type. Expression of ATP2B1 and ATP2B4 were the most abundantly expressed, in accordance with their main role in Ca2+ regulation and we observed all of the PMCAs to have a subcellular punctate localization. Interestingly in hPSCs ATP2B1 and ATP2B3 were highly expressed in a cell cycle specific manner and ATP2B2 and ATP2B4 were highly expressed in a hPSC sub-population. In neural rosettes a strong apical PMCA expression was identified in the luminal region. Lastly, we confirmed all PMCAs to be expressed in mesDA neurons, however at varying levels. Our results reveal that PMCA expression dynamically changes during stem cell differentiation and highlights the diverging needs of cell populations to regulate and properly integrate Ca2+ changes, which can ultimately correspond to changes in specific stem cell transcription states.

The heterogeneity of endothelial cells (ECs) across tissues remains incompletely inventoried. We constructed an atlas of >32,000 single-EC transcriptomes from 11 mouse tissues and identified 78 EC subclusters, including Aqp7+ intestinal capillaries and angiogenic ECs in healthy tissues. ECs from brain/testis, liver/spleen, small intestine/colon, and skeletal muscle/heart pairwise expressed partially overlapping marker genes. Arterial, venous, and lymphatic ECs shared more markers in more tissues than did heterogeneous capillary ECs. ECs from different vascular beds (arteries, capillaries, veins, lymphatics) exhibited transcriptome similarity across tissues, but the tissue (rather than the vessel) type contributed to the EC heterogeneity. Metabolic transcriptome analysis revealed a similar tissue-grouping phenomenon of ECs and heterogeneous metabolic gene signatures in ECs between tissues and between vascular beds within a single tissue in a tissue-type-dependent pattern. The EC atlas taxonomy enabled identification of EC subclusters in public scRNA-seq datasets and provides a powerful discovery tool and resource value.

With a five-year survival rate of 9%, pancreatic ductal adenocarcinoma (PDAC) is the deadliest of all cancers. The rapid mortality makes PDAC difficult to research, and inspires a resolve to create reliable, tractable cellular models for preclinical cancer research. Organoids are increasingly used to model PDAC as they maintain the differentiation status, molecular and genomic signatures of the original tumour. In this paper, we present novel methodologies and experimental approaches to develop PDAC organoids from PDX tumours, and the simultaneous development of matched primary cell lines. Moreover, we also present a method of recapitulating primary cell line cultures to organoids (CLOs). We highlight the usefulness of CLOs as PDAC organoid models, as they maintain similar transcriptomic signatures as their matched patient-derived organoids and patient derived xenografts (PDX)s. These models provide a manageable, expandable in vitro resource for downstream applications such as high throughput screening, functional genomics, and tumour microenvironment studies.

RTEL1 helicase is a component of DNA repair and telomere maintenance machineries. While RTEL1's role in DNA replication is emerging, how RTEL1 preserves genomic stability during replication remains elusive. Here we used a range of proteomic, biochemical, cell, and molecular biology and gene editing approaches to provide further insights into potential role(s) of RTEL1 in DNA replication and genome integrity maintenance. Our results from complementary human cell culture models established that RTEL1 and the Polδ subunit Poldip3 form a complex and are/function mutually dependent in chromatin binding after replication stress. Loss of RTEL1 and Poldip3 leads to marked R-loop accumulation that is confined to sites of active replication, enhances endogenous replication stress, and fuels ensuing genomic instability. The impact of depleting RTEL1 and Poldip3 is epistatic, consistent with our proposed concept of these two proteins operating in a shared pathway involved in DNA replication control under stress conditions. Overall, our data highlight a previously unsuspected role of RTEL1 and Poldip3 in R-loop suppression at genomic regions where transcription and replication intersect, with implications for human diseases including cancer.

The generation of large transgenic animals is impeded by complex cloning, long maturation and gastrulation times. An introduction of multiple gene alterations increases the complexity. We have cloned a transgenic Cas9 minipig to introduce multiple mutations by CRISPR in somatic cells. Transgenic Cas9 pigs were generated by somatic cell nuclear transfer and were backcrossed to Göttingen Minipigs for two generations. Cas9 expression was controlled by FlpO-mediated recombination and was visualized by translation from red to yellow fluorescent protein. In vitro analyses in primary fibroblasts, keratinocytes and lung epithelial cells confirmed the genetic alterations executed by the viral delivery of single guide RNAs (sgRNA) to the target cells. Moreover, multiple gene alterations could be introduced simultaneously in a cell by viral delivery of sgRNAs. Cells with loss of TP53, PTEN and gain-of-function mutation in KRASG12D showed increased proliferation, confirming a transformation of the primary cells. An in vivo activation of Cas9 expression could be induced by viral delivery to the skin. Overall, we have generated a minipig with conditional expression of Cas9, where multiple gene alterations can be introduced to somatic cells by viral delivery of sgRNA. The development of a transgenic Cas9 minipig facilitates the creation of complex pre-clinical models for cancer research.

Advances in molecular profiling have opened up the possibility to map the expression of genes in cells, tissues, and organs in the human body. Here, we combined single-cell transcriptomics analysis with spatial antibody-based protein profiling to create a high-resolution single-cell type map of human tissues. An open access atlas has been launched to allow researchers to explore the expression of human protein-coding genes in 192 individual cell type clusters. An expression specificity classification was performed to determine the number of genes elevated in each cell type, allowing comparisons with bulk transcriptomics data. The analysis highlights distinct expression clusters corresponding to cell types sharing similar functions, both within the same organs and between organs.

The design of CRISPR gRNAs requires accurate on-target efficiency predictions, which demand high-quality gRNA activity data and efficient modeling. To advance, we here report on the generation of on-target gRNA activity data for 10,592 SpCas9 gRNAs. Integrating these with complementary published data, we train a deep learning model, CRISPRon, on 23,902 gRNAs. Compared to existing tools, CRISPRon exhibits significantly higher prediction performances on four test datasets not overlapping with training data used for the development of these tools. Furthermore, we present an interactive gRNA design webserver based on the CRISPRon standalone software, both available via https://rth.dk/resources/crispr/ . CRISPRon advances CRISPR applications by providing more accurate gRNA efficiency predictions than the existing tools.

Recent progress in CRISPR gene editing tools has substantially increased the opportunities for curing devastating genetic diseases. Here we compare in-frame deletion by CRISPR-based non-homologous blunt end joining (NHBEJ), homology-directed repair (HDR), and prime editing (PE, PE2, and PE3)-based correction of two Duchenne Muscular Dystrophy (DMD) loss-of-function mutations (c.5533G>T and c.7893delC). To enable accurate and rapid evaluation of editing efficiency, we generated a genomically integrated synthetic reporter system (VENUS) carrying the DMD mutations. The VENUS contains a modified enhanced green fluorescence protein (EGFP) gene, in which expression was restored upon the CRISPR-mediated correction of DMD loss-of-function mutations. We observed that the highest editing efficiency was achieved by NHBEJ (74-77%), followed by HDR (21-24%) and PE2 (1.5%) in HEK293T VENUS reporter cells. A similar HDR (23%) and PE2 (1.1%) correction efficiency is achieved in fibroblast VENUS cells. With PE3 (PE2 plus nicking gRNA), the c.7893delC correction efficiency was increased 3-fold. Furthermore, an approximately 31% correction efficiency of the endogenous DMD: c.7893delC is achieved in the FACS-enriched HDR-edited VENUS EGFP+ patient fibroblasts. We demonstrated that a highly efficient correction of DMD loss-of-function mutations in patient cells can be achieved by several means of CRISPR gene editing.

A major challenge of CRISPR/Cas9-mediated genome engineering is that not all guide RNAs (gRNAs) cleave the DNA efficiently. Although the heterogeneity of gRNA activity is well recognized, the current understanding of how CRISPR/Cas9 activity is regulated remains incomplete. Here, we identify a sweet spot range of binding free energy change for optimal efficiency which largely explains why gRNAs display changes in efficiency at on- and off-target sites, including why gRNAs can cleave an off-target with higher efficiency than the on-target. Using an energy-based model, we show that local gRNA-DNA interactions resulting from Cas9 "sliding" on overlapping protospacer adjacent motifs (PAMs) profoundly impact gRNA activities. Combining the effects of local sliding for a given PAM context with global off-targets allows us to better identify highly specific, and thus efficient, gRNAs. We validate the effects of local sliding on gRNA efficiency using both public data and in-house data generated by measuring SpCas9 cleavage efficiency at 1024 sites designed to cover all possible combinations of 4-nt PAM and context sequences of 4 gRNAs. Our results provide insights into the mechanisms of Cas9-PAM compatibility and cleavage activation, underlining the importance of accounting for local sliding in gRNA design.

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.