Maternal vitamin B12 deficiency in rats alters DNA methylation in metabolically important genes in their offspring.

Vitamin B12 deficiency is a critical problem worldwide and peri-conceptional deficiency of this vitamin is associated with the risk of complex cardio-metabolic diseases. Nutritional perturbations during these stages of development may lead to changes in the fetal epigenome. Using Wistar rat model system, we have earlier shown that low maternal B12 levels are associated with low birth weight, adiposity, insulin resistance, and increased triglyceride levels in the offspring, which might predispose them to the risk of cardio-metabolic diseases in adulthood. In this study, we have investigated the effects of maternal B12 deficiency on genome-wide DNA methylation profile of the offspring and the effect of rehabilitation of mothers with B12 at conception. We have performed methylated DNA immunoprecipitation sequencing of liver from pups in four groups of Wistar rats: Control (C), B12-restricted (B12R), B12-rehabilitated at conception (B12RC), and B12-rehabilitated at parturition (B12RP). We have analyzed differentially methylated signatures between the three groups as compared to controls. We have identified a total of 214 hypermethylated and 142 hypomethylated regions in the 10 kb upstream region of transcription start site in pups of B12-deficient mothers, which are enriched in genes involved in fatty acid metabolism and mitochondrial transport/metabolism. B12 rehabilitation at conception and parturition is responsible for reversal of methylation status of many of these regions to control levels suggesting a causal association with metabolic phenotypes. Thus, maternal B12 restriction alters DNA methylation of genes involved in important metabolic processes and influences the offspring phenotype, which is reversed by B12 rehabilitation of mothers at conception.

Whole exome sequencing in extended families with autism spectrum disorder implicates four candidate genes.

Autism spectrum disorders (ASDs) are a group of neurodevelopmental disorders, characterized by impairment in communication and social interactions, and by repetitive behaviors. ASDs are highly heritable, and estimates of the number of risk loci range from hundreds to >1000. We considered 7 extended families (size 12-47 individuals), each with ≥3 individuals affected by ASD. All individuals were genotyped with dense SNP panels. A small subset of each family was typed with whole exome sequence (WES). We used a 3-step approach for variant identification. First, we used family-specific parametric linkage analysis of the SNP data to identify regions of interest. Second, we filtered variants in these regions based on frequency and function, obtaining exactly 200 candidates. Third, we compared two approaches to narrowing this list further. We used information from the SNP data to impute exome variant dosages into those without WES. We regressed affected status on variant allele dosage, using pedigree-based kinship matrices to account for relationships. The p value for the test of the null hypothesis that variant allele dosage is unrelated to phenotype was used to indicate strength of evidence supporting the variant. A cutoff of p = 0.05 gave 28 variants. As an alternative third filter, we required Mendelian inheritance in those with WES, resulting in 70 variants. The imputation- and association-based approach was effective. We identified four strong candidate genes for ASD (SEZ6L, HISPPD1, FEZF1, SAMD11), all of which have been previously implicated in other studies, or have a strong biological argument for their relevance.

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.

In vivo protein trapping produces a functional expression codex of the vertebrate proteome.

We describe a conditional in vivo protein-trap mutagenesis system that reveals spatiotemporal protein expression dynamics and can be used to assess gene function in the vertebrate Danio rerio. Integration of pGBT-RP2.1 (RP2), a gene-breaking transposon containing a protein trap, efficiently disrupts gene expression with >97% knockdown of normal transcript amounts and simultaneously reports protein expression for each locus. The mutant alleles are revertible in somatic tissues via Cre recombinase or splice-site-blocking morpholinos and are thus to our knowledge the first systematic conditional mutant alleles outside the mouse model. We report a collection of 350 zebrafish lines that include diverse molecular loci. RP2 integrations reveal the complexity of genomic architecture and gene function in a living organism and can provide information on protein subcellular localization. The RP2 mutagenesis system is a step toward a unified 'codex' of protein expression and direct functional annotation of the vertebrate genome.

Proximity of H2A.Z containing nucleosome to the transcription start site influences gene expression levels in the mammalian liver and brain.

Nucleosome positioning maps of several organisms have shown that Transcription Start Sites (TSSs) are marked by nucleosome depleted regions flanked by strongly positioned nucleosomes. Using genome-wide nucleosome maps and histone variant occupancy in the mouse liver, we show that the majority of genes were associated with a single prominent H2A.Z containing nucleosome in their promoter region. We classified genes into clusters depending on the proximity of H2A.Z to the TSS. The genes with no detectable H2A.Z showed lowest expression level, whereas H2A.Z was positioned closer to the TSS of genes with higher expression levels. We confirmed this relation between the proximity of H2A.Z and expression level in the brain. The proximity of histone variant H2A.Z, but not H3.3 to the TSS, over seven consecutive nucleosomes, was correlated with expression. Further, a nucleosome was positioned over the TSS of silenced genes while it was displaced to expose the TSS in highly expressed genes. Our results suggest that gene expression levels in vivo are determined by accessibility of the TSS and proximity of H2A.Z.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms.

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders; yet, the role of cell type-specific splicing and transcript-isoform diversity during human brain development has not been systematically investigated. In this work, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone and cortical plate regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 distinct isoforms, of which 72.6% were novel (not previously annotated in Gencode version 33), and uncovered a substantial contribution of transcript-isoform diversity-regulated by RNA binding proteins-in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to reprioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

Forward genetic screen using a gene-breaking trap approach identifies a novel role of grin2bb-associated RNA transcript (grin2bbART) in zebrafish heart function.

LncRNA-based control affects cardiac pathophysiologies like myocardial infarction, coronary artery disease, hypertrophy, and myotonic muscular dystrophy. This study used a gene-break transposon (GBT) to screen zebrafish (Danio rerio) for insertional mutagenesis. We identified three insertional mutants where the GBT captured a cardiac gene. One of the adult viable GBT mutants had bradycardia (heart arrhythmia) and enlarged cardiac chambers or hypertrophy; we named it "bigheart." Bigheart mutant insertion maps to grin2bb or N-methyl D-aspartate receptor (NMDAR2B) gene intron 2 in reverse orientation. Rapid amplification of adjacent cDNA ends analysis suggested a new insertion site transcript in the intron 2 of grin2bb. Analysis of the RNA sequencing of wild-type zebrafish heart chambers revealed a possible new transcript at the insertion site. As this putative lncRNA transcript satisfies the canonical signatures, we called this transcript grin2bb associated RNA transcript (grin2bbART). Using in situ hybridization, we confirmed localized grin2bbART expression in the heart, central nervous system, and muscles in the developing embryos and wild-type adult zebrafish atrium and bulbus arteriosus. The bigheart mutant had reduced Grin2bbART expression. We showed that bigheart gene trap insertion excision reversed cardiac-specific arrhythmia and atrial hypertrophy and restored grin2bbART expression. Morpholino-mediated antisense downregulation of grin2bbART in wild-type zebrafish embryos mimicked bigheart mutants; this suggests grin2bbART is linked to bigheart. Cardiovascular tissues use Grin2bb as a calcium-permeable ion channel. Calcium imaging experiments performed on bigheart mutants indicated calcium mishandling in the heart. The bigheart cardiac transcriptome showed differential expression of calcium homeostasis, cardiac remodeling, and contraction genes. Western blot analysis highlighted Camk2d1 and Hdac1 overexpression. We propose that altered calcium activity due to disruption of grin2bbART, a putative lncRNA in bigheart, altered the Camk2d-Hdac pathway, causing heart arrhythmia and hypertrophy in zebrafish.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms.

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders, yet the role of cell-type-specific splicing or transcript-isoform diversity during human brain development has not been systematically investigated. Here, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 unique isoforms, of which 72.6% are novel (unannotated in Gencode-v33), and uncovered a substantial contribution of transcript-isoform diversity, regulated by RNA binding proteins, in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to re-prioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders. One-Sentence Summary: A cell-specific atlas of gene isoform expression helps shape our understanding of brain development and disease. Structured Abstract: INTRODUCTION: The development of the human brain is regulated by precise molecular and genetic mechanisms driving spatio-temporal and cell-type-specific transcript expression programs. Alternative splicing, a major mechanism increasing transcript diversity, is highly prevalent in the human brain, influences many aspects of brain development, and has strong links to neuropsychiatric disorders. Despite this, the cell-type-specific transcript-isoform diversity of the developing human brain has not been systematically investigated.RATIONALE: Understanding splicing patterns and isoform diversity across the developing neocortex has translational relevance and can elucidate genetic risk mechanisms in neurodevelopmental disorders. However, short-read sequencing, the prevalent technology for transcriptome profiling, is not well suited to capturing alternative splicing and isoform diversity. To address this, we employed third-generation long-read sequencing, which enables capture and sequencing of complete individual RNA molecules, to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution.RESULTS: We profiled microdissected GZ and CP regions of post-conception week (PCW) 15-17 human neocortex in bulk and at single-cell resolution across six subjects using high-fidelity long-read sequencing (PacBio IsoSeq). We identified 214,516 unique isoforms, of which 72.6% were novel (unannotated in Gencode), and >7,000 novel exons, expanding the proteome by 92,422 putative proteoforms. We uncovered thousands of isoform switches during cortical neurogenesis predicted to impact RNA regulatory domains or protein structure and implicating previously uncharacterized RNA-binding proteins in cellular identity and neuropsychiatric disease. At the single-cell level, early-stage excitatory neurons exhibited the greatest isoform diversity, and isoform-centric single-cell clustering led to the identification of previously uncharacterized cell states. We systematically assessed the contribution of transcriptomic features, and localized cell and spatio-temporal transcript expression signatures across neuropsychiatric disorders, revealing predominant enrichments in dynamic isoform expression and utilization patterns and that the number and complexity of isoforms per gene is strongly predictive of disease. Leveraging this resource, we re-prioritized thousands of rare de novo risk variants associated with autism spectrum disorders (ASD), intellectual disability (ID), and neurodevelopmental disorders (NDDs), more broadly, to potentially more severe consequences and revealed a larger proportion of cryptic splice variants with the expanded transcriptome annotation provided in this study.CONCLUSION: Our study offers a comprehensive landscape of isoform diversity in the human neocortex during development. This extensive cataloging of novel isoforms and splicing events sheds light on the underlying mechanisms of neurodevelopmental disorders and presents an opportunity to explore rare genetic variants linked to these conditions. The implications of our findings extend beyond fundamental neuroscience, as they provide crucial insights into the molecular basis of developmental brain disorders and pave the way for targeted therapeutic interventions. To facilitate exploration of this dataset we developed an online portal ( https://sciso.gandallab.org/ ).

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.

Family-based exome sequencing and case-control analysis implicate CEP41 as an ASD gene.

Autism Spectrum Disorder (ASD) is a complex neurodevelopmental disorder with a strong genetic component. Although next-generation sequencing (NGS) technologies have been successfully applied to gene identification in de novo ASD, the genetic architecture of familial ASD remains largely unexplored. Our approach, which leverages the high specificity and sensitivity of NGS technology, has focused on rare variants in familial autism. We used NGS exome sequencing in 26 families with distantly related affected individuals to identify genes with private gene disrupting and missense variants of interest (VOI). We found that the genes carrying VOIs were enriched for biological processes related to cell projection organization and neuron development, which is consistent with the neurodevelopmental hypothesis of ASD. For a subset of genes carrying VOIs, we then used targeted NGS sequencing and gene-based variant burden case-control analysis to test for association with ASD. Missense variants in one gene, CEP41, associated significantly with ASD (p = 6.185e-05). Homozygous gene-disrupting variants in CEP41 were initially found to be responsible for recessive Joubert syndrome. Using a zebrafish model, we evaluated the mechanism by which the CEP41 variants might contribute to ASD. We found that CEP41 missense variants affect development of the axonal tract, cranial neural crest migration and social behavior phenotype. Our work demonstrates the involvement of CEP41 heterozygous missense variants in ASD and that biological processes involved in cell projection organization and neuron development are enriched in ASD families we have studied.

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.

Next Generation Sequencing Mitochondrial DNA Analysis in Autism Spectrum Disorder.

Autism is a complex genetic disorder where both de-novo and inherited genetics factors play a role. Next generation sequencing approaches have been extensively used to identify rare variants associated with autism. To date, all such studies were focused on nuclear genome; thereby leaving the role of mitochondrial DNA (mtDNA) variation in autism unexplored. Recently, analytical tools have been developed to evaluate mtDNA in whole-exome data. We have analyzed the mtDNA sequence derived from whole-exome sequencing in 10 multiplex families. In one of the families we have identified two variants of interest in MT-ND5 gene that were previously determined to impair mitochondrial function. In addition in a second family we have identified two VOIs; mtDNA variant in MT-ATP6 and nuclear DNA variant in NDUFS4, where both VOIs are within mitochondrial Respiratory Chain Complex. Our findings provide further support for the role of mitochondria in ASD and confirm that whole-exome sequencing allows for analysis of mtDNA, which sets a stage for further comprehensive genetic investigations of the role of mitochondria in autism. Autism Res 2017, 10: 1338-1343. © 2017 International Society for Autism Research, Wiley Periodicals, Inc.

Comparative whole-genome analysis of clinical isolates reveals characteristic architecture of Mycobacterium tuberculosis pangenome.

The tubercle complex consists of closely related mycobacterium species which appear to be variants of a single species. Comparative genome analysis of different strains could provide useful clues and insights into the genetic diversity of the species. We integrated genome assemblies of 96 strains from Mycobacterium tuberculosis complex (MTBC), which included 8 Indian clinical isolates sequenced and assembled in this study, to understand its pangenome architecture. We predicted genes for all the 96 strains and clustered their respective CDSs into homologous gene clusters (HGCs) to reveal a hard-core, soft-core and accessory genome component of MTBC. The hard-core (HGCs shared amongst 100% of the strains) was comprised of 2,066 gene clusters whereas the soft-core (HGCs shared amongst at least 95% of the strains) comprised of 3,374 gene clusters. The change in the core and accessory genome components when observed as a function of their size revealed that MTBC has an open pangenome. We identified 74 HGCs that were absent from reference strains H37Rv and H37Ra but were present in most of clinical isolates. We report PCR validation on 9 candidate genes depicting 7 genes completely absent from H37Rv and H37Ra whereas 2 genes shared partial homology with them accounting to probable insertion and deletion events. The pangenome approach is a promising tool for studying strain specific genetic differences occurring within species. We also suggest that since selecting appropriate target genes for typing purposes requires the expected target gene be present in all isolates being typed, therefore estimating the core-component of the species becomes a subject of prime importance.

Human 45,X fibroblast transcriptome reveals distinct differentially expressed genes including long noncoding RNAs potentially associated with the pathophysiology of Turner syndrome.

Turner syndrome is a chromosomal abnormality characterized by the absence of whole or part of the X chromosome in females. This X aneuploidy condition is associated with a diverse set of clinical phenotypes such as gonadal dysfunction, short stature, osteoporosis and Type II diabetes mellitus, among others. These phenotypes differ in their severity and penetrance among the affected individuals. Haploinsufficiency for a few X linked genes has been associated with some of these disease phenotypes. RNA sequencing can provide valuable insights to understand molecular mechanism of disease process. In the current study, we have analysed the transcriptome profiles of human untransformed 45,X and 46,XX fibroblast cells and identified differential expression of genes in these two karyotypes. Functional analysis revealed that these differentially expressing genes are associated with bone differentiation, glucose metabolism and gonadal development pathways. We also report differential expression of lincRNAs in X monosomic cells. Our observations provide a basis for evaluation of cellular and molecular mechanism(s) in the establishment of Turner syndrome phenotypes.

Exome sequencing reveals a novel mutation, p.L325H, in the KRT5 gene associated with autosomal dominant Epidermolysis Bullosa Simplex Koebner type in a large family from western India.

We report a large, non-consanguineous family comprising five generations of individuals residing in Gujarat, India affected with localized Epidermolysis Bullosa Simplex (EBS) Koebner type. We analyzed 14 individuals including 9 affected individuals from this family. Exome sequencing in two cases suggested a novel non-synonymous variation, p.L325H, in the KRT5 gene. The present analysis also reports the first causative mutation of EBS Koebner type from India.

The Zebrafish GenomeWiki: a crowdsourcing approach to connect the long tail for zebrafish gene annotation.

A large repertoire of gene-centric data has been generated in the field of zebrafish biology. Although the bulk of these data are available in the public domain, most of them are not readily accessible or available in nonstandard formats. One major challenge is to unify and integrate these widely scattered data sources. We tested the hypothesis that active community participation could be a viable option to address this challenge. We present here our approach to create standards for assimilation and sharing of information and a system of open standards for database intercommunication. We have attempted to address this challenge by creating a community-centric solution for zebrafish gene annotation. The Zebrafish GenomeWiki is a 'wiki'-based resource, which aims to provide an altruistic shared environment for collective annotation of the zebrafish genes. The Zebrafish GenomeWiki has features that enable users to comment, annotate, edit and rate this gene-centric information. The credits for contributions can be tracked through a transparent microattribution system. In contrast to other wikis, the Zebrafish GenomeWiki is a 'structured wiki' or rather a 'semantic wiki'. The Zebrafish GenomeWiki implements a semantically linked data structure, which in the future would be amenable to semantic search. Database URL: http://genome.igib.res.in/twiki.

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.

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.

Draft Genome Sequence of a Clinical Isolate of Multidrug-Resistant Mycobacterium tuberculosis East African Indian Strain OSDD271.

We describe the genome sequencing and analysis of a clinical isolate of Mycobacterium tuberculosis East African Indian (EAI) strain OSDD271 from India.

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 subject's 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.