Reusable Single Cell for Iterative Epigenomic Analyses

Current single-cell epigenome analyses are designed for single use. The cell is discarded after a single use, preventing analysis of multiple epigenetic ...

National Institutes of Health

Frederick National Laboratory for Cancer Research sponsored by the National Cancer Institute

Naked mole rat (MR) Heterocephalus glaber is a rodent model of delayed aging because of its unusually long life span (>28 years). It is also not known to develop cancer. In the current work, tissue imaging by x-ray fluorescence microscopy and direct analyses of trace elements revealed low levels of selenium in the MR liver and kidney, whereas MR and mouse brains had similar selenium levels. This effect was not explained by uniform selenium deficiency because methionine sulfoxide reductase activities were similar in mice and MR. However, glutathione peroxidase activity was an order of magnitude lower in MR liver and kidney than in mouse tissues. In addition, metabolic labeling of MR cells with (75)Se revealed a loss of the abundant glutathione peroxidase 1 (GPx1) band, whereas other selenoproteins were preserved. To characterize the MR selenoproteome, we sequenced its liver transcriptome. Gene reconstruction revealed standard selenoprotein sequences except for GPx1, which had an early stop codon, and SelP, which had low selenocysteine content. When expressed in HEK 293 cells, MR GPx1 was present in low levels, and its expression could be rescued neither by removing the early stop codon nor by replacing its SECIS element. In addition, GPx1 mRNA was present in lower levels in MR liver than in mouse liver. To determine if GPx1 deficiency could account for the reduced selenium content, we analyzed GPx1 knock-out mice and found reduced selenium levels in their livers and kidneys. Thus, MR is characterized by the reduced utilization of selenium due to a specific defect in GPx1 expression.

Reduced utilization of selenium by naked mole rats due to a specific defect in GPx1 expression.

Selenium is an essential trace element in mammals due to its presence in proteins in the form of selenocysteine (Sec). Human genome codes for 25 Sec-containing protein genes, and mouse and rat genomes for 24.

Composition and evolution of the vertebrate and mammalian selenoproteomes.

It is thought that the SelenoCysteine Insertion Sequence (SECIS) element and UGA codon are sufficient for selenocysteine (Sec) insertion. However, we found that UGA supported Sec insertion only at its natural position or in its close proximity in mammalian thioredoxin reductase 1 (TR1). In contrast, Sec could be inserted at any tested position in mammalian TR3. Replacement of the 3'-UTR of TR3 with the corresponding segment of a Euplotes crassus TR restricted Sec insertion into the C-terminal region, whereas the 3'-UTR of TR3 conferred unrestricted Sec insertion into E. crassus TR, in which Sec insertion is normally limited to the C-terminal region. Exchanges of 3'-UTRs between mammalian TR1 and E. crassus TR had no effect, as both proteins restricted Sec insertion. We further found that these effects could be explained by the use of selenoprotein-specific SECIS elements. Examination of Sec insertion into other selenoproteins was consistent with this model. The data indicate that mammals evolved the ability to limit Sec insertion into natural positions within selenoproteins, but do so in a selenoprotein-specific manner, and that this process is controlled by the SECIS element in the 3'-UTR. 

UGA codon position-dependent incorporation of selenocysteine into mammalian selenoproteins.

SelS (Selenoprotein S) is a selenocysteine-containing protein with roles in ER (endoplasmic reticulum) function and inflammation. It has been implicated in ERAD (ER-associated protein degradation), and clinical studies revealed an association of its promoter polymorphism with cytokine levels and human diseases. However, the pathways and interacting proteins that could shed light on pathogenesis of SelS-associated diseases have not been studied systematically. We performed a large-scale affinity isolation of human SelS and its mutant forms and analysed the proteins that interact with them. All previously known SelS targets and nearly two hundred additional proteins were identified that were remarkably enriched for various multiprotein complexes. Subsequent chemical cross-linking experiments identified the specific interacting sites in SelS and its several targets. Most of these interactions involved coiled-coil domains. The data suggest that SelS participates in intracellular membrane transport and maintenance of protein complexes by anchoring them to the ER membrane. 

Selenoprotein S is involved in maintenance and transport of multiprotein complexes.

The human genome contains 25 genes coding for selenocysteine-containing proteins (selenoproteins). These proteins are involved in a variety of functions, most notably redox homeostasis. Selenoprotein enzymes with known functions are designated according to these functions: TXNRD1, TXNRD2, and TXNRD3 (thioredoxin reductases), GPX1, GPX2, GPX3, GPX4, and GPX6 (glutathione peroxidases), DIO1, DIO2, and DIO3 (iodothyronine deiodinases), MSRB1 (methionine sulfoxide reductase B1), and SEPHS2 (selenophosphate synthetase 2). Selenoproteins without known functions have traditionally been denoted by SEL or SEP symbols. However, these symbols are sometimes ambiguous and conflict with the approved nomenclature for several other genes. Therefore, there is a need to implement a rational and coherent nomenclature system for selenoprotein-encoding genes. Our solution is to use the root symbol SELENO followed by a letter. This nomenclature applies to SELENOF (selenoprotein F, the 15-kDa selenoprotein, SEP15), SELENOH (selenoprotein H, SELH, C11orf31), SELENOI (selenoprotein I, SELI, EPT1), SELENOK (selenoprotein K, SELK), SELENOM (selenoprotein M, SELM), SELENON (selenoprotein N, SEPN1, SELN), SELENOO (selenoprotein O, SELO), SELENOP (selenoprotein P, SeP, SEPP1, SELP), SELENOS (selenoprotein S, SELS, SEPS1, VIMP), SELENOT (selenoprotein T, SELT), SELENOV (selenoprotein V, SELV), and SELENOW (selenoprotein W, SELW, SEPW1). This system, approved by the HUGO Gene Nomenclature Committee, also resolves conflicting, missing, and ambiguous designations for selenoprotein genes and is applicable to selenoproteins across vertebrates.

Selenoprotein Gene Nomenclature.

The ribosome can change its reading frame during translation in a process known as programmed ribosomal frameshifting. These rare events are supported by complex mRNA signals. However, we found that the ciliates Euplotes crassus and Euplotes focardii exhibit widespread frameshifting at stop codons. 47 different codons preceding stop signals resulted in either +1 or +2 frameshifts, and +1 frameshifting at AAA was the most frequent. The frameshifts showed unusual plasticity and rapid evolution, and had little influence on translation rates. The proximity of a stop codon to the 3' mRNA end, rather than its occurrence or sequence context, appeared to designate termination. Thus, a 'stop codon' is not a sufficient signal for translation termination, and the default function of stop codons in Euplotes is frameshifting, whereas termination is specific to certain mRNA positions and probably requires additional factors.

Position-dependent termination and widespread obligatory frameshifting in Euplotes translation.

Kaposi's sarcoma-associated herpesvirus (KSHV)-transformed primary effusion lymphoma cell lines contain ∼70 to 150 copies of episomal KSHV genomes per cell and have been widely used for studying the mechanisms of KSHV latency and lytic reactivation. Here, we report the first complete knockout (KO) of viral ORF57 gene from all ∼100 copies of KSHV genome per cell in BCBL-1 cells. This was achieved by a modified CRISPR/Cas9 technology to simultaneously express two guide RNAs (gRNAs) and Cas9 from a single expression vector in transfected cells in combination with multiple rounds of cell selection and single-cell cloning. CRISPR/Cas9-mediated genome engineering induces the targeted gene deletion and inversion  We found the inverted ORF57 gene in the targeted site in the KSHV genome in one of two characterized single cell clones. Knockout of ORF57 from the KSHV genome led to viral genome instability, thereby reducing viral genome copies and expression of viral lytic genes in BCBL-1-derived single-cell clones. The modified CRISPR/Cas9 technology was very efficient in knocking out the ORF57 gene in iSLK/Bac16 and HEK293/Bac36 cells, where each cell contains only a few copies of the KSHV genome. The ORF57 KO genome was stable in iSLK/Bac16 cells, and, upon lytic induction, was partially rescued by ectopic ORF57 to express viral lytic gene ORF59 and produce infectious virions. Together, the technology developed in this study has paved the way to express two separate gRNAs and the Cas9 enzyme simultaneously in the same cell and could be efficiently applied to any genetic alterations from various genomes, including those in extreme high copy numbers. This study provides the first evidence that CRISPR/Cas9 technology can be applied to knock out the ORF57 gene from all ∼100 copies of the KSHV genome in primary effusion lymphoma (PEL) cells by coexpressing two guide RNAs (gRNAs) and Cas9 from a single expression vector in combination with single-cell cloning. The gene knockout efficiency in this system was evaluated rapidly using a direct cell PCR screening. The current CRISPR/Cas9 technology also mediated ORF57 inversion  in the targeted site of the KSHV genome. The successful rescue of viral lytic gene expression and infectious virion production from the ORF57 knockout (KO) genome further reiterates the essential role of ORF57 in KSHV infection and multiplication. This modified technology should be useful for knocking out any viral genes from a genome to dissect functions of individual viral genes in the context of the virus genome and to understand their contributions to viral genetics and the virus life cycle.

CRISPR/Cas9-Mediated Knockout and  Inversion of the ORF57 Gene from All Copies of the Kaposi's Sarcoma-Associated Herpesvirus Genome in BCBL-1 Cells.

Gene expression in individual cells is epigenetically regulated by DNA modifications, histone modifications, transcription factors, and other DNA-binding proteins. It has been shown that multiple histone modifications can predict gene expression and reflect future responses of bulk cells to extracellular cues. However, the predictive ability of epigenomic analysis is still limited for mechanistic research at a single cell level. To overcome this limitation, it would be useful to acquire reliable signals from multiple epigenetic marks in the same single cell. Here, we propose a new approach and a new method for analysis of several components of the epigenome in the same single cell. The new method allows reanalysis of the same single cell. We found that reanalysis of the same single cell is feasible, provides confirmation of the epigenetic signals, and allows application of statistical analysis to identify reproduced reads using data sets generated only from the single cell. Reanalysis of the same single cell is also useful to acquire multiple epigenetic marks from the same single cells. The method can acquire at least five epigenetic marks: H3K27ac, H3K27me3, mediator complex subunit 1, a DNA modification, and a DNA-interacting protein. We can predict active signaling pathways in K562 single cells using the epigenetic data and confirm that the predicted results strongly correlate with actual active signaling pathways identified by RNA-seq results. These results suggest that the new method provides mechanistic insights for cellular phenotypes through multilayered epigenome analysis in the same single cells.

Alexei Lobanov

CCR Collaborative Bioinformatics Resource (CCBR),
Center for Cancer Research,
National Cancer Institute,
CCR Collaborative Bioinformatics Resource (CCBR), Center for Cancer Research, National Cancer Institute

National Institutes of Health

Reusable Single Cell for Iterative Epigenomic Analyses

Alexei Lobanov

.css-o250i3{font-size:18px;font-family:Open Sans,sans-serif;line-height:21.6px;}CCR Collaborative Bioinformatics Resource (CCBR),Center for Cancer Research,National Cancer Institute,CCR Collaborative Bioinformatics Resource (CCBR), Center for Cancer Research, National Cancer Institute

National Institutes of Health

Reusable Single Cell for Iterative Epigenomic Analyses

CCR Collaborative Bioinformatics Resource (CCBR),
Center for Cancer Research,
National Cancer Institute,
CCR Collaborative Bioinformatics Resource (CCBR), Center for Cancer Research, National Cancer Institute