Functional genomic studies of the structure and regulation of eukaryotic transcriptomes

The main focus of this thesis is the use of high-throughput sequencing technologies in functional genomics (in particular in the form of ChIP-seq, chromatin immunoprecipitation coupled with sequencing, and RNA-seq) and the study of the structure and regulation of transcriptomes. Some parts of it are of a more methodological nature while others describe the application of these functional genomic tools to address various biological problems. A significant part of the research presented here was conducted as part of the ENCODE (ENCyclopedia Of DNA Elements) Project.

The first part of the thesis focuses on the structure and diversity of the human transcriptome. Chapter 1 contains an analysis of the diversity of the human polyadenylated transcriptome based on RNA-seq data generated for the ENCODE Project. Chapter 2 presents a simulation-based examination of the performance of some of the most popular computational tools used to assemble and quantify transcriptomes. Chapter 3 includes a study of variation in gene expression, alternative splicing and allelic expression bias on the single-cell level and on a genome-wide scale in human lymphoblastoid cells; it also brings forward a number of critical to the practice of single-cell RNA-seq measurements methodological considerations.

The second part presents several studies applying functional genomic tools to the study of the regulatory biology of organellar genomes, primarily in mammals but also in plants. Chapter 5 contains an analysis of the occupancy of the human mitochondrial genome by TFAM, an important structural and regulatory protein in mitochondria, using ChIP-seq. In Chapter 6, the mitochondrial DNA occupancy of the TFB2M transcriptional regulator, the MTERF termination factor, and the mitochondrial RNA and DNA polymerases is characterized. Chapter 7 consists of an investigation into the curious phenomenon of the physical association of nuclear transcription factors with mitochondrial DNA, based on the diverse collections of transcription factor ChIP-seq datasets generated by the ENCODE, mouseENCODE and modENCODE consortia. In Chapter 8 this line of research is further extended to existing publicly available ChIP-seq datasets in plants and their mitochondrial and plastid genomes.

The third part is dedicated to the analytical and experimental practice of ChIP-seq. As part of the ENCODE Project, a set of metrics for assessing the quality of ChIP-seq experiments was developed, and the results of this activity are presented in Chapter 9. These metrics were later used to carry out a global analysis of ChIP-seq quality in the published literature (Chapter 10). In Chapter 11, the development and initial application of an automated robotic ChIP-seq (in which these metrics also played a major role) is presented.

The fourth part presents the results of some additional projects the author has been involved in, including the study of the role of the Piwi protein in the transcriptional regulation of transposon expression in Drosophila (Chapter 12), and the use of single-cell RNA-seq to characterize the heterogeneity of gene expression during cellular reprogramming (Chapter 13).

The last part of the thesis provides a review of the results of the ENCODE Project and the interpretation of the complexity of the biochemical activity exhibited by mammalian genomes that they have revealed (Chapters 15 and 16), an overview of the expected in the near future technical developments and their impact on the field of functional genomics (Chapter 14), and a discussion of some so far insufficiently explored research areas, the future study of which will, in the opinion of the author, provide deep insights into many fundamental but not yet completely answered questions about the transcriptional biology of eukaryotes and its regulation.

Cell-targeted regulation of gene expression through synthetic RNA devices

The ability to interface with and program cellular function remains a challenging research frontier in biotechnology. Although the emerging field of synthetic biology has recently generated a variety of gene-regulatory strategies based on synthetic RNA molecules, few strategies exist through which to control such regulatory effects in response to specific exogenous or endogenous molecular signals. Here, we present the development of an engineered RNA-based device platform to detect and act on endogenous protein signals, linking these signals to the regulation of genes and thus cellular function.

We describe efforts to develop an RNA-based device framework for regulating endogenous genes in human cells. Previously developed RNA control devices have demonstrated programmable ligand-responsive genetic regulation in diverse cell types, and we attempted to adapt this class of cis-acting control elements to function in trans. We divided the device into two strands that reconstitute activity upon hybridization. Device function was optimized using an in vivo model system, and we found that device sequence is not as flexible as previously reported. After verifying the in vitro activity of our optimized design, we attempted to establish gene regulation in a human cell line using additional elements to direct device stability, structure, and localization. The significant limitations of our platform prevented endogenous gene regulation.

We next describe the development of a protein-responsive RNA-based regulatory platform. Employing various design strategies, we demonstrated functional devices that both up- and downregulate gene expression in response to a heterologous protein in a human cell line. The activity of our platform exceeded that of a similar, small-molecule-responsive platform. We demonstrated the ability of our devices to respond to both cytoplasmic- and nuclear-localized protein, providing insight into the mechanism of action and distinguishing our platform from previously described devices with more restrictive ligand localization requirements. Finally, we demonstrated the versatility of our device platform by developing a regulatory device that responds to an endogenous signaling protein.

The foundational tool we present here possesses unique advantages over previously described RNA-based gene-regulatory platforms. This genetically encoded technology may find future applications in the development of more effective diagnostic tools and targeted molecular therapy strategies.

Studies on 4S and 5S RNA of HeLa cells

SECTION I

Section I is concerned with a partial sequence analysis conducted on 5S RNA from HeLa cells. Analysis of the oligonucleotide pattern after pancreatic ribonuclease digestion of a highly-purified preparation of 5S RNA gave results which were in general agreement with those published for KB cells, both with respect to the identity and the frequency of the partial sequences. However, the presence of a trinucleotide not found in the KB 5S pattern, together with the reproducibly much lower than expected molar yield of the larger oligonucleotides strongly suggested the occurrence of alternate sequences at various sites in the 5S molecules of human cells. The presence of ppGp and pppGp at the 5'-terminus of HeLa 5S RNA was clearly demonstrated. The implications of this finding with regard to the origin of 5S RNA are discussed.

SECTION II

In Section II the proportion of the HeLa cell genome complementary to tRNA was investigated by using RNA- DNA hybridization. The value for saturation of the HeLa DNA by tRNA was found to be 1.1 x 10^-5, which corresponds to about 4900 sites for tRNA per HeLa cell in an exponentially growing culture. Analysis of the nucleotide composition of the hybridized tRNA revealed significant differences from the nucleotide composition of the input tRNA, with the purine to pyrimidine ratio indicating, however, that these differences were not produced by excessive RNase attack of the hybrid. The size of the hybridized tRNA was only moderately smaller than that of the input RNA; the average S value in formaldehyde was 2.7 (corresponding to a length of about 65 nucleotides), suggesting that a relatively small portion near the ends of the hybridized 4S chains had been removed by RNase.

SECTION III

The proportion of the HeLa cell genome complementary to 5S RNA was investigated by using RNA-DNA hybridization. The value for saturation of the HeLa DNA by 5S RNA was found to be 2.3 x 10^-5, which corresponds to about 7,000 sites for 5S RNA per HeLa cell in an exponentially growing culture. Analysis of the nucleotide composition of the hybridized 5S RNA revealed no significant difference from the nucleotide composition of the input RNA. At the RNA to DNA input ratio of 1:1000, the average S value in formaldehyde of the hybridized 5S RNA corresponded to a polynucleotide chain about two-thirds the size of the input RNA.

I. Proton magnetic resonance of polynucleotides and transfer RNA. II. Electron spin relaxation studies of manganese (II) complexes in acetonitrile

Part I. Proton Magnetic Resonance of Polynucleotides and Transfer RNA.

Proton magnetic resonance was used to follow the temperature dependent intramolecular stacking of the bases in the polynucleotides of adenine and cytosine. Analysis of the results on the basis of a two state stacked-unstacked model yielded values of -4.5 kcal/mole and -9.5 kcal/mole for the enthalpies of stacking in polyadenylic and polycytidylic acid, respectively.

The interaction of purine with these molecules was also studied by pmr. Analysis of these results and the comparison of the thermal unstacking of polynucleotides and short chain nucleotides indicates that the bases contained in stacks within the long chain poly nucleotides are, on the average, closer together than the bases contained in stacks in the short chain nucleotides.

Temperature and purine studies were also carried out with an aqueous solution of formylmethionine transfer ribonucleic acid. Comparison of these results with the results of similar experiments with the homopolynucleotides of adenine, cytosine and uracil indicate that the purine is probably intercalating into loop regions of the molecule.

The solvent denaturation of phenylalanine transfer ribonucleic acid was followed by pmr. In a solvent mixture containing 83 volume per cent dimethylsulf oxide and 17 per cent deuterium oxide, the tRNA molecule is rendered quite flexible. It is possible to resolve resonances of protons on the common bases and on certain modified bases.

Part II. Electron Spin Relaxation Studies of Manganese (II) Complexes in Acetonitrile.

The electron paramagnetic resonance spectra of three Mn⁺² complexes, [Mn(CH₃CN)₆]⁺², [MnCl₄]^-2, and [MnBr₄]^-2, in acetonitrile were studied in detail. The objective of this study was to relate changes in the effective spin Hamiltonian parameters and the resonance line widths to the structure of these molecular complexes as well as to dynamical processes in solution.

Of the three systems studied, the results obtained from the [Mn(CH₃CN)₆]⁺² system were the most straight-forward to interpret. Resonance broadening attributable to manganese spin-spin dipolar interactions was observed as the manganese concentration was increased.

In the [MnCl₄]^-2 system, solvent fluctuations and dynamical ion-pairing appear to be significant in determining electron spin relaxation.

In the [MnBr₄]^-2 system, solvent fluctuations, ion-pairing, and Br- ligand exchange provide the principal means of electron spin relaxation. It was also found that the spin relaxation in this system is dependent upon the field strength and is directly related to the manganese concentration. A relaxation theory based on a two state collisional model was developed to account for the observed behavior.

Part I. Regulation of development in the cellular slime mold Dictyostelium discoideum. Part II. Polysomes and RNA synthesis during early development of the surf clam Spisula solidissima

Part I. The cellular slime mold Dictyostelium discoideum is a simple eukaryote which undergoes a multi-cellular developmental process. Single cell myxamoebae divide vegetatively in the presence of a food source. When the food is depleted or removed, the cells aggregate, forming a migrating pseudoplasmodium which differentiates into a fruiting body containing stalk and spore cells. I have shown that during the developmental cycle glycogen phosphorylase, aminopeptidase, and alanine transaminase are developmentally regulated, that is their specific activities increased at a specific time in the developmental cycle. Phosphorylase activity is undetectable in developing cells until mid-aggregation whereupon it increases and reaches a maximum at mid-culmination. Thereafter the enzyme disappears. Actinomycin D and cycloheximide studies as well as studies with morphologically aberrant and temporally deranged mutants indicate that prior RNA and concomitant protein synthesis are necessary for the rise and decrease in activity and support the view that the appearance of the enzyme is regulated at the transcriptional level. Aminopeptidase and alanine transaminase increase 3 fold starting at starvation and reach maximum activity at 18 and 5 hours respectively.

The cellular DNA s of D. discoideum were characterized by CsC1 buoyant density gradient centrifugation and by renaturation kinetics. Whole cell DNA exhibits three bands in CsCl: ρ = 1.676 g/cc (nuclear main band), 1.687 (nuclear satellite), and 1.682 (mitochondrial). Reassociation kinetics at a criterion of Tm -23°C indicates that the nuclear reiterated sequences make up 30% of the genome (Cot_1/2 (pure) 0.28) and the single-copy DNA 70% (Cot_1/2(pure) 70). The complexity of the nuclear genome is 30 x 10⁹ daltons and that of the mitochondrial DNA is 35-40 x 10⁶ daltons (Cot_1/2 0.15). rRNA cistrons constitute 2.2% of nuclear DNA and have a ρ = 1.682.

RNA extracted from 4 stages during developmental cycle of Dictyostelium was hybridized with purified single-copy nuclear DNA. The hybrids had properties indicative of single-copy DNA-RNA hybrids. These studies indicate that there are, during development, qualitative and quantitative changes in the portion of the single-copy of the genome transcribed. Overall, 56% of the genome is represented by transcripts between the amoeba and mid-culmination stages. Some 19% are sequences which are represented at all stages while 37% of the genome consists of stage specific sequences.

Part II. RNA and protein synthesis and polysome formation were studied during early development of the surf clam Spisula solidissima embryos. The oocyte has a small number of polysomes and a low but measurable rate of protein synthesis (leucine-³H incorporation). After fertilization, there is a continual increase in the percentage of ribosomes sedimenting in the polysome region. Newly synthesized RNA (uridine-5-³H incorporation) was found in polysomes as early as the 2-cell stage. During cleavage, the newly formed RNA is associated mainly with the light polysomes.

RNA extracted from polysomes labeled at the 4-cell stage is polydisperse, nonribosomal, and non-4 S. Actinomycin D causes a reduction of about 30% of the polysomes formed between fertilization and the 16-cell stage.

In the early cleavage stages the light polysomes are mostly affected by actinomycin.

Proton magnetic resonance studies of ribonucleic acid complexes. I. Complexes of biological bases and oligonucleotides with RNA. II. Template recognition and the degeneracy of the genetic code

Part I. Complexes of Biological Bases and Oligonucleotides with RNA

The physical nature of complexes of several biological bases and oligonucleotides with single-stranded ribonucleic acids have been studied by high resolution proton magnetic resonance spectroscopy. The importance of various forces in the stabilization of these complexes is also discussed.

Previous work has shown that purine forms an intercalated complex with single-stranded nucleic acids. This complex formation led to severe and stereospecific broadening of the purine resonances. From the field dependence of the linewidths, T₁ measurements of the purine protons and nuclear Overhauser enhancement experiments, the mechanism for the line broadening was ascertained to be dipole-dipole interactions between the purine protons and the ribose protons of the nucleic acid.

The interactions of ethidium bromide (EB) with several RNA residues have been studied. EB forms vertically stacked aggregates with itself as well as with uridine, 3'-uridine monophosphate and 5'-uridine monophosphate and forms an intercalated complex with uridylyl (3' → 5') uridine and polyuridylic acid (poly U). The geometry of EB in the intercalated complex has also been determined.

The effect of chain length of oligo-A-nucleotides on their mode of interaction with poly U in D₂0 at neutral pD have also been studied. Below room temperatures, ApA and ApApA form a rigid triple-stranded complex involving a stoichiometry of one adenine to two uracil bases, presumably via specific adenine-uracil base pairing and cooperative base stacking of the adenine bases. While no evidence was obtained for the interaction of ApA with poly U above room temperature, ApApA exhibited complex formation of a 1:1 nature with poly U by forming Watson-Crick base pairs. The thermodynamics of these systems are discussed.

Part II. Template Recognition and the Degeneracy of the Genetic Code

The interaction of ApApG and poly U was studied as a model system for the codon-anticodon interaction of tRNA and mRNA in vivo. ApApG was shown to interact with poly U below ~20°C. The interaction was of a 1:1 nature which exhibited the Hoogsteen bonding scheme. The three bases of ApApG are in an anti conformation and the guanosine base appears to be in the lactim tautomeric form in the complex.

Due to the inadequacies of previous models for the degeneracy of the genetic code in explaining the observed interactions of ApApG with poly U, the "tautomeric doublet" model is proposed as a possible explanation of the degenerate interactions of tRNA with mRNA during protein synthesis in vivo.

Assembly of viral genomes from metagenomes

Viral infections remain a serious global health issue. Metagenomic approaches are increasingly used in the detection of novel viral pathogens but also to generate complete genomes of uncultivated viruses. In silico identification of complete viral genomes from sequence data would allow rapid phylogenetic characterization of these new viruses. Often, however, complete viral genomes are not recovered, but rather several distinct contigs derived from a single entity are, some of which have no sequence homology to any known proteins. De novo assembly of single viruses from a metagenome is challenging, not only because of the lack of a reference genome, but also because of intrapopulation variation and uneven or insufficient coverage. Here we explored different assembly algorithms, remote homology searches, genome-specific sequence motifs, k-mer frequency ranking, and coverage profile binning to detect and obtain viral target genomes from metagenomes. All methods were tested on 454-generated sequencing datasets containing three recently described RNA viruses with a relatively large genome which were divergent to previously known viruses from the viral families Rhabdoviridae and Coronaviridae. Depending on specific characteristics of the target virus and the metagenomic community, different assembly and in silico gap closure strategies were successful in obtaining near complete viral genomes.

The promoter of cell growth- and RNA protection-associated SND1 gene is activated by endoplasmic reticulum stress in human hepatoma cells

Background: Staphyloccocal nuclease domain-containing protein 1 (SND1) is involved in the regulation of gene expression and RNA protection. While numerous studies have established that SND1 protein expression is modulated by cellular stresses associated with tumor growth, hypoxia, inflammation, heat- shock and oxidative conditions, little is known about the factors responsible for SND1 expression. Here, we have approached this question by analyzing the transcriptional response of human SND1 gene to pharmacological endoplasmic reticulum (ER) stress in liver cancer cells. Results: We provide first evidence that SND1 promoter activity is increased in human liver cancer cells upon exposure to thapsigargin or tunicamycin or by ectopic expression of ATF6, a crucial transcription factor in the unfolded protein response triggered by ER stress. Deletion analysis of the 5'-flanking region of SND1 promoter identified maximal activation in fragment (-934, +221), which contains most of the predicted ER stress response elements in proximal promoter. Quantitative real- time PCR revealed a near 3 fold increase in SND1 mRNA expression by either of the stress- inducers; whereas SND1 protein was maximally upregulated (3.4-fold) in cells exposed to tunicamycin, a protein glycosylation inhibitor. Conclusion: Promoter activity of the cell growth- and RNA-protection associated SND1 gene is up-regulated by ER stress in human hepatoma cells.

Development of nucleic acid vaccines: use of self-amplifying RNA in lipid nanoparticles

Self-amplifying RNA or RNA replicon is a form of nucleic acid-based vaccine derived from either positive-strand or negative-strand RNA viruses. The gene sequences encoding structural proteins in these RNA viruses are replaced by mRNA encoding antigens of interest as well as by RNA polymerase for replication and transcription. This kind of vaccine has been successfully assayed with many different antigens as vaccines candidates, and has been shown to be potent in several animal species, including mice, nonhuman primates, and humans. A key challenge to realizing the broad potential of self-amplifying vaccines is the need for safe and effective delivery methods. Ideally, an RNA nanocarrier should provide protection from blood nucleases and extended blood circulation, which ultimately would increase the possibility of reaching the target tissue. The delivery system must then be internalized by the target cell and, upon receptor-mediated endocytosis, must be able to escape from the endosomal compartment into the cell cytoplasm, where the RNA machinery is located, while avoiding degradation by lysosomal enzymes. Further, delivery systems for systemic administration ought to be well tolerated upon administration. They should be safe, enabling the multiadministration treatment modalities required for improved clinical outcomes and, from a developmental point of view, production of large batches with reproducible specifications is also desirable. In this review, the concept of self-amplifying RNA vaccines and the most promising lipid-based delivery systems are discussed.

Identification of Sex and Female’s Reproductive Stage in Commercial Fish Species through the Quantification of Ribosomal Transcripts in Gonads

The estimation of maturity and sex of fish stocks in European waters is a requirement of the EU Data Collection Framework as part of the policy to improve fisheries management. On the other hand, research on fish biology is increasingly focused in molecular approaches, researchers needing correct identification of fish sex and reproductive stage without necessarily having in house the histological know-how necessary for the task. Taking advantage of the differential gene transcription occurring during fish sex differentiation and gametogenesis, the utility of 5S ribosomal RNA (5S rRNA) and General transcription factor IIIA (gtf3a) in the molecular identification of sex and gametogenic stage was tested in different economically-relevant fish species from the Bay of Biscay. Gonads of 9 fish species (, Atlantic, Atlantic-chub and horse mackerel, blue whiting, bogue, European anchovy, hake and pilchard and megrim), collected from local commercial fishing vessels were histologically sexed and 5S and 18S rRNA concentrations were quantified by capillary electrophoresis to calculate a 5S/18S rRNA index. Degenerate primers permitted cloning and sequencing of gtf3a fragments in 7 of the studied species. 5S rRNA and gtf3a transcript levels, together with 5S/18S rRNA index, distinguished clearly ovaries from testis in all of the studied species. The values were always higher in females than in males. 5S/18S rRNA index values in females were always highest when fish were captured in early phases of ovary development whilst, in later vitellogenic stages, the values decreased significantly. In megrim and European anchovy, where gonads in different oogenesis stages were obtained, the 5S/18S rRNA index identified clearly gametogenic stage. This approach, to the sexing and the quantitative non-subjective identification of the maturity stage of female fish, could have multiple applications in the study of fish stock dynamics, fish reproduction and fecundity and fish biology in general.