Gene and genome duplications in vertebrates: the one-to-four (-to-eight in fish) rule and the evolution of novel gene functions

In Vivo Monitoring of Fabp7 Expression in Transgenic Zebrafish

In zebrafish, like in mammals, radial glial cells (RGCs) can act as neural progenitors during development and regeneration in adults. However, the heterogeneity of glia subpopulations entails the need for different specific markers of zebrafish glia. Currently, fluorescent protein expression mediated by a regulatory element from the glial fibrillary acidic protein (gfap) gene is used as a prominent glia reporter. We now expand this tool by demonstrating that a regulatory element from the mouse Fatty acid binding protein 7 (Fabp7) gene drives reliable expression in fabp7-expressing zebrafish glial cells. By using three different Fabp7 regulatory element-mediated fluorescent protein reporter strains, we reveal in double transgenic zebrafish that progenitor cells expressing fluorescent proteins driven by the Fabp7 regulatory element give rise to radial glia, oligodendrocyte progenitors, and some neuronal precursors. Furthermore, Bergmann glia represent the almost only glial population of the zebrafish cerebellum (besides a few oligodendrocytes), and the radial glia also remain in the mature cerebellum. Fabp7 regulatory element-mediated reporter protein expression in Bergmann glia progenitors suggests their origin from the ventral cerebellar proliferation zone, the ventricular zone, but not from the dorsally positioned upper rhombic lip. These new Fabp7 reporters will be valuable for functional studies during development and regeneration.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Molecular evolution of transcription factors AF4/FMR2 family member (AFF) gene family and the role of lamprey AFF3 in cell proliferation

AF4/FMR2 family member (AFF) proteins are a group of transcriptional regulators that can regulate gene transcription and play an important role in cellular physiological processes such as proliferation and differentiation. The transcriptome data of the lamprey spinal cord injury were analyzed in previous research. We then identified a hub gene, Lr-AFF3, from this dataset. Phylogenetic tree analysis determined the evolutionary relationships of the AFF gene family across different species. In addition, analysis of motifs, domains, and 3D structures further confirmed the conservatism of the AFF gene family. In particular, the gene structure of the AFF3 gene was not conserved, possibly because of intron insertion. It was also found that the neighboring genes of the Lr-AFF3 gene had a higher diversity than that in jawed vertebrates through synteny analysis. The results of the MTT and EdU experiments showed that the C-terminal homology domain (CHD) and N-terminal homology domain (NHD) of Lr-AFF3 promoted cell proliferation. In summary, our research will not only provide new insights into the origin and evolution of the AFF gene family in different species, but also provide new clues for the functions of Lr_AFF3.

Zebrafish Congenital Heart Disease Models: Opportunities and Challenges

Congenital heart defects (CHDs) are common human birth defects. Genetic mutations potentially cause the exhibition of various pathological phenotypes associated with CHDs, occurring alone or as part of certain syndromes. Zebrafish, a model organism with a strong molecular conservation similar to humans, is commonly used in studies on cardiovascular diseases owing to its advantageous features, such as a similarity to human electrophysiology, transparent embryos and larvae for observation, and suitability for forward and reverse genetics technology, to create various economical and easily controlled zebrafish CHD models. In this review, we outline the pros and cons of zebrafish CHD models created by genetic mutations associated with single defects and syndromes and the underlying pathogenic mechanism of CHDs discovered in these models. The challenges of zebrafish CHD models generated through gene editing are also discussed, since the cardiac phenotypes resulting from a single-candidate pathological gene mutation in zebrafish might not mirror the corresponding human phenotypes. The comprehensive review of these zebrafish CHD models will facilitate the understanding of the pathogenic mechanisms of CHDs and offer new opportunities for their treatments and intervention strategies.

Conserved expression of the zebrafish syt4 gene in GABAergic neurons in the cerebellum of adult fishes revealed by mammalian SYT4 immunoreactive-like signals

Synaptotagmin 4 (syt4) belongs to the synaptotagmin protein family, which has 17 and 28 family members in human and zebrafish, respectively. In zebrafish and rodents, syt4 is known to express abundantly in the entire central nervous system in the early developmental stages. In adult rodents, the gene expression shifts to be predominant in the cerebellum, mostly in Purkinje cells, a type of GABAergic neurons. However, there is no report of the expression pattern of syt4 in the adult zebrafish brain. Therefore, we hypothesize that the expression of syt4 is conserved in adult zebrafish and is specific to the GABAergic neurons, likely Purkinje cells, in the cerebellum. To examine the hypothesis, we first show that only one copy of syt4 gene remains in the zebrafish genome, and it is orthologous to the gene in other vertebrates. We further observe mammalian SYT4 antibody immunoreactive-like (mSYT4-ir) signals in several structures in the hindbrain including the medial divisions of the valvula cerebelli and the corpus cerebelli. In addition, our observations indicate the presence of mSYT4-ir signals in GABAergic neurons, most notably in the Purkinje cell layer of the molecular layer in the aforementioned structures. Conversely, mSYT4-ir signals are not observed in glutamatergic or cholinergic neurons. Therefore, we deduce that the syt4 gene in zebrafish exhibits a homologous expression pattern to those of previously studied vertebrate species, which is revealed by the positive immunoreactive-like signals of mammalian SYT4 antibodies.

Zebrafish ApoB-Containing Lipoprotein Metabolism: A Closer Look

Zebrafish have become a powerful model of mammalian lipoprotein metabolism and lipid cell biology. Most key proteins involved in lipid metabolism, including cholesteryl ester transfer protein, are conserved in zebrafish. Consequently, zebrafish exhibit a human-like lipoprotein profile. Zebrafish with mutations in genes linked to human metabolic diseases often mimic the human phenotype. Zebrafish larvae develop rapidly and externally around the maternally deposited yolk. Recent work revealed that any disturbance of lipoprotein formation leads to the accumulation of cytoplasmic lipid droplets and an opaque yolk, providing a visible phenotype to investigate disturbances of the lipoprotein pathway, already leading to discoveries in MTTP (microsomal triglyceride transfer protein) and ApoB (apolipoprotein B). By 5 days of development, the digestive system is functional, making it possible to study fluorescently labeled lipid uptake in the transparent larvae. These and other approaches enabled the first in vivo description of the STAB (stabilin) receptors, showing lipoprotein uptake in endothelial cells. Various zebrafish models have been developed to mimic human diseases by mutating genes known to influence lipoproteins (eg, ldlra, apoC2). This review aims to discuss the most recent research in the zebrafish ApoB-containing lipoprotein and lipid metabolism field. We also summarize new insights into lipid processing within the yolk cell and how changes in lipid flux alter yolk opacity. This curious new finding, coupled with the development of several techniques, can be deployed to identify new players in lipoprotein research directly relevant to human disease.

Effects of serta and sertb knockout on aggression in zebrafish (Danio rerio)

Zebrafish (Danio rerio) are unusual in having two paralogues of the serotonin re-uptake transporter (Sert), slc6a4a (serta) and slc6a4b (sertb), the transporter that serves in serotonin re-uptake from a synapse into the pre-synaptic cell or in serotonin uptake from the extracellular milieu into cells in the peripheral tissues. To address a knowledge gap concerning the specific roles of these paralogues, we used CRISPR/Cas9 technology to generate zebrafish knockout lines predicted to lack functional expression of Serta or Sertb. The consequences of loss-of-function of Serta or Sertb were assessed at the gene expression level, focusing on the serotonergic signalling pathway, and at the behaviour level, focusing on aggression. Whereas serta mRNA was expressed in all tissues examined, with high expression in the heart, gill and brain, only the brain displayed substantial sertb mRNA expression. In both serta-/- and sertb-/- fish, changes in transcript abundances of multiple components of the serotonin signalling pathway were detected, including proteins involved in serotonin synthesis (tph1a, tph1b, tph2, ddc), packaging (vmat2) and degradation (mao), and serotonin receptors (htr1aa, htr1ab). Using a mirror aggression test, serta-/- male but not female fish exhibited greater aggression than wildtype fish. However, both male and female sertb-/- fish displayed less aggression than their wildtype counterparts. These differences in behaviour between serta-/- and sertb-/- individuals hold promise for increasing our understanding of the neurophysiological basis of aggression in zebrafish.

Polyploidization of Indotyphlops braminus: evidence from isoform-sequencing

BACKGROUND

Indotyphlops braminus, the only known triploid parthenogenetic snake, is a compelling species for revealing the mechanism of polyploid emergence in vertebrates.

METHODS

In this study, we applied PacBio isoform sequencing technology to generate the first full-length transcriptome of I. braminus, aiming to improve the understanding of the molecular characteristics of this species.

RESULTS

A total of 51,849 nonredundant full-length transcript assemblies (with an N50 length of 2980 bp) from I. braminus were generated and fully annotated using various gene function databases. Our analysis provides preliminary evidence supporting a recent genome duplication event in I. braminus. Phylogenetic analysis indicated that the divergence of I. braminus subgenomes occurred approximately 11.5 ~ 15 million years ago (Mya). The full-length transcript resource generated as part of this research will facilitate transcriptome analysis and genomic evolution studies in the future.

Functional analysis of germline VANGL2 variants using rescue assays of vangl2 knockout zebrafish

Developmental studies have shown that the evolutionarily conserved Wnt Planar Cell Polarity (PCP) pathway is essential for the development of a diverse range of tissues and organs including the brain, spinal cord, heart and sensory organs, as well as establishment of the left-right body axis. Germline mutations in the highly conserved PCP gene VANGL2 in humans have only been associated with central nervous system malformations, and functional testing to understand variant impact has not been performed. Here we report three new families with missense variants in VANGL2 associated with heterotaxy and congenital heart disease p.(Arg169His), non-syndromic hearing loss p.(Glu465Ala) and congenital heart disease with brain defects p.(Arg135Trp). To test the in vivo impact of these and previously described variants, we have established clinically-relevant assays using mRNA rescue of the vangl2 mutant zebrafish. We show that all variants disrupt Vangl2 function, although to different extents and depending on the developmental process. We also begin to identify that different VANGL2 missense variants may be haploinsufficient and discuss evidence in support of pathogenicity. Together, this study demonstrates that zebrafish present a suitable pipeline to investigate variants of unknown significance and suggests new avenues for investigation of the different developmental contexts of VANGL2 function that are clinically meaningful.

Zebrafish cardiac repolarization does not functionally depend on the expression of the hERG1b-like transcript

Zebrafish provide a translational model of human cardiac function. Their similar cardiac electrophysiology enables screening of human cardiac repolarization disorders, drug arrhythmogenicity, and novel antiarrhythmic therapeutics. However, while zebrafish cardiac repolarization is driven by delayed rectifier potassium channel current (IKr), the relative role of alternate channel transcripts is uncertain. While human ether-a-go-go-related-gene-1a (hERG1a) is the dominant transcript in humans, expression of the functionally distinct alternate transcript, hERG1b, modifies the electrophysiological and pharmacologic IKr phenotype. Studies of zebrafish IKr are frequently translated without consideration for the presence and impact of hERG1b in humans. Here, we performed phylogenetic analyses of all available KCNH genes from Actinopterygii (ray-finned fishes). Our findings confirmed zebrafish cardiac zkcnh6a as the paralog of human hERG1a (hKCNH2a), but also revealed evidence of a hERG1b (hKCNH2b)-like N-terminally truncated gene, zkcnh6b, in zebrafish. zkcnh6b is a teleost-specific variant that resulted from the 3R genome duplication. qRT-PCR showed dominant expression of zkcnh6a in zebrafish atrial and ventricular tissue, with low levels of zkcnh6b. Functional evaluation of zkcnh6b in a heterologous system showed no discernable function under the conditions tested, and no influence on zkcnh6a function during the zebrafish ventricular action potential. Our findings provide the first descriptions of the zkcnh6b gene, and show that, unlike in humans, zebrafish cardiac repolarization does not rely upon co-assembly of zERG1a/zERG1b. Given that hERG1b modifies IKr function and drug binding in humans, our findings highlight the need for consideration when translating hERG variant effects and toxicological screens in zebrafish, which lack a functional hERG1b-equivalent gene.

Protein-protein interaction network module changes associated with the vertebrate fin-to-limb transition

Evolutionary phenotypic transitions, such as the fin-to-limb transition in vertebrates, result from modifications in related proteins and their interactions, often in response to changing environment. Identifying these alterations in protein networks is crucial for a more comprehensive understanding of these transitions. However, previous research has not attempted to compare protein-protein interaction (PPI) networks associated with evolutionary transitions, and most experimental studies concentrate on a limited set of proteins. Therefore, the goal of this work was to develop a network-based platform for investigating the fin-to-limb transition using PPI networks. Quality-enhanced protein networks, constructed by integrating PPI networks with anatomy ontology data, were leveraged to compare protein modules for paired fins (pectoral fin and pelvic fin) of fishes (zebrafish) to those of the paired limbs (forelimb and hindlimb) of mammals (mouse). This also included prediction of novel protein candidates and their validation by enrichment and homology analyses. Hub proteins such as shh and bmp4, which are crucial for module stability, were identified, and their changing roles throughout the transition were examined. Proteins with preserved roles during the fin-to-limb transition were more likely to be hub proteins. This study also addressed hypotheses regarding the role of non-preserved proteins associated with the transition.

Evolution of lysine-specific demethylase 1 and REST corepressor gene families and their molecular interaction

Lysine-specific demethylase 1A (LSD1) binds to the REST corepressor (RCOR) protein family of corepressors to erase transcriptionally active marks on histones. Functional diversity in these complexes depends on the type of RCOR included, which modulates the catalytic activity of the complex. Here, we studied the duplicative history of the RCOR and LSD gene families and analyzed the evolution of their interaction. We found that RCOR genes are the product of the two rounds of whole-genome duplications that occurred early in vertebrate evolution. In contrast, the origin of the LSD genes traces back before to the divergence of animals and plants. Using bioinformatics tools, we show that the RCOR and LSD1 interaction precedes the RCOR repertoire expansion that occurred in the last common ancestor of jawed vertebrates. Overall, we trace LSD1-RCOR complex evolution and propose that animal non-model species offer advantages in addressing questions about the molecular biology of this epigenetic complex.

A chromosome-level reference genome for the common octopus, Octopus vulgaris (Cuvier, 1797)

Cephalopods are emerging animal models and include iconic species for studying the link between genomic innovations and physiological and behavioral complexities. Coleoid cephalopods possess the largest nervous system among invertebrates, both for cell counts and brain-to-body ratio. Octopus vulgaris has been at the center of a long-standing tradition of research into diverse aspects of cephalopod biology, including behavioral and neural plasticity, learning and memory recall, regeneration, and sophisticated cognition. However, no chromosome-scale genome assembly was available for O. vulgaris to aid in functional studies. To fill this gap, we sequenced and assembled a chromosome-scale genome of the common octopus, O. vulgaris. The final assembly spans 2.8 billion basepairs, 99.34% of which are in 30 chromosome-scale scaffolds. Hi-C heatmaps support a karyotype of 1n = 30 chromosomes. Comparisons with other octopus species' genomes show a conserved octopus karyotype and a pattern of local genome rearrangements between species. This new chromosome-scale genome of O. vulgaris will further facilitate research in all aspects of cephalopod biology, including various forms of plasticity and the neural machinery underlying sophisticated cognition, as well as an understanding of cephalopod evolution.

Role of HIF in fish inflammation

The hypoxia-inducing factor (HIF) is a central transcription factor in cellular oxygen sensing and regulation. It is common that the inflammation always appears in many diseases, like infectious diseases in fishes, and the inflammation is often accompanied by hypoxia, as a hallmark of inflammation. Besides coordinating cellular responses to low oxygen, HIF-mediated hypoxia signaling pathway is also crucial for immune responses such as the regulations of innate immune cell phenotype and function, as well as metabolic reprogramming under the inflammation. However, the understanding of the molecular mechanisms by which HIFs regulate the inflammatory response in fish is still very limited. Here, we review the characteristics of HIF as well as its roles in innate immune cells and the infections caused by bacteria and viruses. The regulatory effects of HIF on the metabolic reprogramming of innate immune cells are also discussed and the future research directions are outlooked. This paper will serve as a reference for elucidating the molecular mechanism of HIF regulating inflammation and identifying treatment strategies to target HIF for fish disease.

Unraveling the Unusual Subgenomic Organization in the Neopolyploid Free-Living Flatworm Macrostomum lignano

Whole genome duplication (WGD) is an evolutionary event resulting in a redundancy of genetic material. Different mechanisms of WGD, allo- or autopolyploidization, lead to distinct evolutionary trajectories of newly formed polyploids. Genome studies on such species are important for understanding the early stages of genome evolution. However, assembling neopolyploid is a challenging task due to the presence of 2 homologous (or homeologous) chromosome sets and therefore the existence of the extended paralogous regions in its genome. Post-WGD evolution of polyploids includes cytogenetic diploidization leading to the formation of species, whose polyploid origin might be hidden by disomic inheritance. Earlier we uncovered the hidden polyploid origin of the free-living flatworms of the genus Macrostomum (Macrostomum lignano, M. janickei, and M. mirumnovem). Cytogenetic diploidization in these species is accompanied by intensive chromosomal rearrangements including chromosomes fusions. In this study, we unravel the M. lignano genome organization through generation and sequencing of 2 sublines of the commonly used inbred line of M. lignano (called DV1) differing only in a copy number of the largest chromosome (MLI1). Using nontrivial assembly free comparative analysis of their genomes, we deciphered DNA sequences belonging to MLI1 and validated them by sequencing the pool of microdissected MLI1. Here we presented the uncommon mechanism of genome rediplodization of M. lignano, which consists of (i) presence of 3 subgenomes, which emerged via formation of large fused chromosomes and its variants, and (ii) sustaining their heterozygosity through inter- and intrachromosomal rearrangements.

Genomic organization and transcription of superoxide dismutase genes (sod1, sod2, and sod3b) and response to diazinon toxicity in platyfish (Xiphophorus maculatus) by using SOD enzyme activity

The aim of this study is to determine the effects of 50% of 96 h LC50 (5.25 ppm) diazinon on the expression of superoxide dismutase (SOD) enzyme genes (sod1, sod2, and sod3b) and SOD enzyme activity at the end of 24, 48, 72, and 96 h in platyfish liver and gill tissues. To this end, we determined the tissue-specific distribution of sod1, sod2, and sod3b genes and performed in silico analyses in platyfish (Xiphophorus maculatus). It was determined that malondialdehyde (MDA) level and SOD enzyme activity were increased in the liver [(43.90 EU mg protein-1 (control), 62.45 EU mg protein-1 (24 h), 73.17 EU mg protein-1 (48 h), 82.18 EU mg protein-1 (72 h), 92.93 EU mg protein-1 (96 h)] and gill [(16.44 EU mg protein-1 (control), 33.47 EU mg protein-1 (24 h), 50.38 EU mg protein-1 (48 h), 64.62 EU mg protein-1 (72 h), 74.04 EU mg protein-1 (96 h)] tissues of platyfish exposed to diazinon, while the expression of the sod genes was down-regulated. The tissue-specific distribution of the sod genes varied, with the tissues and the sod genes expression were being predominant in the liver (628.32 in sod1, 637.59 in sod2, 888.5 in sod3b). Thus, the liver was considered a suitable tissue for further gene expression studies. Based on the phylogenetic analyses, platyfish sod genes can be reported to be orthologs of sod/SOD genes from other vertebrates. Identity/similarity analyses supported this determination. Conserved gene synteny proved that there are conserved sod genes in platyfish, zebrafish, and humans.

Evolution of the Spider Homeobox Gene Repertoire by Tandem and Whole Genome Duplication

Gene duplication generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. Duplicated genes can undergo subfunctionalization to partition ancestral functions and/or neofunctionalization to assume a new function. We previously found there had been a whole genome duplication (WGD) in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues originated by WGD versus smaller-scale events such as tandem duplications. Understanding this is a key to determining the contribution of the WGD to arachnopulmonate genome evolution. Here we characterized the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that most duplicated homeobox genes in spiders are consistent with an origin by WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE clusters, in all eight species. Consistently, we observed one copy of each cluster was degenerated in terms of gene content and organization while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalization between the duplicated NK genes in the spider Parasteatoda tepidariorum compared to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the relative contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.

Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes

The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.

Subfunctionalisation of paralogous genes and evolution of differential codon usage preferences: The showcase of polypyrimidine tract binding proteins

Gene paralogs are copies of an ancestral gene that appear after gene or full genome duplication. When two sister gene copies are maintained in the genome, redundancy may release certain evolutionary pressures, allowing one of them to access novel functions. Here, we focused our study on gene paralogs on the evolutionary history of the three polypyrimidine tract binding protein genes (PTBP) and their concurrent evolution of differential codon usage preferences (CUPrefs) in vertebrate species. PTBP1-3 show high identity at the amino acid level (up to 80%) but display strongly different nucleotide composition, divergent CUPrefs and, in humans and in many other vertebrates, distinct tissue-specific expression levels. Our phylogenetic inference results show that the duplication events leading to the three extant PTBP1-3 lineages predate the basal diversification within vertebrates, and genomic context analysis illustrates that local synteny has been well preserved over time for the three paralogs. We identify a distinct evolutionary pattern towards GC3-enriching substitutions in PTBP1, concurrent with enrichment in frequently used codons and with a tissue-wide expression. In contrast, PTBP2s are enriched in AT-ending, rare codons, and display tissue-restricted expression. As a result of this substitution trend, CUPrefs sharply differ between mammalian PTBP1s and the rest of PTBPs. Genomic context analysis suggests that GC3-rich nucleotide composition in PTBP1s is driven by local substitution processes, while the evidence in this direction is thinner for PTBP2-3. An actual lack of co-variation between the observed GC composition of PTBP2-3 and that of the surrounding non-coding genomic environment would raise an interrogation on the origin of CUPrefs, warranting further research on a putative tissue-specific translational selection. Finally, we communicate an intriguing trend for the use of the UUG-Leu codon, which matches the trends of AT-ending codons. Our results are compatible with a scenario in which a combination of directional mutation-selection processes would have differentially shaped CUPrefs of PTBPs in vertebrates: the observed GC-enrichment of PTBP1 in placental mammals may be linked to genomic location and to the strong and broad tissue-expression, while AT-enrichment of PTBP2 and PTBP3 would be associated with rare CUPrefs and thus, possibly to specialized spatio-temporal expression. Our interpretation is coherent with a gene subfunctionalisation process by differential expression regulation associated with the evolution of specific CUPrefs.

Evolution and divergence of teleost adrenergic receptors: why sometimes 'the drugs don't work' in fish

Adrenaline and noradrenaline, released as hormones and/or neurotransmitters, exert diverse physiological functions in vertebrates, and teleost fishes are widely used as model organisms to study adrenergic regulation; however, such investigations often rely on receptor subtype-specific pharmacological agents (agonists and antagonists; see Glossary) developed and validated in mammals. Meanwhile, evolutionary (phylogenetic and comparative genomic) studies have begun to unravel the diversification of adrenergic receptors (ARs) and reveal that whole-genome duplications and pseudogenization events in fishes results in notable distinctions from mammals in their genomic repertoire of ARs, while lineage-specific gene losses within teleosts have generated significant interspecific variability. In this Review, we visit the evolutionary history of ARs (including α1-, α2- and β-ARs) to highlight the prominent interspecific differences in teleosts, as well as between teleosts and other vertebrates. We also show that structural modelling of teleost ARs predicts differences in ligand binding affinity compared with mammalian orthologs. To emphasize the difficulty of studying the roles of different AR subtypes in fish, we collate examples from the literature of fish ARs behaving atypically compared with standard mammalian pharmacology. Thereafter, we focus on specific case studies of the liver, heart and red blood cells, where our understanding of AR expression has benefited from combining pharmacological approaches with molecular genetics. Finally, we briefly discuss the ongoing advances in 'omics' technologies that, alongside classical pharmacology, will provide abundant opportunities to further explore adrenergic signalling in teleosts.

Genome-Wide Identification of Trachinotus ovatus Antimicrobial Peptides and Their Immune Response against Two Pathogen Challenges

Golden pompano, Trachinotus ovatus, as a highly nutritious commercially valuable marine fish, has become one of the preferred species for many fish farmers due to its rapid growth, wide adaptability, and ease of feeding and management. However, with the expansion of aquaculture scale, bacterial and parasitic diseases have also become major threats to the golden pompano industry. This study, based on comparative genomics, shows the possibility of preferential evolution of freshwater fish over marine fish by analyzing the phylogenetic relationships and divergence times of 14 marine fish and freshwater fish. Furthermore, we identified antimicrobial peptide genes from 14 species at the genomic level and found that the number of putative antimicrobial peptides may be related to species evolution. Subsequently, we classified the 341 identified AMPs from golden pompano into 38 categories based on the classification provided by the APD3. Among them, TCP represented the highest proportion, accounting for 23.2% of the total, followed by scolopendin, lectin, chemokine, BPTI, and histone-derived peptides. At the same time, the distribution of AMPs in chromosomes varied with type, and covariance analysis showed the frequency of its repeat events. Enrichment analysis and PPI indicated that AMP was mainly concentrated in pathways associated with disease immunity. In addition, our transcriptomic data measured the expression of putative AMPs of golden pompano in 12 normal tissues, as well as in the liver, spleen, and kidney infected with Streptococcus agalactiae and skin infected with Cryptocaryon irritans. As the infection with S. agalactiae and C. irritans progressed, we observed tissue specificity in the number and types of responsive AMPs. Positive selection of AMP genes may participate in the immune response through the MAPK signaling pathway. The genome-wide identification of antimicrobial peptides in the golden pompano provided a complete database of potential AMPs that can contribute to further understanding the immune mechanisms in pathogens. AMPs were expected to replace traditional antibiotics and be developed into targeted drugs against specific bacterial and parasitic pathogens for more precise and effective treatment to improve aquaculture production.

Parallel evolution of fish bi-modal breathing and expansion of olfactory receptor (OR) genes: toward a universal ORs nomenclature

Olfactory receptors (ORs) play a key role in the prime sensorial perception, being highly relevant for intra/interspecific interactions. ORs are a subgroup of G-protein coupled receptors that exhibit highly complex subgenomes in vertebrates. However, OR repertoires remain poorly studied in fish lineages, precluding finely retracing their origin, evolution, and diversification, especially in the most basal groups. Here, we conduct an exhaustive gene screening upon 43 high-quality fish genomes exhibiting varied gene repertoires (2-583 genes). While the early vertebrates performed gas exchange through gills, we hypothesize that the emergence of new breathing structures (swim bladder and paired lungs) in early osteichthyans may be associated with expansions in the ORs gene families sensitive to airborne molecules. Additionally, we verify that the OR repertoire of moderns actinopterygians has not increased as expected following a whole genome duplication, likely due to regulatory mechanisms compensating the gene load excess. Finally, we identify 25 distinct OR families, allowing us to propose an updated universal nomenclature for the fish ORs.

The Zebrafish Embryo as a Model Organism for Testing mRNA-Based Therapeutics

mRNA-based therapeutics have revolutionized the world of molecular therapy and have proven their potential in the vaccination campaigns for SARS-CoV2 and clinical trials for hereditary disorders. Preclinical studies have mainly focused on in vitro and rodent studies. However, research in rodents is costly and labour intensive, and requires ethical approval for all interventions. Zebrafish embryonic disease models are not always classified as laboratory animals and have been shown to be extremely valuable for high-throughput drug testing. Zebrafish larvae are characterized by their small size, optical transparency and high number of embryos, and are therefore also suited for the study of mRNA-based therapeutics. First, the one-cell stage injection of naked mRNA can be used to assess the effectivity of gene addition in vivo. Second, the intravascular injection in older larvae can be used to assess tissue targeting efficiency of (packaged) mRNA. In this review, we describe how zebrafish can be used as a steppingstone prior to testing mRNA in rodent models. We define the procedures that can be employed for both the one-cell stage and later-stage injections, as well as the appropriate procedures for post-injection follow-up.

Nocturnal melatonin increases glucose uptake via insulin-independent action in the goldfish brain

Melatonin, a neurohormone nocturnally produced by the pineal gland, is known to regulate the circadian rhythm. It has been recently reported that variants of melatonin receptors are associated with an increased risk of hyperglycemia and type 2 diabetes, suggesting that melatonin may be involved in the regulation of glucose homeostasis. Insulin is a key hormone that regulates circulating glucose levels and cellular metabolism after food intake in many tissues, including the brain. Although cells actively uptake glucose even during sleep and without food, little is known regarding the physiological effects of nocturnal melatonin on glucose homeostasis. Therefore, we presume the involvement of melatonin in the diurnal rhythm of glucose metabolism, independent of insulin action after food intake. In the present study, goldfish (Carassius auratus) was used as an animal model, since this species has no insulin-dependent glucose transporter type 4 (GLUT4). We found that in fasted individuals, plasma melatonin levels were significantly higher and insulin levels were significantly lower during the night. Furthermore, glucose uptake in the brain, liver, and muscle tissues also significantly increased at night. After intraperitoneal administration of melatonin, glucose uptake by the brain and liver showed significantly greater increases than in the control group. The administration of melatonin also significantly decreased plasma glucose levels in hyperglycemic goldfish, but failed to alter insulin mRNA expression in Brockmann body and plasma insulin levels. Using an insulin-free medium, we demonstrated that melatonin treatment increased glucose uptake in a dose-dependent manner in primary cell cultures of goldfish brain and liver cells. Moreover, the addition of a melatonin receptor antagonist decreased glucose uptake in hepatocytes, but not in brain cells. Next, treatment with N1-acetyl-5-methoxykynuramine (AMK), a melatonin metabolite in the brain, directly increased glucose uptake in cultured brain cells. Taken together, these findings suggest that melatonin is a possible circadian regulator of glucose homeostasis, whereas insulin acquires its effect on glucose metabolism following food intake.

Fish Models for Exploring Mitochondrial Dysfunction Affecting Neurodegenerative Disorders

Neurodegenerative disorders are characterized by the progressive loss of neuronal structure or function, resulting in memory loss and movement disorders. Although the detailed pathogenic mechanism has not been elucidated, it is thought to be related to the loss of mitochondrial function in the process of aging. Animal models that mimic the pathology of a disease are essential for understanding human diseases. In recent years, small fish have become ideal vertebrate models for human disease due to their high genetic and histological homology to humans, ease of in vivo imaging, and ease of genetic manipulation. In this review, we first outline the impact of mitochondrial dysfunction on the progression of neurodegenerative diseases. Then, we highlight the advantages of small fish as model organisms, and present examples of previous studies regarding mitochondria-related neuronal disorders. Lastly, we discuss the applicability of the turquoise killifish, a unique model for aging research, as a model for neurodegenerative diseases. Small fish models are expected to advance our understanding of the mitochondrial function in vivo, the pathogenesis of neurodegenerative diseases, and be important tools for developing therapies to treat diseases.

Genes enriched in A/T-ending codons are co-regulated and conserved across mammals

Codon usage influences gene expression distinctly depending on the cell context. Yet, the importance of codon bias in the simultaneous turnover of specific groups of protein-coding genes remains to be investigated. Here, we find that genes enriched in A/T-ending codons are expressed more coordinately in general and across tissues and development than those enriched in G/C-ending codons. tRNA abundance measurements indicate that this coordination is linked to the expression changes of tRNA isoacceptors reading A/T-ending codons. Genes with similar codon composition are more likely to be part of the same protein complex, especially for genes with A/T-ending codons. The codon preferences of genes with A/T-ending codons are conserved among mammals and other vertebrates. We suggest that this orchestration contributes to tissue-specific and ontogenetic-specific expression, which can facilitate, for instance, timely protein complex formation.

Genetically engineered zebrafish as models of skeletal development and regeneration

Zebrafish (Danio rerio) are aquatic vertebrates with significant homology to their terrestrial counterparts. While zebrafish have a centuries-long track record in developmental and regenerative biology, their utility has grown exponentially with the onset of modern genetics. This is exemplified in studies focused on skeletal development and repair. Herein, the numerous contributions of zebrafish to our understanding of the basic science of cartilage, bone, tendon/ligament, and other skeletal tissues are described, with a particular focus on applications to development and regeneration. We summarize the genetic strengths that have made the zebrafish a powerful model to understand skeletal biology. We also highlight the large body of existing tools and techniques available to understand skeletal development and repair in the zebrafish and introduce emerging methods that will aid in novel discoveries in skeletal biology. Finally, we review the unique contributions of zebrafish to our understanding of regeneration and highlight diverse routes of repair in different contexts of injury. We conclude that zebrafish will continue to fill a niche of increasing breadth and depth in the study of basic cellular mechanisms of skeletal biology.

CRISPR/Cas9 and FLP-FRT mediated regulatory dissection of the BX-C of Drosophila melanogaster

The homeotic genes or Hox define the anterior-posterior (AP) body axis formation in bilaterians and are often present on the chromosome in an order collinear to their function across the AP axis. However, there are many cases wherein the Hox are not collinear, but their expression pattern is conserved across the AP axis. The expression pattern of Hox is attributed to the cis-regulatory modules (CRMs) consisting of enhancers, initiators, or repressor elements that regulate the genes in a segment-specific manner. In the Drosophila melanogaster Hox complex, the bithorax complex (BX-C) and even the CRMs are organized in an order that is collinear to their function in the thoracic and abdominal segments. In the present study, the regulatorily inert regions were targeted using CRISPR/Cas9 to generate a series of transgenic lines with the insertion of FRT sequences. These FRT lines are repurposed to shuffle the CRMs associated with Abd-B to generate modular deletion, duplication, or inversion of multiple CRMs. The rearrangements yielded entirely novel phenotypes in the fly suggesting the requirement of such complex manipulations to address the significance of higher order arrangement of the CRMs. The functional map and the transgenic flies generated in this study are important resources to decipher the collective ability of multiple regulatory elements in the eukaryotic genome to function as complex modules.

Genome-wide identification and characterization of toll-like receptor 5 (TLR5) in fishes

Toll-like receptors 5 (TLR5), a member of the toll-like receptors (TLRs) family, is a class of pattern recognition receptors (PRRs) that recognize pathogen-associated molecular patterns (PAMPs). It responds to vertebrate recognition of bacterial flagellin and participates in innate immune responses. However, genome-wide identification and characterization of TLR5 in fishes have not been investigated. Here, three TLR5M isotypes (TLR5Ma, TLR5Mb1, and TLR5Mb2) and a TLR5S are all extracted from fish genomes on the basis of phylogenetic and synteny analyses. We confirmed that the non-teleost fishes have one TLR5M gene, as well as additional TLR5 genes (TLR5M and TLR5S) in teleost fishes. In addition, some special teleost fishes possess two to three TLR5 genes, which have undergone the fourth whole-genome duplication (WGD). According to our results, we inferred that the diversity of TLR5 genes in fishes seems to be the result of combinations of WGD and gene loss. Furthermore, TLR5 isoforms displayed differences at the flagellin interaction sites and viral binding sites, and showed lineage-specific, which indicated that TLR5 duplicates may generate functional divergence. Bacterial experiments also supported the idea that CiTLR5Ma and CiTLR5Mb are subfunctionalized to sense bacterial flagellin. In summary, our present comparative genomic survey will benefit for further functional investigations of TLR5 genes in fish.

A Revised Perspective on the Evolution of Troponin I and Troponin T Gene Families in Vertebrates

The troponin (Tn) complex, responsible for the Ca2+ activation of striated muscle, is composed of three interacting protein subunits: TnC, TnI, and TnT, encoded by TNNC, TNNI, and TNNT genes. TNNI and TNNT are sister gene families, and in mammals the three TNNI paralogs (TNNI1, TNNI2, TNNI3), which encode proteins with tissue-specific expression, are each in close genomic proximity with one of the three TNNT paralogs (TNNT2, TNNT3, TNNT1, respectively). It has been widely presumed that all vertebrates broadly possess genes of these same three classes, although earlier work has overlooked jawless fishes (cyclostomes) and cartilaginous fishes (chimeras, rays, and sharks), which are distantly related to other jawed vertebrates. With a new phylogenetic and synteny analysis of a diverse array of vertebrates including these taxonomic groups, we define five distinct TNNI classes (TNNI1-5), with TNNI4 and TNNI5 being only present in non-amniote vertebrates and typically found in tandem, and four classes of TNNT (TNNT1-4). These genes are located in four genomic loci that were generated by the 2R whole-genome duplications. TNNI3, encoding "cardiac TnI" in tetrapods, was independently lost in cartilaginous and ray-finned fishes. Instead, ray-finned fishes predominantly express TNNI1 in the heart. TNNI5 is highly expressed in shark hearts and contains a N-terminal extension similar to that of TNNI3 found in tetrapod hearts. Given that TNNI3 and TNNI5 are distantly related, this supports the hypothesis that the N-terminal extension may be an ancestral feature of vertebrate TNNI and not an innovation unique to TNNI3, as has been commonly believed.

Early calcium and cardiac contraction defects in a model of phospholamban R9C mutation in zebrafish

The phospholamban mutation Arg 9 to Cys (R9C) has been found to cause a dilated cardiomyopathy in humans and in transgenic mice, with ventricular dilation and premature death. Emerging evidence suggests that phospholamban R9C is a loss-of-function mutation with dominant negative effect on SERCA2a activity. We imaged calcium and cardiac contraction simultaneously in 3 and 9 days-post-fertilization (dpf) zebrafish larvae expressing plnb^R9C in the heart to unveil the early pathological pathway that triggers the disease. We generated transgenic zebrafish lines expressing phospholamban wild-type (Tg(myl7:plnb^wt)) and phospholamban R9C (Tg(myl7:plnb^R9C)) in the heart of zebrafish. To measure calcium and cardiac contraction in 3 and 9 dpf larvae, Tg(myl7:plnb^wt) and Tg(myl7:plnb^R9C) fish were outcrossed with a transgenic line expressing the ratiometric fluorescent calcium biosensor mCyRFP1-GCaMP6f. We found that Plnb^R9C raised calcium transient amplitude, induced positive inotropy and lusitropy, and blunted the β-adrenergic response to isoproterenol in 3 dpf larvae. These effects can be attributed to enhanced SERCA2a activity induced by the Plnb^R9C mutation. In contrast, Tg(myl7:plnb^R9C) larvae at 9 dpf exhibited ventricular dilation, systolic dysfunction and negative lusitropy, hallmarks of a dilated cardiomyopathy in humans. Importantly, N-acetyl-L-cysteine rescued this deleterious phenotype, suggesting that reactive oxygen species contribute to the pathological pathway. These results also imply that dysregulation of calcium homeostasis during embryo development contributes to the disease progression at later stages. Our in vivo model in zebrafish allows characterization of pathophysiological mechanisms leading to heart disease, and can be used for screening of potential therapeutical agents.

Functional identification of the zebrafish Interleukin-1 receptor in an embryonic model of Il-1β-induced systemic inflammation

Interleukin-1β (IL-1β) is a potent proinflammatory cytokine that plays a vital role in the innate immune system. To observe the innate immune response in vivo, several transgenic zebrafish lines have been developed to model IL-1β-induced inflammation and to visualize immune cell migration and proliferation in real time. However, our understanding of the IL-1β response in zebrafish is limited due to an incomplete genome annotation and a lack of functional data for the cytokine receptors involved in the inflammatory process. Here, we use a combination of database mining, genetic analyses, and functional assays to identify zebrafish Interleukin-1 receptor, type 1 (Il1r1). We identified putative zebrafish il1r1 candidate genes that encode proteins with predicted structures similar to human IL1R1. To examine functionality of these candidates, we designed highly effective morpholinos to disrupt gene expression in a zebrafish model of embryonic Il-1β-induced systemic inflammation. In this double transgenic model, ubb:Gal4-EcR, uas:il1βmat , the zebrafish ubiquitin b (ubb) promoter drives expression of the modified Gal4 transcription factor fused to the ecdysone receptor (EcR), which in turn drives the tightly-regulated expression and secretion of mature Il-1β only in the presence of the ecdysone analog tebufenozide (Teb). Application of Teb to ubb:Gal4-EcR, uas:il1βmat embryos causes premature death, fin degradation, substantial neutrophil expansion, and generation of reactive oxygen species (ROS). To rescue these deleterious phenotypes, we injected ubb:Gal4-EcR, uas:il1βmat embryos with putative il1r1 morpholinos and found that knockdown of only one candidate gene prevented the adverse effects caused by Il-1β. Mosaic knockout of il1r1 using the CRISPR/Cas9 system phenocopied these results. Taken together, our study identifies the functional zebrafish Il1r1 utilizing a genetic model of Il-1β-induced inflammation and provides valuable new insights to study inflammatory conditions specifically driven by Il-1β or related to Il1r1 function in zebrafish.

Cd248a and Cd248b in zebrafish participate in innate immune responses

CD248, also known as endosialin or tumor endothelial marker 1, is a type I single transmembrane glycoprotein. CD248 has been demonstrated to be upregulated in cancers, tumors and many fibrotic diseases in human and mice, such as liver damage, pulmonary fibrosis, renal fibrosis, arthritis and tumor neovascularization. However, no definite CD248 orthologs in fish have been documented so far. In this study, we report the identification of cd248a and cd248b in the zebrafish. Both the phylogenetic analysis and the conserved synteny strongly suggested that zebrafish cd248a and cd248b are orthologs of the human CD248. Both cd248a and cd248b exhibited similar and dynamic expression pattern in early development, both genes had weak maternal expression, the zygotic transcripts were first seen in anterior somites and head mesenchyme, then shifted to eyes and head mesenchyme, later expanded to branchial arches, and gradually declined with development. The expression profiles of cd248a and cd248b were upregulated upon LPS (Lipopolysaccharide) challenge. Both Cd248a protein and Cd248b protein were localized on the cell membrane and cytoplasm, and overexpression of cd248a and cd248b induced the expression of pro-inflammatory cytokines, in vitro and in vivo. Moreover, deficiency of cd248a or cd248b both downregulated the expression of pro-inflammatory cytokines and upregulated anti-inflammatory cytokine. Additionally, loss of cd248a or cd248b both downregulated the expression of pro-inflammatory cytokines after LPS treatment. Taken together, these results indicated that cd248a and cd248b in zebrafish were involved in immune response and would provide further information to understand functions of Cd248 protein in innate immunity of fish.

Fishing for answers to hemostatic and thrombotic disease: Genome editing in zebrafish

Over the past two decades, the teleost vertebrate Danio rerio (zebrafish) has emerged as a model for hemostasis and thrombosis. At genomic and functional levels, there is a high degree of conservation of the hemostatic system with that of mammals. Numerous features of the fish model offer unique advantages for investigating hemostasis and thrombosis. These include high fecundity, rapid and external development, optical transparency, and extensive functional homology with mammalian hemostasis and thrombosis. Zebrafish are particularly suited to genome-wide mutagenesis experiments for the study of modifier genes. They are also amenable to whole-organism small-molecule screens, a feature that is exceptionally relevant to hemostasis and thrombosis. Zebrafish coagulation factor knockouts that are in utero or neonatal lethal in mammals survive into adulthood before succumbing to hemorrhage or thrombosis, enabling studies not possible in mammals. In this illustrated review, we outline how zebrafish have been employed for the study of hemostasis and thrombosis using modern genome editing techniques, coagulation assays in larvae, and in vivo evaluation of patient-specific variants to infer causality and demonstrate pathogenicity. Zebrafish hemostasis and thrombosis models will continue to serve as a clinically directed basic research tool and powerful alternative to mammals for the development of new diagnostic markers and novel therapeutics for coagulation disorders through high-throughput genetic and small-molecule studies.

Dichloroacetate improves mitochondrial function, physiology, and morphology in FBXL4 disease models

Pathogenic variants in the human F-box and leucine-rich repeat protein 4 (FBXL4) gene result in an autosomal recessive, multisystemic, mitochondrial disorder involving variable mitochondrial depletion and respiratory chain complex deficiencies with lactic acidemia. As no FDA-approved effective therapies for this disease exist, we sought to characterize translational C. elegans and zebrafish animal models, as well as human fibroblasts, to study FBXL4^–/– disease mechanisms and identify preclinical therapeutic leads. Developmental delay, impaired fecundity and neurologic and/or muscular activity, mitochondrial dysfunction, and altered lactate metabolism were identified in fbxl-1(ok3741) C. elegans. Detailed studies of a PDHc activator, dichloroacetate (DCA), in fbxl-1(ok3741)C. elegans demonstrated its beneficial effects on fecundity, neuromotor activity, and mitochondrial function. Validation studies were performed in fbxl4^sa12470 zebrafish larvae and in FBXL4^–/– human fibroblasts; they showed DCA efficacy in preventing brain death, impairment of neurologic and/or muscular function, mitochondrial biochemical dysfunction, and stress-induced morphologic and ultrastructural mitochondrial defects. These data demonstrate that fbxl-1(ok3741) C. elegans and fbxl4^sa12470 zebrafish provide robust translational models to study mechanisms and identify preclinical therapeutic candidates for FBXL4^–/– disease. Furthermore, DCA is a lead therapeutic candidate with therapeutic benefit on diverse aspects of survival, neurologic and/or muscular function, and mitochondrial physiology that warrants rigorous clinical trial study in humans with FBXL4^–/– disease.

Gene expression and functional analysis of Aha1a and Aha1b in stress response in zebrafish

Activator of heat shock protein 90 (hsp90) ATPase (Aha1) is a Hsp90 co-chaperone required for Hsp90 ATPase activation. Aha1 is essential for yeast survival and muscle development in C. elegans under elevated temperature and hsp90-deficeiency induced stress conditions. The roles of Aha1 in vertebrates are poorly understood. Here, we characterized the expression and function of Aha1 in zebrafish. We showed that zebrafish genome contains two aha1 genes, aha1a and aha1b, that show distinct patterns of expression during development. Under the normal physiological conditions, aha1a is primarily expressed in skeletal muscle cells of zebrafish embryos, while aha1b is strongly expressed in the head region. aha1a and aha1b expression increased dramatically in response to heat shock induced stress. In addition, Aha1a-GFP fusion protein exhibited a dynamic translocation in muscle cells in response to heat shock. Moreover, upregulation of aha1 expression was also observed in hsp90a1 knockdown embryos that showed a muscle defect. Genetic studies demonstrated that knockout of aha1a, aha1b or both had no detectable effect on embryonic development, survival, and growth in zebrafish. The aha1a and aha1b mutant embryos showed normal muscle development and stress response in response to heat shock. Single or double aha1a and aha1b mutants could grow into normal reproductive adults with normal skeletal muscle structure and morphology compared with wild type control. Together, data from these studies indicate that Aha1a and Aha1b are involved in stress response. However, they are dispensable in zebrafish embryonic development, growth, and survival.

Loss-of-Function Models of the Metabotropic Glutamate Receptor Genes Grm8a and Grm8b Display Distinct Behavioral Phenotypes in Zebrafish Larvae (Danio rerio)

Members of the family of metabotropic glutamate receptors are involved in the pathomechanism of several disorders of the nervous system. Besides the well-investigated function of dysfunctional glutamate receptor signaling in neurodegenerative diseases, neurodevelopmental disorders (NDD), like autism spectrum disorders (ASD) and attention-deficit and hyperactivity disorder (ADHD) might also be partly caused by disturbed glutamate signaling during development. However, the underlying mechanism of the type III metabotropic glutamate receptor 8 (mGluR8 or GRM8) involvement in neurodevelopment and disease mechanism is largely unknown. Here we show that the expression pattern of the two orthologs of human GRM8, grm8a and grm8b, have evolved partially distinct expression patterns in the brain of zebrafish (Danio rerio), especially at adult stages, suggesting sub-functionalization of these two genes during evolution. Using double in situ hybridization staining in the developing brain we demonstrate that grm8a is expressed in a subset of gad1a-positive cells, pointing towards glutamatergic modulation of GABAergic signaling. Building on this result we generated loss-of-function models of both genes using CRISPR/Cas9. Both mutant lines are viable and display no obvious gross morphological phenotypes making them suitable for further analysis. Initial behavioral characterization revealed distinct phenotypes in larvae. Whereas grm8a mutant animals display reduced swimming velocity, grm8b mutant animals show increased thigmotaxis behavior, suggesting an anxiety-like phenotype. We anticipate that our two novel metabotropic glutamate receptor 8 zebrafish models may contribute to a deeper understanding of its function in normal development and its role in the pathomechanism of disorders of the central nervous system.

Zebrafish Models of Cardiac Disease: From Fortuitous Mutants to Precision Medicine

Heart disease is the leading cause of death worldwide. Despite decades of research, most heart pathologies have limited treatments, and often the only curative approach is heart transplantation. Thus, there is an urgent need to develop new therapeutic approaches for treating cardiac diseases. Animal models that reproduce the human pathophysiology are essential to uncovering the biology of diseases and discovering therapies. Traditionally, mammals have been used as models of cardiac disease, but the cost of generating and maintaining new models is exorbitant, and the studies have very low throughput. In the last decade, the zebrafish has emerged as a tractable model for cardiac diseases, owing to several characteristics that made this animal popular among developmental biologists. Zebrafish fertilization and development are external; embryos can be obtained in high numbers, are cheap and easy to maintain, and can be manipulated to create new genetic models. Moreover, zebrafish exhibit an exceptional ability to regenerate their heart after injury. This review summarizes 25 years of research using the zebrafish to study the heart, from the classical forward screenings to the contemporary methods to model mutations found in patients with cardiac disease. We discuss the advantages and limitations of this model organism and introduce the experimental approaches exploited in zebrafish, including forward and reverse genetics and chemical screenings. Last, we review the models used to induce cardiac injury and essential ideas derived from studying natural regeneration. Studies using zebrafish have the potential to accelerate the discovery of new strategies to treat cardiac diseases.

Stage Specific Transcriptomic Analysis and Database for Zebrafish Oogenesis

Oogenesis produces functional eggs and is essential for fertility, embryonic development, and reproduction. The zebrafish ovary is an excellent model to study oogenesis in vertebrates, and recent studies have identified multiple regulators in oocyte development through forward genetic screens, as well as reverse genetics by CRISPR mutagenesis. However, many developmental steps in oogenesis, in zebrafish and other species, remain poorly understood, and their underlying mechanisms are unknown. Here, we take a genomic approach to systematically uncover biological activities throughout oogenesis. We performed transcriptomic analysis on five stages of oogenesis, from the onset of oocyte differentiation through Stage III, which precedes oocyte maturation. These transcriptomes revealed thousands of differentially expressed genes across stages of oogenesis. We analyzed trends of gene expression dynamics along oogenesis, as well as their expression in pair-wise comparisons between stages. We determined their functionally enriched terms, identifying uniquely characteristic biological activities in each stage. These data identified two prominent developmental phases in oocyte differentiation and traced the accumulation of maternally deposited embryonic regulator transcripts in the developing oocyte. Our analysis provides the first molecular description for oogenesis in zebrafish, which we deposit online as a resource for the community. Further, the presence of multiple gene paralogs in zebrafish, and the exclusive curation by many bioinformatic tools of the single paralogs present in humans, challenge zebrafish genomic analyses. We offer an approach for converting zebrafish gene name nomenclature to the human nomenclature for supporting genomic analyses generally in zebrafish. Altogether, our work provides a valuable resource as a first step to uncover oogenesis mechanisms and candidate regulators and track accumulating transcripts of maternal regulators of embryonic development.

Genome-wide mining of gpx gene family provides new insights into cadmium stress responses in common carp (Cyprinus carpio)

Glutathione peroxidase (Gpx) is an important member of antioxidant enzymes, which can play a vital role in metabolizing reactive oxygen species (ROS) and in maintaining cell homeostasis. In order to study the evolutionary dynamics of gpx gene family in allotetraploid fish species, we identified a total of 14 gpx genes in common carp Cyprinus carpio, while 9 gpx genes were discovered in the diploid progenitor-like species Poropuntius huangchuchieni. Comparative genomic analysis and phylogenetic analysis revealed that the common carp gpx genes had significant expansion and were divided into five distinct subclades. Exon-intron distribution patterns and conserved motif analysis revealed highly conserved evolutionary patterns. Transcript profiles suggested that different gpx genes had specific patterns of regulation during early embryonic development. In adult tissues, gpx genes had a relatively broad expression distribution, most of which were highly expressed in the gills, intestines, and gonads. RT-qPCR studies showed that most gpx genes were downregulated during the initial cd²⁺ treatment stage. Dietary supplementation of Bacillus coagulans at different concentrations (Group 2 of 1.0 × 10⁷ cfu/g, Group 3 of 1.0 × 10⁸ cfu/g, and Group 4 of 1.0 × 10⁹ cfu/g) induced different regulatory responses of gpx subclades. This result suggested that the appropriate concentration of B. coagulans can improve gpx gene expression when exposed to heavy metal cadmium treatment, which may play a vital role in the resistance to oxidative stress and immune responses. This study has expanded our understanding of the functional evolution of the gpx gene family in common carp.

Intraspecific de novo gene birth revealed by presence-absence variant genes in Caenorhabditis elegans

Genes embed their evolutionary history in the form of various alleles. Presence-absence variants (PAVs) are extreme cases of such alleles, where a gene present in one haplotype does not exist in another. Because PAVs may result from either birth or death of a gene, PAV genes and their alternative alleles, if available, can represent a basis for rapid intraspecific gene evolution. Using long-read sequencing technologies, this study traced the possible evolution of PAV genes in the PD1074 and CB4856 C. elegans strains as well as their alternative alleles in 14 other wild strains. We updated the CB4856 genome by filling 18 gaps and identified 46 genes and 7,460 isoforms from both strains not annotated previously. We verified 328 PAV genes, out of which 46 were C. elegans-specific. Among these possible newly born genes, 12 had alternative alleles in other wild strains; in particular, the alternative alleles of three genes showed signatures of active transposons. Alternative alleles of three other genes showed another type of signature reflected in accumulation of small insertions or deletions. Research on gene evolution using both species-specific PAV genes and their alternative alleles may provide new insights into the process of gene evolution.

Genome Structures and Evolution Analysis of Hsp90 Gene Family in Brassica napus Reveal the Possible Roles of Members in Response to Salt Stress and the Infection of Sclerotinia sclerotiorum

Heat shock proteins 90 (Hsp90s) are conserved proteins participating in the responses to heat stress and are found to be involved in different kinds of abiotic and biotic stresses. Brassica napus (B. napus) is an important heteropolyploid crop, producing edible oil. Salt stress is one of the most important hazards to the growth of rape in the world, while Sclerotinia stem rot is one of the most serious diseases, caused by Sclerotinia sclerotiorum (S. sclerotiorum). In this study, the evolution of Hsp90 genes and their responses to these two stresses were elucidated. Bioinformatic analysis through the whole genome of B. napus identified 35 Hsp90 gene family members. Five groups were obtained via phylogenetic analysis with the 35 Hsp genes, Hsps from its two ancestor species Brassica rapa, Brassica oleracea, and AtHsps. Gene structure and conservative motif analysis of these 35 Hsps indicated that the Hsps were relatively conservative in each group. Strong collinearity was also detected between the genomes of Brassica rapa, Brassica oleracea and B. napus, along with identifying syntenic gene pairs of Hsps among the three genomes. In addition, whole genome duplication was discovered as the main reason for the generation of BnHsp gene family. The analysis of cis-acting elements indicated that BnHsp90 might be involved in a variety of abiotic and biotic stress responses. Analysis of the expression pattern indicated that BnHsp90 participates in the responses of B. napus to salt stress and the infection of S. sclerotiorum. Fourteen and nine BnHsp90s were validated to be involved in the defense responses of B. napus against salt stress and S. sclerotiorum, respectively. Our results provide new insights for the roles of BnHsp90s in the responses of B. napus to salt stress and S. sclerotiorum.

Effect of Feed Supplementation with Bacillus coagulans on Nrf Gene Family Expression in Common Carp (Cyprinus carpio) under Long-Term Exposure to Cd2+

Nuclear factor-E2-related factor (Nrf) belongs to the Cap ‘n’ collar basic leucine zipper (CNC-bZIP) family, which plays an important role in the resistance to oxidative stress in the body. In this study, 12 Nrf genes were identified in the common carp genome database. Comparative genomic analysis showed that the Nrf genes of common carp had significant amplification, confirming that the common carp had experienced four genome-wide replication events. Phylogenetic analysis showed that all common carp Nrf clustered with scleractinian fish Nrf, indicating that they were highly conserved during evolution. In addition, tissue distribution results showed that most Nrf genes had a broad tissue distribution but exhibited tissue-specific expression patterns, demonstrating functional differences after WGD events. At 30 and 60 days of Cd2+ stress, most of the Nrf genes showed an increase in expression compared with the control group, indicating that they played a key role in the organism’s response to oxidative stress. To find a suitable concentration of Bacillus coagulans to activate the Nrf genes, we added three different concentrations (2.0 × 107 CFU/g, 2.0 × 108 CFU/g, and 2.0 × 109 CFU/g) of B. coagulans into the feed and defined them as L1, L2, and L3 groups, respectively. We investigated the effect of different concentrations of B. coagulans in the feed on the expression level of Nrf genes in the intestine of common carp under Cd2+ stress at 30 and 60 days. The results showed that, compared with the control/stress group, the expression of different Nrf genes was improved to varying degrees at three concentrations, and the effect of the L2 group (2.0 × 108 CFU/g) was the best. This suggests that the L2 group is the optimum concentration for activating Nrf gene expression when subjected to heavy metal Cd2+ stress and may act as an activation switch with a prominent role in the body’s resistance to oxidative stress and immune response.

Zebrafish as a Model to Study Vascular Elastic Fibers and Associated Pathologies

Many extensible tissues such as skin, lungs, and blood vessels require elasticity to function properly. The recoil of elastic energy stored during a stretching phase is provided by elastic fibers, which are mostly composed of elastin and fibrillin-rich microfibrils. In arteries, the lack of elastic fibers leads to a weakening of the vessel wall with an increased risk to develop cardiovascular defects such as stenosis, aneurysms, and dissections. The development of new therapeutic molecules involves preliminary tests in animal models that recapitulate the disease and whose response to drugs should be as close as possible to that of humans. Due to its superior in vivo imaging possibilities and the broad tool kit for forward and reverse genetics, the zebrafish has become an important model organism to study human pathologies. Moreover, it is particularly adapted to large scale studies, making it an attractive model in particular for the first steps of investigations. In this review, we discuss the relevance of the zebrafish model for the study of elastic fiber-related vascular pathologies. We evidence zebrafish as a compelling alternative to conventional mouse models.

Zebrafish vascular quantification: a tool for quantification of three-dimensional zebrafish cerebrovascular architecture by automated image analysis

Zebrafish transgenic lines and light sheet fluorescence microscopy allow in-depth insights into three-dimensional vascular development in vivo. However, quantification of the zebrafish cerebral vasculature in 3D remains highly challenging. Here, we describe and test an image analysis workflow for 3D quantification of the total or regional zebrafish brain vasculature, called zebrafish vasculature quantification (ZVQ). It provides the first landmark- or object-based vascular inter-sample registration of the zebrafish cerebral vasculature, producing population average maps allowing rapid assessment of intra- and inter-group vascular anatomy. ZVQ also extracts a range of quantitative vascular parameters from a user-specified region of interest, including volume, surface area, density, branching points, length, radius and complexity. Application of ZVQ to 13 experimental conditions, including embryonic development, pharmacological manipulations and morpholino-induced gene knockdown, shows that ZVQ is robust, allows extraction of biologically relevant information and quantification of vascular alteration, and can provide novel insights into vascular biology. To allow dissemination, the code for quantification, a graphical user interface and workflow documentation are provided. Together, ZVQ provides the first open-source quantitative approach to assess the 3D cerebrovascular architecture in zebrafish.

Summary: An image analysis workflow pipeline for 3D quantification of the total or regional zebrafish brain vasculature, called zebrafish vasculature quantification or ZVQ.

Regenerating zebrafish scales express a subset of evolutionary conserved genes involved in human skeletal disease

Background

Scales are mineralised exoskeletal structures that are part of the dermal skeleton. Scales have been mostly lost during evolution of terrestrial vertebrates whilst bony fish have retained a mineralised dermal skeleton in the form of fin rays and scales. Each scale is a mineralised collagen plate that is decorated with both matrix-building and resorbing cells. When removed, an ontogenetic scale is quickly replaced following differentiation of the scale pocket-lining cells that regenerate a scale. Processes promoting de novo matrix formation and mineralisation initiated during scale regeneration are poorly understood. Therefore, we performed transcriptomic analysis to determine gene networks and their pathways involved in dermal scale regeneration.

Results

We defined the transcriptomic profiles of ontogenetic and regenerating scales of zebrafish and identified 604 differentially expressed genes (DEGs). These were enriched for extracellular matrix, ossification, and cell adhesion pathways, but not in enamel or dentin formation processes indicating that scales are reminiscent to bone. Hypergeometric tests involving monogenetic skeletal disorders showed that DEGs were strongly enriched for human orthologues that are mutated in low bone mass and abnormal bone mineralisation diseases (P< 2× 10⁻³). The DEGs were also enriched for human orthologues associated with polygenetic skeletal traits, including height (P< 6× 10⁻⁴), and estimated bone mineral density (eBMD, P< 2× 10⁻⁵). Zebrafish mutants of two human orthologues that were robustly associated with height (COL11A2, P=6× 10⁻²⁴) or eBMD (SPP1, P=6× 10⁻²⁰) showed both exo- and endo- skeletal abnormalities as predicted by our genetic association analyses; col11a2^Y228X/Y228X mutants showed exoskeletal and endoskeletal features consistent with abnormal growth, whereas spp1^P160X/P160X mutants predominantly showed mineralisation defects.

Conclusion

We show that scales have a strong osteogenic expression profile comparable to other elements of the dermal skeleton, enriched in genes that favour collagen matrix growth. Despite the many differences between scale and endoskeletal developmental processes, we also show that zebrafish scales express an evolutionarily conserved sub-population of genes that are relevant to human skeletal disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01209-8.

The Role of Microtubule Associated Serine/Threonine Kinase 3 Variants in Neurodevelopmental Diseases: Genotype-Phenotype Association

Objective: To prove microtubule associated serine/threonine kinase 3 (MAST3) gene is associated with neurodevelopmental diseases (NDD) and the genotype-phenotype correlation.

Methods: Trio exome sequencing (trio ES) was performed on four NDD trios. Bioinformatic analysis was conducted based on large-scale genome sequencing data and human brain transcriptomic data. Further in vivo zebrafish studies were performed.

Results: In our study, we identified four de novo MAST3 variants (NM_015016.1: c.302C &gt; T:p.Ser101Phe; c.311C &gt; T:p.Ser104Leu; c.1543G &gt; A:p.Gly515Ser; and c.1547T &gt; C:p.Leu516Pro) in four patients with developmental and epileptic encephalopathy (DEE) separately. Clinical heterogeneities were observed in patients carrying variants in domain of unknown function (DUF) and serine-threonine kinase (STK) domain separately. Using the published large-scale exome sequencing data, higher CADD scores of missense variants in DUF domain were found in NDD cohort compared with gnomAD database. In addition, we obtained an excess of missense variants in DUF domain when compared autistic spectrum disorder (ASD) cohort with gnomAD database, similarly an excess of missense variants in STK domain when compared DEE cohort with gnomAD database. Based on Brainspan datasets, we showed that MAST3 expression was significantly upregulated in ASD and DEE-related brain regions and was functionally linked with DEE genes. In zebrafish model, abnormal morphology of central nervous system was observed in mast3a/b crispants.

Conclusion: Our results support the possibility that MAST3 is a novel gene associated with NDD which could expand the genetic spectrum for NDD. The genotype-phenotype correlation may contribute to future genetic counseling.

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.