Christoffer Nellåker

Mouse genomic variation and its effect on phenotypes and gene regulation.

We report genome sequences of 17 inbred strains of laboratory mice and identify almost ten times more variants than previously known. We use these genomes to explore the phylogenetic history of the laboratory mouse and to examine the functional consequences of allele-specific variation on transcript abundance, revealing that at least 12% of transcripts show a significant tissue-specific expression bias. By identifying candidate functional variants at 718 quantitative trait loci we show that the molecular nature of functional variants and their position relative to genes vary according to the effect size of the locus. These sequences provide a starting point for a new era in the functional analysis of a key model organism.

Sequence-based characterization of structural variation in the mouse genome.

Structural variation is widespread in mammalian genomes and is an important cause of disease, but just how abundant and important structural variants (SVs) are in shaping phenotypic variation remains unclear. Without knowing how many SVs there are, and how they arise, it is difficult to discover what they do. Combining experimental with automated analyses, we identified 711,920 SVs at 281,243 sites in the genomes of thirteen classical and four wild-derived inbred mouse strains. The majority of SVs are less than 1 kilobase in size and 98% are deletions or insertions. The breakpoints of 160,000 SVs were mapped to base pair resolution, allowing us to infer that insertion of retrotransposons causes more than half of SVs. Yet, despite their prevalence, SVs are less likely than other sequence variants to cause gene expression or quantitative phenotypic variation. We identified 24 SVs that disrupt coding exons, acting as rare variants of large effect on gene function. One-third of the genes so affected have immunological functions.

The RNA-editing enzyme ADAR1 controls innate immune responses to RNA.

Summary The ADAR RNA-editing enzymes deaminate adenosine bases to inosines in cellular RNAs. Aberrant interferon expression occurs in patients in whom ADAR1 mutations cause Aicardi-Goutières syndrome (AGS) or dystonia arising from striatal neurodegeneration. Adar1 mutant mouse embryos show aberrant interferon induction and die by embryonic day E12.5. We demonstrate that Adar1 embryonic lethality is rescued to live birth in Adar1; Mavs double mutants in which the antiviral interferon induction response to cytoplasmic double-stranded RNA (dsRNA) is prevented. Aberrant immune responses in Adar1 mutant mouse embryo fibroblasts are dramatically reduced by restoring the expression of editing-active cytoplasmic ADARs. We propose that inosine in cellular RNA inhibits antiviral inflammatory and interferon responses by altering RLR interactions. Transfecting dsRNA oligonucleotides containing inosine-uracil base pairs into Adar1 mutant mouse embryo fibroblasts reduces the aberrant innate immune response. ADAR1 mutations causing AGS affect the activity of the interferon-inducible cytoplasmic isoform more severely than the nuclear isoform.

High levels of RNA-editing site conservation amongst 15 laboratory mouse strains

BackgroundAdenosine-to-inosine (A-to-I) editing is a site-selective post-transcriptional alteration of double-stranded RNA by ADAR deaminases that is crucial for homeostasis and development. Recently the Mouse Genomes Project generated genome sequences for 17 laboratory mouse strains and rich catalogues of variants. We also generated RNA-seq data from whole brain RNA from 15 of the sequenced strains.ResultsHere we present a computational approach that takes an initial set of transcriptome/genome mismatch sites and filters these calls taking into account systematic biases in alignment, single nucleotide variant calling, and sequencing depth to identify RNA editing sites with high accuracy. We applied this approach to our panel of mouse strain transcriptomes identifying 7,389 editing sites with an estimated false-discovery rate of between 2.9 and 10.5%. The overwhelming majority of these edits were of the A-to-I type, with less than 2.4% not of this class, and only three of these edits could not be explained as alignment artifacts. We validated 24 novel RNA editing sites in coding sequence, including two non-synonymous edits in the Cacna1d gene that fell into the IQ domain portion of the Cav1.2 voltage-gated calcium channel, indicating a potential role for editing in the generation of transcript diversity.ConclusionsWe show that despite over two million years of evolutionary divergence, the sites edited and the level of editing at each site is remarkably consistent across the 15 strains. In the Cds2 gene we find evidence for RNA editing acting to preserve the ancestral transcript sequence despite genomic sequence divergence.

The genomic landscape shaped by selection on transposable elements across 18 mouse strains

BackgroundTransposable element (TE)-derived sequence dominates the landscape of mammalian genomes and can modulate gene function by dysregulating transcription and translation. Our current knowledge of TEs in laboratory mouse strains is limited primarily to those present in the C57BL/6J reference genome, with most mouse TEs being drawn from three distinct classes, namely short interspersed nuclear elements (SINEs), long interspersed nuclear elements (LINEs) and the endogenous retrovirus (ERV) superfamily. Despite their high prevalence, the different genomic and gene properties controlling whether TEs are preferentially purged from, or are retained by, genetic drift or positive selection in mammalian genomes remain poorly defined.ResultsUsing whole genome sequencing data from 13 classical laboratory and 4 wild-derived mouse inbred strains, we developed a comprehensive catalogue of 103,798 polymorphic TE variants. We employ this extensive data set to characterize TE variants across the Mus lineage, and to infer neutral and selective processes that have acted over 2 million years. Our results indicate that the majority of TE variants are introduced though the male germline and that only a minority of TE variants exert detectable changes in gene expression. However, among genes with differential expression across the strains there are twice as many TE variants identified as being putative causal variants as expected.ConclusionsMost TE variants that cause gene expression changes appear to be purged rapidly by purifying selection. Our findings demonstrate that past TE insertions have often been highly deleterious, and help to prioritize TE variants according to their likely contribution to gene expression or phenotype variation.

Genetic heterogeneity in Cornelia de Lange syndrome (CdLS) and CdLS-like phenotypes with observed and predicted levels of mosaicism

Background Cornelia de Lange syndrome (CdLS) is a multisystem disorder with distinctive facial appearance, intellectual disability and growth failure as prominent features. Most individuals with typical CdLS have de novo heterozygous loss-of-function mutations in NIPBL with mosaic individuals representing a significant proportion. Mutations in other cohesin components, SMC1A, SMC3, HDAC8 and RAD21 cause less typical CdLS. Methods We screened 163 affected individuals for coding region mutations in the known genes, 90 for genomic rearrangements, 19 for deep intronic variants in NIPBL and 5 had whole-exome sequencing. Results Pathogenic mutations [including mosaic changes] were identified in: NIPBL 46 [3] (28.2%); SMC1A 5 [1] (3.1%); SMC3 5 [1] (3.1%); HDAC8 6 [0] (3.6%) and RAD21 1 [0] (0.6%). One individual had a de novo 1.3 Mb deletion of 1p36.3. Another had a 520 kb duplication of 12q13.13 encompassing ESPL1, encoding separase, an enzyme that cleaves the cohesin ring. Three de novo mutations were identified in ANKRD11 demonstrating a phenotypic overlap with KBG syndrome. To estimate the number of undetected mosaic cases we used recursive partitioning to identify discriminating features in the NIPBL-positive subgroup. Filtering of the mutation-negative group on these features classified at least 18% as ‘NIPBL-like’. A computer composition of the average face of this NIPBL-like subgroup was also more typical in appearance than that of all others in the mutation-negative group supporting the existence of undetected mosaic cases. Conclusions Future diagnostic testing in ‘mutation-negative’ CdLS thus merits deeper sequencing of multiple DNA samples derived from different tissues.

Elevated levels of human endogenous retrovirus-W transcripts in blood cells from patients with first episode schizophrenia

We previously reported on the differential presence of transcripts related to the human endogenous retrovirus (HERV)‐W family in cerebrospinal fluid and plasma from patients with first‐episode schizophrenia compared with control individuals. Whether this is a consequence of qualitative or quantitative differences in transcription of genomic regions harboring HERV‐W elements is not known. The purpose of the present study was therefore to characterize the transcribed HERV‐W elements in mononuclear cells obtained from 30 patients first hospitalized for schizophrenia‐related psychosis and from 26 healthy control individuals. We observed elevated total levels of HERV‐W gag (2.1‐fold, P < 0.01) but not env transcripts in the cells of patients compared with controls. By using the melting temperatures of the amplicons as a proxy marker for sequence identity, no absolute qualitative differences was detected between the two groups. Mapping of the detected transcripts identified several intronic and intergenic HERV‐W elements transcribed in the cells, including elements previously considered transcriptionally silent. Element‐specific assays revealed elevated levels of intronic transcripts containing HERV‐W gag sequence from the putative gene PTD015 on chromosome 11q13.5 (1.6‐fold, P < 0.05) in the patients compared with the controls. Thus, studies aiming to further understanding of complex human disease such as schizophrenia may need to be extended beyond the strictly protein‐coding fraction of the transcriptome.

A systematic evaluation of expression of HERV-W elements; influence of genomic context, viral structure and orientation.

BackgroundOne member of the W family of human endogenous retroviruses (HERV) appears to have been functionally adopted by the human host. Nevertheless, a highly diversified and regulated transcription from a range of HERV-W elements has been observed in human tissues and cells. Aberrant expression of members of this family has also been associated with human disease such as multiple sclerosis (MS) and schizophrenia. It is not known whether this broad expression of HERV-W elements represents transcriptional leakage or specific transcription initiated from the retroviral promoter in the long terminal repeat (LTR) region. Therefore, potential influences of genomic context, structure and orientation on the expression levels of individual HERV-W elements in normal human tissues were systematically investigated.ResultsWhereas intronic HERV-W elements with a pseudogene structure exhibited a strong anti-sense orientation bias, intronic elements with a proviral structure and solo LTRs did not. Although a highly variable expression across tissues and elements was observed, systematic effects of context, structure and orientation were also observed. Elements located in intronic regions appeared to be expressed at higher levels than elements located in intergenic regions. Intronic elements with proviral structures were expressed at higher levels than those elements bearing hallmarks of processed pseudogenes or solo LTRs. Relative to their corresponding genes, intronic elements integrated on the sense strand appeared to be transcribed at higher levels than those integrated on the anti-sense strand. Moreover, the expression of proviral elements appeared to be independent from that of their corresponding genes.ConclusionsIntronic HERV-W provirus integrations on the sense strand appear to have elicited a weaker negative selection than pseudogene integrations of transcripts from such elements. Our current findings suggest that the previously observed diversified and tissue-specific expression of elements in the HERV-W family is the result of both directed transcription (involving both the LTR and internal sequence) and leaky transcription of HERV-W elements in normal human tissues.

Transcriptional Derepression of the ERVWE1 Locus following Influenza A Virus Infection

ABSTRACT Syncytin-1, a fusogenic protein encoded by a human endogenous retrovirus of the W family (HERV-W) element (ERVWE1), is expressed in the syncytiotrophoblast layer of the placenta. This locus is transcriptionally repressed in adult tissues through promoter CpG methylation and suppressive histone modifications. Whereas syncytin-1 appears to be crucial for the development and functioning of the human placenta, its ectopic expression has been associated with pathological conditions, such as multiple sclerosis and schizophrenia. We previously reported on the transactivation of HERV-W elements, including ERVWE1, during influenza A/WSN/33 virus infection in a range of human cell lines. Here we report the results of quantitative PCR analyses of transcripts encoding syncytin-1 in both cell lines and primary fibroblast cells. We observed that spliced ERVWE1 transcripts and those encoding the transcription factor glial cells missing 1 (GCM1), acting as an enhancer element upstream of ERVWE1, are prominently upregulated in response to influenza A/WSN/33 virus infection in nonplacental cells. Knockdown of GCM1 by small interfering RNA followed by infection suppressed the transactivation of ERVWE1. While the infection had no influence on CpG methylation in the ERVWE1 promoter, chromatin immunoprecipitation assays detected decreased H3K9 trimethylation (H3K9me3) and histone methyltransferase SETDB1 levels along with influenza virus proteins associated with ERVWE1 and other HERV-W loci in infected CCF-STTG1 cells. The present findings suggest that an exogenous influenza virus infection can transactivate ERVWE1 by increasing transcription of GCM1 and reducing H3K9me3 in this region and in other regions harboring HERV-W elements. IMPORTANCE Syncytin-1, a protein encoded by the env gene in the HERV-W locus ERVWE1, appears to be crucial for the development and functioning of the human placenta and is transcriptionally repressed in nonplacental tissues. Nevertheless, its ectopic expression has been associated with pathological conditions, such as multiple sclerosis and schizophrenia. In the present paper, we report findings suggesting that an exogenous influenza A virus infection can transactivate ERVWE1 by increasing the transcription of GCM1 and reducing the repressive histone mark H3K9me3 in this region and in other regions harboring HERV-W elements. These observations have implications of potential relevance for viral pathogenesis and for conditions associated with the aberrant transcription of HERV-W loci.

Rapid turnover of functional sequence in human and other genomes.

The amount of a genomes sequence that is functional has been surprisingly difficult to estimate accurately. This has severely hindered analyses asking whether the amount of functional genomic sequence correlates with organismal complexity. Most studies estimate these amounts by considering nucleotide substitution rates within aligned sequences. These approaches show reduced power to identify sequence that is aligned, functional, and constrained only within narrowly defined phyla. The neutral indel model exploits insertions or deletions (indels) rather than substitutions in predicting functional sequence. Surprisingly, this method indicates that half of all functional sequence is specific to individual eutherian lineages. This review considers the rates at which coding or noncoding and functional or nonfunctional sequence changes among mammalian genomes. In contrast to the slow rate at which protein-coding sequence changes, functional noncoding sequence appears to change or be turned over at rapid rates in mammals.