Giant lungfish genome elucidates the conquest of land by vertebrates.

Lungfishes belong to lobe-fined fish (Sarcopterygii) that, in the Devonian period, 'conquered' the land and ultimately gave rise to all land vertebrates, including humans1-3. Here we determine the chromosome-quality genome of the Australian lungfish (Neoceratodus forsteri), which is known to have the largest genome of any animal. The vast size of this genome, which is about 14× larger than that of humans, is attributable mostly to huge intergenic regions and introns with high repeat content (around 90%), the components of which resemble those of tetrapods (comprising mainly long interspersed nuclear elements) more than they do those of ray-finned fish. The lungfish genome continues to expand independently (its transposable elements are still active), through mechanisms different to those of the enormous genomes of salamanders. The 17 fully assembled lungfish macrochromosomes maintain synteny to other vertebrate chromosomes, and all microchromosomes maintain conserved ancient homology with the ancestral vertebrate karyotype. Our phylogenomic analyses confirm previous reports that lungfish occupy a key evolutionary position as the closest living relatives to tetrapods4,5, underscoring the importance of lungfish for understanding innovations associated with terrestrialization. Lungfish preadaptations to living on land include the gain of limb-like expression in developmental genes such as hoxc13 and sall1 in their lobed fins. Increased rates of evolution and the duplication of genes associated with obligate air-breathing, such as lung surfactants and the expansion of odorant receptor gene families (which encode proteins involved in detecting airborne odours), contribute to the tetrapod-like biology of lungfishes. These findings advance our understanding of this major transition during vertebrate evolution.

Taxonomic and metabolic diversity of Actinomycetota isolated from faeces of a 28,000-year-old mammoth.

Ancient environmental samples, including permafrost soils and frozen animal remains, represent an archive with microbial communities that have barely been explored. This yet unexplored microbial world is a genetic resource that may provide us with new evolutionary insights into recent genomic changes, as well as novel metabolic pathways and chemistry. Here, we describe Actinomycetota Micromonospora, Oerskovia, Saccharopolyspora, Sanguibacter and Streptomyces species were successfully revived and their genome sequences resolved. Surprisingly, the genomes of these bacteria from an ancient source show a large phylogenetic distance to known strains and harbour many novel biosynthetic gene clusters that may well represent uncharacterised biosynthetic potential. Metabolic profiles of the strains display the production of known molecules like antimycin, conglobatin and macrotetrolides, but the majority of the mass features could not be dereplicated. Our work provides insights into Actinomycetota isolated from an ancient source, yielding unexplored genomic information that is not yet present in current databases.

Ultrahigh field diffusion magnetic resonance imaging uncovers intriguing microstructural changes in the adult zebrafish brain caused by Toll-like receptor 2 genomic deletion.

Toll-like receptor 2 (TLR2) belongs to the TLR protein family that plays an important role in the immune and inflammation response system. While TLR2 is predominantly expressed in immune cells, its expression has also been detected in the brain, specifically in microglia and astrocytes. Recent studies indicate that genomic deletion of TLR2 can result in impaired neurobehavioural function. It is currently not clear if the genomic deletion of TLR2 leads to any alterations in the microstructural features of the brain. In the current study, we noninvasively assess microstructural changes in the brain of TLR2-deficient (tlr2-/-) zebrafish using state-of-the art magnetic resonance imaging (MRI) methods at ultrahigh magnetic field strength (17.6 T). A significant increase in cortical thickness and an overall trend towards increased brain volumes were observed in young tlr2-/- zebrafish. An elevated T2 relaxation time and significantly reduced apparent diffusion coefficient (ADC) unveil brain-wide microstructural alterations, potentially indicative of cytotoxic oedema and astrogliosis in the tlr2-/- zebrafish. Multicomponent analysis of the ADC diffusivity signal by the phasor approach shows an increase in the slow ADC component associated with restricted diffusion. Diffusion tensor imaging and diffusion kurtosis imaging analysis revealed diminished diffusivity and enhanced kurtosis in various white matter tracks in tlr2-/- compared with control zebrafish, identifying the microstructural underpinnings associated with compromised white matter integrity and axonal degeneration. Taken together, our findings demonstrate that the genomic deletion of TLR2 results in severe alterations to the microstructural features of the zebrafish brain. This study also highlights the potential of ultrahigh field diffusion MRI techniques in discerning exceptionally fine microstructural details within the small zebrafish brain, offering potential for investigating microstructural changes in zebrafish models of various brain diseases.

Unveiling the Exquisite Microstructural Details in Zebrafish Brain Non-Invasively Using Magnetic Resonance Imaging at 28.2 T.

Zebrafish (Danio rerio) is an important animal model for a wide range of neurodegenerative diseases. However, obtaining the cellular resolution that is essential for studying the zebrafish brain remains challenging as it requires high spatial resolution and signal-to-noise ratios (SNR). In the current study, we present the first MRI results of the zebrafish brain at the state-of-the-art magnetic field strength of 28.2 T. The performance of MRI at 28.2 T was compared to 17.6 T. A 20% improvement in SNR was observed at 28.2 T as compared to 17.6 T. Excellent contrast, resolution, and SNR allowed the identification of several brain structures. The normative T1 and T2 relaxation values were established over different zebrafish brain structures at 28.2 T. To zoom into the white matter structures, we applied diffusion tensor imaging (DTI) and obtained axial, radial, and mean diffusivity, as well as fractional anisotropy, at a very high spatial resolution. Visualisation of white matter structures was achieved by short-track track-density imaging by applying the constrained spherical deconvolution method (stTDI CSD). For the first time, an algorithm for stTDI with multi-shell multi-tissue (msmt) CSD was tested on zebrafish brain data. A significant reduction in false-positive tracks from grey matter signals was observed compared to stTDI with single-shell single-tissue (ssst) CSD. This allowed the non-invasive identification of white matter structures at high resolution and contrast. Our results show that ultra-high field DTI and tractography provide reproducible and quantitative maps of fibre organisation from tiny zebrafish brains, which can be implemented in the future for a mechanistic understanding of disease-related microstructural changes in zebrafish models of various brain diseases.

Thermal Proteome Profiling in Zebrafish Reveals Effects of Napabucasin on Retinoic Acid Metabolism.

Thermal proteome profiling (TPP) allows for the unbiased detection of drug-target protein engagements in vivo. Traditionally, 1 cell type is used for TPP studies, with the risk of missing important differentially expressed target proteins. The use of whole organisms would circumvent this problem. Zebrafish embryos are amenable to such an approach. Here, we used TPP on whole zebrafish embryo lysate to identify protein targets of napabucasin, a compound that may affect signal transducer and activator of transcription 3 (Stat3) signaling through an ill-understood mechanism. In zebrafish embryos, napabucasin induced developmental defects consistent with inhibition of Stat3 signaling. TPP profiling showed no distinct shift in Stat3 upon napabucasin treatment, but effects were detected on the oxidoreductase, Pora, which might explain effects on Stat3 signaling. Interestingly, thermal stability of several aldehyde dehydrogenases was affected. Moreover, napabucasin activated aldehyde dehydrogenase enzymatic activity in vitro. Aldehyde dehydrogenases have crucial roles in retinoic acid metabolism, and functionally, we validated napabucasin-mediated activation of the retinoic acid pathway in zebrafish in vivo. We conclude that TPP profiling in whole zebrafish embryo lysate is feasible and facilitates direct correlation of in vivo effects of small molecule drugs with their protein targets.

Leptin mutation and mycobacterial infection lead non-synergistically to a similar metabolic syndrome.

INTRODUCTION: The leptin signaling pathway plays an important role as a key regulator of glucose homeostasis, metabolism control and systemic inflammatory responses. However, the metabolic effects of leptin on infectious diseases, for example tuberculosis (TB), are still little known. OBJECTIVES: In this study, we aim to investigate the role of leptin on metabolism in the absence and presence of mycobacterial infection in zebrafish larvae and mice. METHODS: Metabolites in entire zebrafish larvae and the blood of mice were studied using high-resolution magic-angle-spinning nuclear magnetic resonance (HR-MAS NMR) spectroscopy and mass spectrometry, respectively. For transcriptome studies of zebrafish larvae, deep RNA sequencing was used. RESULTS: The results show that leptin mutation leads to a similar metabolic syndrome as caused by mycobacterial infection in the two species, characterized by the decrease of 11 amine metabolites. In both species, this metabolic syndrome was not aggravated further when the leptin mutant was infected by mycobacteria. Therefore, we conclude that leptin and mycobacterial infection are both impacting metabolism non-synergistically. In addition, we studied the transcriptomes of lepbibl54 mutant zebrafish larvae and wild type (WT) siblings after mycobacterial infection. These studies showed that mycobacteria induced a very distinct transcriptome signature in the lepbibl54 mutant zebrafish compared to WT sibling control larvae. Furthermore, lepbibl55 Tg (pck1:luc1) zebrafish line was constructed and confirmed this difference in transcriptional responses. CONCLUSIONS: Leptin mutation and TB lead non-synergistically to a similar metabolic syndrome. Moreover, different transcriptomic responses in the lepbibl54  mutant and TB can lead to the similar metabolic end states.

Increased dynamin expression precedes proteinuria in glomerular disease.

Dynamin plays an essential role in maintaining the structure and function of the glomerular filtration barrier. Specifically, dynamin regulates the actin cytoskeleton and the turnover of nephrin in podocytes, and knocking down dynamin expression causes proteinuria. Moreover, promoting dynamin oligomerization with Bis-T-23 restores podocyte function and reduces proteinuria in several animal models of chronic kidney disease. Thus, dynamin is a promising therapeutic target for treating chronic kidney disease. Here, we investigated the pathophysiological role of dynamin under proteinuric circumstances in a rat model and in humans. We found that glomerular Dnm2 and Dnm1 mRNA levels are increased prior to the onset of proteinuria in a rat model of spontaneous proteinuria. Also, in zebrafish embryos, we confirm that knocking down dynamin translation results in proteinuria. Finally, we show that the glomerular expression of dynamin and cathepsin L protein is increased in several human proteinuric kidney diseases. We propose that the increased expression of glomerular dynamin reflects an exhausted attempt to maintain and/or restore integrity of the glomerular filtration barrier. These results confirm that dynamin plays an important role in maintaining the glomerular filtration barrier, and they support the notion that dynamin is a promising therapeutic target in proteinuric kidney disease. © 2018 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Infection and RNA-seq analysis of a zebrafish tlr2 mutant shows a broad function of this toll-like receptor in transcriptional and metabolic control and defense to Mycobacterium marinum infection.

BACKGROUND: The function of Toll-like receptor 2 (TLR2) in host defense against pathogens, especially Mycobacterium tuberculosis (Mtb) is poorly understood. To investigate the role of TLR2 during mycobacterial infection, we analyzed the response of tlr2 zebrafish mutant larvae to infection with Mycobacterium marinum (Mm), a close relative to Mtb, as a model for tuberculosis. We measured infection phenotypes and transcriptome responses using RNA deep sequencing in mutant and control larvae. RESULTS: tlr2 mutant embryos at 2 dpf do not show differences in numbers of macrophages and neutrophils compared to control embryos. However, we found substantial changes in gene expression in these mutants, particularly in metabolic pathways, when compared with the heterozygote tlr2+/- control. At 4 days after Mm infection, the total bacterial burden and the presence of extracellular bacteria were higher in tlr2-/- larvae than in tlr2+/-, or tlr2+/+ larvae, whereas granuloma numbers were reduced, showing a function of Tlr2 in zebrafish host defense. RNAseq analysis of infected tlr2-/- versus tlr2+/- shows that the number of up-regulated and down-regulated genes in response to infection was greatly diminished in tlr2 mutants by at least 2 fold and 10 fold, respectively. Analysis of the transcriptome data and qPCR validation shows that Mm infection of tlr2 mutants leads to decreased mRNA levels of genes involved in inflammation and immune responses, including il1b, tnfb, cxcl11aa/ac, fosl1a, and cebpb. Furthermore, RNAseq analyses revealed that the expression of genes for Maf family transcription factors, vitamin D receptors, and Dicps proteins is altered in tlr2 mutants with or without infection. In addition, the data indicate a function of Tlr2 in the control of induction of cytokines and chemokines, such as the CXCR3-CXCL11 signaling axis. CONCLUSION: The transcriptome and infection burden analyses show a function of Tlr2 as a protective factor against mycobacteria. Transcriptome analysis revealed tlr2-specific pathways involved in Mm infection, which are related to responses to Mtb infection in human macrophages. Considering its dominant function in control of transcriptional processes that govern defense responses and metabolism, the TLR2 protein can be expected to be also of importance for other infectious diseases and interactions with the microbiome.

Glomerular permeability is not affected by heparan sulfate glycosaminoglycan deficiency in zebrafish embryos.

Proteinuria develops when specific components in the glomerular filtration barrier have impaired function. Although the precise components involved in maintaining this barrier have not been fully identified, heparan sulfate proteoglycans are believed to play an essential role in maintaining glomerular filtration. Although in situ studies have shown that a loss of heparan sulfate glycosaminoglycans increases the permeability of the glomerular filtration barrier, recent studies using experimental models have shown that podocyte-specific deletion of heparan sulfate glycosaminoglycan assembly does not lead to proteinuria. However, tubular reabsorption of leaked proteins might have masked an increase in glomerular permeability in these models. Furthermore, not only podocytes but also glomerular endothelial cells are involved in heparan sulfate synthesis in the glomerular filtration barrier. Therefore, we investigated the effect of a global heparan sulfate glycosaminoglycan deficiency on glomerular permeability. We used a zebrafish embryo model carrying a homozygous germline mutation in the ext2 gene. Glomerular permeability was assessed with a quantitative dextran tracer injection method. In this model, we accounted for tubular reabsorption. Loss of anionic sites in the glomerular basement membrane was measured using polyethyleneimine staining. Although mutant animals had significantly fewer negatively charged areas in the glomerular basement membrane, glomerular permeability was unaffected. Moreover, heparan sulfate glycosaminoglycan-deficient embryos had morphologically intact podocyte foot processes. Glomerular filtration remains fully functional despite a global reduction of heparan sulfate.

Mechanistic and Quantitative Understanding of Pharmacokinetics in Zebrafish Larvae through Nanoscale Blood Sampling and Metabolite Modeling of Paracetamol.

Zebrafish larvae are increasingly used for pharmacological research, but internal drug exposure is often not measured. Understanding pharmacokinetics is necessary for reliable translation of pharmacological results to higher vertebrates, including humans. Quantification of drug clearance and distribution requires measurements of blood concentrations. Additionally, measuring drug metabolites is of importance to understand clearance in this model organism mechanistically. We therefore mechanistically studied and quantified pharmacokinetics in zebrafish larvae, and compared this to higher vertebrates, using paracetamol (acetaminophen) as a paradigm compound. A method was developed to sample blood from zebrafish larvae 5 days post fertilization. Blood concentrations of paracetamol and its major metabolites, paracetamol-glucuronide and paracetamol-sulfate, were measured. Blood concentration data were combined with measured amounts in larval homogenates and excreted amounts and simultaneously analyzed through nonlinear mixed-effects modeling, quantifying absolute clearance and distribution volume. Blood sampling from zebrafish larvae was most successful from the posterior cardinal vein, with a median volume (interquartile range) of 1.12 nl (0.676-1.66 nl) per blood sample. Samples were pooled (n = 15-35) to reach measurable levels. Paracetamol blood concentrations at steady state were only 10% of the external paracetamol concentration. Paracetamol-sulfate was the major metabolite, and its formation was quantified using a time-dependent metabolic formation rate. Absolute clearance and distribution volume correlated well with reported values in higher vertebrates, including humans. Based on blood concentrations and advanced data analysis, the mechanistic and quantitative understanding of paracetamol pharmacokinetics in zebrafish larvae has been established. This will improve the translational value of this vertebrate model organism in drug discovery and development. SIGNIFICANCE STATEMENT: In early phases of drug development, new compounds are increasingly screened in zebrafish larvae, but the internal drug exposure is often not taken into consideration. We developed innovative experimental and computational methods, including a blood-sampling technique, to measure the paradigm drug paracetamol (acetaminophen) and its major metabolites and quantify pharmacokinetics (absorption, distribution, elimination) in zebrafish larvae of 5 days post fertilization with a total volume of only 300 nl. These parameter values were scaled to higher vertebrates, including humans.

Expression of distinct maternal and somatic 5.8S, 18S, and 28S rRNA types during zebrafish development.

There is mounting evidence that the ribosome is not a static translation machinery, but a cell-specific, adaptive system. Ribosomal variations have mostly been studied at the protein level, even though the essential transcriptional functions are primarily performed by rRNAs. At the RNA level, oocyte-specific 5S rRNAs are long known for Xenopus. Recently, we described for zebrafish a similar system in which the sole maternal-type 5S rRNA present in eggs is replaced completely during embryonic development by a somatic-type. Here, we report the discovery of an analogous system for the 45S rDNA elements: 5.8S, 18S, and 28S. The maternal-type 5.8S, 18S, and 28S rRNA sequences differ substantially from those of the somatic-type, plus the maternal-type rRNAs are also replaced by the somatic-type rRNAs during embryogenesis. We discuss the structural and functional implications of the observed sequence differences with respect to the translational functions of the 5.8S, 18S, and 28S rRNA elements. Finally, in silico evidence suggests that expansion segments (ES) in 18S rRNA, previously implicated in ribosome-mRNA interaction, may have a preference for interacting with specific mRNA genes. Taken together, our findings indicate that two distinct types of ribosomes exist in zebrafish during development, each likely conducting the translation machinery in a unique way.

Linking maternal and somatic 5S rRNA types with different sequence-specific non-LTR retrotransposons.

5S rRNA is a ribosomal core component, transcribed from many gene copies organized in genomic repeats. Some eukaryotic species have two 5S rRNA types defined by their predominant expression in oogenesis or adult tissue. Our next-generation sequencing study on zebrafish egg, embryo, and adult tissue identified maternal-type 5S rRNA that is exclusively accumulated during oogenesis, replaced throughout the embryogenesis by a somatic-type, and thus virtually absent in adult somatic tissue. The maternal-type 5S rDNA contains several thousands of gene copies on chromosome 4 in tandem repeats with small intergenic regions, whereas the somatic-type is present in only 12 gene copies on chromosome 18 with large intergenic regions. The nine-nucleotide variation between the two 5S rRNA types likely affects TFIII binding and riboprotein L5 binding, probably leading to storage of maternal-type rRNA. Remarkably, these sequence differences are located exactly at the sequence-specific target site for genome integration by the 5S rRNA-specific Mutsu retrotransposon family. Thus, we could define maternal- and somatic-type MutsuDr subfamilies. Furthermore, we identified four additional maternal-type and two new somatic-type MutsuDr subfamilies, each with their own target sequence. This target-site specificity, frequently intact maternal-type retrotransposon elements, plus specific presence of Mutsu retrotransposon RNA and piRNA in egg and adult tissue, suggest an involvement of retrotransposons in achieving the differential copy number of the two types of 5S rDNA loci.

Biological clock function is linked to proactive and reactive personality types.

BACKGROUND: Many physiological processes in our body are controlled by the biological clock and show circadian rhythmicity. It is generally accepted that a robust rhythm is a prerequisite for optimal functioning and that a lack of rhythmicity can contribute to the pathogenesis of various diseases. Here, we tested in a heterogeneous laboratory zebrafish population whether and how variation in the rhythmicity of the biological clock is associated with the coping styles of individual animals, as assessed in a behavioural assay to reliably measure this along a continuum between proactive and reactive extremes. RESULTS: Using RNA sequencing on brain samples, we demonstrated a prominent difference in the expression level of genes involved in the biological clock between proactive and reactive individuals. Subsequently, we tested whether this correlation between gene expression and coping style was due to a consistent change in the level of clock gene expression or to a phase shift or to altered amplitude of the circadian rhythm of gene expression. Our data show a remarkable individual variation in amplitude of the clock gene expression rhythms, which was also reflected in the fluctuating concentrations of melatonin and cortisol, and locomotor activity. This variation in rhythmicity showed a strong correlation with the coping style of the individual, ranging from robust rhythms with large amplitudes in proactive fish to a complete absence of rhythmicity in reactive fish. The rhythmicity of the proactive fish decreased when challenged with constant light conditions whereas the rhythmicity of reactive individuals was not altered. CONCLUSION: These results shed new light on the role of the biological clock by demonstrating that large variation in circadian rhythmicity of individuals may occur within populations. The observed correlation between coping style and circadian rhythmicity suggests that the level of rhythmicity forms an integral part of proactive or reactive coping styles.

COMICS: Cartoon Visualization of Omics Data in Spatial Context Using Anatomical Ontologies.

COMICS is an interactive and open-access web platform for integration and visualization of molecular expression data in anatomograms of zebrafish, carp, and mouse model systems. Anatomical ontologies are used to map omics data across experiments and between an experiment and a particular visualization in a data-dependent manner. COMICS is built on top of several existing resources. Zebrafish and mouse anatomical ontologies with their controlled vocabulary (CV) and defined hierarchy are used with the ontoCAT R package to aggregate data for comparison and visualization. Libraries from the QGIS geographical information system are used with the R packages "maps" and "maptools" to visualize and interact with molecular expression data in anatomical drawings of the model systems. COMICS allows users to upload their own data from omics experiments, using any gene or protein nomenclature they wish, as long as CV terms are used to define anatomical regions or developmental stages. Common nomenclatures such as the ZFIN gene names and UniProt accessions are provided additional support. COMICS can be used to generate publication-quality visualizations of gene and protein expression across experiments. Unlike previous tools that have used anatomical ontologies to interpret imaging data in several animal models, including zebrafish, COMICS is designed to take spatially resolved data generated by dissection or fractionation and display this data in visually clear anatomical representations rather than large data tables. COMICS is optimized for ease-of-use, with a minimalistic web interface and automatic selection of the appropriate visual representation depending on the input data.

Efferocytosis and extrusion of leukocytes determine the progression of early mycobacterial pathogenesis.

Macrophages and neutrophils are the first responders to invading pathogens and contribute strongly to the host defense against intracellular pathogens. The collective interplay and dynamic interactions between these leukocytes are to a large extent not understood. In the present study, we have investigated their role using a combination of confocal laser-scanning and electron microscopy in a zebrafish model for tuberculosis, a local Mycobacterium marinum infection in the tissue of the larval tail fin. Our results show that neutrophils are efficient in phagocytosis of mycobacteria and that they contribute largely to their dissemination. Macrophages appear to play a major role in efferocytosis, phagocytosis of dead cells that contain bacterial content. Phagocytic cells with large bacterial aggregates are formed that can be extruded out of the tissue after cell death. Alternatively, these excessively infected cells can undergo necrosis leading to immediate recruitment of surrounding leukocytes and subsequent phagocytosis of released bacteria. Our data show that these necrotic burst events result in progression of the infection, whereas extrusion abates the infection.

Identifying small RNAs derived from maternal- and somatic-type rRNAs in zebrafish development.

rRNAs are non-coding RNAs present in all prokaryotes and eukaryotes. In eukaryotes there are four rRNAs: 18S, 5.8S, 28S, originating from a common precursor (45S), and 5S. We have recently discovered the existence of two distinct developmental types of rRNA: a maternal-type, present in eggs and a somatic-type, expressed in adult tissues. Lately, next-generation sequencing has allowed the discovery of new small-RNAs deriving from longer non-coding RNAs, including small-RNAs from rRNAs (srRNAs). Here, we systemically investigated srRNAs of maternal- or somatic-type 18S, 5.8S, 28S, with small-RNAseq from many zebrafish developmental stages. We identified new srRNAs for each rRNA. For 5.8S, we found srRNA consisting of the 5' or 3' halves, with only the latter having different sequence for the maternal- and somatic-types. For 18S, we discovered 21 nt srRNA from the 5' end of the 18S rRNA with a striking resemblance to microRNAs; as it is likely processed from a stem-loop precursor and present in human and mouse Argonaute-complexed small-RNA. For 28S, an abundant 80 nt srRNA from the 3' end of the 28S rRNA was found. The expression levels during embryogenesis of these srRNA indicate they are not generated from rRNA degradation and might have a role in the zebrafish development.

Transcriptome dynamics in early zebrafish embryogenesis determined by high-resolution time course analysis of 180 successive, individual zebrafish embryos.

BACKGROUND: Recently, much progress has been made in the field of gene-expression in early embryogenesis. However, the dynamic behaviour of transcriptomes in individual embryos has hardly been studied yet and the time points at which pools of embryos are collected are usually still quite far apart. Here, we present a high-resolution gene-expression time series with 180 individual zebrafish embryos, obtained from nine different spawns, developmentally ordered and profiled from late blastula to mid-gastrula stage. On average one embryo per minute was analysed. The focus was on identification and description of the transcriptome dynamics of the expressed genes in this embryonic stage, rather than to biologically interpret profiles in cellular processes and pathways. RESULTS: In the late blastula to mid-gastrula stage, we found 6,734 genes being expressed with low variability and rather gradual changes. Ten types of dynamic behaviour were defined, such as genes with continuously increasing or decreasing expression, and all expressed genes were grouped into these types. Also, the exact expression starting and stopping points of several hundred genes during this developmental period could be pinpointed. Although the resolution of the experiment was so high, that we were able to clearly identify four known oscillating genes, no genes were observed with a peaking expression. Additionally, several genes showed expression at two or three distinct levels that strongly related to the spawn an embryo originated from. CONCLUSION: Our unique experimental set-up of whole-transcriptome analysis of 180 individual embryos, provided an unparalleled in-depth insight into the dynamics of early zebrafish embryogenesis. The existence of a tightly regulated embryonic transcriptome program, even between individuals from different spawns is shown. We have made the expression profile of all genes available for domain experts. The fact that we were able to separate the different spawns by their gene-expression variance over all expressed genes, underlines the importance of spawn specificity, as well as the unexpectedly tight gene-expression regulation in early zebrafish embryogenesis.

Improving small RNA-seq by using a synthetic spike-in set for size-range quality control together with a set for data normalization.

There is an increasing interest in complementing RNA-seq experiments with small-RNA (sRNA) expression data to obtain a comprehensive view of a transcriptome. Currently, two main experimental challenges concerning sRNA-seq exist: how to check the size distribution of isolated sRNAs, given the sensitive size-selection steps in the protocol; and how to normalize data between samples, given the low complexity of sRNA types. We here present two separate sets of synthetic RNA spike-ins for monitoring size-selection and for performing data normalization in sRNA-seq. The size-range quality control (SRQC) spike-in set, consisting of 11 oligoribonucleotides (10-70 nucleotides), was tested by intentionally altering the size-selection protocol and verified via several comparative experiments. We demonstrate that the SRQC set is useful to reproducibly track down biases in the size-selection in sRNA-seq. The external reference for data-normalization (ERDN) spike-in set, consisting of 19 oligoribonucleotides, was developed for sample-to-sample normalization in differential-expression analysis of sRNA-seq data. Testing and applying the ERDN set showed that it can reproducibly detect differential expression over a dynamic range of 2(18). Hence, biological variation in sRNA composition and content between samples is preserved while technical variation is effectively minimized. Together, both spike-in sets can significantly improve the technical reproducibility of sRNA-seq.

Contrasted innate responses to two viruses in zebrafish: insights into the ancestral repertoire of vertebrate IFN-stimulated genes.

Ease of imaging and abundance of genetic tools make the zebrafish an attractive model host to understand host-pathogen interactions. However, basic knowledge regarding the identity of genes involved in antiviral immune responses is still lagging in this species. We conducted a microarray analysis of the larval zebrafish response to two models of RNA virus infections with very different outcomes. Chikungunya virus (CHIKV) induces a rapid and protective IFN response. Infection with infectious hematopoietic necrosis virus is lethal and is associated with a delayed and inefficient IFN response. A typical signature of IFN-stimulated genes (ISGs) was observed with both viruses, but was stronger for CHIKV. We further compared the zebrafish and human ISG repertoires and made a genomic and phylogenic characterization of the main gene families. We describe a core set of well-induced ISGs conserved across vertebrates, as well as multigenic families diversified independently in each taxon. The conservation of ISGs involved in antiviral signaling indicates conservation of the main feedback loops in these pathways. Whole-mount in situ hybridization of selected transcripts in infected larvae revealed a typical pattern of expression for ISGs in the liver, gut, and blood vessels with both viruses. We further show that some inflammatory genes were additionally induced through IFN-independent pathways by infectious hematopoietic necrosis virus and not by CHIKV. This study provides a useful reference set for the analysis of host-virus interactions in zebrafish and highlights the differences between protective and nonprotective antiviral innate responses.

Changes in ovarian gene expression profiles and plasma hormone levels in maturing European eel (Anguilla anguilla); Biomarkers for broodstock selection.

Complete sexual maturation of European eels (Anguilla anguilla) in captivity can only be achieved via injections with gonadotropins. For female eels this procedure takes 4-6months and the response ranges from "unresponsive" to final maturation and ovulation. Reproductive success could be significantly increased via early selection of responders based on predictive markers and minimally invasive sampling methods. To get a better understanding of the genetic background of ovarian maturation of the European eel we performed a pilot deep-sequencing transcriptome analysis of ovarian tissue derived from a yellow eel, a prepubertal silver eel and a post-spawning matured eel. Two key players in steroidogenesis were strongly correlated with advanced sexual maturation, namely P450c17 and liver receptor homolog-1, suggesting that blood plasma steroids might qualify as minimally invasive markers for early detection of responders. Since the predictive value of plasma sex steroid levels for final maturation of the European eel had not yet been carefully examined, we performed an extensive artificial maturation trial. Farmed silver eels were treated with pituitary extracts and sampled at multiple time intervals. Expression of steroidogenesis-related genes in ovarian tissue of responding and non-responding eels after four weekly injections with pituitary extract was compared using a custom-built microarray and RNAseq. Increased expression of 17ß-hsd1 was strongly linked to sexual maturation. Blood plasma levels of sex steroids were measured using ELISAs. We show that a 2.5-fold increase in blood-plasma estradiol level after 4 weekly pituitary extract injections is a strong predictor of final sexual maturation of female European eel.