Transcriptome analysis indicates dominant effects on ribosome and mitochondrial function of a premature termination codon mutation in the zebrafish gene psen2

Differential allelic representation (DAR) identifies candidate eQTLs and improves transcriptome analysis

In comparisons between mutant and wild-type genotypes, transcriptome analysis can reveal the direct impacts of a mutation, together with the homeostatic responses of the biological system. Recent studies have highlighted that, when the effects of homozygosity for recessive mutations are studied in non-isogenic backgrounds, genes located proximal to the mutation on the same chromosome often appear over-represented among those genes identified as differentially expressed (DE). One hypothesis suggests that DE genes chromosomally linked to a mutation may not reflect functional responses to the mutation but, instead, result from an unequal distribution of expression quantitative trait loci (eQTLs) between sample groups of mutant or wild-type genotypes. This is problematic because eQTL expression differences are difficult to distinguish from genes that are DE due to functional responses to a mutation. Here we show that chromosomally co-located differentially expressed genes (CC-DEGs) are also observed in analyses of dominant mutations in heterozygotes. We define a method and a metric to quantify, in RNA-sequencing data, localised differential allelic representation (DAR) between those sample groups subjected to differential expression analysis. We show how the DAR metric can predict regions prone to eQTL-driven differential expression, and how it can improve functional enrichment analyses through gene exclusion or weighting-based approaches. Advantageously, this improved ability to identify probable eQTLs also reveals examples of CC-DEGs that are likely to be functionally related to a mutant phenotype. This supports a long-standing prediction that selection for advantageous linkage disequilibrium influences chromosome evolution. By comparing the genomes of zebrafish (Danio rerio) and medaka (Oryzias latipes), a teleost with a conserved ancestral karyotype, we find possible examples of chromosomal aggregation of CC-DEGs during evolution of the zebrafish lineage. Our method for DAR analysis requires only RNA-sequencing data, facilitating its application across new and existing datasets.

The Use of Zebrafish in Transcriptome Analysis of the Early Effects of Mutations Causing Early Onset Familial Alzheimer's Disease and Other Inherited Neurodegenerative Conditions

The degree to which non-human animals can be used to model Alzheimer's disease is a contentious issue, particularly as there is still widespread disagreement regarding the pathogenesis of this neurodegenerative dementia. The currently popular transgenic models are based on artificial expression of genes mutated in early onset forms of familial Alzheimer's disease (EOfAD). Uncertainty regarding the veracity of these models led us to focus on heterozygous, single mutations of endogenous genes (knock-in models) as these most closely resemble the genetic state of humans with EOfAD, and so incorporate the fewest assumptions regarding pathological mechanism. We have generated a number of lines of zebrafish bearing EOfAD-like and non-EOfAD-like mutations in genes equivalent to human PSEN1, PSEN2, and SORL1. To analyze the young adult brain transcriptomes of these mutants, we exploited the ability of zebrafish to produce very large families of simultaneous siblings composed of a variety of genotypes and raised in a uniform environment. This "intra-family" analysis strategy greatly reduced genetic and environmental "noise" thereby allowing detection of subtle changes in gene sets after bulk RNA sequencing of entire brains. Changes to oxidative phosphorylation were predicted for all EOfAD-like mutations in the three genes studied. Here we describe some of the analytical lessons learned in our program combining zebrafish genome editing with transcriptomics to understand the molecular pathologies of neurodegenerative disease.

Differential allelic representation (DAR) identifies candidate eQTLs and improves transcriptome analysis

In comparisons between mutant and wild-type genotypes, transcriptome analysis can reveal the direct impacts of a mutation, together with the homeostatic responses of the biological system. Recent studies have highlighted that, when homozygous mutations are studied in non-isogenic backgrounds, genes from the same chromosome as a mutation often appear over-represented among differentially expressed (DE) genes. One hypothesis suggests that DE genes chromosomally linked to a mutation may not reflect true biological responses to the mutation but, instead, result from differences in representation of expression quantitative trait loci (eQTLs) between sample groups selected on the basis of mutant or wild-type genotype. This is problematic when inclusion of spurious DE genes in a functional enrichment study results in incorrect inferences of mutation effect. Here we show that chromosomally co-located differentially expressed genes (CC-DEGs) can also be observed in analyses of dominant mutations in heterozygotes. We define a method and a metric to quantify, in RNA-sequencing data, localised differential allelic representation (DAR) between groups of samples subject to differential expression analysis. We show how the DAR metric can predict regions prone to eQTL-driven differential expression, and how it can improve functional enrichment analyses through gene exclusion or weighting of gene-level rankings. Advantageously, this improved ability to identify probable eQTLs also reveals examples of CC-DEGs that are likely to be functionally related to a mutant phenotype. This supports a long-standing prediction that selection for advantageous linkage disequilibrium influences chromosome evolution. By comparing the genomes of zebrafish (Danio rerio) and medaka (Oryzias latipes), a teleost with a conserved ancestral karyotype, we find possible examples of chromosomal aggregation of CC-DEGs during evolution of the zebrafish lineage. The DAR metric provides a solid foundation for addressing the eQTL issue in new and existing datasets because it relies solely on RNA-sequencing data.

Brain transcriptomes of zebrafish and mouse Alzheimer's disease knock-in models imply early disrupted energy metabolism

Energy production is the most fundamentally important cellular activity supporting all other functions, particularly in highly active organs, such as brains. Here, we summarise transcriptome analyses of young adult (pre-disease) brains from a collection of 11 early-onset familial Alzheimer's disease (EOFAD)-like and non-EOFAD-like mutations in three zebrafish genes. The one cellular activity consistently predicted as affected by only the EOFAD-like mutations is oxidative phosphorylation, which produces most of the energy of the brain. All the mutations were predicted to affect protein synthesis. We extended our analysis to knock-in mouse models of APOE alleles and found the same effect for the late onset Alzheimer's disease risk allele ε4. Our results support a common molecular basis for the initiation of the pathological processes leading to both early and late onset forms of Alzheimer's disease, and illustrate the utility of zebrafish and knock-in single EOFAD mutation models for understanding the causes of this disease.

Summary: Young adult zebrafish mutants and a mouse model of a genetic variant promoting early- and late-onset Alzheimer's disease, respectively, share changes in brain gene expression, indicating disturbance of oxidative phosphorylation.

PRESENILIN 1 Mutations Causing Early-Onset Familial Alzheimer's Disease or Familial Acne Inversa Differ in Their Effects on Genes Facilitating Energy Metabolism and Signal Transduction

BACKGROUND

The most common cause of early-onset familial Alzheimer's disease (EOfAD) is mutations in PRESENILIN 1 (PSEN1) allowing production of mRNAs encoding full-length, but mutant, proteins. In contrast, a single known frameshift mutation in PSEN1 causes familial acne inversa (fAI) without EOfAD. The molecular consequences of heterozygosity for these mutation types, and how they cause completely different diseases, remains largely unexplored.

OBJECTIVE

To analyze brain transcriptomes of young adult zebrafish to identify similarities and differences in the effects of heterozygosity for psen1 mutations causing EOfAD or fAI.

METHODS

RNA sequencing was performed on mRNA isolated from the brains of a single family of 6-month-old zebrafish siblings either wild type or possessing a single, heterozygous EOfAD-like or fAI-like mutation in their endogenous psen1 gene.

RESULTS

Both mutations downregulate genes encoding ribosomal subunits, and upregulate genes involved in inflammation. Genes involved in energy metabolism appeared significantly affected only by the EOfAD-like mutation, while genes involved in Notch, Wnt and neurotrophin signaling pathways appeared significantly affected only by the fAI-like mutation. However, investigation of direct transcriptional targets of Notch signaling revealed possible increases in γ-secretase activity due to heterozygosity for either psen1 mutation. Transcriptional adaptation due to the fAI-like frameshift mutation was evident.

CONCLUSION

We observed both similar and contrasting effects on brain transcriptomes of the heterozygous EOfAD-like and fAI-like mutations. The contrasting effects may illuminate how these mutation types cause distinct diseases.

In-Frame and Frameshift Mutations in Zebrafish Presenilin 2 Affect Different Cellular Functions in Young Adult Brains

BACKGROUND

Mutations in PRESENILIN 2 (PSEN2) cause early onset familial Alzheimer's disease (EOfAD) but their mode of action remains elusive. One consistent observation for all PRESENILIN gene mutations causing EOfAD is that a transcript is produced with a reading frame terminated by the normal stop codon-the "reading frame preservation rule". Mutations that do not obey this rule do not cause the disease. The reasons for this are debated.

OBJECTIVE

To predict cellular functions affected by heterozygosity for a frameshift, or a reading frame-preserving mutation in zebrafish psen2 using bioinformatic techniques.

METHODS

A frameshift mutation (psen2 N140fs ) and a reading frame-preserving (in-frame) mutation (psen2 T141 _ L142delinsMISLISV ) were previously isolated during genome editing directed at the N140 codon of zebrafish psen2 (equivalent to N141 of human PSEN2). We mated a pair of fish heterozygous for each mutation to generate a family of siblings including wild type and heterozygous mutant genotypes. Transcriptomes from young adult (6 months) brains of these genotypes were analyzed.

RESULTS

The in-frame mutation uniquely caused subtle, but statistically significant, changes to expression of genes involved in oxidative phosphorylation, long-term potentiation and the cell cycle. The frameshift mutation uniquely affected genes involved in Notch and MAPK signaling, extracellular matrix receptor interactions and focal adhesion. Both mutations affected ribosomal protein gene expression but in opposite directions.

CONCLUSION

A frameshift and an in-frame mutation at the same position in zebrafish psen2 cause discrete effects. Changes in oxidative phosphorylation, long-term potentiation and the cell cycle may promote EOfAD pathogenesis in humans.