Zebrafish transposable elements show extensive diversification in age, genomic distribution, and developmental expression.

There is considerable interest in understanding the effect of transposable elements (TEs) on embryonic development. Studies in humans and mice are limited by the difficulty of working with mammalian embryos and by the relative scarcity of active TEs in these organisms. The zebrafish is an outstanding model for the study of vertebrate development, and over half of its genome consists of diverse TEs. However, zebrafish TEs remain poorly characterized. Here we describe the demography and genomic distribution of zebrafish TEs and their expression throughout embryogenesis using bulk and single-cell RNA sequencing data. These results reveal a highly dynamic genomic ecosystem comprising nearly 2000 distinct TE families, which vary in copy number by four orders of magnitude and span a wide range of ages. Longer retroelements tend to be retained in intergenic regions, whereas short interspersed nuclear elements (SINEs) and DNA transposons are more frequently found nearby or within genes. Locus-specific mapping of TE expression reveals extensive TE transcription during development. Although two-thirds of TE transcripts are likely driven by nearby gene promoters, we still observe stage- and tissue-specific expression patterns in self-regulated TEs. Long terminal repeat (LTR) retroelements are most transcriptionally active immediately following zygotic genome activation, whereas DNA transposons are enriched among transcripts expressed in later stages of development. Single-cell analysis reveals several endogenous retroviruses expressed in specific somatic cell lineages. Overall, our study provides a valuable resource for using zebrafish as a model to study the impact of TEs on vertebrate development.

Specific photoreceptor cell fate pathways are differentially altered in NR2E3-associated diseases.

Mutations in NR2E3, a gene encoding an orphan nuclear transcription factor, cause two retinal dystrophies with a distinct phenotype, but the precise role of NR2E3 in rod and cone transcriptional networks remains unclear. To dissect NR2E3 function, we performed scRNA-seq in the retinas of wildtype and two different Nr2e3 mouse models that show phenotypes similar to patients carrying NR2E3 mutations. Our results reveal that rod and cone populations are not homogeneous and can be separated into different sub-classes. We identify a previously unreported cone pathway that generates hybrid cones co-expressing both cone- and rod-related genes. In mutant retinas, this hybrid cone subpopulation is more abundant and includes a subpopulation of rods transitioning towards a cone cell fate. Hybrid photoreceptors with high misexpression of cone- and rod-related genes are prone to regulated necrosis. Overall, our results shed light on the role of NR2E3 in modulating photoreceptor differentiation towards cone and rod fates and explain how different mutations in NR2E3 lead to distinct visual disorders in humans.

A developmentally programmed splicing failure contributes to DNA damage response attenuation during mammalian zygotic genome activation.

Transition from maternal to embryonic transcriptional control is crucial for embryogenesis. However, alternative splicing regulation during this process remains understudied. Using transcriptomic data from human, mouse, and cow preimplantation development, we show that the stage of zygotic genome activation (ZGA) exhibits the highest levels of exon skipping diversity reported for any cell or tissue type. Much of this exon skipping is temporary, leads to disruptive noncanonical isoforms, and occurs in genes enriched for DNA damage response in the three species. Two core spliceosomal components, Snrpb and Snrpd2, regulate these patterns. These genes have low maternal expression at ZGA and increase sharply thereafter. Microinjection of Snrpb/d2 messenger RNA into mouse zygotes reduces the levels of exon skipping at ZGA and leads to increased p53-mediated DNA damage response. We propose that mammalian embryos undergo an evolutionarily conserved, developmentally programmed splicing failure at ZGA that contributes to the attenuation of cellular responses to DNA damage.

Repression of endogenous retroviruses prevents antiviral immune response and is required for mammary gland development.

The role of heterochromatin in cell fate specification during development is unclear. We demonstrate that loss of the lysine 9 of histone H3 (H3K9) methyltransferase G9a in the mammary epithelium results in de novo chromatin opening, aberrant formation of the mammary ductal tree, impaired stem cell potential, disrupted intraductal polarity, and loss of tissue function. G9a loss derepresses long terminal repeat (LTR) retroviral sequences (predominantly the ERVK family). Transcriptionally activated endogenous retroviruses generate double-stranded DNA (dsDNA) that triggers an antiviral innate immune response, and knockdown of the cytosolic dsDNA sensor Aim2 in G9a knockout (G9acKO) mammary epithelium rescues mammary ductal invasion. Mammary stem cell transplantation into immunocompromised or G9acKO-conditioned hosts shows partial dependence of the G9acKO mammary morphological defects on the inflammatory milieu of the host mammary fat pad. Thus, altering the chromatin accessibility of retroviral elements disrupts mammary gland development and stem cell activity through both cell-autonomous and non-autonomous mechanisms.

Heterochromatin establishment during early mammalian development is regulated by pericentromeric RNA and characterized by non-repressive H3K9me3.

Following fertilization in mammals, the gametes are reprogrammed to create a totipotent zygote, a process that involves de novo establishment of chromatin domains. A major feature occurring during preimplantation development is the dramatic remodelling of constitutive heterochromatin, although the functional relevance of this is unknown. Here, we show that heterochromatin establishment relies on the stepwise expression and regulated activity of SUV39H enzymes. Enforcing precocious acquisition of constitutive heterochromatin results in compromised development and epigenetic reprogramming, which demonstrates that heterochromatin remodelling is essential for natural reprogramming at fertilization. We find that de novo H3K9 trimethylation (H3K9me3) in the paternal pronucleus after fertilization is catalysed by SUV39H2 and that pericentromeric RNAs inhibit SUV39H2 activity and reduce H3K9me3. De novo H3K9me3 is initially non-repressive for gene expression, but instead bookmarks promoters for compaction. Overall, we uncover the functional importance for the restricted transmission of constitutive heterochromatin during reprogramming and a non-repressive role for H3K9me3.

Natural and laboratory mutations in kuzbanian are associated with zinc stress phenotypes in Drosophila melanogaster.

Organisms must cope with altered environmental conditions such as high concentrations of heavy metals. Stress response to heavy metals is mediated by the metal-responsive transcription factor 1 (MTF-1), which is conserved from Drosophila to humans. MTF-1 binds to metal response elements (MREs) and changes the expression of target genes. kuzbanian (kuz), a metalloendopeptidase that activates the evolutionary conserved Notch signaling pathway, has been identified as an MTF-1 target gene. We have previously identified a putatively adaptive transposable element in the Drosophila melanogaster genome, named FBti0019170, inserted in a kuz intron. In this work, we investigated whether a laboratory mutant stock overexpressing kuz is associated with zinc stress phenotypes. We found that both embryos and adult flies overexpressing kuz are more tolerant to zinc compared with wild-type flies. On the other hand, we found that the effect of FBti0019170 on zinc stress tolerance depends on developmental stage and genetic background. Moreover, in the majority of the genetic backgrounds analyzed, FBti0019170 has a deleterious effect in unpolluted environments in pre-adult stages. These results highlight the complexity of natural mutations and suggest that besides laboratory mutations, natural mutations should be studied in order to accurately characterize gene function and evolution.