Search

Epilepsy will affect nearly 3% of people at some point during their lifetime. Previous copy number variants (CNVs) studies of epilepsy have used array-based technology and were restricted to the detection of large or exonic events. In contrast, whole-genome sequencing (WGS) has the potential to more comprehensively profile CNVs but existing analytic methods suffer from limited accuracy. We show that this is in part due to the non-uniformity of read coverage, even after intra-sample normalization. To improve on this, we developed PopSV, an algorithm that uses multiple samples to control for technical variation and enables the robust detection of CNVs. Using WGS and PopSV, we performed a comprehensive characterization of CNVs in 198 individuals affected with epilepsy and 301 controls. For both large and small variants, we found an enrichment of rare exonic events in epilepsy patients, especially in genes with predicted loss-of-function intolerance. Notably, this genome-wide survey also revealed an enrichment of rare non-coding CNVs near previously known epilepsy genes. This enrichment was strongest for non-coding CNVs located within 100 Kbp of an epilepsy gene and in regions associated with changes in the gene expression, such as expression QTLs or DNase I hypersensitive sites. Finally, we report on 21 potentially damaging events that could be associated with known or new candidate epilepsy genes. Our results suggest that comprehensive sequence-based profiling of CNVs could help explain a larger fraction of epilepsy cases., Author summary Epilepsy is a common neurological disorder affecting around 3% of the population. In some cases, epilepsy is caused by brain trauma or other brain anomalies but there are often no clear causes. Genetic factors have been associated with epilepsy in the past such as rare genetic variations found by linkage studies as well as common genetic variations found by genome-wide association studies and large copy-number variants. We sequenced the genome of ∼200 epilepsy patients and ∼300 healthy controls and compared the distribution of deletion (loss of a copy) and duplication (additional copy) of genomic regions. Thanks to the sequencing technology and a new method that takes advantage of the large sample size, we could compare the distribution of small copy-number variants between epilepsy patients and controls. Overall, we found that small variants are also associated with epilepsy. Indeed, the genome of epilepsy patients had more exonic copy-number variants, especially when rare or affecting genes with predicted loss-of-function intolerance. Focusing on regions around genes that have been previously associated with epilepsy, we also found more non-coding variants in epilepsy patients, especially deletions or variants in regulatory regions. Finally, we provide a list of 21 regions in which we found likely pathogenic variants.

International audience; To assess the role of rare copy number variations in Alzheimer's disease (AD), we conducted a case-control study using whole-exome sequencing data from 522 early-onset cases and 584 controls. The most recurrent rearrangement was a 17q21.31 microduplication, overlapping the CRHR1, MAPT, STH and KANSL1 genes that was found in four cases, including one de novo rearrangement, and was absent in controls. The increased MAPT gene dosage led to a 1.6-1.9-fold expression of the MAPT messenger RNA. Clinical signs, neuroimaging and cerebrospinal fluid biomarker profiles were consistent with an AD diagnosis in MAPT duplication carriers. However, amyloid positon emission tomography (PET) imaging, performed in three patients, was negative. Analysis of an additional case with neuropathological examination confirmed that the MAPT duplication causes a complex tauopathy, including prominent neurofibrillary tangle pathology in the medial temporal lobe without amyloid-β deposits. 17q21.31 duplication is the genetic basis of a novel entity marked by prominent tauopathy, leading to early-onset dementia with an AD clinical phenotype. This entity could account for a proportion of probable AD cases with negative amyloid PET imaging recently identified in large clinical series.

DNA methylation studies have enabled researchers to understand methylation patterns and their regulatory roles in biological processes and disease. However, only a limited number of statistical approaches have been developed to provide formal quantitative analysis. Specifically, a few available methods do identify differentially methylated CpG (DMC) sites or regions (DMR), but they suffer from limitations that arise mostly due to challenges inherent in bisulfite sequencing data. These challenges include: (1) that read-depths vary considerably among genomic positions and are often low; (2) both methylation and autocorrelation patterns change as regions change; and (3) CpG sites are distributed unevenly. Furthermore, there are several methodological limitations: almost none of these tools is capable of comparing multiple groups and/or working with missing values, and only a few allow continuous or multiple covariates. The last of these is of great interest among researchers, as the goal is often to find which regions of the genome are associated with several exposures and traits. To tackle these issues, we have developed an efficient DMC identification method based on Hidden Markov Models (HMMs) called "DMCHMM" which is a three-step approach (model selection, prediction, testing) aiming to address the aforementioned drawbacks. Our proposed method is different from other HMM methods since it profiles methylation of each sample separately, hence exploiting inter-CpG autocorrelation within samples, and it is more flexible than previous approaches by allowing multiple hidden states. Using simulations, we show that DMCHMM has the best performance among several competing methods. An analysis of cell-separated blood methylation profiles is also provided.

Androgen receptor (AR) signaling reprograms cellular metabolism to support prostate cancer (PCa) growth and survival. Another key regulator of cellular metabolism is mTOR, a kinase found in diverse protein complexes and cellular localizations, including the nucleus. However, whether nuclear mTOR plays a role in PCa progression and participates in direct transcriptional cross-talk with the AR is unknown. Here, via the intersection of gene expression, genomic, and metabolic studies, we reveal the existence of a nuclear mTOR-AR transcriptional axis integral to the metabolic rewiring of PCa cells. Androgens reprogram mTOR-chromatin associations in an AR-dependent manner in which activation of mTOR-dependent metabolic gene networks is essential for androgen-induced aerobic glycolysis and mitochondrial respiration. In models of castration-resistant PCa cells, mTOR was capable of transcriptionally regulating metabolic gene programs in the absence of androgens, highlighting a potential novel castration resistance mechanism to sustain cell metabolism even without a functional AR. Remarkably, we demonstrate that increased mTOR nuclear localization is indicative of poor prognosis in patients, with the highest levels detected in castration-resistant PCa tumors and metastases. Identification of a functional mTOR targeted multigene signature robustly discriminates between normal prostate tissues, primary tumors, and hormone refractory metastatic samples but is also predictive of cancer recurrence. This study thus underscores a paradigm shift from AR to nuclear mTOR as being the master transcriptional regulator of metabolism in PCa.

The primary cilium is an outward projecting antenna-like organelle with an important role in bone mechanotransduction. The capacity to sense mechanical stimuli can affect important cellular and molecular aspects of bone tissue. Idiopathic scoliosis (IS) is a complex pediatric disease of unknown cause, defined by abnormal spinal curvatures. We demonstrate significant elongation of primary cilia in IS patient bone cells. In response to mechanical stimulation, these IS cells differentially express osteogenic factors, mechanosensitive genes, and signaling genes. Considering that numerous ciliary genes are associated with a scoliosis phenotype, among ciliopathies and knockout animal models, we expected IS patients to have an accumulation of rare variants in ciliary genes. Instead, our SKAT-O analysis of whole exomes showed an enrichment among IS patients for rare variants in genes with a role in cellular mechanotransduction. Our data indicates defective cilia in IS bone cells, which may be linked to heterogeneous gene variants pertaining to cellular mechanotransduction.

Recent genomic studies of sporadic clear cell renal cell carcinoma (ccRCC) have uncovered novel driver genes and pathways. Given the unequal incidence rates among men and women (male:female incidence ratio approaches 2:1), we compared the genome-wide distribution of the chromosomal abnormalities in both sexes. We observed a higher frequency for the somatic recurrent chromosomal copy number variations (CNVs) of autosomes in male subjects, whereas somatic loss of chromosome X was detected exclusively in female patients (17.1%). Furthermore, somatic loss of chromosome Y (LOY) was detected in about 40% of male subjects, while mosaic LOY was detected in DNA isolated from peripheral blood in 9.6% of them, and was the only recurrent CNV in constitutional DNA samples. LOY in constitutional DNA, but not in tumor DNA was associated with older age. Amongst Y-linked genes that were downregulated due to LOY, KDM5D and KDM6C epigenetic modifiers have functionally-similar X-linked homologs whose deficiency is involved in ccRCC progression. Our findings establish somatic LOY as a highly recurrent genetic defect in ccRCC that leads to downregulation of hitherto unsuspected epigenetic factors, and suggest that different mechanisms may underlie the somatic and mosaic LOY observed in tumors and peripheral blood, respectively.

Background The functional impact of genetic variation has been extensively surveyed, revealing that genetic changes correlated to phenotypes lie mostly in non-coding genomic regions. Studies have linked allele-specific genetic changes to gene expression, DNA methylation, and histone marks but these investigations have only been carried out in a limited set of samples. Results We describe a large-scale coordinated study of allelic and non-allelic effects on DNA methylation, histone mark deposition, and gene expression, detecting the interrelations between epigenetic and functional features at unprecedented resolution. We use information from whole genome and targeted bisulfite sequencing from 910 samples to perform genotype-dependent analyses of allele-specific methylation (ASM) and non-allelic methylation (mQTL). In addition, we introduce a novel genotype-independent test to detect methylation imbalance between chromosomes. Of the ~2.2 million CpGs tested for ASM, mQTL, and genotype-independent effects, we identify ~32% as being genetically regulated (ASM or mQTL) and ~14% as being putatively epigenetically regulated. We also show that epigenetically driven effects are strongly enriched in repressed regions and near transcription start sites, whereas the genetically regulated CpGs are enriched in enhancers. Known imprinted regions are enriched among epigenetically regulated loci, but we also observe several novel genomic regions (e.g., HOX genes) as being epigenetically regulated. Finally, we use our ASM datasets for functional interpretation of disease-associated loci and show the advantage of utilizing naïve T cells for understanding autoimmune diseases. Conclusions Our rich catalogue of haploid methylomes across multiple tissues will allow validation of epigenome association studies and exploration of new biological models for allelic exclusion in the human genome. Electronic supplementary material The online version of this article (doi:10.1186/s13059-017-1173-7) contains supplementary material, which is available to authorized users.

International audience; We performed whole-exome and whole-genome sequencing in 927 late-onset Alzheimer disease (LOAD) cases, 852 early-onset AD (EOAD) cases, and 1273 controls from France. We assessed the evidence for gene-based association of rare variants with AD in 6 genes for which an association with such variants was previously claimed. When aggregating protein-truncating and missense-predicted damaging variants, we found exome-wide significant association between EOAD risk~and rare variants in SORL1, TREM2, and ABCA7. No exome-wide significant signal was obtained in the LOAD sample, and significance of the order of 10$^{-6}$ was observed in the whole AD group for TREM2. Our study confirms previous gene-level results for TREM2, SORL1, and ABCA7 and provides a~clearer insight into the classes of rare variants involved. Despite different effect sizes and~varying~cumulative minor allele frequencies, the rare protein-truncating and missense-predicted~damaging variants in TREM2, SORL1, and ABCA7 contribute similarly to the heritability of EOAD and explain between 1.1% and 1.5% of EOAD heritability each, compared with 9.12% for APOE $\varepsilon$4.

Transposable elements (TEs) have contributed greatly to the regulatory landscape of many cell types, including to human embryonic stem cells. In this chapter we summarize recent studies that have shown that TE insertions have silenced nearby genes, introduced transcription factor-binding sites, as well as produced novel transcripts including those that are translated into viral particles and those that are functional as noncoding RNA. We focus on presenting results on human and other mammalian pluripotent stem cells because of their importance in studying human health and evolution. Finally we also discuss the applications, implications, and questions that are raised from this newfound knowledge.

Item does not contain fulltext Developmental and epileptic encephalopathy (DEE) is a group of conditions characterized by the co-occurrence of epilepsy and intellectual disability (ID), typically with developmental plateauing or regression associated with frequent epileptiform activity. The cause of DEE remains unknown in the majority of cases. We performed whole-genome sequencing (WGS) in 197 individuals with unexplained DEE and pharmaco-resistant seizures and in their unaffected parents. We focused our attention on de novo mutations (DNMs) and identified candidate genes containing such variants. We sought to identify additional subjects with DNMs in these genes by performing targeted sequencing in another series of individuals with DEE and by mining various sequencing datasets. We also performed meta-analyses to document enrichment of DNMs in candidate genes by leveraging our WGS dataset with those of several DEE and ID series. By combining these strategies, we were able to provide a causal link between DEE and the following genes: NTRK2, GABRB2, CLTC, DHDDS, NUS1, RAB11A, GABBR2, and SNAP25. Overall, we established a molecular diagnosis in 63/197 (32%) individuals in our WGS series. The main cause of DEE in these individuals was de novo point mutations (53/63 solved cases), followed by inherited mutations (6/63 solved cases) and de novo CNVs (4/63 solved cases). De novo missense variants explained a larger proportion of individuals in our series than in other series that were primarily ascertained because of ID. Moreover, these DNMs were more frequently recurrent than those identified in ID series. These observations indicate that the genetic landscape of DEE might be different from that of ID without epilepsy.

We have previously reported on a unique patient in whom homozygosity for a mutation at IRF8 (IRF8(K108E)) causes a severe immunodeficiency. Laboratory evaluation revealed a highly unusual myeloid compartment, remarkable for the complete absence of CD141 and CD161 monocytes, absence of CD11c1 conventional dendritic cells (DCs) and CD11c1/CD1231 plasmacytoid DCs, and striking granulocytic hyperplasia. The patient initially presented with severe disseminated mycobacterial and mucocutaneous fungal infections and was ultimately cured by cord blood transplant. Sequencing RNA from the IRF8(K108E) patient's primary blood cells prior to transplant shows not only depletion of IRF8-bound and IRF8-regulated transcriptional targets, in keeping with the distorted composition of the myeloid compartment, but also a paucity of transcripts associated with activated CD41 and CD81 T lymphocytes. This suggests that T cells reared in the absence of a functional antigen-presenting compartment in IRF8(K108E) are anergic. Biochemical characterization of the IRF8(K108E) mutant in vitro shows that loss of the positively charged side chain at K108 causes loss of nuclear localization and loss of transcriptional activity, which is concomitant with decreased protein stability, increased ubiquitination, increased small ubiquitin-like modification, and enhanced proteasomal degradation. These findings provide functional insight into the molecular basis of immunodeficiency associated with loss of IRF8.

Infectious diseases are responsible for over 25% of deaths globally, but many more individuals are exposed to deadly pathogens. The outcome of infection results from a set of diverse factors including pathogen virulence factors, the environment, and the genetic make-up of the host. The completion of the human reference genome sequence in 2004 along with technological advances have tremendously accelerated and renovated the tools to study the genetic etiology of infectious diseases in humans and its best characterized mammalian model, the mouse. Advancements in mouse genomic resources have accelerated genome-wide functional approaches, such as gene-driven and phenotype-driven mutagenesis, bringing to the fore the use of mouse models that reproduce accurately many aspects of the pathogenesis of human infectious diseases. Treatment with the mutagen N-ethyl-N-nitrosourea (ENU) has become the most popular phenotype-driven approach. Our team and others have employed mouse ENU mutagenesis to identify host genes that directly impact susceptibility to pathogens of global significance. In this review, we first describe the strategies and tools used in mouse genetics to understand immunity to infection with special emphasis on chemical mutagenesis of the mouse germ-line together with current strategies to efficiently identify functional mutations using next generation sequencing. Then, we highlight illustrative examples of genes, proteins, and cellular signatures that have been revealed by ENU screens and have been shown to be involved in susceptibility or resistance to infectious diseases caused by parasites, bacteria, and viruses.

Studies attempting to describe fish community structure in shallow riverine environments typically rely on electrofishing and/or visual (snorkelling) surveys, but few have addressed the relative efficiencies of these two methods at estimating fish density and biomass across wide ranges of geography, taxonomy and life history stages. Multiple paired electrofishing and visual surveys were conducted in 18 temperate Canadian rivers in order to obtain community-wide density and biomass estimates from both methods. Partial canonical multivariate analyses were applied to the paired fish community matrices comparing the results of both surveying methods at the taxonomic levels of family, genus and species, as well as size classes within families and species, to assess the particular effectiveness of each sampling method. Although electrofishing estimates of family and species richness were generally greater, snorkelling surveys tended to generate higher density and biomass estimates for different size classes of many salmonid and cyprinid species. Moreover, mean river biomass estimates derived from visual surveying matched those obtained from our best mean river biomass estimates arising from the two methods combined. This study provides empirical evidence that electrofishing and visual survey methods generate different types of information when assessing fish community structure at the family level or by size classes. Our results provide ample background information for determining the most accurate sampling method for a particular fish community assemblage, which is fundamental to fisheries management and research. Copyright © 2014 John Wiley & Sons, Ltd.

Diffuse intrinsic pontine glioma (DIPG) is a fatal brain cancer that arises in the brainstem of children, with no effective treatment and near 100% fatality. The failure of most therapies can be attributed to the delicate location of these tumors and to the selection of therapies on the basis of assumptions that DIPGs are molecularly similar to adult disease. Recent studies have unraveled the unique genetic makeup of this brain cancer, with nearly 80% found to harbor a p.Lys27Met histone H3.3 or p.Lys27Met histone H3.1 alteration. However, DIPGs are still thought of as one disease, with limited understanding of the genetic drivers of these tumors. To understand what drives DIPGs, we integrated whole-genome sequencing with methylation, expression and copy number profiling, discovering that DIPGs comprise three molecularly distinct subgroups (H3-K27M, silent and MYCN) and uncovering a new recurrent activating mutation affecting the activin receptor gene ACVR1 in 20% of DIPGs. Mutations in ACVR1 were constitutively activating, leading to SMAD phosphorylation and increased expression of the downstream activin signaling targets ID1 and ID2. Our results highlight distinct molecular subgroups and novel therapeutic targets for this incurable pediatric cancer.

Understanding the mechanisms driving lineage-specific evolution in both primates and rodents has been hindered by the lack of sister clades with a similar phylogenetic structure having high-quality genome assemblies. Here, we have created chromosome-level assemblies of the Mus caroli and Mus pahari genomes. Together with the Mus musculus and Rattus norvegicus genomes, this set of rodent genomes is similar in divergence times to the Hominidae (human-chimpanzee-gorilla-orangutan). By comparing the evolutionary dynamics between the Muridae and Hominidae, we identified punctate events of chromosome reshuffling that shaped the ancestral karyotype of Mus musculus and Mus caroli between 3 and 6 million yr ago, but that are absent in the Hominidae. Hominidae show between four- and sevenfold lower rates of nucleotide change and feature turnover in both neutral and functional sequences, suggesting an underlying coherence to the Muridae acceleration. Our system of matched, high-quality genome assemblies revealed how specific classes of repeats can play lineage-specific roles in related species. Recent LINE activity has remodeled protein-coding loci to a greater extent across the Muridae than the Hominidae, with functional consequences at the species level such as reproductive isolation. Furthermore, we charted a Muridae-specific retrotransposon expansion at unprecedented resolution, revealing how a single nucleotide mutation transformed a specific SINE element into an active CTCF binding site carrier specifically in Mus caroli, which resulted in thousands of novel, species-specific CTCF binding sites. Our results show that the comparison of matched phylogenetic sets of genomes will be an increasingly powerful strategy for understanding mammalian biology.

Following publication of the original article [1], the authors reported an error in Additional file 1.

SummarymTOR and ERRα are key regulators of common metabolic processes, including lipid homeostasis. However, it is currently unknown whether these factors cooperate in the control of metabolism. ChIP-sequencing analyses of mouse liver reveal that mTOR occupies regulatory regions of genes on a genome-wide scale including enrichment at genes shared with ERRα that are involved in the TCA cycle and lipid biosynthesis. Genetic ablation of ERRα and rapamycin treatment, alone or in combination, alter the expression of these genes and induce the accumulation of TCA metabolites. As a consequence, both genetic and pharmacological inhibition of ERRα activity exacerbates hepatic hyperlipidemia observed in rapamycin-treated mice. We further show that mTOR regulates ERRα activity through ubiquitin-mediated degradation via transcriptional control of the ubiquitin-proteasome pathway. Our work expands the role of mTOR action in metabolism and highlights the existence of a potent mTOR/ERRα regulatory axis with significant clinical impact.

We evaluated the effects of temporal variation of fish density estimates on the explanatory power of habitat use models. Fish density estimates were obtained using visual surveys (10 visits) in eighteen 100 m reaches over a 7-week period. Physical attributes of reaches were estimated. Field data were used to develop a simulation domain (10 000 reaches) that reflected the spatio-temporal variability of fish density estimates and physical attributes. Simulations indicated that for a sampling effort of approximately 200 surveys, the number of reaches surveyed (25 to 200) and the number of surveys per reach (1 to 8) affected the adjusted R2 of models by 5% to 42%. The established practice of sampling a maximized number of reaches once did not appear necessarily optimal for developing habitat use models. Analysis of temporal coefficients of variation suggests that species within the same family may require a similar survey design. Hence, for salmonids, it may be more appropriate to sample more reaches once, and for cyprinids, it may be more optimal to repeatedly sample fewer reaches. Resume : Nous avons evalue l'effet de variations temporelles de la densite de poissons sur le pouvoir explicatif de modeles d'utilisation de l'habitat et les variables explicatives incluses dans ces modeles. Les densites d'especes de poissons ont ete estimees 10 fois dans dix-huit troncons de 100 m durant 7 semaines de l'ete 2007. Les caracteristiques physiques ont ete documentees. Les donnees de terrain ont ete utilisees pour developper un domaine de simulations (10 000 troncons) refletant la variabilite spatiale et temporelle des estimes de densites de poissons et la variabilite spatiale des caracteristiques des habitats. Des simulations ont indique que, pour un effort constant d'approximativement 200 inventaires, le nombre de troncons echan- tillonnes (25 a ` 200) et le nombre d'inventaires par troncon (1 a ` 8) ont affecte le R2 ajuste des modeles par 5 a 42 %. Nos resultats suggerent que la strategie d'echantillonnage optimale differe selon la famille d'appartenance de l'espece. Pour certaines (sal- monides), il semble optimal d'echantillonner une seule fois un nombre maximal de sites. Pour d'autres (cyprinides), il serait preferable d'echantillonner de facon repetee un nombre diminue de sites.

The transcriptome is composed of different types of RNA molecules including mRNAs, tRNAs, rRNAs, and other noncoding RNAs that are found inside a cell at a given time. Analyzing transcriptome patterns can shed light on the functional state of the cell as well as on the dynamics of cellular behavior associated with genomic and environmental changes. Likewise, transcriptome analysis has been a major help in solving biological issues and understanding the molecular basis of many diseases including human cancers. Specifically, since targeted and whole genome sequencing studies are becoming more common in identifying the driving factors of cancer, a comprehensive and high-resolution analysis of the transcriptome, as provided by RNA-Sequencing (RNA-Seq), plays a key role in investigating the functional relevance of the identified genomic aberrations. Here, we describe experimental procedures of RNA-Seq and downstream data processing and analysis, with a focus on the identification of abnormally expressed transcripts and genes.

Motivation Many peak detection algorithms have been proposed for ChIP-seq data analysis, but it is not obvious which algorithm and what parameters are optimal for any given dataset. In contrast, regions with and without obvious peaks can be easily labeled by visual inspection of aligned read counts in a genome browser. We propose a supervised machine learning approach for ChIP-seq data analysis, using labels that encode qualitative judgments about which genomic regions contain or do not contain peaks. The main idea is to manually label a small subset of the genome, and then learn a model that makes consistent peak predictions on the rest of the genome. Results We created 7 new histone mark datasets with 12 826 visually determined labels, and analyzed 3 existing transcription factor datasets. We observed that default peak detection parameters yield high false positive rates, which can be reduced by learning parameters using a relatively small training set of labeled data from the same experiment type. We also observed that labels from different people are highly consistent. Overall, these data indicate that our supervised labeling method is useful for quantitatively training and testing peak detection algorithms. Availability and Implementation Labeled histone mark data http://cbio.ensmp.fr/~thocking/chip-seq-chunk-db/, R package to compute the label error of predicted peaks https://github.com/tdhock/PeakError Supplementary information Supplementary data are available at Bioinformatics online.

International audience; Causative variants in APP, PSEN1 or PSEN2 account for a majority of cases of autosomal dominant early-onset Alzheimer disease (ADEOAD, onset before 65 years). Variant detection rates in other EOAD patients, that is, with family history of late-onset AD (LOAD) (and no incidence of EOAD) and sporadic cases might be much lower. We analyzed the genomes from 264 patients using whole-exome sequencing (WES) with high depth of coverage: 90 EOAD patients with family history of LOAD and no incidence of EOAD in the family and 174 patients with sporadic AD starting between 51 and 65 years. We found three PSEN1 and one PSEN2 causative, probably or possibly causative variants in four patients (1.5%). Given the absence of PSEN1, PSEN2 and APP causative variants, we investigated whether these 260 patients might be burdened with protein-modifying variants in 20 genes that were previously shown to cause other types of dementia when mutated. For this analysis, we included an additional set of 160 patients who were previously shown to be free of causative variants in PSEN1, PSEN2 and APP: 107 ADEOAD patients and 53 sporadic EOAD patients with an age of onset before 51 years. In these 420 patients, we detected no variant that might modify the function of the 20 dementia-causing genes. We conclude that EOAD patients with family history of LOAD and no incidence of EOAD in the family or patients with sporadic AD starting between 51 and 65 years have a low variant-detection rate in AD genes.

Transposable elements (TEs) have recently been shown to have many regulatory roles within the genome. In this chapter, we will examine two in silico methods for analyzing TEs and identifying families that may have acquired such functions. The first method will look at how the overrepresentation of a repeat family in a set of genomic features can be discovered. The example situation of OCT4 binding sites originating from LTR7 TE sequences will be used to show how this method could be applied. The second method will describe how to determine if a TE family exhibits a cell type-specific expression pattern. As an example, we will look at the expression of HERV-H, an endogenous retrovirus known to act as an lncRNA in embryonic stem cells. We will use this example to demonstrate how RNA-seq data can be used to compare cell type expression of repeats.

International audience; Investigating 5-methylcytosine (5mC) has led to many hypotheses regarding molecular mechanism underlying human diseases and disorders. Many of these studies, however, utilize bisulfite conversion alone, which cannot distinguish 5mC from its recently discovered oxidative product, 5-hydroxymethylcytosine (5hmC). Furthermore, previous array-based technologies do not have the necessary probes to adequately investigate both modifications simultaneously. In this manuscript, we used technical replicates of DNA from human brain, human blood, and human saliva, in combination with oxidative bisul-fite conversion and Illumina's Infinium MethylationEPIC array, to analyze 5mC and 5hmC at more than 650 000 and 450 000 relevant loci, respectively, in the human genome. We show the presence of loci with detectable 5mC and 5hmC to be equally distributed across chromosomes and genomic features, while also being present in genomic regions with transcriptional regulatory properties. We also describe 2528 5hmC sites common across tissue types that show a strong association with immune-related functions. Lastly, in human brain, we show that 5hmC accounts for one-third of the total signal from bisulfite-converted data. As such, not only do our results confirm the efficacy and sensitivity of pairing oxidative bisulfite conversion and the EPIC array to detect 5mC and 5hmC in all three tissue types, but they also highlight the importance of dissociating 5hmC from 5mC in future studies related to cytosine modifications.

We recently reported that atypical teratoid rhabdoid tumors (ATRTs) comprise at least two transcriptional subtypes with different clinical outcomes; however, the mechanisms underlying therapeutic heterogeneity remained unclear. In this study, we analyzed 191 primary ATRTs and 10 ATRT cell lines to define the genomic and epigenomic landscape of ATRTs and identify subgroup-specific therapeutic targets. We found ATRTs segregated into three epigenetic subgroups with distinct genomic profiles, SMARCB1 genotypes, and chromatin landscape that correlated with differential cellular responses to a panel of signaling and epigenetic inhibitors. Significantly, we discovered that differential methylation of a PDGFRB-associated enhancer confers specific sensitivity of group 2 ATRT cells to dasatinib and nilotinib, and suggest that these are promising therapies for this highly lethal ATRT subtype.

Transposable elements (TEs) have recently been shown to have many regulatory roles within the genome. In this chapter, we will examine two in silico methods for analyzing TEs and identifying families that may have acquired such functions. The first method will look at how the overrepresentation of a repeat family in a set of genomic features can be discovered. The example situation of OCT4 binding sites originating from LTR7 TE sequences will be used to show how this method could be applied. The second method will describe how to determine if a TE family exhibits a cell type-specific expression pattern. As an example, we will look at the expression of HERV-H, an endogenous retrovirus known to act as an lncRNA in embryonic stem cells. We will use this example to demonstrate how RNA-seq data can be used to compare cell type expression of repeats.

Significance Embryonic stem cells (ESCs) maintain a low translation rate; therefore control of mRNA translation is critical for preserving their stemness. We identified a hitherto unstudied transcription factor, Yin-yang 2 (YY2), which is translationaly regulated and controls self-renewal and differentiation of mouse ESCs (mESCs). Although YY2 is essential for mESC self-renewal, increased YY2 expression directs differentiation of mESCs toward cardiovascular lineages. Examination of the Yy2 5′-UTR revealed a multilayer regulatory mechanism through which YY2 expression is dictated by the combined actions of the splicing regulator, Polypyrimidine tract-binding protein 1 (PTBP1), and the translation inhibitors, Eukaryotic initiation factor 4E-binding proteins (4E-BPs). YY2 directly controls the expression of several pluripotency and development-related genes. This study describes a synchronized network of alternative splicing and mRNA translation in controlling self-renewal and differentiation.

The transcriptome is composed of different types of RNA molecules including mRNAs, tRNAs, rRNAs, and other noncoding RNAs that are found inside a cell at a given time. Analyzing transcriptome patterns can shed light on the functional state of the cell as well as on the dynamics of cellular behavior associated with genomic and environmental changes. Likewise, transcriptome analysis has been a major help in solving biological issues and understanding the molecular basis of many diseases including human cancers. Specifically, since targeted and whole genome sequencing studies are becoming more common in identifying the driving factors of cancer, a comprehensive and high-resolution analysis of the transcriptome, as provided by RNA-Sequencing (RNA-Seq), plays a key role in investigating the functional relevance of the identified genomic aberrations. Here, we describe experimental procedures of RNA-Seq and downstream data processing and analysis, with a focus on the identification of abnormally expressed transcripts and genes.

Pediatric glioblastomas (GBM) including diffuse intrinsic pontine gliomas (DIPG) are devastating brain tumors with no effective therapy. Here, we investigated clinical and biological impacts of histone H3.3 mutations. Forty-two DIPGs were tested for H3.3 mutations. Wild-type versus mutated (K27M-H3.3) subgroups were compared for HIST1H3B, IDH, ATRX and TP53 mutations, copy number alterations and clinical outcome. K27M-H3.3 occurred in 71 %, TP53 mutations in 77 % and ATRX mutations in 9 % of DIPGs. ATRX mutations were more frequent in older children (p

Human endogenous retrovirus subfamily H (HERVH) is a class of transposable elements expressed preferentially in human embryonic stem cells (hESCs). Here, we report that the long terminal repeats of HERVH function as enhancers and that HERVH is a nuclear long noncoding RNA required to maintain hESC identity. Furthermore, HERVH is associated with OCT4, coactivators and Mediator subunits. Together, these results uncover a new role of species-specific transposable elements in hESCs. ASTAR (Agency for Sci., Tech. and Research, S’pore)

Intratumoral and interpatient heterogeneity are characteristics of glioblastoma and constitute important challenges in overcoming treatment resistance and developing new, more effective therapies. Using single-cell RNA sequencing, we characterized 60 933 cells from 4 developing fetal brains and 8 glioblastomas. By using fetal brain development as a road map, we show a tri-lineage (astrocytic, oligodendrocytic, and neuronal) hierarchical organization in all glioblastomas. In each patient, a population of progenitor cancer cells was found at the apex of this hierarchy. These cells were enriched in our patient-derived glioma stem cell samples, and, like progenitors in the developing brain, were the main dividing cell population within the cancer. Using expression signatures obtained from single-cell RNA-sequencing, we isolated progenitor cancer cells and compared them to other glioblastoma cell types. We showed the progenitors are the most resistant to chemotherapy and the most tumorigenic in mouse xenograft models. This newly found conserved developmental organization points to the cell of origin and suggests new therapeutic approaches.

Emerging evidence has illustrated the importance of epigenomic reprogramming in cancer, with altered post-translational modifications of histones contributing to the pathogenesis of many tumor types, including breast cancer. However, the contributions of histone modifiers to breast cancer progression are unclear, and how these processes vary between molecular subtypes has yet to be adequately addressed. In this study, we demonstrate that genetic or pharmacologic targeting of the epigenetic modifier Ezh2 dramatically hinders metastatic behavior in a mouse model of breast cancer reflective of the luminal B subtype. We also employ patient-derived xenografts from different intrinsic subtypes of breast cancer to confirm that pharmacologic inhibition of Ezh2 activity hinders metastasis specifically in the luminal B subtype. We further define a molecular mechanism intrinsic to the luminal B subtype whereby EZH2 maintains H3K27me3-mediated repression of the FOXC1 gene, thereby inactivating a FOXC1-driven, antiinvasive transcriptional program. We demonstrate that higher FOXC1 levels are predictive of favorable outcome specifically in luminal B breast cancer patients and establish the use of EZH2 methyltransferase inhibitors as a viable strategy to block metastasis in luminal B breast cancer, where options for targeted therapy are currently limited. Citation Format: Alison Hirukawa, Harvey Smith, Dongmei Zou, Paul Savage, Radia Johnson, Guillaume Bourque, Vincent Giguere, Mark Basik, Cathy Dufour, Morag Park, William Muller. Targeting EZH2 reactivates a breast cancer subtype-specific antimetastatic transcriptional program [abstract]. In: Proceedings of the AACR Special Conference: Advances in Breast Cancer Research; 2017 Oct 7-10; Hollywood, CA. Philadelphia (PA): AACR; Mol Cancer Res 2018;16(8_Suppl):Abstract nr B46.

Intratumoral genetic heterogeneity may relate to aggressiveness and treatment-resistance. However, the genomic abnormalities underlie this heterogeneity remain poorly characterized. To better describe the picture of tumour heterogeneity, single cell approaches to accommodate larger number of cells are highly in demand. Here, we present results of single cell DNA sequencing by using a novel microfluidic droplet system. We applied this approach to discover copy number profiles in two common and therapeutically challenging cancers, isocitrate dehydrogenase wild-type glioblastoma (GBM) and breast cancer. The GBM cells were extracted from patients operated on site and cultured for 2 months in restricted media to isolate glioma stem cells. The breast cancer cells were extracted from patient-derived xenografts (PDXs). A novel microfluidic droplet based platform (10x GenomicsTM) was used to perform sensitive copy number profiling of thousands of single cells simultaneously. Briefly, the platform integrates single cell encapsulation, cell lysis, genomic amplification, and barcoding with oligonucleotides. The barcoded genomic DNA fragments from individual cells were pooled and converted into libraries that are compatible with Illumina sequencers. Thus far, we sequenced one glioma stem cell sample and two breast cancer PDXs. Uniquely mapped reads were counted in variable-bins with GC bias corrected by lowess fit. The Circular Binary Segmentation was performed and the integer copy number was estimated by least square fitting. The single cell copy number profiles were clustered by hierarchical clustering. We achieved approximately 2% and 6% coverage of the genome per cell for the glioblastoma and the breast tumors, respectively. The uniformity of coverage was comparable to other whole genome amplification and single cell multiplex approaches, while requiring 700k-2.1mi reads/single cell with the median absolute deviation of pair-wise differences in read counts between neighboring bins between 0.17-0.21 vs 0.16)1. We identified two major clonal subpopulations in the glioblastoma sample. One subclone had chr1q amplification and the other subclone had a 9.8Mb deletion on chr4. The clonal CNVs observed were validated using single-cell RNA-seq-based CNV analysis, and the frequent CNVs using fluorescent in situ hybridization. We also identified the HER2 copy number profile within a basal like triple-negative breast cancer as 40% of 2.5N and 25% of 3N. This was previously validated by fluorescence in situ hybridization2. By applying this new microfluidic system, our work will provide more insights into genomic heterogeneity for malignant properties of cancers. 1. Navin, N. et al. Tumour evolution inferred by single-cell sequencing. Nature 472, 90-94 (2011). 2. Savage P. et al. A Targetable EGFR-Dependent Tumor-Initiating Program in Breast Cancer. Cell Rep. 21, 1140-1149 (2017). Citation Format: Rui Li, Charles Couturier, Paul Savage, Jean Monlong, Guillaume Bourque, Kevin Petrecca, Morag Park, Jiannis Ragoussis. Sensitive single cell copy number profiling using a novel microfluidic droplet based platform [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2177.

Detection of new genomic control elements is critical in understanding transcriptional regulatory networks in their entirety. We studied the genome-wide binding locations of three key regulatory proteins (POU5F1, also known as OCT4; NANOG; and CTCF) in human and mouse embryonic stem cells. In contrast to CTCF, we found that the binding profiles of OCT4 and NANOG are markedly different, with only approximately 5% of the regions being homologously occupied. We show that transposable elements contributed up to 25% of the bound sites in humans and mice and have wired new genes into the core regulatory network of embryonic stem cells. These data indicate that species-specific transposable elements have substantially altered the transcriptional circuitry of pluripotent stem cells.

Genomes are organized into high-level three-dimensional structures, and DNA elements separated by long genomic distances can in principle interact functionally. Many transcription factors bind to regulatory DNA elements distant from gene promoters. Although distal binding sites have been shown to regulate transcription by long-range chromatin interactions at a few loci, chromatin interactions and their impact on transcription regulation have not been investigated in a genome-wide manner. Here we describe the development of a new strategy, chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) for the de novo detection of global chromatin interactions, with which we have comprehensively mapped the chromatin interaction network bound by oestrogen receptor α (ER-α) in the human genome. We found that most high-confidence remote ER-α-binding sites are anchored at gene promoters through long-range chromatin interactions, suggesting that ER-α functions by extensive chromatin looping to bring genes together for coordinated transcriptional regulation. We propose that chromatin interactions constitute a primary mechanism for regulating transcription in mammalian genomes. © 2009 Macmillan Publishers Limited. All rights reserved.

Nuclear signaling by estrogens rapidly induces the global recruitment of estrogen receptors (ERs) to thousands of highly specific locations in the genome. Here, we have examined whether ER binding sites that are located distal from the transcription start sites of estrogen target genes are functionally relevant. Similar to ER binding sites near the proximal promoter region, ER binding sites located at distal locations are occupied by ERs after estrogen stimulation. And, like proximal bound ERs, ERs occupied at distal sites can recruit coactivators and the RNA polymerase transcription machinery and mediate specific structural changes to chromatin. Furthermore, ERs occupied at the distal sites are capable of communicating with ERs bound at the promoter region, possibly via long range chromosome looping. In functional analysis, disruption of the response elements in the distal ER binding sites abrogated ER binding and significantly reduced transcriptional response. Finally, sequence comparison of the response elements at the distal sites suggests a high level of conservation across different species. Together, our data indicate that distal ER binding sites are bona fide transcriptional enhancers that are involved in long range chromosomal interaction, transcription complex formation, and distinct structural modifications of chromatin across large genomic spans.

Genome rearrangement events, including inversions and translocations, are frequently observed across related microbial species, but the impact of such events on functional diversity is unclear. To clarify this relationship, we compared 4 members of the Gram-negative Burkholderia family (Burkholderia pseudomallei, Burkholderia mallei, Burkholderia thailandensis, and Burkholderia cenocepacia) and identified a core set of 2,590 orthologs present in all 4 species (metagenes). The metagenes were organized into 255 synteny blocks whose relative order has been altered by a predicted minimum of 242 genome rearrangement events. Functionally, metagenes within individual synteny blocks were often related. The molecular divergence of metagenes adjacent to synteny breakpoints (boundary metagenes) was significantly greater compared with metagenes within blocks, suggesting an association between breakpoint locations and local fine-scale nucleotide alterations. This phenomenon, referred to as boundary element associated divergence, was also observed in Pseudomonas and Shigella, suggesting that this is a common phenomenon in prokaryotes. We also observed preferential localization of species-specific genes and insertion sequence element to synteny breakpoints in Burkholderia. Our results suggest that in prokaryotes, genome rearrangements may influence functional diversity through the enhanced divergence of boundary genes and the creation of foci for acquiring and deleting species-specific genes.

Identification of unconventional functional features such as fusion transcripts is a challenging task in the effort to annotate all functional DNA elements in the human genome. Paired-End diTag (PET) analysis possesses a unique capability to accurately and efficiently characterize the two ends of DNA fragments, which may have either normal or unusual compositions. This unique nature of PET analysis makes it an ideal tool for uncovering unconventional features residing in the human genome. Using the PET approach for comprehensive transcriptome analysis, we were able to identify fusion transcripts derived from genome rearrangements and actively expressed retrotransposed pseudogenes, which would be difficult to capture by other means. Here, we demonstrate this unique capability through the analysis of 865,000 individual transcripts in two types of cancer cells. In addition to the characterization of a large number of differentially expressed alternative 5′ and 3′ transcript variants and novel transcriptional units, we identified 70 fusion transcript candidates in this study. One was validated as the product of a fusion gene between BCAS4 and BCAS3 resulting from an amplification followed by a translocation event between the two loci, chr20q13 and chr17q23. Through an examination of PETs that mapped to multiple genomic locations, we identified 4055 retrotransposed loci in the human genome, of which at least three were found to be transcriptionally active. The PET mapping strategy presented here promises to be a useful tool in annotating the human genome, especially aberrations in human cancer genomes.

Objective:To study the association between ABCA7 rare coding variants and Alzheimer disease (AD) in a case-control setting.Methods:We conducted a whole exome analysis among 484 French patients with early-onset AD and 590 ethnically matched controls.Results:After collapsing rare variants (minor allele frequency ≤1%), we detected an enrichment of ABCA7 loss of function (LOF) and predicted damaging missense variants in cases (odds ratio [OR] 3.40, 95% confidence interval [CI] 1.68–7.35, p = 0.0002). Performing a meta-analysis with previously published data, we found that in a combined sample of 1,256 patients and 1,347 controls from France and Belgium, the OR was 2.81 (95% CI 1.89–4.20, p = 3.60 × 10−7).Conclusions:These results confirm that ABCA7 LOF variants are enriched in patients with AD and extend this finding to predicted damaging missense variants.

Despite the initial benefits of treating HER2-amplified breast cancer patients with the tyrosine kinase inhibitor lapatinib, resistance inevitably develops. Here we report that lapatinib induces the degradation of the nuclear receptor ERRα, a master regulator of cellular metabolism, and that the expression of ERRα is restored in lapatinib-resistant breast cancer cells through reactivation of mTOR signalling. Re-expression of ERRα in resistant cells triggers metabolic adaptations favouring mitochondrial energy metabolism through increased glutamine metabolism, as well as ROS detoxification required for cell survival under therapeutic stress conditions. An ERRα inverse agonist counteracts these metabolic adaptations and overcomes lapatinib resistance in a HER2-induced mammary tumour mouse model. This work reveals a molecular mechanism by which ERRα-induced metabolic reprogramming promotes survival of lapatinib-resistant cancer cells and demonstrates the potential of ERRα inhibition as an effective adjuvant therapy in poor outcome HER2-positive breast cancer., Despite initial benefits in treating HER2-positive breast cancer patients with lapatinib, resistance is prevalent. Here the authors show that lapatinib resistance can be ascribed to mTOR-mediated re-activation of ERRα and to the consequent induction of a metabolic adaptation.

Background CpG methylation variation is involved in human trait formation and disease susceptibility. Analyses within populations have been biased towards CpG-dense regions through the application of targeted arrays. We generate whole-genome bisulfite sequencing data for approximately 30 adipose and blood samples from monozygotic and dizygotic twins for the characterization of non-genetic and genetic effects at single-site resolution. Results Purely invariable CpGs display a bimodal distribution with enrichment of unmethylated CpGs and depletion of fully methylated CpGs in promoter and enhancer regions. Population-variable CpGs account for approximately 15–20 % of total CpGs per tissue, are enriched in enhancer-associated regions and depleted in promoters, and single nucleotide polymorphisms at CpGs are a frequent confounder of extreme methylation variation. Differential methylation is primarily non-genetic in origin, with non-shared environment accounting for most of the variance. These non-genetic effects are mainly tissue-specific. Tobacco smoking is associated with differential methylation in blood with no evidence of this exposure impacting cell counts. Opposite to non-genetic effects, genetic effects of CpG methylation are shared across tissues and thus limit inter-tissue epigenetic drift. CpH methylation is rare, and shows similar characteristics of variation patterns as CpGs. Conclusions Our study highlights the utility of low pass whole-genome bisulfite sequencing in identifying methylome variation beyond promoter regions, and suggests that targeting the population dynamic methylome of tissues requires assessment of understudied intergenic CpGs distal to gene promoters to reveal the full extent of inter-individual variation. Electronic supplementary material The online version of this article (doi:10.1186/s13059-015-0856-1) contains supplementary material, which is available to authorized users.

Dietary folate is a major source of methyl groups required for DNA methylation, an epigenetic modification that is actively maintained and remodeled during spermatogenesis. While high-dose folic acid supplementation (up to 10 times the daily recommended dose) has been shown to improve sperm parameters in infertile men, the effects of supplementation on the sperm epigenome are unknown. To assess the impact of 6 months of high-dose folic acid supplementation on the sperm epigenome, we studied 30 men with idiopathic infertility. Blood folate concentrations increased significantly after supplementation with no significant improvements in sperm parameters. Methylation levels of the differentially methylated regions of several imprinted loci (H19, DLK1/GTL2, MEST, SNRPN, PLAGL1, KCNQ1OT1) were normal both before and after supplementation. Reduced representation bisulfite sequencing (RRBS) revealed a significant global loss of methylation across different regions of the sperm genome. The most marked loss of DNA methylation was found in sperm from patients homozygous for the methylenetetrahydrofolate reductase (MTHFR) C677T polymorphism, a common polymorphism in a key enzyme required for folate metabolism. RRBS analysis also showed that most of the differentially methylated tiles were located in DNA repeats, low CpG-density and intergenic regions. Ingenuity Pathway Analysis revealed that methylation of promoter regions was altered in several genes involved in cancer and neurobehavioral disorders including CBFA2T3, PTPN6, COL18A1, ALDH2, UBE4B, ERBB2, GABRB3, CNTNAP4 and NIPA1. Our data reveal alterations of the human sperm epigenome associated with high-dose folic acid supplementation, effects that were exacerbated by a common polymorphism in MTHFR.

Study design A replication association study that used genomic data generated from French-Canadian case and control cohorts. Objectives To determine whether the 53 single nucleotide polymorphisms (SNPs) that were previously associated with spinal deformity progression in an American Caucasian cohort are similarly associated in French-Canadian population. Summary of background data It is widely accepted that genetic factors contribute to adolescent idiopathic scoliosis. The identification of genetic variants associated with the predisposition or progression of curvature could facilitate diagnostic/prognostic tool development. Although 53 SNPs have been associated with spinal curve progression in Caucasian cohorts in the United States, these associations were not replicated in a large Japanese population study, arguing that such a discrepancy could be explained by ethnicity, thus raising the importance of a replication study in an independent Caucasian population of European descent. Methods Genomic data were collected from the French-Canadian population, using the Illumina HumanOmni 2.5M BeadChip. Fifty-two SNPs, tested in ScoliScore or in high linkage disequilibrium with SNPs in the test, were selected to assess their association with scoliosis generally, and with spinal curve progression. One SNP in ScoliScore, rs16909285, could not be evaluated in our Genome-Wide association study. Results None of the SNPs used in ScoliScore were associated with adolescent idiopathic scoliosis curve progression or curve occurrence in French-Canadian population. We evaluated 52 SNPs in severe patients by comparing risk allele frequencies with those in nonsevere patients and with those in control individuals. There was no significant difference between the severe group and the nonsevere group or between the severe group and the control group. Conclusion Although the 52 SNPs studied here were previously associated with curve progression in an American population of European descent, we found no association in French-Canadian patients with adolescent idiopathic scoliosis. This second replication cohort suggests that the lack of association of these SNPs in a Japanese cohort is not due to ethnicity. Level of evidence 4.

Genome-wide demethylation and remethylation of DNA during early embryogenesis is essential for development. Imprinted germline differentially methylated domains (gDMDs) established by sex-specific methylation in either male or female germ cells, must escape these dynamic changes and sustain precise inheritance of both methylated and unmethylated parental alleles. To identify other, gDMD-like sequences with the same epigenetic inheritance properties, we used a modified embryonic stem (ES) cell line that emulates the early embryonic demethylation and remethylation waves. Transient DNMT1 suppression revealed gDMD-like sequences requiring continuous DNMT1 activity to sustain a highly methylated state. Remethylation of these sequences was also compromised in vivo in a mouse model of transient DNMT1 loss in the preimplantation embryo. These novel regions, possessing heritable epigenetic features similar to imprinted-gDMDs are required for normal physiological and developmental processes and when disrupted are associated with disorders such as cancer and autism spectrum disorders. This study presents new perspectives on DNA methylation heritability during early embryo development that extend beyond conventional imprinted-gDMDs.

Large-scale epigenome mapping by the NIH Roadmap Epigenomics Project, the ENCODE Consortium and the International Human Epigenome Consortium (IHEC) produces genome-wide DNA methylation data at one base-pair resolution. We examine how such data can be made open-access while balancing appropriate interpretation and genomic privacy. We propose guidelines for data release that both reduce ambiguity in the interpretation of open-access data and limit immediate access to genetic variation data that are made available through controlled access. Electronic supplementary material The online version of this article (doi:10.1186/s13059-015-0723-0) contains supplementary material, which is available to authorized users.

DNA derived from transposable elements (TEs) constitutes large parts of the genomes of complex eukaryotes, with major impacts not only on genomic research but also on how organisms evolve and function. Although a variety of methods and tools have been developed to detect and annotate TEs, there are as yet no standard benchmarks—that is, no standard way to measure or compare their accuracy. This lack of accuracy assessment calls into question conclusions from a wide range of research that depends explicitly or implicitly on TE annotation. In the absence of standard benchmarks, toolmakers are impeded in improving their tools, annotators cannot properly assess which tools might best suit their needs, and downstream researchers cannot judge how accuracy limitations might impact their studies. We therefore propose that the TE research community create and adopt standard TE annotation benchmarks, and we call for other researchers to join the authors in making this long-overdue effort a success., This work was funded by a Bioinformatics and Computational Biology grant to MB and TB from Genome Canada, Genome Québec, and the Canadian Institutes for Health Research.

Oct4 and Nanog are transcription factors required to maintain the pluripotency and self-renewal of embryonic stem (ES) cells. Using the chromatin immunoprecipitation paired-end ditags method, we mapped the binding sites of these factors in the mouse ES cell genome. We identified 1,083 and 3,006 high-confidence binding sites for Oct4 and Nanog, respectively. Comparative location analyses indicated that Oct4 and Nanog overlap substantially in their targets, and they are bound to genes in different configurations. Using de novo motif discovery algorithms, we defined the cis-acting elements mediating their respective binding to genomic sites. By integrating RNA interference-mediated depletion of Oct4 and Nanog with microarray expression profiling, we demonstrated that these factors can activate or suppress transcription. We further showed that common core downstream targets are important to keep ES cells from differentiating. The emerging picture is one in which Oct4 and Nanog control a cascade of pathways that are intricately connected to govern pluripotency, self-renewal, genome surveillance and cell fate determination.

The MDS1 and ecotropic viral integration site 1 (EVI1) complex locus (MECOM) gene encodes several transcription factor variants including MDS1-EVI1, EVI1 and EVI1Δ324. Although MDS1-EVI1 has been associated with tumor-suppressing activity, EVI1 is a known oncogene in various cancers, whose expression is associated with poor patient survival. Although EVI1Δ324 is co-transcribed with EVI1, its activity in cancer cells is not fully understood. Previous reports described that unlike EVI1, EVI1Δ324 protein cannot transform fibroblasts because of its disrupted N-terminal zinc finger (ZNF) domain. To better understand EVI1Δ324 biology and function, we obtained genome-wide binding occupancies and expression data in ovarian cancer cells. We characterized its DNA-binding sites, binding motif and target genes. Comparative analyses with previous study show that EVI1 and EVI1Δ324 share similar transcriptional activities linked to their common C-terminus ZNF domain. They bind to an E-twenty-six family (ETS)-like motif, target to a large extent the same genes and cooperate with AP1 transcription factor. EVI1Δ324-occupied genes were 70.7% similar to EVI1-bound genes. More strikingly, EVI1 and EVI1Δ324 differentially expressed genes were 99.87% identical, indicating comparable transcriptional regulatory functions. Consistently with gene ontologies linked to these target genes, EVI1Δ324 expression in HeLa cells could enhance anchorage-independent growth, such as EVI1, showing that EVI1Δ324 expression also lead to pro-oncogenic effects. The main specific feature of EVI1 variant is its N-terminus ZNF domain that binds DNA through GATA-like motif. We found that most GATA-like EVI1 chromatin immunoprecipitation sequencing peaks are far from genes and are not involved in transcriptional regulation. These genomic regions were enriched in simple sequence repeats and displayed high meiotic recombination rates. Overall, our genomics analyses uncovered common and specific features of two major MECOM isoforms. Their influence on transcription and downstream cell proliferation was comparable. However, EVI1-specific GATA-like binding sites, from its N-terminus ZNF domain, associated with high recombination rates, suggesting possible additional oncogenic potential for EVI1 in modulating genomic stability.

We present a method for gene network inference and revision based on time-series data. Gene networks are modeled using linear differential equations and a generalized stepwise multiple linear regression procedure is used to recover the interaction coefficients. Our system is designed for the recovery of gene interactions concurrently in many gene regulatory networks related by a tree or a more general graph. We show how this comparative framework can facilitate the recovery of the networks and improve the quality of the solutions inferred.

Recent analysis of genome rearrangements in human and mouse genomes revealed evidence for more rearrangements than thought previously and shed light on previously unknown features of mammalian evolution, like breakpoint reuse and numerous microrearrangements. However, two-way analysis cannot reveal the genomic architecture of ancestral mammals or assign rearrangement events to different lineages. Thus, the “original synteny” problem introduced by Nadeau and Sankoff previously, remains unsolved, as at least three mammalian genomes are required to derive the ancestral mammalian karyotype. We show that availability of the rat genome allows one to reconstruct a putative genomic architecture of the ancestral murid rodent genome. This reconstruction suggests that this ancestral genome retained many previously postulated chromosome associations in the placental ancestor and reveals others that were beyond the resolution of cytogenetic, radiation hybrid mapping, and chromosome painting techniques. Three-way analysis of rearrangements leads to a reliable reconstruction of the genomic architecture of specific regions in the murid ancestor, including the X chromosome, and for the first time allows one to assign major rearrangement events to one of human, mouse, and rat lineages. Our analysis implies that the rate of rearrangements is much higher in murid rodents than in the human lineage and confirms the existence of rearrangement hot-spots in all three lineages.