Regulation of Coding and Non-coding Genes: New insights obtained through analysis of high-throughput sequencing data

The genetic code of a cell is kept in its DNA. However, a vast number of functions of a cell are carried out by proteins. Through gene expression the genetic code can be expressed and give rise to proteins. The expression of genes into proteins follows two steps: transcription of DNA into RNA and translation of (messenger-)RNA into proteins. This thesis concerns these two steps in gene expression, which are very different from a molecular point-ofview

An iterative workflow for mining the human intestinal metaproteome

Peer reviewe

Genome-wide profiling of p53-regulated enhancer RNAs uncovers a subset of enhancers controlled by a lncRNA

p53 binds enhancers to regulate key target genes. Here, we globally mapped p53-regulated enhancers by looking at enhancer RNA (eRNA) production. Intriguingly, while many p53-induced enhancers contained p53-binding sites, most did not. As long non-coding RNAs(lncRNAs) are prominent regulators of chromatin dynamics, we hypothesized that p53-induced lncRNAs contribute to the activation of enhancers by p53. Among p53-induced lncRNAs, we identified LED and demonstrate that its suppression attenuates p53 function. Chromatin-binding and eRNA expression analyses show that LED associates with and activates strong enhancers. One prominent target of LED was located at an enhancer region within CDKN1A gene, a potent p53-responsive cell cycle inhibitor. LED knockdown reduces CDKN1A enhancer induction and activity, and cell cycle arrest following p53 activation. Finally, promoter-associated hypermethylation analysis shows silencing of LED in human tumours. Thus, our study identifies a new layer of complexity in the p53 pathway and suggests its dysregulation in cancer

Alcohol-derived DNA crosslinks are repaired by two distinct mechanisms

Acetaldehyde is a highly reactive, DNA-damaging metabolite that is produced upon alcohol consumption1. Impaired detoxification of acetaldehyde is common in the Asian population, and is associated with alcohol-related cancers1,2. Cells are protected against acetaldehyde-induced damage by DNA crosslink repair, which when impaired causes Fanconi anaemia (FA), a disease resulting in failure to produce blood cells and a predisposition to cancer3,4. The combined inactivation of acetaldehyde detoxification and the FA pathway induces mutation, accelerates malignancies and causes the rapid attrition of blood stem cells5,6,7. However, the nature of the DNA damage induced by acetaldehyde and how this is repaired remains a key question. Here we generate acetaldehyde-induced DNA interstrand crosslinks and determine their repair mechanism in Xenopus egg extracts. We find that two replication-coupled pathways repair these lesions. The first is the FA pathway, which operates using excision—analogous to the mechanism used to repair the interstrand crosslinks caused by the chemotherapeutic agent cisplatin. However, the repair of acetaldehyde-induced crosslinks results in increased mutation frequency and an altered mutational spectrum compared with the repair of cisplatin-induced crosslinks. The second repair mechanism requires replication fork convergence, but does not involve DNA incisions—instead the acetaldehyde crosslink itself is broken. The Y-family DNA polymerase REV1 completes repair of the crosslink, culminating in a distinct mutational spectrum. These results define the repair pathways of DNA interstrand crosslinks caused by an endogenous and alcohol-derived metabolite, and identify an excision-independent mechanism

Simultaneous quantification of protein-DNA contacts and transcriptomes in single cells.

Protein-DNA interactions are critical to the regulation of gene expression, but it remains challenging to define how cell-to-cell heterogeneity in protein-DNA binding influences gene expression variability. Here we report a method for the simultaneous quantification of protein-DNA contacts by combining single-cell DNA adenine methyltransferase identification (DamID) with messenger RNA sequencing of the same cell (scDam&T-seq). We apply scDam&T-seq to reveal how genome-lamina contacts or chromatin accessibility correlate with gene expression in individual cells. Furthermore, we provide single-cell genome-wide interaction data on a polycomb-group protein, RING1B, and the associated transcriptome. Our results show that scDam&T-seq is sensitive enough to distinguish mouse embryonic stem cells cultured under different conditions and their different chromatin landscapes. Our method will enable the analysis of protein-mediated mechanisms that regulate cell-type-specific transcriptional programs in heterogeneous tissues

Using mitoribosomal profiling to investigate human mitochondrial translation [version 2; referees: 2 approved]

Background: Gene expression in human mitochondria has various idiosyncratic features. One of these was recently revealed as the unprecedented recruitment of a mitochondrially-encoded tRNA as a structural component of the large mitoribosomal subunit. In porcine particles this is mt-tRNAPhe whilst in humans it is mt-tRNAVal. We have previously shown that when a mutation in mt-tRNAVal causes very low steady state levels, there is preferential recruitment of mt-tRNAPhe. We have investigated whether this altered mitoribosome affects intra-organellar protein synthesis. Methods: By using mitoribosomal profiling we have revealed aspects of mitoribosome behaviour with its template mt-mRNA under both normal conditions as well as those where the mitoribosome has incorporated mt-tRNAPhe. Results: Analysis of the mitoribosome residency on transcripts under control conditions reveals that although mitochondria employ only 22 mt-tRNAs for protein synthesis, the use of non-canonical wobble base pairs at codon position 3 does not cause any measurable difference in mitoribosome occupancy irrespective of the codon. Comparison of the profile of aberrant mt-tRNAPhe containing mitoribosomes with those of controls that integrate mt-tRNAVal revealed that the impaired translation seen in the latter was not due to stalling on triplets encoding either of these amino acids. The alterations in mitoribosome interactions with start codons was not directly attributable to the either the use of non-cognate initiation codons or the presence or absence of 5’ leader sequences, except in the two bicistronic RNA units, RNA7 and RNA14 where the initiation sites are internal. Conclusions: These data report the power of mitoribosomal profiling in helping to understand the subtleties of mammalian mitochondrial protein synthesis. Analysis of profiles from the mutant mt-tRNAVal cell line suggest that despite mt-tRNAPhe being preferred in the porcine mitoribosome, its integration into the human counterpart results in a suboptimal structure that modifies its interaction with mt-mRNAs

Using mitoribosomal profiling to investigate human mitochondrial translation [version 1; referees: 2 approved]

Background: Gene expression in human mitochondria has various idiosyncratic features. One of these was recently revealed as the unprecedented recruitment of a mitochondrially-encoded tRNA as a structural component of the large mitoribosomal subunit. In porcine particles this is mt-tRNAPhe whilst in humans it is mt-tRNAVal. We have previously shown that when a mutation in mt-tRNAVal causes very low steady state levels, there is preferential recruitment of mt-tRNAPhe. We have investigated whether this altered mitoribosome affects intra-organellar protein synthesis. Methods: By using mitoribosomal profiling we have revealed aspects of mitoribosome behaviour with its template mt-mRNA under both normal conditions as well as those where the mitoribosome has incorporated mt-tRNAPhe. Results: Analysis of the mitoribosome residency on transcripts under control conditions reveals that although mitochondria employ only 22 mt-tRNAs for protein synthesis, the use of non-canonical wobble base pairs at codon position 3 does not cause any measurable difference in mitoribosome occupancy irrespective of the codon. Comparison of the profile of aberrant mt-tRNAPhe containing mitoribosomes with those of controls that integrate mt-tRNAVal revealed that the impaired translation seen in the latter was not due to stalling on triplets encoding either of these amino acids. The alterations in mitoribosome interactions with start codons was not directly attributable to the either the use of non-cognate initiation codons or the presence or absence of 5’ leader sequences, except in the two bicistronic RNA units, RNA7 and RNA14 where the initiation sites are internal. Conclusions: These data report the power of mitoribosomal profiling in helping to understand the subtleties of mammalian mitochondrial protein synthesis. Analysis of profiles from the mutant mt-tRNAVal cell line suggest that despite mt-tRNAPhe being preferred in the porcine mitoribosome, its integration into the human counterpart results in a suboptimal structure that modifies its interaction with mt-mRNAs

Simultaneous quantification of protein-DNA interactions and transcriptomes in single cells with scDam&T-seq.

Protein-DNA interactions are essential for establishing cell type-specific chromatin architecture and gene expression. We recently developed scDam&T-seq, a multi-omics method that can simultaneously quantify protein-DNA interactions and the transcriptome in single cells. The method effectively combines two existing methods: DNA adenine methyltransferase identification (DamID) and CEL-Seq2. DamID works through the tethering of a protein of interest (POI) to the Escherichia coli DNA adenine methyltransferase (Dam). Upon expression of this fusion protein, DNA in proximity to the POI is methylated by Dam and can be selectively digested and amplified. CEL-Seq2, in contrast, makes use of poly-dT primers to reverse transcribe mRNA, followed by linear amplification through in vitro transcription. scDam&T-seq is the first technique capable of providing a combined readout of protein-DNA contact and transcription from single-cell samples. Once suitable cell lines have been established, the protocol can be completed in 5 d, with a throughput of hundreds to thousands of cells. The processing of raw sequencing data takes an additional 1-2 d. Our method can be used to understand the transcriptional changes a cell undergoes upon the DNA binding of a POI. It can be performed in any laboratory with access to FACS, robotic and high-throughput-sequencing facilities