Ashok Patowary

Detection of altered global DNA methylation in coronary artery disease patients.

Epigenetic modifications, especially alteration in DNA methylation, are increasingly being recognized as a key factor in the pathogenesis of complex disorders, including atherosclerosis. However, there are limited data on the epigenetic changes in the coronary artery disease (CAD) patients. In the present study we evaluated the methylation status of genomic DNA from peripheral lymphocytes in a cohort of 287 individuals: 137 angiographically confirmed CAD patients and 150 controls. The differential susceptibility of genomic DNA to methylation-sensitive restriction enzymes was utilized to assess the methylation status of the genome. We observed that the genomic DNA methylation in CAD patients is significantly higher than in controls (p < 0.05). Since elevated homocysteine levels are known to be an independent risk factor for CAD and a key modulator of macromolecular methylation, we investigated the probable correlation between plasma homocysteine levels and global DNA methylation. We observed a significant positive correlation of global DNA methylation with plasma homocysteine levels in CAD patients (p = 0.001). Further, within a higher range of serum homocysteine levels (>/=12-50 muM), global DNA methylation was significantly higher in CAD patients than in controls. The alteration in genomic DNA methylation associated with cardiovascular disease per se appears to be further accentuated by higher homocysteine levels.

High resolution methylome map of rat indicates role of intragenic DNA methylation in identification of coding region.

DNA methylation is crucial for gene regulation and maintenance of genomic stability. Rat has been a key model system in understanding mammalian systemic physiology, however detailed rat methylome remains uncharacterized till date. Here, we present the first high resolution methylome of rat liver generated using Methylated DNA immunoprecipitation and high throughput sequencing (MeDIP-Seq) approach. We observed that within the DNA/RNA repeat elements, simple repeats harbor the highest degree of methylation. Promoter hypomethylation and exon hypermethylation were common features in both RefSeq genes and expressed genes (as evaluated by proteomic approach). We also found that although CpG islands were generally hypomethylated, about 6% of them were methylated and a large proportion (37%) of methylated islands fell within the exons. Notably, we obeserved significant differences in methylation of terminal exons (UTRs); methylation being more pronounced in coding/partially coding exons compared to the non-coding exons. Further, events like alternate exon splicing (cassette exon) and intron retentions were marked by DNA methylation and these regions are retained in the final transcript. Thus, we suggest that DNA methylation could play a crucial role in marking coding regions thereby regulating alternative splicing. Apart from generating the first high resolution methylome map of rat liver tissue, the present study provides several critical insights into methylome organization and extends our understanding of interplay between epigenome, gene expression and genome stability.

Dynamic expression of long non-coding RNAs (lncRNAs) in adult zebrafish.

Long non-coding RNAs (lncRNA) represent an assorted class of transcripts having little or no protein coding capacity and have recently gained importance for their function as regulators of gene expression. Molecular studies on lncRNA have uncovered multifaceted interactions with protein coding genes. It has been suggested that lncRNAs are an additional layer of regulatory switches involved in gene regulation during development and disease. LncRNAs expressing in specific tissues or cell types during adult stages can have potential roles in form, function, maintenance and repair of tissues and organs. We used RNA sequencing followed by computational analysis to identify tissue restricted lncRNA transcript signatures from five different tissues of adult zebrafish. The present study reports 442 predicted lncRNA transcripts from adult zebrafish tissues out of which 419 were novel lncRNA transcripts. Of these, 77 lncRNAs show predominant tissue restricted expression across the five major tissues investigated. Adult zebrafish brain expressed the largest number of tissue restricted lncRNA transcripts followed by cardiovascular tissue. We also validated the tissue restricted expression of a subset of lncRNAs using independent methods. Our data constitute a useful genomic resource towards understanding the expression of lncRNAs in various tissues in adult zebrafish. Our study is thus a starting point and opens a way towards discovering new molecular interactions of gene expression within the specific adult tissues in the context of maintenance of organ form and function.

miR-34 is maternally inherited in Drosophila melanogaster and Danio rerio

MicroRNAs (miRNAs) are small, endogenous, regulatory RNA molecules that can bind to partially complementary regions on target messenger RNAs and impede their expression or translation. We rationalized that miRNAs, being localized to the cytoplasm, will be maternally inherited during fertilization and may play a role in early development. Although Dicer is known to be essential for the transition from single-celled zygote to two-cell embryo, a direct role for miRNAs has not yet been demonstrated. We identified miRNAs with targets in zygotically expressed transcripts in Drosophila using a combination of transcriptome analysis and miRNA target prediction. We experimentally established that Drosophila miRNA dme-miR-34, the fly homologue of the cancer-related mammalian miRNA miR-34, involved in somatic-cell reprogramming and having critical role in early neuronal differentiation, is present in Drosophila embryos before initiation of zygotic transcription. We also show that the Drosophila miR-34 is dependent on maternal Dicer-1 for its expression in oocytes. Further, we show that miR-34 is also abundant in unfertilized oocytes of zebrafish. Its temporal expression profile during early development showed abundant expression in unfertilized oocytes that gradually decreased by 5 days post-fertilization (dpf). We find that knocking down the maternal, but not the zygotic, miR-34 led to developmental defects in the neuronal system during early embryonic development in zebrafish. Here, we report for the first time, the maternal inheritance of an miRNA involved in development of the neuronal system in a vertebrate model system.

Reverse genetics screen in zebrafish identifies a role of miR-142a-3p in vascular development and integrity.

MicroRNAs are a well-studied class of non-coding RNA and are known to regulate developmental processes in eukaryotes. Their role in key biological processes such as vasculature development has attracted interest. However, a comprehensive understanding of molecular regulation of angiogenesis and vascular integrity during development remains less explored. Here we identified miRNAs involved in the development and maintenance of vasculature in zebrafish embryos using a reverse genetics approach. Using a combination of bioinformatics predictions and literature based evidences we mined over 701 Human and 329 Zebrafish miRNAs to derive a list of 29 miRNAs targeting vascular specific genes in zebrafish. We shortlisted eight miRNAs and investigated their potential role in regulating vascular development in zebrafish transgenic model. In this screen we identified three miRNAs, namely miR-1, miR-144 and miR-142a-3p that have the potential to influence vascular development in zebrafish. We show that miR-142a-3p mediates vascular integrity and developmental angiogenesis in vivo. Overexpression of miR-142a-3p results in loss of vascular integrity, hemorrhage and vascular remodeling during zebrafish embryonic development, while loss of function of miR-142a-3p causes abnormal vascular remodeling. MiR-142a-3p functions in part by directly repressing cdh5 (VE-cadherin). The vascular abnormalities that results from modulation of miR-142a-3p are reminiscent of cdh5 perturbation in zebrafish embryos. We also demonstrate that the action of miR-142a on cdh5 is potentially regulated by Lmo2, an important transcription factor, known for its role in vasculature development. The miR142a-3p mediated control of cdh5 constitutes an additional layer of regulation for maintaining vascular integrity and developmental angiogenesis. These findings have implications in development, wound repair and tumor growth.

A Sequence-Based Variation Map of Zebrafish

Zebrafish (Danio rerio) is a popular vertebrate model organism largely deployed using outbred laboratory animals. The nonisogenic nature of the zebrafish as a model system offers the opportunity to understand natural variations and their effect in modulating phenotype. In an effort to better characterize the range of natural variation in this model system and to complement the zebrafish reference genome project, the whole genome sequence of a wild zebrafish at 39-fold genome coverage was determined. Comparative analysis with the zebrafish reference genome revealed approximately 5.2 million single nucleotide variations and over 1.6 million insertion-deletion variations. This dataset thus represents a new catalog of genetic variations in the zebrafish genome. Further analysis revealed selective enrichment for variations in genes involved in immune function and response to the environment, suggesting genome-level adaptations to environmental niches. We also show that human disease gene orthologs in the sequenced wild zebrafish genome show a lower ratio of nonsynonymous to synonymous single nucleotide variations.

Silencing c-MYC expression by targeting quadruplex in P1 promoter using locked nucleic acid trap.

The nuclease hypersensitive element of P1 promoter in c-MYC gene harbors a potential of unusual structure called quadruplex, which is involved in molecular recognition and function. This Hoogsteen bonded structure is in dynamic equilibrium with the usual Watson-Crick duplex structure, and these competing secondary structures undergo interconversion for execution of their respective biological roles. Herein, we investigate the sensitivity of the c-MYC quadruplex-duplex equilibrium by employing a locked nucleic acid (LNA) modified complementary strand as a pharmacological agent. Our biophysical experiments indicate that the c-MYC quadruplex under physiological conditions is stable and dominates the quadruplex-WC duplex equilibrium in both sodium and potassium buffers. This equilibrium is perturbed upon introducing the LNA modified complementary strand, which demonstrates efficient invasion of stable c-MYC quadruplex and duplex formation in contrast to the unmodified complementary strand. Our data indicate that LNA modifications confer increased thermodynamic stability to the duplex and thus favor the predominance of the duplex population over that of the quadruplex. Further, we demonstrate that this perturbation of equilibrium by a pharmacological agent results in altered gene expression. Our in vivo experiment performed using the LNA modified complementary strand suggests the influence of the quadruplex-duplex structural switch in the modulation of gene expression. We believe that this exploratory approach utilizing the selectivity and specificity of Watson-Crick base pairing of LNA bases would allow the modulation of quadruplex regulated gene expression.

Systematic Pharmacogenomics Analysis of a Malay Whole Genome: Proof of Concept for Personalized Medicine

Background With a higher throughput and lower cost in sequencing, second generation sequencing technology has immense potential for translation into clinical practice and in the realization of pharmacogenomics based patient care. The systematic analysis of whole genome sequences to assess patient to patient variability in pharmacokinetics and pharmacodynamics responses towards drugs would be the next step in future medicine in line with the vision of personalizing medicine. Methods Genomic DNA obtained from a 55 years old, self-declared healthy, anonymous male of Malay descent was sequenced. The subjects mother died of lung cancer and the father had a history of schizophrenia and deceased at the age of 65 years old. A systematic, intuitive computational workflow/pipeline integrating custom algorithm in tandem with large datasets of variant annotations and gene functions for genetic variations with pharmacogenomics impact was developed. A comprehensive pathway map of drug transport, metabolism and action was used as a template to map non-synonymous variations with potential functional consequences. Principal Findings Over 3 million known variations and 100,898 novel variations in the Malay genome were identified. Further in-depth pharmacogenetics analysis revealed a total of 607 unique variants in 563 proteins, with the eventual identification of 4 drug transport genes, 2 drug metabolizing enzyme genes and 33 target genes harboring deleterious SNVs involved in pharmacological pathways, which could have a potential role in clinical settings. Conclusions The current study successfully unravels the potential of personal genome sequencing in understanding the functionally relevant variations with potential influence on drug transport, metabolism and differential therapeutic outcomes. These will be essential for realizing personalized medicine through the use of comprehensive computational pipeline for systematic data mining and analysis.

Systematic analysis and functional annotation of variations in the genome of an Indian individual

Whole genome sequencing of personal genomes has revealed a large repertoire of genomic variations and has provided a rich template for identification of common and rare variants in genomes in addition to understanding the genetic basis of diseases. The widespread application of personal genome sequencing in clinical settings for predictive and preventive medicine has been limited due to the lack of comprehensive computational analysis pipelines. We have used next‐generation sequencing technology to sequence the whole genome of a self‐declared healthy male of Indian origin. We have generated around 28X of the reference human genome with over 99% coverage. Analysis revealed over 3 million single nucleotide variations and about 490,000 small insertion–deletion events including several novel variants. Using this dataset as a template, we designed a comprehensive computational analysis pipeline for the systematic analysis and annotation of functionally relevant variants in the genome. This study follows a systematic and intuitive data analysis workflow to annotate genome variations and its potential functional effects. Moreover, we integrate predictive analysis of pharmacogenomic traits with emphasis on drugs for which pharmacogenomic testing has been recommended. This study thus provides the template for genome‐scale analysis of personal genomes for personalized medicine. Hum Mutat 33:1133–1140, 2012.

FishMap: a community resource for zebrafish genomics.

An enormous amount of information on a genomics scale is available for zebrafish (Danio rerio), which is a well-studied model organism for human diseases. However, a majority of this annotation is scattered in obscure data sources. There have been limited efforts to present it on a unified and integrated platform, which would help to understand the biological processes in this organism better. FishMap is a unified and centralized resource for storage, retrieval, and display of genomic information of zebrafish. The datasets have been methodically collected from various resources and supplementary information of publications and mapped to the zebrafish genome. The data are organized into nine major sections, which include comparative genomics, mapping and sequencing, gene and gene predictions, expression and regulation, and variation and repeats. A number of unique sections have been incorporated, which include tracks on noncoding gene annotation, location of retrovirus/transposon integrations in the genome, and their flanking genomic sequences and novel transcripts. The datasets are linked to related data sources. FishMap is built on the Gbrowse, which is a part of the Generic Model Organism Database Consortium Project. The resource also features a Web-based BLAST server for sequence homology search and a gene ID converter that would enable users to sift through different interchangeable gene annotation identifier systems. The database is amenable to programmatic access through the Distributed Annotation System as well as BioMoby protocols, thus making it a central community resource that can be integrated with existing data mining and analysis workflows. We hope that FishMap would be an integral resource for community participation in zebrafish genomics. The resource is freely available at http://miracle.igib.res.in/fishmap, or at http://fishmap.igib.res.in.