Personal omics profiling reveals dynamic molecular and medical phenotypes.

Personalized medicine is expected to benefit from combining genomic information with regular monitoring of physiological states by multiple high-throughput methods. Here, we present an integrative personal omics profile (iPOP), an analysis that combines genomic, transcriptomic, proteomic, metabolomic, and autoantibody profiles from a single individual over a 14 month period. Our iPOP analysis revealed various medical risks, including type 2 diabetes. It also uncovered extensive, dynamic changes in diverse molecular components and biological pathways across healthy and diseased conditions. Extremely high-coverage genomic and transcriptomic data, which provide the basis of our iPOP, revealed extensive heteroallelic changes during healthy and diseased states and an unexpected RNA editing mechanism. This study demonstrates that longitudinal iPOP can be used to interpret healthy and diseased states by connecting genomic information with additional dynamic omics activity.

Spatiotemporal DNA methylome dynamics of the developing mouse fetus.

Cytosine DNA methylation is essential for mammalian development but understanding of its spatiotemporal distribution in the developing embryo remains limited1,2. Here, as part of the mouse Encyclopedia of DNA Elements (ENCODE) project, we profiled 168 methylomes from 12 mouse tissues or organs at 9 developmental stages from embryogenesis to adulthood. We identified 1,808,810 genomic regions that showed variations in CG methylation by comparing the methylomes of different tissues or organs from different developmental stages. These DNA elements predominantly lose CG methylation during fetal development, whereas the trend is reversed after birth. During late stages of fetal development, non-CG methylation accumulated within the bodies of key developmental transcription factor genes, coinciding with their transcriptional repression. Integration of genome-wide DNA methylation, histone modification and chromatin accessibility data enabled us to predict 461,141 putative developmental tissue-specific enhancers, the human orthologues of which were enriched for disease-associated genetic variants. These spatiotemporal epigenome maps provide a resource for studies of gene regulation during tissue or organ progression, and a starting point for investigating regulatory elements that are involved in human developmental disorders.

Human body epigenome maps reveal noncanonical DNA methylation variation.

Understanding the diversity of human tissues is fundamental to disease and requires linking genetic information, which is identical in most of an individual's cells, with epigenetic mechanisms that could have tissue-specific roles. Surveys of DNA methylation in human tissues have established a complex landscape including both tissue-specific and invariant methylation patterns. Here we report high coverage methylomes that catalogue cytosine methylation in all contexts for the major human organ systems, integrated with matched transcriptomes and genomic sequence. By combining these diverse data types with each individuals' phased genome, we identified widespread tissue-specific differential CG methylation (mCG), partially methylated domains, allele-specific methylation and transcription, and the unexpected presence of non-CG methylation (mCH) in almost all human tissues. mCH correlated with tissue-specific functions, and using this mark, we made novel predictions of genes that escape X-chromosome inactivation in specific tissues. Overall, DNA methylation in several genomic contexts varies substantially among human tissues.

Integrative analysis of haplotype-resolved epigenomes across human tissues.

Allelic differences between the two homologous chromosomes can affect the propensity of inheritance in humans; however, the extent of such differences in the human genome has yet to be fully explored. Here we delineate allelic chromatin modifications and transcriptomes among a broad set of human tissues, enabled by a chromosome-spanning haplotype reconstruction strategy. The resulting large collection of haplotype-resolved epigenomic maps reveals extensive allelic biases in both chromatin state and transcription, which show considerable variation across tissues and between individuals, and allow us to investigate cis-regulatory relationships between genes and their control sequences. Analyses of histone modification maps also uncover intriguing characteristics of cis-regulatory elements and tissue-restricted activities of repetitive elements. The rich data sets described here will enhance our understanding of the mechanisms by which cis-regulatory elements control gene expression programs.

Abnormalities in human pluripotent cells due to reprogramming mechanisms.

Human pluripotent stem cells hold potential for regenerative medicine, but available cell types have significant limitations. Although embryonic stem cells (ES cells) from in vitro fertilized embryos (IVF ES cells) represent the 'gold standard', they are allogeneic to patients. Autologous induced pluripotent stem cells (iPS cells) are prone to epigenetic and transcriptional aberrations. To determine whether such abnormalities are intrinsic to somatic cell reprogramming or secondary to the reprogramming method, genetically matched sets of human IVF ES cells, iPS cells and nuclear transfer ES cells (NT ES cells) derived by somatic cell nuclear transfer (SCNT) were subjected to genome-wide analyses. Both NT ES cells and iPS cells derived from the same somatic cells contained comparable numbers of de novo copy number variations. In contrast, DNA methylation and transcriptome profiles of NT ES cells corresponded closely to those of IVF ES cells, whereas iPS cells differed and retained residual DNA methylation patterns typical of parental somatic cells. Thus, human somatic cells can be faithfully reprogrammed to pluripotency by SCNT and are therefore ideal for cell replacement therapies.

Architecture of the human regulatory network derived from ENCODE data.

Transcription factors bind in a combinatorial fashion to specify the on-and-off states of genes; the ensemble of these binding events forms a regulatory network, constituting the wiring diagram for a cell. To examine the principles of the human transcriptional regulatory network, we determined the genomic binding information of 119 transcription-related factors in over 450 distinct experiments. We found the combinatorial, co-association of transcription factors to be highly context specific: distinct combinations of factors bind at specific genomic locations. In particular, there are significant differences in the binding proximal and distal to genes. We organized all the transcription factor binding into a hierarchy and integrated it with other genomic information (for example, microRNA regulation), forming a dense meta-network. Factors at different levels have different properties; for instance, top-level transcription factors more strongly influence expression and middle-level ones co-regulate targets to mitigate information-flow bottlenecks. Moreover, these co-regulations give rise to many enriched network motifs (for example, noise-buffering feed-forward loops). Finally, more connected network components are under stronger selection and exhibit a greater degree of allele-specific activity (that is, differential binding to the two parental alleles). The regulatory information obtained in this study will be crucial for interpreting personal genome sequences and understanding basic principles of human biology and disease.

Annotation of functional variation in personal genomes using RegulomeDB.

As the sequencing of healthy and disease genomes becomes more commonplace, detailed annotation provides interpretation for individual variation responsible for normal and disease phenotypes. Current approaches focus on direct changes in protein coding genes, particularly nonsynonymous mutations that directly affect the gene product. However, most individual variation occurs outside of genes and, indeed, most markers generated from genome-wide association studies (GWAS) identify variants outside of coding segments. Identification of potential regulatory changes that perturb these sites will lead to a better localization of truly functional variants and interpretation of their effects. We have developed a novel approach and database, RegulomeDB, which guides interpretation of regulatory variants in the human genome. RegulomeDB includes high-throughput, experimental data sets from ENCODE and other sources, as well as computational predictions and manual annotations to identify putative regulatory potential and identify functional variants. These data sources are combined into a powerful tool that scores variants to help separate functional variants from a large pool and provides a small set of putative sites with testable hypotheses as to their function. We demonstrate the applicability of this tool to the annotation of noncoding variants from 69 full sequenced genomes as well as that of a personal genome, where thousands of functionally associated variants were identified. Moreover, we demonstrate a GWAS where the database is able to quickly identify the known associated functional variant and provide a hypothesis as to its function. Overall, we expect this approach and resource to be valuable for the annotation of human genome sequences.

Tcf7 is an important regulator of the switch of self-renewal and differentiation in a multipotential hematopoietic cell line.

A critical problem in biology is understanding how cells choose between self-renewal and differentiation. To generate a comprehensive view of the mechanisms controlling early hematopoietic precursor self-renewal and differentiation, we used systems-based approaches and murine EML multipotential hematopoietic precursor cells as a primary model. EML cells give rise to a mixture of self-renewing Lin-SCA+CD34+ cells and partially differentiated non-renewing Lin-SCA-CD34- cells in a cell autonomous fashion. We identified and validated the HMG box protein TCF7 as a regulator in this self-renewal/differentiation switch that operates in the absence of autocrine Wnt signaling. We found that Tcf7 is the most down-regulated transcription factor when CD34+ cells switch into CD34- cells, using RNA-Seq. We subsequently identified the target genes bound by TCF7, using ChIP-Seq. We show that TCF7 and RUNX1 (AML1) bind to each other's promoter regions and that TCF7 is necessary for the production of the short isoforms, but not the long isoforms of RUNX1, suggesting that TCF7 and the short isoforms of RUNX1 function coordinately in regulation. Tcf7 knock-down experiments and Gene Set Enrichment Analyses suggest that TCF7 plays a dual role in promoting the expression of genes characteristic of self-renewing CD34+ cells while repressing genes activated in partially differentiated CD34- state. Finally a network of up-regulated transcription factors of CD34+ cells was constructed. Factors that control hematopoietic stem cell (HSC) establishment and development, cell growth, and multipotency were identified. These studies in EML cells demonstrate fundamental cell-intrinsic properties of the switch between self-renewal and differentiation, and yield valuable insights for manipulating HSCs and other differentiating systems.

High-throughput sequencing for biology and medicine.

Advances in genome sequencing have progressed at a rapid pace, with increased throughput accompanied by plunging costs. But these advances go far beyond faster and cheaper. High-throughput sequencing technologies are now routinely being applied to a wide range of important topics in biology and medicine, often allowing researchers to address important biological questions that were not possible before. In this review, we discuss these innovative new approaches-including ever finer analyses of transcriptome dynamics, genome structure and genomic variation-and provide an overview of the new insights into complex biological systems catalyzed by these technologies. We also assess the impact of genotyping, genome sequencing and personal omics profiling on medical applications, including diagnosis and disease monitoring. Finally, we review recent developments in single-cell sequencing, and conclude with a discussion of possible future advances and obstacles for sequencing in biology and health.

Human Immune Cell Epigenomic Signatures in Response to Infectious Diseases and Chemical Exposures.

Variations in DNA methylation patterns in human tissues have been linked to various environmental exposures and infections. Here, we identified the DNA methylation signatures associated with multiple exposures in nine major immune cell types derived from peripheral blood mononuclear cells (PBMCs) at single-cell resolution. We performed methylome sequencing on 111,180 immune cells obtained from 112 individuals who were exposed to different viruses, bacteria, or chemicals. Our analysis revealed 790,662 differentially methylated regions (DMRs) associated with these exposures, which are mostly individual CpG sites. Additionally, we integrated methylation and ATAC-seq data from same samples and found strong correlations between the two modalities. However, the epigenomic remodeling in these two modalities are complementary. Finally, we identified the minimum set of DMRs that can predict exposures. Overall, our study provides the first comprehensive dataset of single immune cell methylation profiles, along with unique methylation biomarkers for various biological and chemical exposures.

Consensus miRNA expression profiles derived from interplatform normalization of microarray data.

Eukaryotic gene expression is controlled at the post-transcriptional level by small noncoding RNAs called microRNAs (miRNA). miRNAs play important roles during early development and participate in gene regulatory circuits in the cell. Different high-throughput expression analysis methods including microarrays, bead-based detection, and small RNA cloning have been applied to quantitatively detect miRNAs in various tissues, cell types, and biological conditions. High-throughput expression data was collected from public repositories and processed to create a database of miRNA expression profiles. Several commonly used normalization methods were compared to identify suitable methods for cross-platform comparison of high-throughput miRNA expression data. The database provides interlaboratory and interplatform validated reference expression levels for miRNAs. The normalized expression profiles were validated by querying for well-established features of miRNA expression. Firstly, expression profiles of several tissue-specific miRNAs showed good agreement between the database and previously reported profiles. We have also identified a set of miRNAs that are constitutively expressed across mammalian tissues. Secondly, we used the database to compare the expression patterns of miRNAs belonging to the let-7 family, where the divergence in expression patterns implies that they may have diversified functionally. Lastly, we compared expression profiles of intronic and clustered miRNAs. Expression profiles of intronic miRNAs and clustered miRNAs showed either very good, or in certain cases, very poor correlation with the host gene. Interplatform comparison of miRNA expression profiles thus provides a resource of consensus expression profiles that can be used in the future for studying miRNA function and regulation.

Posttranscriptional regulation of interleukin-10 expression by hsa-miR-106a.

IL-10 is a key regulator of the immune system that critically determines health and disease. Its expression is finely tuned both at the transcriptional and posttranscriptional levels. Although the importance of posttranscriptional regulation of IL-10 has been previously shown, understanding the underlying mechanisms is still in its infancy. In this study, using a combination of bioinformatics and molecular approaches, we report that microRNA (hsa-miR-106a) regulates IL-10 expression. The hsa-miR-106a binding site in the 3' UTR of IL10 has been identified by site-directed mutagenesis studies. Also, the involvement of transcription factors, Sp1 and Egr1, in the regulation of hsa-miR-106a expression and concomitant decrease in the IL-10 expression, has also been demonstrated. In summary, our results showed that IL-10 expression may be regulated by miR-106a, which is in turn transcriptionally regulated by Egr1 and Sp1.

Inhibition of neuronal cholesterol biosynthesis with lovastatin leads to impaired synaptic vesicle release even in the presence of lipoproteins or geranylgeraniol.

Cholesterol is highly enriched in the brain, and plays a key role in synapse formation and function. The brain does not derive cholesterol from the circulation; instead, the majority of cholesterol is made in glia and secreted in form of lipoproteins. Neurons can synthesize cholesterol, but the extent of neuronal cholesterol biosynthesis in the adult brain is unknown. Cholesterol biosynthesis inhibitors of the statin family are widely used to lower circulating cholesterol and cardiovascular risk. Lipophilic statins can cross the blood brain barrier and inhibit brain cholesterol biosynthesis with possible consequences for synaptic cholesterol homeostasis. We have investigated the effects of lovastatin on synapse maturation and synaptic vesicle release. Treatment of primary hippocampal neurons with low levels of lovastatin for one week reduced synapse density and impaired synaptic vesicle release. Neither lipoproteins nor geranylgeraniol fully counteracted the lovastatin-induced decrease of synaptic vesicle exocytosis, even when cholesterol depletion was prevented. In contrast, restoration of neuronal cholesterol synthesis with mevalonate prevented defects in vesicle exocytosis without fully normalizing neuronal cholesterol content. These results raise the possibility that chronic exposure of neurons to lipophilic statins may affect synaptic transmission, and indicate that hippocampal neurons need a certain level of endogenous cholesterol biosynthesis.

Incomplete penetrance and variable expressivity: is there a microRNA connection?

Incomplete penetrance and variable expressivity are non-Mendelian phenomena resulting in the lack of correlation between genotype and phenotype. Not withstanding the diversity in mechanisms, differential expression of homologous alleles within cells manifests as variations in penetrance and expressivity of mutations between individuals of the same genotype. These phenomena are seen most often in dominantly inherited diseases, implying that they are sensitive to concentration of the gene product. In this framework and the advances in understanding the role of microRNA (miRNA) in fine-tuning gene expression at translational level, we propose miRNA-mediated regulation as a mechanism for incomplete penetrance and variable expressivity. The presence of miRNA binding sites at 3' UTR, co-expression of target gene-miRNA pairs for genes showing incomplete penetrance and variable expressivity derived from available data lend support to our hypothesis. Single nucleotide polymorphisms in the miRNA target site facilitate the implied differential targeting of the transcripts from homologous alleles.

Inhibition of translation in living eukaryotic cells by an RNA G-quadruplex motif.

Guanine-rich sequences can adopt intramolecular four-stranded structures, called G-quadruplexes. These motifs have been intensively investigated on the DNA level, but their overall biological relevance remains elusive. Only recently has research concerning the function of G-quadruplexes in RNAs commenced. Here, we demonstrate for the first time, that an RNA G-quadruplex structure inhibits translation in vivo in eukaryotic cells. We investigated the function of a highly conserved, thermodynamically stable RNA G-quadruplex in the 5'-UTR of the mRNA of the human Zic-1 zinc-finger protein. Using dual luciferase reporter assay, we demonstrate that the Zic-1 RNA G-quadruplex represses protein synthesis inside eukaryotic cells. Quantitative RT-PCR assays confirmed that the reduction of protein synthesis is due to regulation of the translation process and not a consequence of reduced transcription. Western blot analysis revealed that expression of Zic-1 is strongly reduced by a 73 nucleotides-long fragment of the UTR containing the G-quadruplex motif. These structures might add to the more recently discovered elements in untranslated regions of mRNAs that regulate their translation.

dbSMR: a novel resource of genome-wide SNPs affecting microRNA mediated regulation.

BACKGROUND: MicroRNAs (miRNAs) regulate several biological processes through post-transcriptional gene silencing. The efficiency of binding of miRNAs to target transcripts depends on the sequence as well as intramolecular structure of the transcript. Single Nucleotide Polymorphisms (SNPs) can contribute to alterations in the structure of regions flanking them, thereby influencing the accessibility for miRNA binding. DESCRIPTION: The entire human genome was analyzed for SNPs in and around predicted miRNA target sites. Polymorphisms within 200 nucleotides that could alter the intramolecular structure at the target site, thereby altering regulation were annotated. Collated information was ported in a MySQL database with a user-friendly interface accessible through the URL: (http://miracle.igib.res.in/dbSMR). CONCLUSION: The database has a user-friendly interface where the information can be queried using either the gene name, microRNA name, polymorphism ID or transcript ID. Combination queries using 'AND' or 'OR' is also possible along with specifying the degree of change of intramolecular bonding with and without the polymorphism. Such a resource would enable researchers address questions like the role of regulatory SNPs in the 3' UTRs and population specific regulatory modulations in the context of microRNA targets.

Human cellular microRNA hsa-miR-29a interferes with viral nef protein expression and HIV-1 replication.

BACKGROUND: Cellular miRNAs play an important role in the regulation of gene expression in eukaryotes. Recently, miRNAs have also been shown to be able to target and inhibit viral gene expression. Computational predictions revealed earlier that the HIV-1 genome includes regions that may be potentially targeted by human miRNAs. Here we report the functionality of predicted miR-29a target site in the HIV-1 nef gene. RESULTS: We find that the human miRNAs hsa-miR-29a and 29b are expressed in human peripheral blood mononuclear cells. Expression of a luciferase reporter bearing the nef miR-29a target site was decreased compared to the luciferase construct without the target site. Locked nucleic acid modified anti-miRNAs targeted against hsa-miR-29a and 29b specifically reversed the inhibitory effect mediated by cellular miRNAs on the target site. Ectopic expression of the miRNA results in repression of the target Nef protein and reduction of virus levels. CONCLUSION: Our results show that the cellular miRNA hsa-miR29a downregulates the expression of Nef protein and interferes with HIV-1 replication.

Quadfinder: server for identification and analysis of quadruplex-forming motifs in nucleotide sequences.

G-quadruplex secondary structures, which play a structural role in repetitive DNA such as telomeres, may also play a functional role at other genomic locations as targetable regulatory elements which control gene expression. The recent interest in application of quadruplexes in biological systems prompted us to develop a tool for the identification and analysis of quadruplex-forming nucleotide sequences especially in the RNA. Here we present Quadfinder, an online server for prediction and bioinformatics of uni-molecular quadruplex-forming nucleotide sequences. The server is designed to be user-friendly and needs minimal intervention by the user, while providing flexibility of defining the variants of the motif. The server is freely available at URL http://miracle.igib.res.in/quadfinder/.

Functional Human Oocytes Generated by Transfer of Polar Body Genomes.

Oocyte defects lie at the heart of some forms of infertility and could potentially be addressed therapeutically by alternative routes for oocyte formation. Here, we describe the generation of functional human oocytes following nuclear transfer of first polar body (PB1) genomes from metaphase II (MII) oocytes into enucleated donor MII cytoplasm (PBNT). The reconstructed oocytes supported the formation of de novo meiotic spindles and, after fertilization with sperm, meiosis completion and formation of normal diploid zygotes. While PBNT zygotes developed to blastocysts less frequently (42%) than controls (75%), genome-wide genetic, epigenetic, and transcriptional analyses of PBNT and control ESCs indicated comparable numbers of structural variations and markedly similar DNA methylation and transcriptome profiles. We conclude that rescue of PB1 genetic material via introduction into donor cytoplasm may offer a source of oocytes for infertility treatment or mitochondrial replacement therapy for mtDNA disease.