Long-read sequencing of the zebrafish genome reorganizes genomic architecture

Skeletal Muscle Aging: Lessons From Teleosts

Aging is the greatest risk factor for a multitude of age-related diseases including sarcopenia-the loss of skeletal muscle mass and strength-which occurs at remarkable rates each year. There is an unmet need not only to understand the mechanisms that drive sarcopenia but also to identify novel therapeutic strategies. Given the ease and affordability of husbandry, along with advances in genomics, genome editing technologies, and imaging capabilities, teleost models are increasingly used for aging and sarcopenia research. Here, we explain how teleost species such as zebrafish, African turquoise killifish, and medaka recapitulate many of the classical hallmarks of sarcopenia, and discuss the various dietary, pharmacological, and genetic approaches that have been used in teleosts to understand the mechanistic basis of sarcopenia.

Lampbrush chromosomes of Danio rerio

Danio rerio, commonly known as zebrafish, is an established model organism for the developmental and cell biology studies. Although significant progress has been made in the analysis of the D. rerio genome, cytogenetic studies face challenges due to the unclear identification of chromosomes. Here, we present a novel approach to the study of the D. rerio karyotype, focusing on the analysis of lampbrush chromosomes isolated from growing oocytes. Lampbrush chromosomes, existing during diplotene, serve as a powerful tool for high-resolution mapping and transcription analysis due to their profound decondensation and remarkable lateral loops decorated by RNA polymerases and ribonucleoprotein (RNP) matrix. In D. rerio, lampbrush chromosomes are about 20 times longer than corresponding metaphase chromosomes. We found that the lampbrush chromosome stage karyotype of D. rerio is generally undifferentiated, except for several bivalents bearing distinct marker structures, including loops with complex RNP matrix and locus-associated nuclear bodies. Locus-associated nuclear bodies were enriched for coilin and snRNAs; the loci where they formed presumably correspond to the histone gene clusters. Further, we observed the accumulation of splicing factors in giant terminal RNP aggregates on one bivalent. DAPI staining of Danio rerio lampbrush chromosomes revealed large and small chromomeres non-uniformly distributed along the axis. For example, D. rerio lampbrush chromosome 4, comprising the sex-determining region, is divided into two halves-with small chromomeres bearing long lateral loops and with large dense chromomeres bearing no or very tiny lateral loops. As centromeres were not distinguishable, we identified centromeric regions in all bivalents by FISH mapping of pericentromeric RFAL1, RFAL2, and RFAM tandem repeats. Through a combination of morphological analysis, immunostaining of marker structures, and centromere mapping, we developed cytological maps of D. rerio lampbrush chromosomes. Finally, by RNA FISH we revealed transcripts of pericentromeric and telomeric tandem repeats at the lampbrush chromosome stage.

The dynamics and functional impact of tRNA repertoires during early embryogenesis in zebrafish

Embryogenesis entails dramatic shifts in mRNA translation and turnover that reprogram gene expression during cellular proliferation and differentiation. Codon identity modulates mRNA stability during early vertebrate embryogenesis, but how the composition of tRNA pools is matched to translational demand is unknown. By quantitative profiling of tRNA repertoires in zebrafish embryos during the maternal-to-zygotic transition, we show that zygotic tRNA repertoires are established after the onset of gastrulation, succeeding the major wave of zygotic mRNA transcription. Maternal and zygotic tRNA pools are distinct, but their reprogramming does not result in a better match to the codon content of the zygotic transcriptome. Instead, we find that an increase in global translation at gastrulation sensitizes decoding rates to tRNA supply, thus destabilizing maternal mRNAs enriched in slowly translated codons. Translational activation and zygotic tRNA expression temporally coincide with an increase of TORC1 activity at gastrulation, which phosphorylates and inactivates the RNA polymerase III repressor Maf1a/b. Our data indicate that a switch in global translation, rather than tRNA reprogramming, determines the onset of codon-dependent maternal mRNA decay during zebrafish embryogenesis.

Astyanax mexicanus surface and cavefish chromosome-scale assemblies for trait variation discovery

The ability of organisms to adapt to sudden extreme environmental changes produces some of the most drastic examples of rapid phenotypic evolution. The Mexican Tetra, Astyanax mexicanus, is abundant in the surface waters of northeastern Mexico, but repeated colonizations of cave environments have resulted in the independent evolution of troglomorphic phenotypes in several populations. Here, we present three chromosome-scale assemblies of this species, for one surface and two cave populations, enabling the first whole-genome comparisons between independently evolved cave populations to evaluate the genetic basis for the evolution of adaptation to the cave environment. Our assemblies represent the highest quality of sequence completeness with predicted protein-coding and noncoding gene metrics far surpassing prior resources and, to our knowledge, all long-read assembled teleost genomes, including zebrafish. Whole-genome synteny alignments show highly conserved gene order among cave forms in contrast to a higher number of chromosomal rearrangements when compared with other phylogenetically close or distant teleost species. By phylogenetically assessing gene orthology across distant branches of amniotes, we discover gene orthogroups unique to A. mexicanus. When compared with a representative surface fish genome, we find a rich amount of structural sequence diversity, defined here as the number and size of insertions and deletions as well as expanding and contracting repeats across cave forms. These new more complete genomic resources ensure higher trait resolution for comparative, functional, developmental, and genetic studies of drastic trait differences within a species.

High Diversity of Long Terminal Repeat Retrotransposons in Compact Vertebrate Genomes: Insights from Genomes of Tetraodontiformes

This study aimed to investigate the evolutionary profile (including diversity, activity, and abundance) of retrotransposons (RTNs) with long terminal repeats (LTRs) in ten species of Tetraodontiformes. These species, Arothron firmamentum, Lagocephalus sceleratus, Pao palembangensis, Takifugu bimaculatus, Takifugu flavidus, Takifugu ocellatus, Takifugu rubripes, Tetraodon nigroviridis, Mola mola, and Thamnaconus septentrionalis, are known for having the smallest genomes among vertebrates. Data mining revealed a high diversity and wide distribution of LTR retrotransposons (LTR-RTNs) in these compact vertebrate genomes, with varying abundances among species. A total of 819 full-length LTR-RTN sequences were identified across these genomes, categorized into nine families belonging to four different superfamilies: ERV (Orthoretrovirinae and Epsilon retrovirus), Copia, BEL-PAO, and Gypsy (Gmr, Mag, V-clade, CsRN1, and Barthez). The Gypsy superfamily exhibited the highest diversity. LTR family distribution varied among species, with Takifugu bimaculatus, Takifugu flavidus, Takifugu ocellatus, and Takifugu rubripes having the highest richness of LTR families and sequences. Additionally, evidence of recent invasions was observed in specific tetraodontiform genomes, suggesting potential transposition activity. This study provides insights into the evolution of LTR retrotransposons in Tetraodontiformes, enhancing our understanding of their impact on the structure and evolution of host genomes.

Recent insights into crosstalk between genetic parasites and their host genome

The bulk of higher order organismal genomes is comprised of transposable element (TE) copies, i.e. genetic parasites. The host-parasite relation is multi-faceted, varying across genomic region (genic versus intergenic), life-cycle stages, tissue-type and of course in health versus pathological state. The reach of functional genomics though, in investigating genotype-to-phenotype relations, has been limited when TEs are involved. The aim of this review is to highlight recent progress made in understanding how TE origin biochemical activity interacts with the central dogma stages of the host genome. Such interaction can also bring about modulation of the immune context and this could have important repercussions in disease state where immunity has a role to play. Thus, the review is to instigate ideas and action points around identifying evolutionary adaptations that the host genome and the genetic parasite have evolved and why they could be relevant.

Astyanax mexicanus surface and cavefish chromosome-scale assemblies for trait variation discovery

The ability of organisms to adapt to sudden extreme environmental changes produces some of the most drastic examples of rapid phenotypic evolution. The Mexican Tetra, Astyanax mexicanus, is abundant in the surface waters of northeastern Mexico, but repeated colonizations of cave environments have resulted in the independent evolution of troglomorphic phenotypes in several populations. Here, we present three chromosome-scale assemblies of this species, for one surface and two cave populations, enabling the first whole-genome comparisons between independently evolved cave populations to evaluate the genetic basis for the evolution of adaptation to the cave environment. Our assemblies represent the highest quality of sequence completeness with predicted protein-coding and non-coding gene metrics far surpassing prior resources and, to our knowledge, all long-read assembled teleost genomes, including zebrafish. Whole genome synteny alignments show highly conserved gene order among cave forms in contrast to a higher number of chromosomal rearrangements when compared to other phylogenetically close or distant teleost species. By phylogenetically assessing gene orthology across distant branches of amniotes, we discover gene orthogroups unique to A. mexicanus. When compared to a representative surface fish genome, we find a rich amount of structural sequence diversity, defined here as the number and size of insertions and deletions as well as expanding and contracting repeats across cave forms. These new more complete genomic resources ensure higher trait resolution for comparative, functional, developmental, and genetic studies of drastic trait differences within a species.

The revised reference genome of the leopard gecko (Eublepharis macularius) provides insight into the considerations of genome phasing and assembly

Genomic resources across squamate reptiles (lizards and snakes) have lagged behind other vertebrate systems and high-quality reference genomes remain scarce. Of the 23 chromosome-scale reference genomes across the order, only 12 of the ~60 squamate families are represented. Within geckos (infraorder Gekkota), a species-rich clade of lizards, chromosome-level genomes are exceptionally sparse representing only two of the seven extant families. Using the latest advances in genome sequencing and assembly methods, we generated one of the highest-quality squamate genomes to date for the leopard gecko, Eublepharis macularius (Eublepharidae). We compared this assembly to the previous, short-read only, E. macularius reference genome published in 2016 and examined potential factors within the assembly influencing contiguity of genome assemblies using PacBio HiFi data. Briefly, the read N50 of the PacBio HiFi reads generated for this study was equal to the contig N50 of the previous E. macularius reference genome at 20.4 kilobases. The HiFi reads were assembled into a total of 132 contigs, which was further scaffolded using HiC data into 75 total sequences representing all 19 chromosomes. We identified 9 of the 19 chromosomal scaffolds were assembled as a near-single contig, whereas the other 10 chromosomes were each scaffolded together from multiple contigs. We qualitatively identified that the percent repeat content within a chromosome broadly affects its assembly contiguity prior to scaffolding. This genome assembly signifies a new age for squamate genomics where high-quality reference genomes rivaling some of the best vertebrate genome assemblies can be generated for a fraction of previous cost estimates. This new E. macularius reference assembly is available on NCBI at JAOPLA010000000.

Zebrafish regulatory genomic resources for disease modelling and regeneration

In the past decades, the zebrafish has become a disease model with increasing popularity owing to its advantages that include fast development, easy genetic manipulation, simplicity for imaging, and sharing conserved disease-associated genes and pathways with those of human. In parallel, studies of disease mechanisms are increasingly focusing on non-coding mutations, which require genome annotation maps of regulatory elements, such as enhancers and promoters. In line with this, genomic resources for zebrafish research are expanding, producing a variety of genomic data that help in defining regulatory elements and their conservation between zebrafish and humans. Here, we discuss recent developments in generating functional annotation maps for regulatory elements of the zebrafish genome and how this can be applied to human diseases. We highlight community-driven developments, such as DANIO-CODE, in generating a centralised and standardised catalogue of zebrafish genomics data and functional annotations; consider the advantages and limitations of current annotation maps; and offer considerations for interpreting and integrating existing maps with comparative genomics tools. We also discuss the need for developing standardised genomics protocols and bioinformatic pipelines and provide suggestions for the development of analysis and visualisation tools that will integrate various multiomic bulk sequencing data together with fast-expanding data on single-cell methods, such as single-cell assay for transposase-accessible chromatin with sequencing. Such integration tools are essential to exploit the multiomic chromatin characterisation offered by bulk genomics together with the cell-type resolution offered by emerging single-cell methods. Together, these advances will build an expansive toolkit for interrogating the mechanisms of human disease in zebrafish.

Sequencing and assembling bear genomes: the bare necessities

Unique genetic adaptations are present in bears of every species across the world. From (nearly) shutting down important organs during hibernation to preventing harm from lifestyles that could easily cause metabolic diseases in humans, bears may hold the answer to various human ailments. However, only a few of these unique traits are currently being investigated at the molecular level, partly because of the lack of necessary tools. One of these tools is well-annotated genome assemblies from the different, extant bear species. These reference genomes are needed to allow us to identify differences in genetic variants, isoforms, gene expression, and genomic features such as transposons and identify those that are associated with biomedical-relevant traits. In this review we assess the current state of the genome assemblies of the eight different bear species, discuss current gaps, and the future benefits these reference genomes may have in informing human biomedical applications, while at the same time improving bear conservation efforts.

ONT-Based Alternative Assemblies Impact on the Annotations of Unique versus Repetitive Features in the Genome of a Romanian Strain of Drosophila melanogaster

To date, different strategies of whole-genome sequencing (WGS) have been developed in order to understand the genome structure and functions. However, the analysis of genomic sequences obtained from natural populations is challenging and the biological interpretation of sequencing data remains the main issue. The MinION device developed by Oxford Nanopore Technologies (ONT) is able to generate long reads with minimal costs and time requirements. These valuable assets qualify it as a suitable method for performing WGS, especially in small laboratories. The long reads resulted using this sequencing approach can cover large structural variants and repetitive sequences commonly present in the genomes of eukaryotes. Using MinION, we performed two WGS assessments of a Romanian local strain of Drosophila melanogaster, referred to as Horezu_LaPeri (Horezu). In total, 1,317,857 reads with a size of 8.9 gigabytes (Gb) were generated. Canu and Flye de novo assembly tools were employed to obtain four distinct assemblies with both unfiltered and filtered reads, achieving maximum reference genome coverages of 94.8% (Canu) and 91.4% (Flye). In order to test the quality of these assemblies, we performed a two-step evaluation. Firstly, we considered the BUSCO scores and inquired for a supplemental set of genes using BLAST. Subsequently, we appraised the total content of natural transposons (NTs) relative to the reference genome (ISO1 strain) and mapped the mdg1 retroelement as a resolution assayer. Our results reveal that filtered data provide only slightly enhanced results when considering genes identification, but the use of unfiltered data had a consistent positive impact on the global evaluation of the NTs content. Our comparative studies also revealed differences between Flye and Canu assemblies regarding the annotation of unique versus repetitive genomic features. In our hands, Flye proved to be moderately better for gene identification, while Canu clearly outperformed Flye for NTs analysis. Data concerning the NTs content were compared to those obtained with ONT for the D. melanogaster ISO1 strain, revealing that our strategy conducted to better results. Additionally, the parameters of our ONT reads and assemblies are similar to those reported for ONT experiments performed on various model organisms, revealing that our assembly data are appropriate for a proficient annotation of the Horezu genome.

Ongoing endeavors to detect mobilization of transposable elements

Transposable elements (TEs) are DNA sequences capable of mobilization from one location to another in the genome. Since the discovery of ‘Dissociation (Dc) locus’ by Barbara McClintock in maize (1), mounting evidence in the era of genomics indicates that a significant fraction of most eukaryotic genomes is composed of TE sequences, involving in various aspects of biological processes such as development, physiology, diseases and evolution. Although technical advances in genomics have discovered numerous functional impacts of TE across species, our understanding of TEs is still ongoing process due to challenges resulted from complexity and abundance of TEs in the genome. In this mini-review, we briefly summarize biology of TEs and their impacts on the host genome, emphasizing importance of understanding TE landscape in the genome. Then, we introduce recent endeavors especially in vivo retrotransposition assays and long read sequencing technology for identifying de novo insertions/TE polymorphism, which will broaden our knowledge of extraordinary relationship between genomic cohabitants and their host.

Ongoing endeavors to detect mobilization of transposable elements

Transposable elements (TEs) are DNA sequences capable of mobilization from one location to another in the genome. Since the discovery of 'Dissociation (Dc) locus' by Barbara McClintock in maize (1), mounting evidence in the era of genomics indicates that a significant fraction of most eukaryotic genomes is composed of TE sequences, involving in various aspects of biological processes such as development, physiology, diseases and evolution. Although technical advances in genomics have discovered numerous functional impacts of TE across species, our understanding of TEs is still ongoing process due to challenges resulted from complexity and abundance of TEs in the genome. In this mini-review, we briefly summarize biology of TEs and their impacts on the host genome, emphasizing importance of understanding TE landscape in the genome. Then, we introduce recent endeavors especially in vivo retrotransposition assays and long read sequencing technology for identifying de novo insertions/TE polymorphism, which will broaden our knowledge of extraordinary relationship between genomic cohabitants and their host. [BMB Reports 2022; 55(7): 305-315].