Comparison of sequencing-based methods to profile DNA methylation and identification of monoallelic epigenetic modifications.

Analysis of DNA methylation patterns relies increasingly on sequencing-based profiling methods. The four most frequently used sequencing-based technologies are the bisulfite-based methods MethylC-seq and reduced representation bisulfite sequencing (RRBS), and the enrichment-based techniques methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylated DNA binding domain sequencing (MBD-seq). We applied all four methods to biological replicates of human embryonic stem cells to assess their genome-wide CpG coverage, resolution, cost, concordance and the influence of CpG density and genomic context. The methylation levels assessed by the two bisulfite methods were concordant (their difference did not exceed a given threshold) for 82% for CpGs and 99% of the non-CpG cytosines. Using binary methylation calls, the two enrichment methods were 99% concordant and regions assessed by all four methods were 97% concordant. We combined MeDIP-seq with methylation-sensitive restriction enzyme (MRE-seq) sequencing for comprehensive methylome coverage at lower cost. This, along with RNA-seq and ChIP-seq of the ES cells enabled us to detect regions with allele-specific epigenetic states, identifying most known imprinted regions and new loci with monoallelic epigenetic marks and monoallelic expression

Computational Analysis of Whole-Genome Differential Allelic Expression Data in Human

Allelic imbalance (AI) is a phenomenon where the two alleles of a given gene are expressed at different levels in a given cell, either because of epigenetic inactivation of one of the two alleles, or because of genetic variation in regulatory regions. Recently, Bing et al. have described the use of genotyping arrays to assay AI at a high resolution (∼750,000 SNPs across the autosomes). In this paper, we investigate computational approaches to analyze this data and identify genomic regions with AI in an unbiased and robust statistical manner. We propose two families of approaches: (i) a statistical approach based on z-score computations, and (ii) a family of machine learning approaches based on Hidden Markov Models. Each method is evaluated using previously published experimental data sets as well as with permutation testing. When applied to whole genome data from 53 HapMap samples, our approaches reveal that allelic imbalance is widespread (most expressed genes show evidence of AI in at least one of our 53 samples) and that most AI regions in a given individual are also found in at least a few other individuals. While many AI regions identified in the genome correspond to known protein-coding transcripts, others overlap with recently discovered long non-coding RNAs. We also observe that genomic regions with AI not only include complete transcripts with consistent differential expression levels, but also more complex patterns of allelic expression such as alternative promoters and alternative 3′ end. The approaches developed not only shed light on the incidence and mechanisms of allelic expression, but will also help towards mapping the genetic causes of allelic expression and identify cases where this variation may be linked to diseases

Meta-Analysis for Genome-Wide Association Study Identifies Multiple Variants at the BIN1 Locus Associated with Late-Onset Alzheimer's Disease

Recent GWAS studies focused on uncovering novel genetic loci related to AD have revealed associations with variants near CLU, CR1, PICALM and BIN1. In this study, we conducted a genome-wide association study in an independent set of 1034 cases and 1186 controls using the Illumina genotyping platforms. By coupling our data with available GWAS datasets from the ADNI and GenADA, we replicated the original associations in both PICALM (rs3851179) and CR1 (rs3818361). The PICALM variant seems to be non-significant after we adjusted for APOE e4 status. We further tested our top markers in 751 independent cases and 751 matched controls. Besides the markers close to the APOE locus, a marker (rs12989701) upstream of BIN1 locus was replicated and the combined analysis reached genome-wide significance level (p = 5E-08). We combined our data with the published Harold et al. study and meta-analysis with all available 6521 cases and 10360 controls at the BIN1 locus revealed two significant variants (rs12989701, p = 1.32E-10 and rs744373, p = 3.16E-10) in limited linkage disequilibrium (r2 = 0.05) with each other. The independent contribution of both SNPs was supported by haplotype conditional analysis. We also conducted multivariate analysis in canonical pathways and identified a consistent signal in the downstream pathways targeted by Gleevec (P = 0.004 in Pfizer; P = 0.028 in ADNI and P = 0.04 in GenADA). We further tested variants in CLU, PICALM, BIN1 and CR1 for association with disease progression in 597 AD patients where longitudinal cognitive measures are sufficient. Both the PICALM and CLU variants showed nominal significant association with cognitive decline as measured by change in Clinical Dementia Rating-sum of boxes (CDR-SB) score from the baseline but did not pass multiple-test correction. Future experiments will help us better understand potential roles of these genetic loci in AD pathology

Parp1 Localizes within the Dnmt1 Promoter and Protects Its Unmethylated State by Its Enzymatic Activity

Aberrant hypermethylation of CpG islands in housekeeping gene promoters and widespread genome hypomethylation are typical events occurring in cancer cells. The molecular mechanisms behind these cancer-related changes in DNA methylation patterns are not well understood. Two questions are particularly important: (i) how are CpG islands protected from methylation in normal cells, and how is this protection compromised in cancer cells, and (ii) how does the genome-wide demethylation in cancer cells occur. The latter question is especially intriguing since so far no DNA demethylase enzyme has been found.Our data show that the absence of ADP-ribose polymers (PARs), caused by ectopic over-expression of poly(ADP-ribose) glycohydrolase (PARG) in L929 mouse fibroblast cells leads to aberrant methylation of the CpG island in the promoter of the Dnmt1 gene, which in turn shuts down its transcription. The transcriptional silencing of Dnmt1 may be responsible for the widespread passive hypomethylation of genomic DNA which we detect on the example of pericentromeric repeat sequences. Chromatin immunoprecipitation results show that in normal cells the Dnmt1 promoter is occupied by poly(ADP-ribosyl)ated Parp1, suggesting that PARylated Parp1 plays a role in protecting the promoter from methylation.In conclusion, the genome methylation pattern following PARG over-expression mirrors the pattern characteristic of cancer cells, supporting our idea that the right balance between Parp/Parg activities maintains the DNA methylation patterns in normal cells. The finding that in normal cells Parp1 and ADP-ribose polymers localize on the Dnmt1 promoter raises the possibility that PARylated Parp1 marks those sequences in the genome that must remain unmethylated and protects them from methylation, thus playing a role in the epigenetic regulation of gene expression

Assembly and dynamics of the bacteriophage T4 homologous recombination machinery

Homologous recombination (HR), a process involving the physical exchange of strands between homologous or nearly homologous DNA molecules, is critical for maintaining the genetic diversity and genome stability of species. Bacteriophage T4 is one of the classic systems for studies of homologous recombination. T4 uses HR for high-frequency genetic exchanges, for homology-directed DNA repair (HDR) processes including DNA double-strand break repair, and for the initiation of DNA replication (RDR). T4 recombination proteins are expressed at high levels during T4 infection in E. coli, and share strong sequence, structural, and/or functional conservation with their counterparts in cellular organisms. Biochemical studies of T4 recombination have provided key insights on DNA strand exchange mechanisms, on the structure and function of recombination proteins, and on the coordination of recombination and DNA synthesis activities during RDR and HDR. Recent years have seen the development of detailed biochemical models for the assembly and dynamics of presynaptic filaments in the T4 recombination system, for the atomic structure of T4 UvsX recombinase, and for the roles of DNA helicases in T4 recombination. The goal of this chapter is to review these recent advances and their implications for HR and HDR mechanisms in all organisms

A Functional Proteomic Method for Biomarker Discovery

The sequencing of the human genome holds out the hope for personalized medicine, but it is clear that analysis of DNA or RNA content alone is not sufficient to understand most disease processes. Proteomic strategies that allow unbiased identification of proteins and their post-transcriptional and -translation modifications are an essential complement to genomic strategies. However, the enormity of the proteome and limitations in proteomic methods make it difficult to determine the targets that are particularly relevant to human disease. Methods are therefore needed that allow rational identification of targets based on function and relevance to disease. Screening methodologies such as phage display, SELEX, and small-molecule combinatorial chemistry have been widely used to discover specific ligands for cells or tissues of interest, such as tumors. Those ligands can be used in turn as affinity probes to identify their cognate molecular targets when they are not known in advance. Here we report an easy, robust and generally applicable approach in which phage particles bearing cell- or tissue-specific peptides serve directly as the affinity probes for their molecular targets. For proof of principle, the method successfully identified molecular binding partners, three of them novel, for 15 peptides specific for pancreatic cancer

Discovery of Novel MicroRNAs in Rat Kidney Using Next Generation Sequencing and Microarray Validation

MicroRNAs (miRNAs) are small non-coding RNAs that regulate a variety of biological processes. The latest version of the miRBase database (Release 18) includes 1,157 mouse and 680 rat mature miRNAs. Only one new rat mature miRNA was added to the rat miRNA database from version 16 to version 18 of miRBase, suggesting that many rat miRNAs remain to be discovered. Given the importance of rat as a model organism, discovery of the completed set of rat miRNAs is necessary for understanding rat miRNA regulation. In this study, next generation sequencing (NGS), microarray analysis and bioinformatics technologies were applied to discover novel miRNAs in rat kidneys. MiRanalyzer was utilized to analyze the sequences of the small RNAs generated from NGS analysis of rat kidney samples. Hundreds of novel miRNA candidates were examined according to the mappings of their reads to the rat genome, presence of sequences that can form a miRNA hairpin structure around the mapped locations, Dicer cleavage patterns, and the levels of their expression determined by both NGS and microarray analyses. Nine novel rat hairpin precursor miRNAs (pre-miRNA) were discovered with high confidence. Five of the novel pre-miRNAs are also reported in other species while four of them are rat specific. In summary, 9 novel pre-miRNAs (14 novel mature miRNAs) were identified via combination of NGS, microarray and bioinformatics high-throughput technologies

The Evolution of Epigenetic Regulators CTCF and BORIS/CTCFL in Amniotes

CTCF is an essential, ubiquitously expressed DNA-binding protein responsible for insulator function, nuclear architecture, and transcriptional control within vertebrates. The gene CTCF was proposed to have duplicated in early mammals, giving rise to a paralogue called “brother of regulator of imprinted sites” (BORIS or CTCFL) with DNA binding capabilities similar to CTCF, but testis-specific expression in humans and mice. CTCF and BORIS have opposite regulatory effects on human cancer-testis genes, the anti-apoptotic BAG1 gene, the insulin-like growth factor 2/H19 imprint control region (IGF2/H19 ICR), and show mutually exclusive expression in humans and mice, suggesting that they are antagonistic epigenetic regulators. We discovered orthologues of BORIS in at least two reptilian species and found traces of its sequence in the chicken genome, implying that the duplication giving rise to BORIS occurred much earlier than previously thought. We analysed the expression of CTCF and BORIS in a range of amniotes by conventional and quantitative PCR. BORIS, as well as CTCF, was found widely expressed in monotremes (platypus) and reptiles (bearded dragon), suggesting redundancy or cooperation between these genes in a common amniote ancestor. However, we discovered that BORIS expression was gonad-specific in marsupials (tammar wallaby) and eutherians (cattle), implying that a functional change occurred in BORIS during the early evolution of therian mammals. Since therians show imprinting of IGF2 but other vertebrate taxa do not, we speculate that CTCF and BORIS evolved specialised functions along with the evolution of imprinting at this and other loci, coinciding with the restriction of BORIS expression to the germline and potential antagonism with CTCF

Hypothermia for moderate or severe neonatal encephalopathy in low-income and middle-income countries (HELIX): a randomised controlled trial in India, Sri Lanka, and Bangladesh

Background: Although therapeutic hypothermia reduces death or disability after neonatal encephalopathy in high-income countries, its safety and efficacy in low-income and middle-income countries is unclear. We aimed to examine whether therapeutic hypothermia alongside optimal supportive intensive care reduces death or moderate or severe disability after neonatal encephalopathy in south Asia. Methods: We did a multicountry open-label, randomised controlled trial in seven tertiary neonatal intensive care units in India, Sri Lanka, and Bangladesh. We enrolled infants born at or after 36 weeks of gestation with moderate or severe neonatal encephalopathy and a need for continued resuscitation at 5 min of age or an Apgar score of less than 6 at 5 min of age (for babies born in a hospital), or both, or an absence of crying by 5 min of age (for babies born at home). Using a web-based randomisation system, we allocated infants into a group receiving whole body hypothermia (33·5°C) for 72 h using a servo-controlled cooling device, or to usual care (control group), within 6 h of birth. All recruiting sites had facilities for invasive ventilation, cardiovascular support, and access to 3 Tesla MRI scanners and spectroscopy. Masking of the intervention was not possible, but those involved in the magnetic resonance biomarker analysis and neurodevelopmental outcome assessments were masked to the allocation. The primary outcome was a combined endpoint of death or moderate or severe disability at 18–22 months, assessed by the Bayley Scales of Infant and Toddler Development (third edition) and a detailed neurological examination. Analysis was by intention to treat. This trial is registered with ClinicalTrials.gov, NCT02387385. Findings: We screened 2296 infants between Aug 15, 2015, and Feb 15, 2019, of whom 576 infants were eligible for inclusion. After exclusions, we recruited 408 eligible infants and we assigned 202 to the hypothermia group and 206 to the control group. Primary outcome data were available for 195 (97%) of the 202 infants in the hypothermia group and 199 (97%) of the 206 control group infants. 98 (50%) infants in the hypothermia group and 94 (47%) infants in the control group died or had a moderate or severe disability (risk ratio 1·06; 95% CI 0·87–1·30; p=0·55). 84 infants (42%) in the hypothermia group and 63 (31%; p=0·022) infants in the control group died, of whom 72 (36%) and 49 (24%; p=0·0087) died during neonatal hospitalisation. Five serious adverse events were reported: three in the hypothermia group (one hospital readmission relating to pneumonia, one septic arthritis, and one suspected venous thrombosis), and two in the control group (one related to desaturations during MRI and other because of endotracheal tube displacement during transport for MRI). No adverse events were considered causally related to the study intervention. Interpretation: Therapeutic hypothermia did not reduce the combined outcome of death or disability at 18 months after neonatal encephalopathy in low-income and middle-income countries, but significantly increased death alone. Therapeutic hypothermia should not be offered as treatment for neonatal encephalopathy in low-income and middle-income countries, even when tertiary neonatal intensive care facilities are available. Funding: National Institute for Health Research, Garfield Weston Foundation, and Bill & Melinda Gates Foundation. Translations: For the Hindi, Malayalam, Telugu, Kannada, Singhalese, Tamil, Marathi and Bangla translations of the abstract see Supplementary Materials section