Allele-specific gene expression can underlie altered transcript abundance in zebrafish mutants

Genetic context of transgene insertion can influence neurodevelopment in zebrafish

The Gal4/UAS system is used across model organisms to overexpress target genes in precise cell types and relies on generating transgenic Gal4 driver lines. In zebrafish, the Tg(elavl3:KalTA4) (HuC) Gal4 line drives robust expression in neurons. We observed an increased prevalence of swim bladder defects in Tg(elavl3:KalTA4) zebrafish larvae compared to wildtype siblings, which prompted us to investigate whether transgenic larvae display additional neurobehavioral phenotypes. Tg(elavl3:KalTA4) larvae showed alterations in brain activity, brain morphology, and behavior, including increased hindbrain size and reduced activity of the cerebellum. Bulk RNA-seq analysis revealed dysregulation of the transcriptome and suggested an increased ratio of neuronal progenitor cells compared to differentiated neurons. To understand whether these phenotypes derive from Gal4 toxicity or from positional effects related to transgenesis, we used economical low-pass whole genome sequencing to map the Tol2-mediated insertion site to chromosome eight. Reduced expression of the neighboring gene gadd45ga, a known cell cycle regulator, is consistent with increased proliferation and suggests a role for positional effects. Challenges with creating alternative pan-neuronal lines include the length of the elavl3 promoter (over 8 kb) and random insertion using traditional transgenesis methods. To facilitate the generation of alternative lines, we cloned five neuronal promoters (atp6v0cb, smaller elavl3, rtn1a, sncb, and stmn1b) ranging from 1.7 kb to 4.3 kb and created KalTA4 lines using Tol2 and the phiC31 integrase-based pIGLET system. Our study highlights the importance of using appropriate genetic controls and interrogating potential positional effects in new transgenic lines.

Zebrafish are resilient to the loss of major diacylglycerol acyltransferase enzymes

In zebrafish, maternally deposited yolk is the source of nutrients for embryogenesis prior to digestive system maturation. Yolk nutrients are processed and secreted to the growing organism by an extra-embryonic tissue, the yolk syncytial layer (YSL). The export of lipids from the YSL occurs through the production of triacylglycerol-rich lipoproteins. Here we report that mutations in the triacylglycerol synthesis enzyme, diacylglycerol acyltransferase-2 (Dgat2), cause yolk sac opacity due to aberrant accumulation of cytoplasmic lipid droplets in the YSL. Although triacylglycerol synthesis continues, it is not properly coupled to lipoprotein production as dgat2 mutants produce fewer, smaller, ApoB-containing lipoproteins. Unlike DGAT2-null mice, which are lipopenic and die soon after birth, zebrafish dgat2 mutants are viable, fertile, and exhibit normal mass and adiposity. Residual Dgat activity cannot be explained by the activity of other known Dgat isoenzymes, as dgat1a;dgat1b;dgat2 triple mutants continue to produce YSL lipid droplets and remain viable as adults. Further, the newly identified diacylglycerol acyltransferase, Tmem68, is also not responsible for the residual triacylglycerol synthesis activity. Unlike overexpression of Dgat1a and Dgat1b, monoacylglycerol acyltransferase-3 (Mogat3b) overexpression does not rescue yolk opacity, suggesting it does not possess Dgat activity in the YSL. However, mogat3b;dgat2 double mutants exhibit increased yolk opacity and often have structural alterations of the yolk extension. Quadruple mogat3b;dgat1a;dgat1b;dgat2 mutants either have severely reduced viability and stunted growth or do not survive past 3 days post fertilization, depending on the dgat2 mutant allele present. Our study highlights the remarkable ability of vertebrates to synthesize triacylglycerol through multiple biosynthetic pathways.

Shared and unique consequences of Joubert Syndrome gene dysfunction on the zebrafish central nervous system

Joubert Syndrome (JBTS) is a neurodevelopmental ciliopathy defined by a highly specific midbrain-hindbrain malformation, variably associated with additional neurological features. JBTS displays prominent genetic heterogeneity with >40 causative genes that encode proteins localising to the primary cilium, a sensory organelle that is essential for transduction of signalling pathways during neurodevelopment, among other vital functions. JBTS proteins localise to distinct ciliary subcompartments, suggesting diverse functions in cilium biology. Currently, there is no unifying pathomechanism to explain how dysfunction of such diverse primary cilia-related proteins results in such a highly specific brain abnormality. To identify the shared consequence of JBTS gene dysfunction, we carried out transcriptomic analysis using zebrafish mutants for the JBTS-causative genes cc2d2aw38, cep290fh297, inpp5ezh506, talpid3i264 and togaram1zh510 and the Bardet-Biedl syndrome-causative gene bbs1k742. We identified no commonly dysregulated signalling pathways in these mutants and yet all mutants displayed an enrichment of altered gene sets related to central nervous system function. We found that JBTS mutants have altered primary cilia throughout the brain but do not display abnormal brain morphology. Nonetheless, behavioural analyses revealed reduced locomotion and loss of postural control which, together with the transcriptomic results, hint at underlying abnormalities in neuronal activity and/or neuronal circuit function. These zebrafish models therefore offer the unique opportunity to study the role of primary cilia in neuronal function beyond early patterning, proliferation and differentiation.

Knockout, Knockdown, and the Schrödinger Paradox: Genetic Immunity to Phenotypic Recapitulation in Zebrafish

Previous research has highlighted significant phenotypic discrepancies between knockout and knockdown approaches in zebrafish, raising concerns about the reliability of these methods. However, our study suggests that these differences are not as pronounced as was once believed. By carefully examining the roles of maternal and zygotic gene contributions, we demonstrate that these factors significantly influence phenotypic outcomes, often accounting for the observed discrepancies. Our findings emphasize that morpholinos, despite their potential off-target effects, can be effective tools when used with rigorous controls. We introduce the concept of graded maternal contribution, which explains how the uneven distribution of maternal mRNA and proteins during gametogenesis impacts phenotypic variability. Our research categorizes genes into three types-susceptible, immune, and "Schrödinger" (conditional)-based on their phenotypic expression and interaction with genetic compensation mechanisms. This distinction provides new insights into the paradoxical outcomes observed in genetic studies. Ultimately, our work underscores the importance of considering both maternal and zygotic contributions, alongside rigorous experimental controls, to accurately interpret gene function and the mechanisms underlying disease. This study advocates for the continued use of morpholinos in conjunction with advanced genetic tools like CRISPR/Cas9, stressing the need for a meticulous experimental design to optimize the utility of zebrafish in genetic research and therapeutic development.

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.

Comprehensive analyses of 435 goat transcriptomes provides insight into male reproduction

A systematic analysis of genes related to reproduction is crucial for obtaining a comprehensive understanding of the molecular mechanisms that underlie male reproductive traits in mammals. Here, we utilized 435 goat transcriptome datasets to unveil the testicular tissue-specific genes (TSGs), allele-specific expression (ASE) genes and their uncharacterized transcriptional features related to male goat reproduction. Results showed a total of 1790 TSGs were identified in goat testis, which was the most among all tissues. GO enrichment analyses suggested that testicular TSGs were mainly involved in spermatogenesis, multicellular organism development, spermatid development, and flagellated sperm motility. Subsequently, a total of 95 highly conserved TSGs (HCTSGs), 508 middle conserved TSGs (MCTSGs) and 42 no conserved TSGs (NCTSGs) were identified in goat testis. GO enrichment analyses suggested that the HCTSGs and MCTSGs has a more important association with male reproduction than NCTSGs. Additionally, we identified 644 ASE genes, including 88 tissue-specific ASE (TS-ASE) genes (e.g., FSIP2, TDRD9). GO enrichment analyses indicated that both ASE genes and TS-ASE genes were associated with goat male reproduction. Overall, this study revealed an extensive gene set involved in the regulation of male goat reproduction and their dynamic transcription patterns. Data reported here provide valuable insights for a further improvement of the economic benefits of goats as well as future treatments for male infertility.

The effects of mutations on gene expression and alternative splicing

Understanding the relationship between mutations and their genomic and phenotypic consequences has been a longstanding goal of evolutionary biology. However, few studies have investigated the impact of mutations on gene expression and alternative splicing on the genome-wide scale. In this study, we aim to bridge this knowledge gap by utilizing whole-genome sequencing data and RNA sequencing data from 16 obligately parthenogenetic Daphnia mutant lines to investigate the effects of ethyl methanesulfonate-induced mutations on gene expression and alternative splicing. Using rigorous analyses of mutations, expression changes and alternative splicing, we show that trans-effects are the major contributor to the variance in gene expression and alternative splicing between the wild-type and mutant lines, whereas cis mutations only affected a limited number of genes and do not always alter gene expression. Moreover, we show that there is a significant association between differentially expressed genes and exonic mutations, indicating that exonic mutations are an important driver of altered gene expression.

The EN-TEx resource of multi-tissue personal epigenomes & variant-impact models

Understanding how genetic variants impact molecular phenotypes is a key goal of functional genomics, currently hindered by reliance on a single haploid reference genome. Here, we present the EN-TEx resource of 1,635 open-access datasets from four donors (∼30 tissues × ∼15 assays). The datasets are mapped to matched, diploid genomes with long-read phasing and structural variants, instantiating a catalog of >1 million allele-specific loci. These loci exhibit coordinated activity along haplotypes and are less conserved than corresponding, non-allele-specific ones. Surprisingly, a deep-learning transformer model can predict the allele-specific activity based only on local nucleotide-sequence context, highlighting the importance of transcription-factor-binding motifs particularly sensitive to variants. Furthermore, combining EN-TEx with existing genome annotations reveals strong associations between allele-specific and GWAS loci. It also enables models for transferring known eQTLs to difficult-to-profile tissues (e.g., from skin to heart). Overall, EN-TEx provides rich data and generalizable models for more accurate personal functional genomics.