Genomic features underlie the co-option of SVA transposons as cis-regulatory elements in human pluripotent stem cells

Domestication of transposable elements (TEs) into functional cis-regulatory elements is a widespread phenomenon. However, the mechanisms behind why some TEs are co-opted as functional enhancers while others are not are underappreciated. SINE-VNTR-Alus (SVAs) are the youngest group of transposons in the human genome, where ~3,700 copies are annotated, nearly half of which are human-specific. Many studies indicate that SVAs are among the most frequently co-opted TEs in human gene regulation, but the mechanisms underlying such processes have not yet been thoroughly investigated. Here, we leveraged CRISPR-interference (CRISPRi), computational and functional genomics to elucidate the genomic features that underlie SVA domestication into human stem-cell gene regulation. We found that ~750 SVAs are co-opted as functional cis-regulatory elements in human induced pluripotent stem cells. These SVAs are significantly closer to genes and harbor more transcription factor binding sites than non-co-opted SVAs. We show that a long DNA motif composed of flanking YY1/2 and OCT4 binding sites is enriched in the co-opted SVAs and that these two transcription factors bind consecutively on the TE sequence. We used CRISPRi to epigenetically repress active SVAs in stem cell-like NCCIT cells. Epigenetic perturbation of active SVAs strongly attenuated YY1/OCT4 binding and influenced neighboring gene expression. Ultimately, SVA repression resulted in ~3,000 differentially expressed genes, 131 of which were the nearest gene to an annotated SVA. In summary, we demonstrated that SVAs modulate human gene expression, and uncovered that location and sequence composition contribute to SVA domestication into gene regulatory networks.

Epigenetics of Mitochondria-Associated Genes in Striated Muscle

Striated muscle has especially large energy demands. We identified 97 genes preferentially expressed in skeletal muscle and heart, but not in aorta, and found significant enrichment for mitochondrial associations among them. We compared the epigenomic and transcriptomic profiles of the 27 genes associated with striated muscle and mitochondria. Many showed strong correlations between their tissue-specific transcription levels, and their tissue-specific promoter, enhancer, or open chromatin as well as their DNA hypomethylation. Their striated muscle-specific enhancer chromatin was inside, upstream, or downstream of the gene, throughout much of the gene as a super-enhancer (CKMT2, SLC25A4, and ACO2), or even overlapping a neighboring gene (COX6A2, COX7A1, and COQ10A). Surprisingly, the 3′ end of the 1.38 Mb PRKN (PARK2) gene (involved in mitophagy and linked to juvenile Parkinson’s disease) displayed skeletal muscle/myoblast-specific enhancer chromatin, a myoblast-specific antisense RNA, as well as brain-specific enhancer chromatin. We also found novel tissue-specific RNAs in brain and embryonic stem cells within PPARGC1A (PGC-1α), which encodes a master transcriptional coregulator for mitochondrial formation and metabolism. The tissue specificity of this gene’s four alternative promoters, including a muscle-associated promoter, correlated with nearby enhancer chromatin and open chromatin. Our in-depth epigenetic examination of these genes revealed previously undescribed tissue-specific enhancer chromatin, intragenic promoters, regions of DNA hypomethylation, and intragenic noncoding RNAs that give new insights into transcription control for this medically important set of genes.

GntR-like SCO3932 Protein Provides a Link between Actinomycete Integrative and Conjugative Elements and Secondary Metabolism

Streptomyces bacteria produce a plethora of secondary metabolites including the majority of medically important antibiotics. The onset of secondary metabolism is correlated with morphological differentiation and controlled by a complex regulatory network involving numerous regulatory proteins. Control over these pathways at the molecular level has a medical and industrial importance. Here we describe a GntR-like DNA binding transcription factor SCO3932, encoded within an actinomycete integrative and conjugative element, which is involved in the secondary metabolite biosynthesis regulation. Affinity chromatography, electrophoresis mobility shift assay, footprinting and chromatin immunoprecipitation experiments revealed, both in vitro and in vivo, SCO3932 binding capability to its own promoter region shared with the neighboring gene SCO3933, as well as promoters of polyketide metabolite genes, such as cpkD, a coelimycin biosynthetic gene, and actII-orf4—an activator of actinorhodin biosynthesis. Increased activity of SCO3932 target promoters, as a result of SCO3932 overproduction, indicates an activatory role of this protein in Streptomyces coelicolor A3(2) metabolite synthesis pathways.

Phage-encoded ribosomal protein S21 expression is linked to late stage phage replication

The ribosomal protein S21 (bS21) gene has been detected in diverse viruses with a large range of genome sizes, yet its in situ expression and potential significance have not been investigated. Here, we report five closely related clades of bacteriophages (phages) represented by 47 genomes (8 curated to completion and up to 331 kbp in length) that encode a bS21 gene. The bS21 gene is on the reverse strand within a conserved region that encodes the large terminase, major capsid protein, prohead protease, portal vertex proteins and some hypothetical proteins. These phages are predicted to infect Bacteroidetes species that inhabit a range of depths in freshwater lakes. Transcriptionally active bS21-encoding phages were sampled in the late-stage of replication, when core structural genes, bS21 and a neighboring gene of unknown function were highly expressed. Thus, our analyses suggest that bS21, which is involved in translation initiation, substitutes into the Bacteroidetes ribosomes and selects for phage transcripts during the late-stage replication when large-scale phage protein production is required for assembly of phage particles.

A standard knockout procedure alters expression of adjacent loci at the translational level

The Saccharomyces cerevisiae gene deletion collection is widely used for functional gene annotation and genetic interaction analyses. However, the standard G418-resistance cassette used to produce knockout mutants delivers strong regulatory elements into the target genetic loci. To date, its side effects on the expression of neighboring genes have never been systematically assessed. Here, using ribosome profiling data, RT-qPCR, and reporter expression, we investigated perturbations induced by the KanMX module. Our analysis revealed significant alterations in the transcription efficiency of neighboring genes and, more importantly, severe impairment of their mRNA translation, leading to changes in protein abundance. In the ‘head-to-head’ orientation of the deleted and neighboring genes, knockout often led to a shift of the transcription start site of the latter, introducing new uAUG codon(s) into the expanded 5′ untranslated region (5′ UTR). In the ‘tail-to-tail’ arrangement, knockout led to activation of alternative polyadenylation signals in the neighboring gene, thus altering its 3′ UTR. These events may explain the so-called neighboring gene effect (NGE), i.e. false genetic interactions of the deleted genes. We estimate that in as much as ∼1/5 of knockout strains the expression of neighboring genes may be substantially (>2-fold) deregulated at the level of translation.

A standard knockout procedure alters expression of adjacent loci at the translational level

The S. cerevisiae gene deletion collection is widely used for functional gene annotation and genetic interaction analyses. However, the standard G418-resistance cassette used to produce knockout mutants delivers strong regulatory elements into the target genetic loci. To date, its side effects on the expression of neighboring genes have never been systematically assessed. Here, using ribosome profiling data, RT-qPCR, and reporter expression, we investigated perturbations induced by the KanMX module. Our analysis revealed significant alterations in the transcription efficiency of neighboring genes and, more importantly, severe impairment of their mRNA translation, leading to changes in protein abundance. In the “head-to-head” orientation of the neighbor and the deleted gene, knockout often led to a shift of the transcription start site of the neighboring gene, introducing new uAUG codon(s) into the expanded 5’ untranslated region (5’ UTR). In the “tail-to-tail” arrangement, knockout led to activation of alternative polyadenylation signals in the neighboring gene, thus altering its 3’ UTR. These events may explain the so-called neighboring gene effect (NGE), i.e. false genetic interactions of the deleted genes. We estimate that in as much as ∼1/5 of knockout strains the expression of neighboring genes may be substantially (>2-fold) deregulated at the level of translation.

Targeting RB1 Loss in Cancers

Retinoblastoma protein 1 (RB1) is encoded by a tumor suppressor gene that was discovered more than 30 years ago. Almost all mitogenic signals promote cell cycle progression by braking on the function of RB1 protein through mono- and subsequent hyper-phosphorylation mediated by cyclin-CDK complexes. The loss of RB1 function drives tumorigenesis in limited types of malignancies including retinoblastoma and small cell lung cancer. In a majority of human cancers, RB1 function is suppressed during tumor progression through various mechanisms. The latter gives rise to the acquisition of various phenotypes that confer malignant progression. The RB1-targeted molecules involved in such phenotypic changes are good quarries for cancer therapy. Indeed, a variety of novel therapies have been proposed to target RB1 loss. In particular, the inhibition of a number of mitotic kinases appeared to be synthetic lethal with RB1 deficiency. A recent study focusing on a neighboring gene that is often collaterally deleted together with RB1 revealed a pharmacologically targetable vulnerability in RB1-deficient cancers. Here we summarize current understanding on possible therapeutic approaches targeting functional or genomic aberration of RB1 in cancers.

Associations between forensic loci and neighboring gene expression levels may compromise medical privacy

A set of 20 short tandem repeats (STRs) is used by the United States criminal justice system to identify suspects, and to maintain a database of genetic profiles for individuals who have been previously convicted or arrested. Some of these STRs were identified in the 1990s, with a preference for markers in putative gene deserts to avoid forensic profiles revealing protected medical information. We revisit that assumption, investigating whether forensic genetic profiles reveal information about gene expression variation, or potential medical information. We find six significant correlations (FDR = 0.23) between the forensic STRs and the expression levels of neighboring genes in lymphoblastoid cell lines. We explore possible mechanisms for these associations, with evidence compatible with forensic STRs causing expression variation, or being in LD with a causal locus in three cases, and weaker or potentially spurious associations in the other three cases. Together, these results suggest that forensic genetic loci may reveal expression level and, perhaps, medical information.

Transposable element activity in the transcriptomic analysis of mouse pancreatic tumors

Transposable elements (TEs) are middle-repeated DNA sequences that can move along chromosomes using internal coding and regulatory regions. By their ability to move and because they are repeated, TEs can promote mutations. Especially they can alter the expression pattern of neighboring genes and have been shown to be involved in the mammalian regulatory network evolution. Human and mouse share more than 95% of their genomes and are affected by comparable diseases, which makes the mouse a perfect model in cancer research. However not much investigation concerning the mouse TE content has been made on this topics. In human cancer condition, a global activation of TEs can been observed which may ask the question of their impact on neighboring gene functioning. In this work, we used RNA sequences of highly aggressive pancreatic tumors from mouse to analyze the gene and TE deregulation happening in this condition compared to pancreas from healthy animals. Our results show that several TE families are deregulated and that the presence of TEs is associated with the expression divergence of genes in the tumor condition. These results illustrate the potential role of TEs in the global deregulation at work in the cancer cells.

Identification and elimination of genomic regions irrelevant for magnetosome biosynthesis by large-scale deletion in Magnetospirillum gryphiswaldense

Background Magnetosome formation in the alphaproteobacterium Magnetospirillum gryphiswaldense is controlled by more than 30 known mam and mms genes clustered within a large genomic region, the ‘magnetosome island’ (MAI), which also harbors numerous mobile genetic elements, repeats, and genetic junk. Because of the inherent genetic instability of the MAI caused by neighboring gene content, the elimination of these regions and their substitution by a compact, minimal magnetosome expression cassette would be important for future analysis and engineering. In addition, the role of the MAI boundaries and adjacent regions are still unclear, and recent studies indicated that further auxiliary determinants for magnetosome biosynthesis are encoded outside the MAI. However, techniques for large-scale genome editing of magnetic bacteria are still limited, and the full complement of genes controlling magnetosome formation has remained uncertain. Results Here we demonstrate that an allelic replacement method based on homologous recombination can be applied for large-scale genome editing in M. gryphiswaldense. By analysis of 24 deletion mutants covering about 167 kb of non-redundant genome content, we identified genes and regions inside and outside the MAI irrelevant for magnetosome biosynthesis. A contiguous stretch of ~ 100 kb, including the scattered mam and mms6 operons, could be functionally substituted by a compact and contiguous ~ 38 kb cassette comprising all essential biosynthetic gene clusters, but devoid of interspersing irrelevant or problematic gene content. Conclusions Our results further delineate the genetic complement for magnetosome biosynthesis and will be useful for future large-scale genome editing and genetic engineering of magnetosome biosynthesis.