Expressed repeat elements improve RT-qPCR normalization across a wide range of zebrafish gene expression studies

Adolescent exposure to tris(1,3-dichloro-2-propyl) phosphate (TCPP) induces reproductive toxicity in zebrafish through hypothalamic-pituitary-gonadal axis disruption

Tris(1,3-dichloro-2-propyl) phosphate (TCPP), a prevalent organophosphorus flame retardant in aquatic environments, has raised significant concerns regarding its ecological risks. This study aims to explore the impacts of TCPP on the reproductive functions of zebrafish and delineate its gender-related toxic mechanisms. By assessing the effects on zebrafish of 10 mg/L TCPP exposure from 30 to 120 days post-fertilization (dpf), we thoroughly evaluated the reproductive capability and endocrine system alterations. Our findings indicated that TCPP exposure disrupted gender differentiation in zebrafish and markedly impaired their reproductive capacity, resulting in decreased egg laying and offspring development quality. Histological analyses of gonadal tissues showed an abnormal increase in immature oocytes in females and a reduction in mature sperm count and spermatogonial structure integrity in males, collectively leading to compromised embryo quality. Additionally, molecular docking results indicated that TCPP showed a strong affinity for estrogen receptors, and TCPP-treated zebrafish exhibited imbalanced sex hormones and increased estrogen receptor expression. Alterations in genes associated with the hypothalamic-pituitary-gonadal (HPG) axis and activation of the steroidogenesis pathway suggested that TCPP targets the HPG axis to regulate sex hormone homeostasis. Tamoxifen (TAM), as a competitive inhibitor of estrogen, exhibited a biphasic effect, as evidenced by the counteraction of TCPP-induced effects in both male and female zebrafish after TAM addition. Overall, our study underscored the gender-dependent reproductive toxicity of TCPP exposure in zebrafish, characterized by diminished reproductive capacity and hormonal disturbances, likely due to interference in the HPG axis and steroidogenesis pathways. These findings emphasize the critical need to consider gender differences in chemical risk assessments for ecosystems and highlight the importance of understanding the mechanisms underlying the effects of chemical pollutants on the reproductive health of aquatic species.

Comparative impact of pristine and aged microplastics with triclosan on lipid metabolism in larval zebrafish: Unveiling the regulatory role of miR-217

The prevalence of microplastics (MPs), especially aged particles, interacting with contaminants like triclosan (TCS), raises concerns about their toxicological effects on aquatic life. This study focused on the impact of aged polyamide (APA) MPs and TCS on zebrafish lipid metabolism. APA MPs, with rougher surfaces and lower hydrophobicity, exhibited reduced TCS adsorption than unaged polyamide (PA) MPs. Co-exposure to PA/APA MPs and TCS resulted in higher TCS accumulation in zebrafish larvae, notably more with PA than APA. Larvae exposed to PA + TCS exhibited greater oxidative stress, disrupted lipid metabolism, and altered insulin pathway genes than those exposed to TCS. However, these negative effects were lessened in the APA + TCS group. Through miRNA-seq and miR-217 microinjection, it was revealed that PA + TCS co-exposure upregulated miR-217, linked to lipid metabolic disorders in zebrafish. Moreover, molecular docking showed stable interactions formed between PA, TCS, and the insulin signaling protein Pik3r2. This study demonstrated that PA and TCS co-exposure significantly inhibited the insulin signaling in zebrafish, triggering lipid metabolism dysregulation mediated by miR-217 upregulation, while APA and TCS co-exposure alleviated these disruptions. This research underscored the ecological and toxicological risks of aged MPs and pollutants in aquatic environments, providing crucial insights into the wider implications of MPs pollution.

Optogenetic induction of chronic glucocorticoid exposure in early-life leads to blunted stress-response in larval zebrafish

Early life stress (ELS) exposure alters stress susceptibility in later life and affects vulnerability to stress-related disorders, but how ELS changes the long-lasting responsiveness of the stress system is not well understood. Zebrafish provides an opportunity to study conserved mechanisms underlying the development and function of the stress response that is regulated largely by the neuroendocrine hypothalamus-pituitary-adrenal/interrenal (HPA/I) axis, with glucocorticoids (GC) as the final effector. In this study, we established a method to chronically elevate endogenous GC levels during early life in larval zebrafish. To this end, we employed an optogenetic actuator, beggiatoa photoactivated adenylyl cyclase, specifically expressed in the interrenal cells of zebrafish and demonstrate that its chronic activation leads to hypercortisolaemia and dampens the acute-stress evoked cortisol levels, across a variety of stressor modalities during early life. This blunting of stress-response was conserved in ontogeny at a later developmental stage. Furthermore, we observe a strong reduction of proopiomelanocortin (pomc)-expression in the pituitary as well as upregulation of fkbp5 gene expression. Going forward, we propose that this model can be leveraged to tease apart the mechanisms underlying developmental programming of the HPA/I axis by early-life GC exposure and its implications for vulnerability and resilience to stress in adulthood.

Syntaxin 18 Defects in Human and Zebrafish Unravel Key Roles in Early Cartilage and Bone Development

SNARE proteins comprise a conserved protein family responsible for catalyzing membrane fusion during vesicle traffic. Syntaxin18 (STX18) is a poorly characterized endoplasmic reticulum (ER)-resident t-SNARE. Recently, together with TANGO1 and SLY1, its involvement was shown in ER to Golgi transport of collagen II during chondrogenesis. We report a fetus with a severe osteochondrodysplasia in whom we identified a homozygous substitution of the highly conserved p.Arg10 to Pro of STX18. CRISPR/Cas9-mediated Stx18 deficiency in zebrafish reveals a crucial role for Stx18 in cartilage and bone development. Furthermore, increased expression of multiple components of the Stx18 SNARE complex and of COPI and COPII proteins suggests that Stx18 deficiency impairs antero- and retrograde vesicular transport in the crispant stx18 zebrafish. Taken together, our studies highlight a new candidate gene for a recessive form of osteochondrodysplasia, thereby possibly broadening the SNAREopathy phenotypic spectrum and opening new doors toward future research avenues. © 2023 American Society for Bone and Mineral Research (ASBMR).

Incorporating Primer Amplification Efficiencies in Quantitative Reverse Transcription Polymerase Chain Reaction Experiments; Considerations for Differential Gene Expression Analyses in Zebrafish

Quantitative reverse transcription polymerase chain reaction (RT-qPCR) is commonly used to measure the mRNA expression of target genes in zebrafish. Gene expression values from RT-qPCR are typically reported as relative fold-changes, and relative quantification of RT-qPCR data incorporates primer amplification efficiency values for each target gene. We describe the influence of the primer amplification efficiency analysis method on RT-qPCR gene expression fold-change calculations. This report describes (1) a sample analysis demonstrating incorporation of primer amplification efficiency into RT-qPCR analysis for comparing gene expression of a gene of interest between two groups when normalized to multiple reference genes, (2) the influence of differences in primer amplification efficiencies between measured genes on gene expression differences calculated from theoretical delta-Cq (dCq) values, and (3) an empirical comparison of the influence of three methods of defining primer amplification efficiency in gene expression analyses (delta-delta-Cq [ddCq], standard curve, LinRegPCR) using mRNA measurements of a set of genes in zebrafish embryonic development. Given the need to account for the influence of primer amplification efficiency along with the simplicity of using software programs (LinRegPCR) to measure primer amplification efficiency from RT-qPCR data, we encourage using empirical measurements of primer amplification efficiency for RT-qPCR analysis of differential gene expression in zebrafish.

Cell differentiation and matrix organization are differentially affected during bone formation in osteogenesis imperfecta zebrafish models with different genetic defects impacting collagen type I structure

Osteogenesis imperfecta (OI) is a family of rare heritable skeletal disorders associated with dominant mutations in the collagen type I encoding genes and recessive defects in proteins involved in collagen type I synthesis and processing and in osteoblast differentiation and activity. Historically, it was believed that the OI bone phenotype was only caused by abnormal collagen type I fibrils in the extracellular matrix, but more recently it became clear that the altered bone cell homeostasis, due to mutant collagen retention, plays a relevant role in modulating disease severity in most of the OI forms and it is correlated to impaired bone cell differentiation. Despite in vitro evidence, in vivo data are missing. To better understand the physiopathology of OI, we used two zebrafish models: Chihuahua (Chi/+), carrying a dominant p.G736D substitution in the α1 chain of collagen type I, and the recessive p3h1-/-, lacking prolyl 3-hydroxylase (P3h1) enzyme. Both models share the delay of collagen type I folding, resulting in its overmodification and partial intracellular retention. The regeneration of the bony caudal fin of Chi/+ and p3h1-/- was employed to investigate the impact of abnormal collagen synthesis on bone cell differentiation. Reduced regenerative ability was evident in both models, but it was associated to impaired osteoblast differentiation and osteoblastogenesis/adipogenesis switch only in Chi/+. On the contrary, reduced osteoclast number and activity were found in both models during regeneration. The dominant OI model showed a more detrimental effect in the extracellular matrix organization. Interestingly, the chemical chaperone 4-phenylbutyrate (4-PBA), known to reduce cellular stress and increase collagen secretion, improved bone formation only in p3h1-/- by favoring caudal fin growth without affecting bone cell markers expression. Taken together, our in vivo data proved the negative impact of structurally abnormal collagen type I on bone formation but revealed a gene mutation-specific effect on bone cell differentiation and matrix organization in OI. These, together with the distinct ability to respond to the chaperone treatment, underline the need for precision medicine approaches to properly treat the disease.

Combined effects of rearing temperature regime (thermocycle vs. constant temperature) during early development and thermal treatment on Nile tilapia (Oreochromis niloticus) sex differentiation

In nature, water temperature experiences daily variations known as thermocycles. Temperature is the main environmental factor that influences sex determination in most teleost fish. The purpose of this study was to examine the effects of rearing temperature (thermocycle (TC) vs. constant (CTE)) on development and a posterior thermal shock throughout the period of sex differentiation of Nile tilapia (Oreochromis niloticus). Embryos and larvae were kept under two temperature regimes: TC of 31 °C:25 °C day:night vs. CTE of 28 °C from 0 to 11 dpf. After this period, the larvae from each group were subjected to either heat treatment (HT, 36 °C for 12 days) or kept under the same rearing temperatures until 23 dpf (Control, C). Then all the groups remained at constant temperature until 270 dpf, when blood and gonads were collected. Larval samples were used to examine the expression of genes related to male (amh, ara, sox9a, dmrt1a) and female (cyp19a1a, foxl2, era) sexual differentiation. In juveniles, sex was characterized by histology, the gonadal expression of the genes involved in the sex steroid synthesis was analyzed by qPCR, and plasma testosterone (T) and estradiol (E₂) levels were analyzed by ELISA. In larvae, daily TCs increased the survival rate against HT and up-regulated the expression of ovarian differentiation genes. In juveniles, TC + C induced a higher proportion of females and higher cyp19a1a expression compared to CTE + C. HT induced changes in the CTE group by up-regulating testicular differentiation genes and down-regulating female promoting genes, which did not occur in the TC group. Juveniles from TC + C group presented a higher proportion of females with higher E₂ and cyp19a1a than CTE + HT. Fish from the CTE + HT group showed a higher percentage of males with highest T and amh. These findings indicate that daily TCs during larval development promote ovarian differentiation and diminish the masculinizing effects of HT.

A tapt1 knock-out zebrafish line with aberrant lens development and impaired vision models human early-onset cataract

Bi-allelic mutations in the gene coding for human trans-membrane anterior-posterior transformation protein 1 (TAPT1) result in a broad phenotypic spectrum, ranging from syndromic disease with severe skeletal and congenital abnormalities to isolated early-onset cataract. We present here the first patient with a frameshift mutation in the TAPT1 gene, resulting in both bilateral early-onset cataract and skeletal abnormalities, in addition to several dysmorphic features, in this way further expanding the phenotypic spectrum associated with TAPT1 mutations. A tapt1a/tapt1b double knock-out (KO) zebrafish model generated by CRISPR/Cas9 gene editing revealed an early larval phenotype with eye malformations, loss of vision, increased photokinetics and hyperpigmentation, without visible skeletal involvement. Ultrastructural analysis of the eyes showed a smaller condensed lens, loss of integrity of the lens capsule with formation of a secondary lens and hyperplasia of the cells in the ganglion and inner plexiform layers of the retina. Transcriptomic analysis pointed to an impaired lens development with aberrant expression of many of the crystallin and other lens-specific genes. Furthermore, the phototransduction and visual perception pathways were found to be significantly disturbed. Differences in light perception are likely the cause of the increased dark photokinetics and generalized hyperpigmentation observed in this zebrafish model. In conclusion, this study validates TAPT1 as a new gene for early-onset cataract and sheds light on its ultrastructural and molecular characteristics.

Zebrafish Tric-b is required for skeletal development and bone cells differentiation

INTRODUCTION

Trimeric intracellular potassium channels TRIC-A and -B are endoplasmic reticulum (ER) integral membrane proteins, involved in the regulation of calcium release mediated by ryanodine (RyRs) and inositol 1,4,5-trisphosphate (IP₃Rs) receptors, respectively. While TRIC-A is mainly expressed in excitable cells, TRIC-B is ubiquitously distributed at moderate level. TRIC-B deficiency causes a dysregulation of calcium flux from the ER, which impacts on multiple collagen specific chaperones and modifying enzymatic activity, leading to a rare form of osteogenesis imperfecta (OI Type XIV). The relevance of TRIC-B on cell homeostasis and the molecular mechanism behind the disease are still unknown.

RESULTS

In this study, we exploited zebrafish to elucidate the role of TRIC-B in skeletal tissue. We demonstrated, for the first time, that tmem38a and tmem38b genes encoding Tric-a and -b, respectively are expressed at early developmental stages in zebrafish, but only the latter has a maternal expression. Two zebrafish mutants for tmem38b were generated by CRISPR/Cas9, one carrying an out of frame mutation introducing a premature stop codon (tmem38b^-/- ) and one with an in frame deletion that removes the highly conserved KEV domain (tmem38b^{Δ120-7/Δ120-7} ). In both models collagen type I is under-modified and partially intracellularly retained in the endoplasmic reticulum, as described in individuals affected by OI type XIV. Tmem38b^-/- showed a mild skeletal phenotype at the late larval and juvenile stages of development whereas tmem38b^{Δ120-7/Δ120-7} bone outcome was limited to a reduced vertebral length at 21 dpf. A caudal fin regeneration study pointed towards impaired activity of osteoblasts and osteoclasts associated with mineralization impairment.

DISCUSSION

Our data support the requirement of Tric-b during early development and for bone cell differentiation.

Hypoplasia of dopaminergic neurons by hypoxia-induced neurotoxicity is associated with disrupted swimming development of larval zebrafish

Hypoxic injury to the developing brain increases the risk of permanent behavioral deficits, but the precise mechanisms of hypoxic injury to the developing nervous system are poorly understood. In this study, we characterized the effects of developmental hypoxia (1% pO₂ from 24 to 48 h post-fertilization, hpf) on diencephalic dopaminergic (DA) neurons in larval zebrafish and the consequences on the development of swimming behavior. Hypoxia reduced the number of diencephalic DA neurons at 48 hpf. Returning zebrafish larvae to normoxia after the hypoxia (i.e., hypoxia-recovery, HR) induced reactive oxygen species (ROS) accumulation. Real-time qPCR results showed that HR caused upregulation of proapoptotic genes, including p53 and caspase3, suggesting the potential for ROS-induced cell death. With HR, we also found an increase in TUNEL-positive DA neurons, a persistent reduction in the number of diencephalic DA neurons, and disrupted swimming development and behavior. Interestingly, post-hypoxia (HR) with the antioxidant N-acetylcysteine partially restored the number of DA neurons and spontaneous swimming behavior, demonstrating potential recovery from hypoxic injury. The present study provides new insights for understanding the mechanisms responsible for motor disability due to developmental hypoxic injury.

AP-3 complex delta subunit gene, ap3d1, regulates melanogenesis and melanophore survival via autophagy in zebrafish (Danio rerio)

Zebrafish are an emerging model organism to study the syndromic albinism disorder, Hermansky–Pudlak syndrome (HPS), due to visible pigment development at 24 hours postfertilization, and conserved melanogenesis mechanisms. We describe crasher, a novel HPS type 10 (HPS10) zebrafish model, with a mutation in AP‐3 complex subunit delta gene, ap3d1. Exon 14 of ap3d1 is overexpressed in crasher mutants, while the expression of ap3d1 as a whole is reduced. ap3d1 knockout in *AB zebrafish recapitulates the mutant crasher phenotype. We show ap3d1 loss‐of‐function mutations cause significant expression changes in the melanogenesis genes, dopachrome tautomerase (dct) and tyrosinase‐related protein 1b (tyrp1b), but not tyrosinase (tyr). Last, Generally Applicable Gene‐set Enrichment (GAGE) analysis suggests autophagy pathway genes are upregulated together in crasher. Treatment with autophagy‐inhibitor, bafilomycin A1, significantly decreases melanophore number in crasher, suggesting ap3d1 promotes melanophore survival by limiting excessive autophagy. crasher is a valuable model to explore the regulation of melanogenesis gene expression and pigmentation disease.

Structure, evolution and expression of zebrafish cartilage oligomeric matrix protein (COMP, TSP5). CRISPR-Cas mutants show a dominant phenotype in myosepta

COMP (Cartilage Oligomeric Matrix Protein), also named thrombospondin-5, is a member of the thrombospondin family of extracellular matrix proteins. It is of clinical relevance, as in humans mutations in COMP lead to chondrodysplasias. The gene encoding zebrafish Comp is located on chromosome 11 in synteny with its mammalian orthologs. Zebrafish Comp has a domain structure identical to that of tetrapod COMP and shares 74% sequence similarity with murine COMP. Zebrafish comp is expressed from 5 hours post fertilization (hpf) on, while the protein is first detectable in somites of 11 hpf embryos. During development and in adults comp is strongly expressed in myosepta, craniofacial tendon and ligaments, around ribs and vertebra, but not in its name-giving tissue cartilage. As in mammals, zebrafish Comp forms pentamers. It is easily extracted from 5 days post fertilization (dpf) whole zebrafish. The lack of Comp expression in zebrafish cartilage implies that its cartilage function evolved recently in tetrapods. The expression in tendon and myosepta may indicate a more fundamental function, as in evolutionary distant Drosophila muscle-specific adhesion to tendon cells requires thrombospondin. A sequence encoding a calcium binding motif within the first TSP type-3 repeat of zebrafish Comp was targeted by CRISPR-Cas. The heterozygous and homozygous mutant Comp zebrafish displayed a patchy irregular Comp staining in 3 dpf myosepta, indicating a dominant phenotype. Electron microscopy revealed that the endoplasmic reticulum of myosepta fibroblasts is not affected in homozygous fish. The disorganized extracellular matrix may indicate that this mutation rather interferes with extracellular matrix assembly, similar to what is seen in a subgroup of chondrodysplasia patients. The early expression and easy detection of mutant Comp in zebrafish points to the potential of using the zebrafish model for large scale screening of small molecules that can improve secretion or function of disease-associated COMP mutants.

Lrp5 Mutant and Crispant Zebrafish Faithfully Model Human Osteoporosis, Establishing the Zebrafish as a Platform for CRISPR-Based Functional Screening of Osteoporosis Candidate Genes

Genomewide association studies (GWAS) have improved our understanding of the genetic architecture of common complex diseases such as osteoporosis. Nevertheless, to attribute functional skeletal contributions of candidate genes to osteoporosis-related traits, there is a need for efficient and cost-effective in vivo functional testing. This can be achieved through CRISPR-based reverse genetic screens, where phenotyping is traditionally performed in stable germline knockout (KO) mutants. Recently it was shown that first-generation (F0) mosaic mutant zebrafish (so-called crispants) recapitulate the phenotype of germline KOs. To demonstrate feasibility of functional validation of osteoporosis candidate genes through crispant screening, we compared a crispant to a stable KO zebrafish model for the lrp5 gene. In humans, recessive loss-of-function mutations in LRP5, a co-receptor in the Wnt signaling pathway, cause osteoporosis-pseudoglioma syndrome. In addition, several GWAS studies identified LRP5 as a major risk locus for osteoporosis-related phenotypes. In this study, we showed that early stage lrp5 KO larvae display decreased notochord mineralization and malformations of the head cartilage. Quantitative micro-computed tomography (micro-CT) scanning and mass-spectrometry element analysis of the adult skeleton revealed decreased vertebral bone volume and bone mineralization, hallmark features of osteoporosis. Furthermore, regenerating fin tissue displayed reduced Wnt signaling activity in lrp5 KO adults. We next compared lrp5 mutants with crispants. Next-generation sequencing analysis of adult crispant tissue revealed a mean out-of-frame mutation rate of 76%, resulting in strongly reduced levels of Lrp5 protein. These crispants generally showed a milder but nonetheless highly comparable skeletal phenotype and a similarly reduced Wnt pathway response compared with lrp5 KO mutants. In conclusion, we show through faithful modeling of LRP5-related primary osteoporosis that crispant screening in zebrafish is a promising approach for rapid functional screening of osteoporosis candidate genes. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Bi-allelic premature truncating variants in LTBP1 cause cutis laxa syndrome

Latent transforming growth factor β (TGFβ)-binding proteins (LTBPs) are microfibril-associated proteins essential for anchoring TGFβ in the extracellular matrix (ECM) as well as for correct assembly of ECM components. Variants in LTBP2, LTBP3, and LTBP4 have been identified in several autosomal recessive Mendelian disorders with skeletal abnormalities with or without impaired development of elastin-rich tissues. Thus far, the human phenotype associated with LTBP1 deficiency has remained enigmatic. In this study, we report homozygous premature truncating LTBP1 variants in eight affected individuals from four unrelated consanguineous families. Affected individuals present with connective tissue features (cutis laxa and inguinal hernia), craniofacial dysmorphology, variable heart defects, and prominent skeletal features (craniosynostosis, short stature, brachydactyly, and syndactyly). In vitro studies on proband-derived dermal fibroblasts indicate distinct molecular mechanisms depending on the position of the variant in LTBP1. C-terminal variants lead to an altered LTBP1 loosely anchored in the microfibrillar network and cause increased ECM deposition in cultured fibroblasts associated with excessive TGFβ growth factor activation and signaling. In contrast, N-terminal truncation results in a loss of LTBP1 that does not alter TGFβ levels or ECM assembly. In vivo validation with two independent zebrafish lines carrying mutations in ltbp1 induce abnormal collagen fibrillogenesis in skin and intervertebral ligaments and ectopic bone formation on the vertebrae. In addition, one of the mutant zebrafish lines shows voluminous and hypo-mineralized vertebrae. Overall, our findings in humans and zebrafish show that LTBP1 function is crucial for skin and bone ECM assembly and homeostasis.

Zebrafish mbnl mutants model physical and molecular phenotypes of myotonic dystrophy

The muscleblind RNA-binding proteins (MBNL1, MBNL2 and MBNL3) are highly conserved across vertebrates and are important regulators of RNA alternative splicing. Loss of MBNL protein function through sequestration by CUG or CCUG RNA repeats is largely responsible for the phenotypes of the human genetic disorder myotonic dystrophy (DM). We generated the first stable zebrafish (Danio rerio) models of DM-associated MBNL loss of function through mutation of the three zebrafish mbnl genes. In contrast to mouse models, zebrafish double and triple homozygous mbnl mutants were viable to adulthood. Zebrafish mbnl mutants displayed disease-relevant physical phenotypes including decreased body size and impaired movement. They also exhibited widespread alternative splicing changes, including the misregulation of many DM-relevant exons. Physical and molecular phenotypes were more severe in compound mbnl mutants than in single mbnl mutants, suggesting partially redundant functions of Mbnl proteins. The high fecundity and larval optical transparency of this complete series of zebrafish mbnl mutants will make them useful for studying DM-related phenotypes and how individual Mbnl proteins contribute to them, and for testing potential therapeutics. This article has an associated First Person interview with the first author of the paper.

Identification and Characterization of Zebrafish Tlr4 Coreceptor Md-2

The zebrafish (Danio rerio) is a powerful model organism for studies of the innate immune system. One apparent difference between human and zebrafish innate immunity is the cellular machinery for LPS sensing. In amniotes, the protein complex formed by TLR4 and myeloid differentiation factor 2 (Tlr4/Md-2) recognizes the bacterial molecule LPS and triggers an inflammatory response. It is believed that zebrafish have neither Md-2 nor Tlr4; Md-2 has not been identified outside of amniotes, whereas the zebrafish tlr4 genes appear to be paralogs, not orthologs, of amniote TLR4s We revisited these conclusions. We identified a zebrafish gene encoding Md-2, ly96 Using single-cell RNA sequencing, we found that ly96 is transcribed in cells that also transcribe genes diagnostic for innate immune cells, including the zebrafish tlr4-like genes. In larval zebrafish, ly96 is expressed in a small number of macrophage-like cells. In a functional assay, zebrafish Md-2 and Tlr4ba form a complex that activates NF-κB signaling in response to LPS. In larval zebrafish ly96 loss-of-function mutations perturbed LPS-induced cytokine production but gave little protection against LPS toxicity. Finally, by analyzing the genomic context of tlr4 genes in 11 jawed vertebrates, we found that tlr4 arose prior to the divergence of teleosts and tetrapods. Thus, an LPS-sensitive Tlr4/Md-2 complex is likely an ancestral feature shared by mammals and zebrafish, rather than a de novo invention on the tetrapod lineage. We hypothesize that zebrafish retain an ancestral, low-sensitivity Tlr4/Md-2 complex that confers LPS responsiveness to a specific subset of innate immune cells.

Rearing temperature conditions (constant vs. thermocycle) affect daily rhythms of thermal tolerance and sensing in zebrafish

In the wild, the environment does not remain constant, but periodically oscillates so that temperature rises in the daytime and drops at night, which generates a daily thermocycle. The effects of thermocycles on thermal tolerance have been previously described in fish. However, the impact of thermocycles on daytime-dependent thermal responses and daily rhythms of temperature tolerance and sensing expression mechanisms remain poorly understood. This study investigates the effects of two rearing conditions: constant (26 °C, C) versus a daily thermocycle (28 °C in the daytime; 24 °C at night, T) on the thermal tolerance response in zebrafish. Thermal tolerance (mortality) was assessed in 4dpf (days post fertilization) zebrafish larvae after acute heat shock (39 °C for 1 h) at two time points: middle of the light phase (ML) or middle of the dark phase (MD). Thermal stress responses were evaluated in adult zebrafish after a 37 °C challenge for 1 h at ML or MD to examine the expression of the heat-shock protein (HSP) (hsp70, hsp90ab1, grp94, hsp90aa1, hspb1, hsp47, cirbp) and transient receptor potential (TRP) channels (trpv4, trpm4a, trpm2, trpa1b) in the brain. Finally, the daily rhythms of gene expression of HSPs and TRPs were measured every 4 h for 24 h. The results revealed the larval mortality rates and the expression induction of most HSPs in adult zebrafish brain reached the highest values in fish reared under constant temperature and subjected to thermal shock at MD. The expression of most HSPs and TRPs was mainly synchronized to the light/dark (LD) cycle, regardless of the temperature regime. Most HSPs involved in hyperthermic challenges displayed diurnal rhythms with their acrophases in phase with warm-sensing thermoTRPs acrophases. The cold-sensing trpa1b peaked in the second half of the light period and slightly shifted toward the dark phase anticipating the acrophase of cirpb, which is involved in hypothermic challenges. These findings indicated that: a) thermal shocks are best tolerated in the daytime; b) the implementation of daily thermocycles during larval development reduces mortality and stress-cellular expression of HSPs to an acute thermal stress at MD; c) daily rhythms need to be considered when discussing physiological responses of thermal sensing and thermotolerance in zebrafish.

Comprehensive evaluation of candidate reference genes for quantitative real-time PCR-based analysis in Caucasian clover

The forage species Caucasian clover (Trifolium ambiguum M. Bieb.), a groundcover plant, is resistant to both cold and drought. However, reference genes for qRT-PCR-based analysis of Caucasian clover are lacking. In this study, 12 reference genes were selected on the basis of transcriptomic data. These genes were used to determine the most stably expressed genes in various organs of Caucasian clover under cold, salt and drought stress for qRT-PCR-based analysis. Reference gene stability was analyzed by geNorm, NormFinder, BestKeeper, the ∆Ct method and RefFinder. Under salt stress, RCD1 and PPIL3 were the most stable reference genes in the leaves, and NLI1 and RCD1 were the most stable references genes in the roots. Under low-temperature stress, APA and EFTu-GTP were the most stable reference genes in the leaves, and the RCD1 and NLI2 genes were highly stable in the roots. Under 10% PEG-6000 stress, NLI1 and NLI2 were highly stable in the leaves, and RCD1 and PPIL3 were the most stable in the roots. Overall, RCD1 and NLI2 were the most stable reference genes in organs under normal conditions and across all samples. The most and least stable reference genes were validated by assessing their appropriateness for normalization via WRKY genes.

Hypomorphic zebrafish models mimic the musculoskeletal phenotype of β4GalT7-deficient Ehlers-Danlos syndrome

β4GalT7 is a transmembrane Golgi enzyme, encoded by B4GALT7, that plays a pivotal role in the proteoglycan linker region formation during proteoglycan biosynthesis. Defects in this enzyme give rise to a rare autosomal recessive form of Ehlers-Danlos syndrome (EDS), currently known as 'spondylodysplastic EDS (spEDS-B4GALT7)'. This EDS subtype is mainly characterized by short stature, hypotonia and skeletal abnormalities, thereby illustrating its pleiotropic importance during human development. Insights into the pathogenic mechanisms underlying this disabling disease are very limited, in part due to the lack of a relevant in vivo model. As the majority of mutations identified in patients with spEDS-B4GALT7 are hypomorphic, we generated zebrafish models with partial loss of B4galt7 function, including different knockdown (morphant) and mosaic knockout (crispant) b4galt7 zebrafish models and studied the morphologic, functional and molecular aspects in embryonic and larval stages. Morphant and crispant zebrafish show highly similar morphological abnormalities in early development including a small, round head, bowed pectoral fins, short body-axis and mild developmental delay. Several craniofacial cartilage and bone structures are absent or strongly misshapen. In addition, the total amount of sulfated glycosaminoglycans is significantly diminished and particularly heparan and chondroitin sulfate proteoglycan levels are greatly reduced. We also show impaired cartilage patterning and loss of chondrocyte organization in a cartilage-specific Tg(Col2a1aBAC:mcherry) zebrafish reporter line. The occurrence of the same abnormalities in the different models confirms these are specifically caused by B4galt7 deficiency. A disturbed actin pattern, along with a lack of muscle tone, was only noted in morphants in which translation of b4galt7 was blocked. In conclusion, we generated the first viable animal models for spEDS-B4GALT7, and show that in early development the human spEDS-B4GALT7 phenotype is faithfully mimicked in these zebrafish models. Our findings underscore a key role for β4GalT7 in early development of cartilage, bone and muscle. These models will lead to a better understanding of spEDS-B4GALT7 and can be used in future efforts focusing on therapeutic applications.

Selection of suitable reference genes for gene expression studies in myxosporean (Myxozoa, Cnidaria) parasites

Myxozoans (Cnidaria: Myxozoa) are an extremely diversified group of endoparasites some of which are causative agents of serious diseases in fish. New methods involving gene expression studies have emerged over the last years to better understand and control myxozoan diseases. Quantitative RT-PCR is the most extensively used approach for gene expression studies. However, the accuracy of the results depends on the normalization of the data to reference genes. We studied the expression of eight commonly used reference genes, adenosylhomocysteinase (AHC1), beta actin (ACTB), eukaryotic translation elongation factor 2 (EF2), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), hypoxanthine-guanine phosphoribosyltransferase 1 (HPRT1), DNA-directed RNA polymerase II (RPB2), 18S ribosomal RNA (18S), 28S ribosomal RNA (28S) across different developmental stages of three myxozoan species, Sphaerospora molnari, Myxobolus cerebralis and Ceratonova shasta, representing the three major myxozoan linages from the largest class Myxosporea. The stable reference genes were identified using four algorithms: geNorm, NormFinder, Bestkeeper and ΔCq method. Additionally, we analyzed transcriptomic data from S. molnari proliferative and spore-forming stages to compare the relative amount of expressed transcripts with the most stable reference genes suggested by RT-qPCR. Our results revealed that GAPDH and EF2 are the most uniformly expressed genes across the different developmental stages of the studied myxozoan species.

Purification of high-quality RNA from a small number of fluorescence activated cell sorted zebrafish cells for RNA sequencing purposes

Background

Transgenic zebrafish lines with the expression of a fluorescent reporter under the control of a cell-type specific promoter, enable transcriptome analysis of FACS sorted cell populations. RNA quality and yield are key determinant factors for accurate expression profiling. Limited cell number and FACS induced cellular stress make RNA isolation of sorted zebrafish cells a delicate process. We aimed to optimize a workflow to extract sufficient amounts of high-quality RNA from a limited number of FACS sorted cells from Tg(fli1a:GFP) zebrafish embryos, which can be used for accurate gene expression analysis.

Results

We evaluated two suitable RNA isolation kits (the RNAqueous micro and the RNeasy plus micro kit) and determined that sorting cells directly into lysis buffer is a critical step for success. For low cell numbers, this ensures direct cell lysis, protects RNA from degradation and results in a higher RNA quality and yield. We showed that this works well up to 0.5× dilution of the lysis buffer with sorted cells. In our sort settings, this corresponded to 30,000 and 75,000 cells for the RNAqueous micro kit and RNeasy plus micro kit respectively. Sorting more cells dilutes the lysis buffer too much and requires the use of a collection buffer. We also demonstrated that an additional genomic DNA removal step after RNA isolation is required to completely clear the RNA from any contaminating genomic DNA. For cDNA synthesis and library preparation, we combined SmartSeq v4 full length cDNA library amplification, Nextera XT tagmentation and sample barcoding. Using this workflow, we were able to generate highly reproducible RNA sequencing results.

Conclusions

The presented optimized workflow enables to generate high quality RNA and allows accurate transcriptome profiling of small populations of sorted zebrafish cells.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5608-2) contains supplementary material, which is available to authorized users.

The validation of Short Interspersed Nuclear Elements (SINEs) as a RT-qPCR normalization strategy in a rodent model for temporal lobe epilepsy

Background

In gene expression studies via RT-qPCR many conclusions are inferred by using reference genes. However, it is generally known that also reference genes could be differentially expressed between various tissue types, experimental conditions and animal models. An increasing amount of studies have been performed to validate the stability of reference genes. In this study, two rodent-specific Short Interspersed Nuclear Elements (SINEs), which are located throughout the transcriptome, were validated and assessed against nine reference genes in a model of Temporal Lobe Epilepsy (TLE). Two different brain regions (i.e. hippocampus and cortex) and two different disease stages (i.e. acute phase and chronic phase) of the systemic kainic acid rat model for TLE were analyzed by performing expression analyses with the geNorm and NormFinder algorithms. Finally, we performed a rank aggregation analysis and validated the reference genes and the rodent-specific SINEs (i.e. B elements) individually via Gfap gene expression.

Results

GeNorm ranked Hprt1, Pgk1 and Ywhaz as the most stable genes in the acute phase, while Gusb and B2m were ranked as the most unstable, being significantly upregulated. The two B elements were ranked as most stable for both brain regions in the chronic phase by geNorm. In contrast, NormFinder ranked the B1 element only once as second best in cortical tissue for the chronic phase. Interestingly, using only one of the two algorithms would have led to skewed conclusions. Finally, the rank aggregation method indicated the use of the B1 element as the best option to normalize target genes, independent of the disease progression and brain region. This result was supported by the expression profile of Gfap.

Conclusion

In this study, we demonstrate the potential of implementing SINEs -notably the B1 element- as a stable normalization factor in a rodent model of TLE, independent of brain region or disease progression.

Modeling Tuberculosis in Mycobacterium marinum Infected Adult Zebrafish

Mycobacterium tuberculosis is currently the deadliest human pathogen causing 1.7 million deaths and 10.4 million infections every year. Exposure to this bacterium causes a wide disease spectrum in humans ranging from a sterilized infection to an actively progressing deadly disease. The most common form is the latent tuberculosis, which is asymptomatic, but has the potential to reactivate into a fulminant disease. Adult zebrafish and its natural pathogen Mycobacterium marinum have recently proven to be an applicable model to study the wide disease spectrum of tuberculosis. Importantly, spontaneous latency and reactivation as well as adaptive immune responses in the context of mycobacterial infection can be studied in this model. In this article, we describe methods for the experimental infection of adult zebrafish, the collection of internal organs for the extraction of nucleic acids for the measurement of mycobacterial loads and host immune responses by quantitative PCR. The in-house-developed, M. marinum-specific qPCR assay is more sensitive than the traditional plating methods as it also detects DNA from non-dividing, dormant or recently dead mycobacteria. As both DNA and RNA are extracted from the same individual, it is possible to study the relationships between the diseased state, and the host and pathogen gene-expression. The adult zebrafish model for tuberculosis thus presents itself as a highly applicable, non-mammalian in vivo system to study host-pathogen interactions.

Evaluation of reference genes for reverse transcription-quantitative PCR assays in organs of zebrafish exposed to glyphosate-based herbicide, Roundup

Roundup is a glyphosate-based herbicide (GBH) widely used in agriculture and may cause toxic effects in non-target organisms. Model organisms, as zebrafish, and analysis of gene expression by reverse transcription-quantitative PCR (RT-qPCR) could be used to better understand the Roundup toxicity. A prerequisite for RT-qPCR is the availability of appropriate reference genes; however, they have not been described for Roundup-exposed fish. The aim of this study was to evaluate the expression stability of six reference genes (rpl8, β-act, gapdh, b2m, ef1α, hprt1) and one expressed repetitive element (hatn10) in organs of males (brain, gill, testis) and females (ovary) of zebrafish exposed to Roundup WG at three concentrations (0.065, 0.65 and 6.5 mg N-(phosphonomethyl) glycine/l) for 7 days. Genes were ranked by geNorm, NormFinder, BestKeeper, Delta C t and RefFinder, and their best combinations were determined by geNorm and NormFinder programs. The two most stable ranked genes were specific to each organ: gill (β-act; rpl8); brain (rpl8; β-act); testis (ef1α; gapdh); and ovary (rpl8; hprt1). The cat transcript level was used to evaluate the effect of normalization with these reference genes. These are the first suitable reference genes described for the analysis of gene expression in organs of Roundup-exposed zebrafish, and will allow investigations of the molecular mechanisms of Roundup toxicity.

Interactions between the circadian clock and TGF-β signaling pathway in zebrafish

Background

TGF-β signaling is a cellular pathway that functions in most cells and has been shown to play a role in multiple processes, such as the immune response, cell differentiation and proliferation. Recent evidence suggests a possible interaction between TGF-β signaling and the molecular circadian oscillator. The current study aims to characterize this interaction in the zebrafish at the molecular and behavioral levels, taking advantage of the early development of a functional circadian clock and the availability of light-entrainable clock-containing cell lines.

Results

Smad3a, a TGF-β signaling-related gene, exhibited a circadian expression pattern throughout the brain of zebrafish larvae. Both pharmacological inhibition and indirect activation of TGF-β signaling in zebrafish Pac-2 cells caused a concentration dependent disruption of rhythmic promoter activity of the core clock gene Per1b. Inhibition of TGF-β signaling in intact zebrafish larvae caused a phase delay in the rhythmic expression of Per1b mRNA. TGF-β inhibition also reversibly disrupted, phase delayed and increased the period of circadian rhythms of locomotor activity in zebrafish larvae.

Conclusions

The current research provides evidence for an interaction between the TGF-β signaling pathway and the circadian clock system at the molecular and behavioral levels, and points to the importance of TGF-β signaling for normal circadian clock function. Future examination of this interaction should contribute to a better understanding of its underlying mechanisms and its influence on a variety of cellular processes including the cell cycle, with possible implications for cancer development and progression.

Expressed repetitive elements are broadly applicable reference targets for normalization of reverse transcription-qPCR data in mice

Reverse transcription quantitative PCR (RT-qPCR) is the gold standard method for gene expression analysis on mRNA level. To remove experimental variation, expression levels of the gene of interest are typically normalized to the expression level of stably expressed endogenous reference genes. Identifying suitable reference genes and determining the optimal number of reference genes should precede each quantification study. Popular reference genes are not necessarily stably expressed in the examined conditions, possibly leading to inaccurate results. Stably and universally expressed repetitive elements (ERE) have previously been shown to be an excellent alternative for normalization using classic reference genes in human and zebrafish samples. Here, we confirm that in mouse tissues, EREs are broadly applicable reference targets for RT-qPCR normalization, provided that the RNA samples undergo a thorough DNase treatment. We identified Orr1a0, Rltr2aiap, and Rltr13a3 as the most stably expressed mouse EREs across six different experimental conditions. Therefore, we propose this set of ERE reference targets as good candidates for normalization of RT-qPCR data in a plethora of conditions. The identification of widely applicable stable mouse RT-qPCR reference targets for normalization has great potential to facilitate future murine gene expression studies and improve the validity of RT-qPCR data.

Identification of transcripts potentially involved in neural tube closure using RNA sequencing

Anencephaly is a fatal human neural tube defect (NTD) in which the anterior neural tube remains open. Zebrafish embryos with reduced Nodal signaling display an open anterior neural tube phenotype that is analogous to anencephaly. Previous work from our laboratory suggests that Nodal signaling acts through induction of the head mesendoderm and mesoderm. Head mesendoderm/mesoderm then, through an unknown mechanism, promotes formation of the polarized neuroepithelium that is capable of undergoing the movements required for closure. We compared the transcriptome of embryos treated with a Nodal signaling inhibitor at sphere stage, which causes NTDs, to embryos treated at 30% epiboly, which does not cause NTDs. This screen identified over 3,000 transcripts with potential roles in anterior neurulation. Expression of several genes encoding components of tight and adherens junctions was significantly reduced, supporting the model that Nodal signaling regulates formation of the neuroepithelium. mRNAs involved in Wnt, FGF, and BMP signaling were also differentially expressed, suggesting these pathways might regulate anterior neurulation. In support of this, we found that pharmacological inhibition of FGF-receptor function causes an open anterior NTD as well as loss of mesodermal derivatives. This suggests that Nodal and FGF signaling both promote anterior neurulation through induction of head mesoderm.

Priming of innate antimycobacterial immunity by heat-killed Listeria monocytogenes induces sterilizing response in the adult zebrafish tuberculosis model

Mycobacterium tuberculosis remains one of the most problematic infectious agents, owing to its highly developed mechanisms to evade host immune responses combined with the increasing emergence of antibiotic resistance. Host-directed therapies aiming to optimize immune responses to improve bacterial eradication or to limit excessive inflammation are a new strategy for the treatment of tuberculosis. In this study, we have established a zebrafish-Mycobacterium marinum natural host-pathogen model system to study induced protective immune responses in mycobacterial infection. We show that priming adult zebrafish with heat-killed Listeria monocytogenes (HKLm) at 1 day prior to M. marinum infection leads to significantly decreased mycobacterial loads in the infected zebrafish. Using rag1^−/− fish, we show that the protective immunity conferred by HKLm priming can be induced through innate immunity alone. At 24 h post-infection, HKLm priming leads to a significant increase in the expression levels of macrophage-expressed gene 1 (mpeg1), tumor necrosis factor α (tnfa) and nitric oxide synthase 2b (nos2b), whereas superoxide dismutase 2 (sod2) expression is downregulated, implying that HKLm priming increases the number of macrophages and boosts intracellular killing mechanisms. The protective effects of HKLm are abolished when the injected material is pretreated with nucleases or proteinase K. Importantly, HKLm priming significantly increases the frequency of clearance of M. marinum infection by evoking sterilizing immunity (25 vs 3.7%, P=0.0021). In this study, immune priming is successfully used to induce sterilizing immunity against mycobacterial infection. This model provides a promising new platform for elucidating the mechanisms underlying sterilizing immunity and to develop host-directed treatment or prevention strategies against tuberculosis.

This article has an associated First Person interview with the first author of the paper.

Summary: Heat-killed Listeria monocytogenes induces immune responses that lead to increased clearance of mycobacterial infection in the adult zebrafish tuberculosis model via innate immune mechanisms.

Evolutionary recruitment of flexible Esrp-dependent splicing programs into diverse embryonic morphogenetic processes

Epithelial-mesenchymal interactions are crucial for the development of numerous animal structures. Thus, unraveling how molecular tools are recruited in different lineages to control interplays between these tissues is key to understanding morphogenetic evolution. Here, we study Esrp genes, which regulate extensive splicing programs and are essential for mammalian organogenesis. We find that Esrp homologs have been independently recruited for the development of multiple structures across deuterostomes. Although Esrp is involved in a wide variety of ontogenetic processes, our results suggest ancient roles in non-neural ectoderm and regulating specific mesenchymal-to-epithelial transitions in deuterostome ancestors. However, consistent with the extensive rewiring of Esrp-dependent splicing programs between phyla, most developmental defects observed in vertebrate mutants are related to other types of morphogenetic processes. This is likely connected to the origin of an event in Fgfr, which was recruited as an Esrp target in stem chordates and subsequently co-opted into the development of many novel traits in vertebrates.

Protective Effects of Otophylloside N on Pentylenetetrazol-Induced Neuronal Injury In vitro and In vivo

Approximately 30% of epileptic patients worldwide are medically unable to control their seizures. In addition, repeated epileptic seizures generally lead to neural damage. Pentylenetetrazol (PTZ) is a clinical circulatory and respiratory stimulant that is experimentally used to mimic epileptic convulsion in epilepsy research. Here, we systematically explore the neuroprotective effects of a pure compound isolated from Cynanchum otophyllum Schneid (Qingyangshen), Otophylloside N (OtoN), against PTZ-induced neuronal injury. We used three models: in vitro primary cortical neurons, in vivo mice, and in vivo zebrafish. Our results revealed that OtoN treatment may attenuate PTZ-induced morphology changes, cell death, LDH efflux in embryonic neuronal cells of C57BL/6J mice, and convulsive behavior in zebrafish. Additionally, our Western blot and RT-PCR results demonstrated that OtoN may attenuate PTZ-induced apoptosis and neuronal activation in neuronal cells, mice, and zebrafish. OtoN may reduce PTZ-induced cleavage of poly ADP-ribose polymerase and upregulation of the Bax/Bcl-2 ratio and decrease the expression level of c-Fos. This study is the first investigation of the neuroprotective effects of OtoN, which might be developed as a novel antiepileptic drug.

Daily Rhythms of the Expression of Key Genes Involved in Steroidogenesis and Gonadal Function in Zebrafish

Fish present daily and seasonal rhythms in spawning and plasmatic levels of steroids that control reproduction. However, the existence of the rhythms of expression of the genes that underlie the endocrine mechanisms responsible for processes such as steroidogenesis and reproduction in fish have still been poorly explored to date. Here we investigated the daily pattern of the expression of key genes involved in sex steroid production that ultimately set the sex ratio in fish. Adult zebrafish were maintained under a 12:12 h light-dark cycle at a constant temperature of 27°C and were sampled every 4 h during a 24-hour cycle. The expression of key genes in the gonads and brains of female and male individuals were analyzed. In gonads, the expression of aromatase (cyp19a1a, ovarian aromatase) and the antimüllerian hormone (amh, testis) was rhythmic, with almost opposite acrophases: ZT 5:13 h (in the light phase) and ZT 15:39 h (at night), respectively. The expression of foxl2 (forkhead box L2) was also rhythmic in the ovary (acrophase located at ZT 5:02 h) and the expression of dmrt1 (doublesex and mab-3-related transcription factor 1) was rhythmic in testes (acrophase at ZT 18:36 h). In the brain, cyp19a1b (brain aromatase) and cyp11b (11beta-hydroxylase) presented daily differences, especially in males, where the expression peaked at night. These results provide the first evidence for marked time-of-the-day-dependent differences in the expression of the genes involved in sex ratio control, which should be considered when investigating processes such as reproduction, sex differentiation and steroidogenesis in fish.

RT-qPCR gene expression analysis in zebrafish: Preanalytical precautions and use of expressed repetitive elements for normalization

Gene expression analysis is increasingly important in many fields of biological research. Understanding patterns of expressed genes is assumed to provide insight into complex regulatory networks and can lead to the identification of genes relevant to specific biological processes, including disease. Among different techniques, reverse transcription quantitative polymerase chain reaction (RT-qPCR) is currently regarded as the gold standard for targeted quantification of RNA gene expression, especially because of its high sensitivity, specificity, accuracy, and precision, and also because of its practical simplicity and processing speed. However, different critical factors can influence the outcome of RT-qPCR studies, including isolation of RNA, reverse transcription to cDNA, and data analysis. These factors need to be addressed in order to obtain biologically meaningful results. In this chapter, we describe how RT-qPCR can be used in a reliable way to successfully study differential gene expression in zebrafish. Hereby, we especially focus on how expressed repetitive elements can be employed as reference targets in zebrafish RT-qPCR studies and how they can further improve the quality of the data.

Identification of Novel Reference Genes Suitable for qRT-PCR Normalization with Respect to the Zebrafish Developmental Stage

Reference genes used in normalizing qRT-PCR data are critical for the accuracy of gene expression analysis. However, many traditional reference genes used in zebrafish early development are not appropriate because of their variable expression levels during embryogenesis. In the present study, we used our previous RNA-Seq dataset to identify novel reference genes suitable for gene expression analysis during zebrafish early developmental stages. We first selected 197 most stably expressed genes from an RNA-Seq dataset (29,291 genes in total), according to the ratio of their maximum to minimum RPKM values. Among the 197 genes, 4 genes with moderate expression levels and the least variation throughout 9 developmental stages were identified as candidate reference genes. Using four independent statistical algorithms (delta-CT, geNorm, BestKeeper and NormFinder), the stability of qRT-PCR expression of these candidates was then evaluated and compared to that of actb1 and actb2, two commonly used zebrafish reference genes. Stability rankings showed that two genes, namely mobk13 (mob4) and lsm12b, were more stable than actb1 and actb2 in most cases. To further test the suitability of mobk13 and lsm12b as novel reference genes, they were used to normalize three well-studied target genes. The results showed that mobk13 and lsm12b were more suitable than actb1 and actb2 with respect to zebrafish early development. We recommend mobk13 and lsm12b as new optimal reference genes for zebrafish qRT-PCR analysis during embryogenesis and early larval stages.

Zebrafish Collagen Type I: Molecular and Biochemical Characterization of the Major Structural Protein in Bone and Skin

Over the last years the zebrafish imposed itself as a powerful model to study skeletal diseases, but a limit to its use is the poor characterization of collagen type I, the most abundant protein in bone and skin. In tetrapods collagen type I is a trimer mainly composed of two α1 chains and one α2 chain, encoded by COL1A1 and COL1A2 genes, respectively. In contrast, in zebrafish three type I collagen genes exist, col1a1a, col1a1b and col1a2 coding for α1(I), α3(I) and α2(I) chains. During embryonic and larval development the three collagen type I genes showed a similar spatio-temporal expression pattern, indicating their co-regulation and interdependence at these stages. In both embryonic and adult tissues, the presence of the three α(I) chains was demonstrated, although in embryos α1(I) was present in two distinct glycosylated states, suggesting a developmental-specific collagen composition. Even though in adult bone, skin and scales equal amounts of α1(I), α3(I) and α2(I) chains are present, the presented data suggest a tissue-specific stoichiometry and/or post-translational modification status for collagen type I. In conclusion, this data will be useful to properly interpret results and insights gained from zebrafish models of skeletal diseases.

The Light Wavelength Affects the Ontogeny of Clock Gene Expression and Activity Rhythms in Zebrafish Larvae

Light plays a key role in synchronizing rhythms and setting the phase of early development. However, to date, little is known about the impact of light wavelengths during the ontogeny of the molecular clock and the behavioural rhythmicity. The aim of this research was to determine the effect of light of different wavelengths (white, blue and red) on the onset of locomotor activity and clock gene (per1b, per2, clock1, bmal1 and dbp) expression rhythms. For this purpose, 4 groups of zebrafish embryo/larvae were raised from 0 to 7 days post-fertilization (dpf) under the following lighting conditions: three groups maintained under light:dark (LD) cycles with white (full visible spectrum, LDW), blue (LDB), or red light (LDR), and one group raised under constant darkness (DD). The results showed that lighting conditions influenced activity rhythms. Larvae were arrhythmic under DD, while under LD cycles they developed wavelength-dependent daily activity rhythms which appeared earlier under LDB (4 dpf) than under LDW or LDR (5 dpf). The results also revealed that development and lighting conditions influenced clock gene expression. While clock1 rhythmic expression appeared in all lighting conditions at 7 dpf, per1b, per2 and dbp showed daily variations already at 3 dpf. Curiously, bmal1 showed consistent rhythmic expression from embryonic stage (0 dpf). Summarizing, the data revealed that daily rhythms appeared earlier in the larvae reared under LDB than in those reared under LDW and LDR. These results emphasize the importance of lighting conditions and wavelengths during early development for the ontogeny of daily rhythms of gene expression and how these rhythms are reflected on the behavioural rhythmicity of zebrafish larvae.