Piwil2 (Mili) sustains neurogenesis and prevents cellular senescence in the postnatal hippocampus.

Adult neural progenitor cells (aNPCs) ensure lifelong neurogenesis in the mammalian hippocampus. Proper regulation of aNPC fate has thus important implications for brain plasticity and healthy aging. Piwi proteins and the small noncoding RNAs interacting with them (piRNAs) have been proposed to control memory and anxiety, but the mechanism remains elusive. Here, we show that Piwil2 (Mili) is essential for proper neurogenesis in the postnatal mouse hippocampus. RNA sequencing of aNPCs and their differentiated progeny reveal that Mili and piRNAs are dynamically expressed in neurogenesis. Depletion of Mili and piRNAs in the adult hippocampus impairs aNPC differentiation toward a neural fate, induces senescence, and generates reactive glia. Transcripts modulated upon Mili depletion bear sequences complementary or homologous to piRNAs and include repetitive elements and mRNAs encoding essential proteins for proper neurogenesis. Our results provide evidence of a critical role for Mili in maintaining fitness and proper fate of aNPCs, underpinning a possible involvement of the piRNA pathway in brain plasticity and successful aging.

Foxg1 regulates translation of neocortical neuronal genes, including the main NMDA receptor subunit gene, Grin1.

BACKGROUND: Mainly known as a transcription factor patterning the rostral brain and governing its histogenesis, FOXG1 has been also detected outside the nucleus; however, biological meaning of that has been only partially clarified. RESULTS: Prompted by FOXG1 expression in cytoplasm of pallial neurons, we investigated its implication in translational control. We documented the impact of FOXG1 on ribosomal recruitment of Grin1-mRNA, encoding for the main subunit of NMDA receptor. Next, we showed that FOXG1 increases GRIN1 protein level by enhancing the translation of its mRNA, while not increasing its stability. Molecular mechanisms underlying this activity included FOXG1 interaction with EIF4E and, possibly, Grin1-mRNA. Besides, we found that, within murine neocortical cultures, de novo synthesis of GRIN1 undergoes a prominent and reversible, homeostatic regulation and FOXG1 is instrumental to that. Finally, by integrated analysis of multiple omic data, we inferred that FOXG1 is implicated in translational control of hundreds of neuronal genes, modulating ribosome engagement and progression. In a few selected cases, we experimentally verified such inference. CONCLUSIONS: These findings point to FOXG1 as a key effector, potentially crucial to multi-scale temporal tuning of neocortical pyramid activity, an issue with profound physiological and neuropathological implications.

Paternal starvation affects metabolic gene expression during zebrafish offspring development and lifelong fitness.

Dietary restriction in the form of fasting is a putative key to a healthier and longer life, but these benefits may come at a trade-off with reproductive fitness and may affect the following generation(s). The potential inter- and transgenerational effects of long-term fasting and starvation are particularly poorly understood in vertebrates when they originate from the paternal line. We utilised the externally fertilising zebrafish amenable to a split-egg clutch design to explore the male-specific effects of fasting/starvation on fertility and fitness of offspring independently of maternal contribution. Eighteen days of fasting resulted in reduced fertility in exposed males. While average offspring survival was not affected, we detected increased larval growth rate in F1 offspring from starved males and more malformed embryos at 24 h post-fertilisation in F2 offspring produced by F1 offspring from starved males. Comparing the transcriptomes of F1 embryos sired by starved and fed fathers revealed robust and reproducible increased expression of muscle composition genes but lower expression of lipid metabolism and lysosome genes in embryos from starved fathers. A large proportion of these genes showed enrichment in the yolk syncytial layer suggesting gene regulatory responses associated with metabolism of nutrients through paternal effects on extra-embryonic tissues which are loaded with maternal factors. We compared the embryo transcriptomes to published adult transcriptome datasets and found comparable repressive effects of starvation on metabolism-associated genes. These similarities suggest a physiologically relevant, directed and potentially adaptive response transmitted by the father, independently from the offspring's nutritional state, which was defined by the mother.

TEspeX: consensus-specific quantification of transposable element expression preventing biases from exonized fragments.

SUMMARY: Transposable elements (TEs) play key roles in crucial biological pathways. Therefore, several tools enabling the quantification of their expression were recently developed. However, many of the existing tools lack the capability to distinguish between the transcription of autonomously expressed TEs and TE fragments embedded in canonical coding/non-coding non-TE transcripts. Consequently, an apparent change in the expression of a given TE may simply reflect the variation in the expression of the transcripts containing TE-derived sequences. To overcome this issue, we have developed TEspeX, a pipeline for the quantification of TE expression at the consensus level. TEspeX uses Illumina RNA-seq short reads to quantify TE expression avoiding counting reads deriving from inactive TE fragments embedded in canonical transcripts. AVAILABILITY AND IMPLEMENTATION: The tool is implemented in python3, distributed under the GNU General Public License (GPL) and available on Github at https://github.com/fansalon/TEspeX (Zenodo URL: https://doi.org/10.5281/zenodo.6800331). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Evolutionary genomics of the cold-adapted diatom Fragilariopsis cylindrus.

The Southern Ocean houses a diverse and productive community of organisms. Unicellular eukaryotic diatoms are the main primary producers in this environment, where photosynthesis is limited by low concentrations of dissolved iron and large seasonal fluctuations in light, temperature and the extent of sea ice. How diatoms have adapted to this extreme environment is largely unknown. Here we present insights into the genome evolution of a cold-adapted diatom from the Southern Ocean, Fragilariopsis cylindrus, based on a comparison with temperate diatoms. We find that approximately 24.7 per cent of the diploid F. cylindrus genome consists of genetic loci with alleles that are highly divergent (15.1 megabases of the total genome size of 61.1 megabases). These divergent alleles were differentially expressed across environmental conditions, including darkness, low iron, freezing, elevated temperature and increased CO2. Alleles with the largest ratio of non-synonymous to synonymous nucleotide substitutions also show the most pronounced condition-dependent expression, suggesting a correlation between diversifying selection and allelic differentiation. Divergent alleles may be involved in adaptation to environmental fluctuations in the Southern Ocean.

Identification of LINE retrotransposons and long non-coding RNAs expressed in the octopus brain.

BACKGROUND: Transposable elements (TEs) widely contribute to the evolution of genomes allowing genomic innovations, generating germinal and somatic heterogeneity, and giving birth to long non-coding RNAs (lncRNAs). These features have been associated to the evolution, functioning, and complexity of the nervous system at such a level that somatic retrotransposition of long interspersed element (LINE) L1 has been proposed to be associated to human cognition. Among invertebrates, octopuses are fascinating animals whose nervous system reaches a high level of complexity achieving sophisticated cognitive abilities. The sequencing of the genome of the Octopus bimaculoides revealed a striking expansion of TEs which were proposed to have contributed to the evolution of its complex nervous system. We recently found a similar expansion also in the genome of Octopus vulgaris. However, a specific search for the existence and the transcription of full-length transpositionally competent TEs has not been performed in this genus. RESULTS: Here, we report the identification of LINE elements competent for retrotransposition in Octopus vulgaris and Octopus bimaculoides and show evidence suggesting that they might be transcribed and determine germline and somatic polymorphisms especially in the brain. Transcription and translation measured for one of these elements resulted in specific signals in neurons belonging to areas associated with behavioral plasticity. We also report the transcription of thousands of lncRNAs and the pervasive inclusion of TE fragments in the transcriptomes of both Octopus species, further testifying the crucial activity of TEs in the evolution of the octopus genomes. CONCLUSIONS: The neural transcriptome of the octopus shows the transcription of thousands of putative lncRNAs and of a full-length LINE element belonging to the RTE class. We speculate that a convergent evolutionary process involving retrotransposons activity in the brain has been important for the evolution of sophisticated cognitive abilities in this genus.

Supporting Cells of the Human Olfactory Epithelium Co-Express the Lipid Scramblase TMEM16F and ACE2 and May Cause Smell Loss by SARS-CoV-2 Spike-Induced Syncytia.

BACKGROUND/AIMS: Quantitative and qualitative alterations in the sense of smell are well established symptoms of COVID-19. Some reports have shown that non-neuronal supporting (also named sustentacular) cells of the human olfactory epithelium co-express ACE2 and TMPRSS2 necessary for SARS-CoV-2 infection. In COVID-19, syncytia were found in many tissues but were not investigated in the olfactory epithelium. Some studies have shown that syncytia in some tissues are formed when SARS-CoV-2 Spike expressed at the surface of an infected cell binds to ACE2 on another cell, followed by activation of the scramblase TMEM16F (also named ANO6) which exposes phosphatidylserine to the external side of the membrane. Furthermore, niclosamide, an approved antihelminthic drug, inhibits Spike-induced syncytia by blocking TMEM16F activity. The aim of this study was to investigate if proteins involved in Spike-induced syncytia formation, i.e., ACE2 and TMEM16F, are expressed in the human olfactory epithelium. METHODS: We analysed a publicly available single-cell RNA-seq dataset from human nasal epithelium and performed immunohistochemistry in human nasal tissues from biopsies. RESULTS: We found that ACE2 and TMEM16F are co-expressed both at RNA and protein levels in non-neuronal supporting cells of the human olfactory epithelium. CONCLUSION: Our results provide the first evidence that TMEM16F is expressed in human olfactory supporting cells and indicate that syncytia formation, that could be blocked by niclosamide, is one of the pathogenic mechanisms worth investigating in COVID-19 smell loss.

Type I interferon activation in RAS-associated autoimmune leukoproliferative disease (RALD).

RAS-associated autoimmune leukoproliferative disease (RALD) is a rare immune dysregulation syndrome caused by somatic gain-of-function mutations of either NRAS or KRAS gene in hematopoietic cells. We describe a 27-year-old patient presenting at 5 months of age with recurrent infections and generalized lymphadenopathy who developed a complex multi-organ autoimmune syndrome with hypogammaglobulinemia, partially controlled with oral steroids, hydroxichloroquine, mofetil mycophenolate and IVIG prophylaxis. Activation of type I interferon pathway was observed in peripheral blood. Since 18 years of age, the patient developed regenerative nodular hyperplasia of the liver evolving into hepatopulmonary syndrome. Whole exome sequencing analysis of the peripheral blood DNA showed the NRAS p.Gly13Asp mutation validated as somatic. Our report highlights the possibility of detecting somatic NRAS gene mutations in patients with inflammatory immune dysregulation and type I interferon activation.

Prolyl 3-Hydroxylase 2 Is a Molecular Player of Angiogenesis.

Prolyl 3-hydroxylase 2 (P3H2) catalyzes the post-translational formation of 3-hydroxyproline on collagens, mainly on type IV. Its activity has never been directly associated to angiogenesis. Here, we identified P3H2 gene through a deep-sequencing transcriptome analysis of human umbilical vein endothelial cells (HUVECs) stimulated with vascular endothelial growth factor A (VEGF-A). Differently from many previous studies we carried out the stimulation not on starved HUVECs, but on cells grown to maintain the best condition for their in vitro survival and propagation. We showed that P3H2 is induced by VEGF-A in two primary human endothelial cell lines and that its transcription is modulated by VEGF-A/VEGF receptor 2 (VEGFR-2) signaling pathway through p38 mitogen-activated protein kinase (MAPK). Then, we demonstrated that P3H2, through its activity on type IV Collagen, is essential for angiogenesis properties of endothelial cells in vitro by performing experiments of gain- and loss-of-function. Immunofluorescence studies showed that the overexpression of P3H2 induced a more condensed status of Collagen IV, accompanied by an alignment of the cells along the Collagen IV bundles, so towards an evident pro-angiogenic status. Finally, we found that P3H2 knockdown prevents pathological angiogenesis in vivo, in the model of laser-induced choroid neovascularization. Together these findings reveal that P3H2 is a new molecular player involved in new vessels formation and could be considered as a potential target for anti-angiogenesis therapy.

A human minisatellite hosts an alternative transcription start site for NPRL3 driving its expression in a repeat number-dependent manner.

Minisatellites, also called variable number of tandem repeats (VNTRs), are a class of repetitive elements that may affect gene expression at multiple levels and have been correlated to disease. Their identification and role as expression quantitative trait loci (eQTL) have been limited by their absence in comparative genomic hybridization and single nucleotide polymorphisms arrays. By taking advantage of cap analysis of gene expression (CAGE), we describe a new example of a minisatellite hosting a transcription start site (TSS) which expression is dependent on the repeat number. It is located in the third intron of the gene nitrogen permease regulator like protein 3 (NPRL3). NPRL3 is a component of the GAP activity toward rags 1 protein complex that inhibits mammalian target of rapamycin complex 1 (mTORC1) activity and it is found mutated in familial focal cortical dysplasia and familial focal epilepsy. CAGE tags represent an alternative TSS identifying TAGNPRL3 messenger RNAs (mRNAs). TAGNPRL3 is expressed in red blood cells both at mRNA and protein levels, it interacts with its protein partner NPRL2 and its overexpression inhibits cell proliferation. This study provides an example of a minisatellite that is both a TSS and an eQTL as well as identifies a new VNTR that may modify mTORC1 activity.

Meta-Omics Reveals Genetic Flexibility of Diatom Nitrogen Transporters in Response to Environmental Changes.

Diatoms (Bacillariophyta), one of the most abundant and diverse groups of marine phytoplankton, respond rapidly to the supply of new nutrients, often out-competing other phytoplankton. Herein, we integrated analyses of the evolution, distribution, and expression modulation of two gene families involved in diatom nitrogen uptake (DiAMT1 and DiNRT2), in order to infer the main drivers of divergence in a key functional trait of phytoplankton. Our results suggest that major steps in the evolution of the two gene families reflected key events triggering diatom radiation and diversification. Their expression is modulated in the contemporary ocean by seawater temperature, nitrate, and iron concentrations. Moreover, the differences in diversity and expression of these gene families throughout the water column hint at a possible link with bacterial activity. This study represents a proof-of-concept of how a holistic approach may shed light on the functional biology of organisms in their natural environment.

The RNA-binding protein ILF3 binds to transposable element sequences in SINEUP lncRNAs.

Transposable elements (TEs) compose about half of the mammalian genome and, as embedded sequences, up to 40% of long noncoding RNA (lncRNA) transcripts. Embedded TEs may represent functional domains within lncRNAs, providing a structured RNA platform for protein interaction. Here we show the interactome profile of the mouse inverted short interspersed nuclear element (SINE) of subfamily B2 (invSINEB2) alone and embedded in antisense (AS) ubiquitin C-terminal hydrolase L1 (Uchl1), an lncRNA that is AS to Uchl1 gene. AS Uchl1 is the representative member of a functional class of AS lncRNAs, named SINEUPs, in which the invSINEB2 acts as effector domain (ED)-enhancing translation of sense protein-coding mRNAs. By using RNA-interacting domainome technology, we identify the IL enhancer-binding factor 3 (ILF3) as a protein partner of AS Uchl1 RNA. We determine that this interaction is mediated by the RNA-binding motif 2 of ILF3 and the invSINEB2. Furthermore, we show that ILF3 is able to bind a free right Arthrobacter luteus (Alu) monomer sequence, the embedded TE acting as ED in human SINEUPs. Bioinformatic analysis of Encyclopedia of DNA Elements-enhanced cross-linking immunoprecipitation data reveals that ILF3 binds transcribed human SINE sequences at transcriptome-wide levels. We then demonstrate that the embedded TEs modulate AS Uchl1 RNA nuclear localization to an extent moderately influenced by ILF3. This work unveils the existence of a specific interaction between embedded TEs and an RNA-binding protein, strengthening the model of TEs as functional modules in lncRNAs.-Fasolo, F., Patrucco, L., Volpe, M., Bon, C., Peano, C., Mignone, F., Carninci, P., Persichetti, F., Santoro, C., Zucchelli, S., Sblattero, D., Sanges, R., Cotella, D., Gustincich, S. The RNA-binding protein ILF3 binds to transposable element sequences in SINEUP lncRNAs.

De Novo Transcriptome Assembly and Gene Expression Profiling of the Copepod Calanus helgolandicus Feeding on the PUA-Producing Diatom Skeletonema marinoi.

Diatoms are the dominant component of the marine phytoplankton. Several diatoms produce secondary metabolites, namely oxylipins, with teratogenic effects on their main predators, crustacean copepods. Our study reports the de novo assembled transcriptome of the calanoid copepod Calanus helgolandicus feeding on the oxylipin-producing diatom Skeletonema marinoi. Differential expression analysis was also performed between copepod females exposed to the diatom and the control flagellate Prorocentrum minimum, which does not produce oxylipins. Our results showed that transcripts involved in carbohydrate, amino acid, folate and methionine metabolism, embryogenesis, and response to stimulus were differentially expressed in the two conditions. Expression of 27 selected genes belonging to these functional categories was also analyzed by RT-qPCR in C. helgolandicus females exposed to a mixed solution of the oxylipins heptadienal and octadienal at the concentration of 10 µM, 15 µM, and 20 µM. The results confirmed differential expression analysis, with up-regulation of genes involved in stress response and down-regulation of genes associated with folate and methionine metabolism, embryogenesis, and signaling. Overall, we offer new insights on the mechanism of action of oxylipins on maternally-induced embryo abnormality. Our results may also help identify biomarker genes associated with diatom-related reproductive failure in the natural copepod population at sea.

Exploratory analysis of transposable elements expression in the C. elegans early embryo.

BACKGROUND: Transposable Elements (TE) are mobile sequences that make up large portions of eukaryote genomes. The functions they play within the complex cellular architecture are still not clearly understood, but it is becoming evident that TE have a role in several physiological and pathological processes. In particular, it has been shown that TE transcription is necessary for the correct development of mice embryos and that their expression is able to finely modulate transcription of coding and non-coding genes. Moreover, their activity in the central nervous system (CNS) and other tissues has been correlated with the creation of somatic mosaicisms and with pathologies such as neurodevelopmental and neurodegenerative diseases as well as cancers. RESULTS: We analyzed TE expression among different cell types of the Caenorhabditis elegans (C. elegans) early embryo asking if, where and when TE are expressed and whether their expression is correlated with genes playing a role in early embryo development. To answer these questions, we took advantage of a public C. elegans embryonic single-cell RNA-seq (sc-RNAseq) dataset and developed a bioinformatics pipeline able to quantify reads mapping specifically against TE, avoiding counting reads mapping on TE fragments embedded in coding/non-coding transcripts. Our results suggest that i) canonical TE expression analysis tools, which do not discard reads mapping on TE fragments embedded in annotated transcripts, may over-estimate TE expression levels, ii) Long Terminal Repeats (LTR) elements are mostly expressed in undifferentiated cells and might play a role in pluripotency maintenance and activation of the innate immune response, iii) non-LTR are expressed in differentiated cells, in particular in neurons and nervous system-associated tissues, and iv) DNA TE are homogenously expressed throughout the C. elegans early embryo development. CONCLUSIONS: TE expression appears finely modulated in the C. elegans early embryo and different TE classes are expressed in different cell types and stages, suggesting that TE might play diverse functions during early embryo development.

Impact of polymorphic transposable elements on transcription in lymphoblastoid cell lines from public data.

BACKGROUND: Transposable elements (TEs) are DNA sequences able to mobilize themselves and to increase their copy-number in the host genome. In the past, they have been considered mainly selfish DNA without evident functions. Nevertheless, currently they are believed to have been extensively involved in the evolution of primate genomes, especially from a regulatory perspective. Due to their recent activity they are also one of the primary sources of structural variants (SVs) in the human genome. By taking advantage of sequencing technologies and bioinformatics tools, recent surveys uncovered specific TE structural variants (TEVs) that gave rise to polymorphisms in human populations. When combined with RNA-seq data this information provides the opportunity to study the potential impact of TEs on gene expression in human. RESULTS: In this work, we assessed the effects of the presence of specific TEs in cis on the expression of flanking genes by producing associations between polymorphic TEs and flanking gene expression levels in human lymphoblastoid cell lines. By using public data from the 1000 Genome Project and the Geuvadis consortium, we exploited an expression quantitative trait loci (eQTL) approach integrated with additional bioinformatics data mining analyses. We uncovered human loci enriched for common, less common and rare TEVs and identified 323 significant TEV-cis-eQTL associations. SINE-R/VNTR/Alus (SVAs) resulted the TE class with the strongest effects on gene expression. We also unveiled differential functional enrichments on genes associated to TEVs, genes associated to TEV-cis-eQTLs and genes associated to the genomic regions mostly enriched in TEV-cis-eQTLs highlighting, at multiple levels, the impact of TEVs on the host genome. Finally, we also identified polymorphic TEs putatively embedded in transcriptional units, proposing a novel mechanism in which TEVs may mediate individual-specific traits. CONCLUSION: We contributed to unveiling the effect of polymorphic TEs on transcription in lymphoblastoid cell lines.

Identification of sex determination genes and their evolution in Phlebotominae sand flies (Diptera, Nematocera).

BACKGROUND: Phlebotomine sand flies (Diptera, Nematocera) are important vectors of several pathogens, including Leishmania parasites, causing serious diseases of humans and dogs. Despite their importance as disease vectors, most aspects of sand fly biology remain unknown including the molecular basis of their reproduction and sex determination, aspects also relevant for the development of novel vector control strategies. RESULTS: Using comparative genomics/transcriptomics data mining and transcriptional profiling, we identified the sex determining genes in phlebotomine sand flies and proposed the first model for the sex determination cascade of these insects. For all the genes identified, we produced manually curated gene models, developmental gene expression profile and performed evolutionary molecular analysis. We identified and characterized, for the first time in a Nematocera species, the transformer (tra) homolog which exhibits both conserved and novel features. The analysis of the tra locus in sand flies and its expression pattern suggest that this gene is able to autoregulate its own splicing, as observed in the fruit fly Ceratitis capitata and several other insect species. CONCLUSIONS: Our results permit to fill the gap about sex determination in sand flies, contribute to a better understanding of this developmental pathway in Nematocera and open the way for the identification of sex determining orthologs in other species of this important Diptera sub-order. Furthermore, the sex determination genes identified in our work also provide the opportunity of future biotechnological applications to control natural population of sand flies, reducing their impact on public health.

Correction to: Identification of the meiotic toolkit in diatoms and exploration of meiosis-specific SPO11 and RAD51 homologs in the sexual species Pseudo-nitzschia multistriata and Seminavis robusta.

Following the publication of this article [1], the authors reported that the link to Additional file 11 linked to the wrong set of data. The correct supplementary data is provided in this Correction article (Additional file 11).

ClusterScan: simple and generalistic identification of genomic clusters.

Summary: Studies on gene clusters proved to be an excellent source of information to understand genomes evolution and identifying specific metabolic pathways or gene families. Improvements in sequencing methods have resulted in a large increase of sequenced genomes for which cluster annotation could be performed and standardized. Currently available programs are developed to search for specific cluster types and none of them is suitable for a broad range of user-based choices. We have developed ClusterScan which allows identifying clusters of any kind of feature simply based on their genomic coordinates and user-defined categorical annotations. Availability and implementation: The tool is written in Python, distributed under the GNU General Public License (GPL) and available on Github at http://bit.ly/ClusterScan or as Docker image at sangeslab/clusterscan: latest. It is supported through a mailing-list on http://bit.ly/ClusterScanSupport. Supplementary information: Supplementary data are available at Bioinformatics online.

Diatom flagellar genes and their expression during sexual reproduction in Leptocylindrus danicus.

BACKGROUND: Flagella have been lost in the vegetative phase of the diatom life cycle, but they are still present in male gametes of centric species, thereby representing a hallmark of sexual reproduction. This process, besides maintaining and creating new genetic diversity, in diatoms is also fundamental to restore the maximum cell size following its reduction during vegetative division. Nevertheless, sexual reproduction has been demonstrated in a limited number of diatom species, while our understanding of its different phases and of their genetic control is scarce. RESULTS: In the transcriptome of Leptocylindrus danicus, a centric diatom widespread in the world's seas, we identified 22 transcripts related to the flagella development and confirmed synchronous overexpression of 6 flagellum-related genes during the male gamete formation process. These transcripts were mostly absent in the closely related species L. aporus, which does not have sexual reproduction. Among the 22 transcripts, L. danicus showed proteins that belong to the Intra Flagellar Transport (IFT) subcomplex B as well as IFT-A proteins, the latter previously thought to be absent in diatoms. The presence of flagellum-related proteins was also traced in the transcriptomes of several other centric species. Finally, phylogenetic reconstruction of the IFT172 and IFT88 proteins showed that their sequences are conserved across protist species and have evolved similarly to other phylogenetic marker genes. CONCLUSION: Our analysis describes for the first time the diatom flagellar gene set, which appears to be more complete and functional than previously reported based on the genome sequence of the model centric diatom, Thalassiosira pseudonana. This first recognition of the whole set of diatom flagellar genes and of their activation pattern paves the way to a wider recognition of the relevance of sexual reproduction in individual species and in the natural environment.

Targeted transgene integration overcomes variability of position effects in zebrafish.

Zebrafish transgenesis is increasingly popular owing to the optical transparency and external development of embryos, which provide a scalable vertebrate model for in vivo experimentation. The ability to express transgenes in a tightly controlled spatio-temporal pattern is an important prerequisite for exploitation of zebrafish in a wide range of biomedical applications. However, conventional transgenesis methods are plagued by position effects: the regulatory environment of genomic integration sites leads to variation of expression patterns of transgenes driven by engineered cis-regulatory modules. This limitation represents a bottleneck when studying the precise function of cis-regulatory modules and their subtle variants or when various effector proteins are to be expressed for labelling and manipulation of defined sets of cells. Here, we provide evidence for the efficient elimination of variability of position effects by developing a PhiC31 integrase-based targeting method. To detect targeted integration events, a simple phenotype scoring of colour change in the lens of larvae is used. We compared PhiC31-based integration and Tol2 transgenesis in the analysis of the activity of a novel conserved enhancer from the developmentally regulated neural-specific esrrga gene. Reporter expression was highly variable among independent lines generated with Tol2, whereas all lines generated with PhiC31 into a single integration site displayed nearly identical, enhancer-specific reporter expression in brain nuclei. Moreover, we demonstrate that a modified integrase system can also be used for the detection of enhancer activity in transient transgenesis. These results demonstrate the power of the PhiC31-based transgene integration for the annotation and fine analysis of transcriptional regulatory elements and it promises to be a generally desirable tool for a range of applications, which rely on highly reproducible patterns of transgene activity in zebrafish.