Comparative analysis of transposable elements highlights mobilome diversity and evolution in vertebrates

Genomic landscape of repetitive DNAs in Neotropical electric fishes

The Neotropical electric fishes of the order Gymnotiformes, known for their unique electrogenic and electrosensory systems, provide exceptional models for investigating ecological and evolutionary questions. While their phylogenetic and cytogenetic diversity is well documented, information about the diversity and evolution of repetitive DNA in gymnotiform genomes remains limited. To understand how repetitive DNA has shaped genome evolution in this group, we conducted bioinformatic analyses on raw sequencing data and genome assemblies from multiple species representing three major families. Our analysis revealed that closely related species share similar patterns of repetitive DNA composition, with Gymnotidae and Apteronotidae exhibiting the highest (16.5%) and lowest (9.4%) proportions of repetitive DNA, respectively. We identified 40 satellite DNA families, with five exclusive to Gymnotidae, twenty-two to Hypopomidae, and six shared between Apteronotidae species. Only one satellite DNA (NEFSat3-18) was conserved across all analyzed species, suggesting rapid evolutionary turnover of these sequences. The evolutionary dynamics of transposable elements varied among families, with Gymnotidae showing recent expansion of LINE elements and DNA transposons, while Apteronotidae displayed more ancient patterns of transposon activity. Analysis of transposable element landscapes revealed that all species experienced at least one burst of transposition during their evolution. Our findings demonstrate that repetitive DNA diversification parallels the evolutionary history of gymnotiform species. We also present a comprehensive dataset of repetitive DNA sequences that can be used as cytogenomic markers in comparative and evolutionary studies.

The landscape of transposable element distribution in the genome of Neotropical fish Apareiodon sp. (Characiformes: Parodontidae)

Transposable elements (TEs) are widely present in eukaryotic genomes, where they can contribute to genome size and functional modifications. As new genomes are sequenced and annotated, more studies can be conducted regarding TE content, distribution, and genome evolution. TEs are extensively diversified in fish genomes resulting in an important role in genome and chromosome evolution. However, curated TE libraries are still scarce in non-model organisms, making it difficult to evaluate TE's impact on genomic modifications thoroughly. Here, we aimed to obtain a curated TE library from the neotropical fish Apareiodon sp. genome. The prospection and curation of the TE library resulted in 244 families from 18 superfamilies of DNA transposons and retrotransposons, which comprise about 10% of the genome, with most insertions fitting in one or a few families. A greater diversity of retrotransposon families is present, especially for Ty3 superfamily. Despite the greater number of retrotransposon families, DNA transposons are the most abundant in the genome, with 37% of all TE insertions belonging to the Tc1-Mariner superfamily. Complete TE copies were observed for almost all superfamilies, with most of the sequences on the Tc1-Mariner group. DNA transposons and SINEs presented older insertions in the genome, followed by LINEs and LTR retrotransposons. TE genome density is highest in the cs25 scaffold, and enriched for Helitron elements. With these data, allied to previous studies on chromosome evolution, we suggest that cs25 bears the W chromosome specific region of the Apareiodon sp. genome, with the presence of significant amount of Helitron insertions.

Evolutionary Dynamics of Transposable Element Activity and Regulation in the Apennine Yellow-Bellied Toad (Bombina pachypus)

Transposable elements (TEs) expansion and accumulation represent one of the main drivers of genomic gigantism. Different host genome silencing mechanisms have evolved to counteract TE amplification, leading to a genomic arms race between them. Nevertheless, the evolutionary relationship between TEs and host genome silencing pathways remains poorly understood. Here, we investigate the activity of TEs and TE-silencing mechanisms in somatic and germline tissues of Bombina pachypus, a 10 Gb anuran genome. Our findings reveal a higher activity of TEs in the gonads compared to the brain, with retrotransposons as the most active class in both gonads (∼15% increased expression compared to the brain) and DNA transposons showing a 2-fold higher activity in ovaries. However, analysis of differentially expressed TEs between male and female gonads revealed a greater number of overexpressed TEs in testes (231 vs. 169), with maximum fold changes up to 22 in testes versus 8 in ovaries. This suggests a more permissive environment for TE expression in male gonads. Accordingly, increased activity of TE-silencing pathways was observed in ovaries compared to testes, with the KRAB-ZFP complex showing not only the highest overall expression levels but also a distinct ovary-specific expression pattern. Summarizing, while the higher TE activity in the male gonad may result from the lower efficiency of the KRAB-ZFP complex, the elevated activity of KRAB-ZFPs in ovaries, along with growing evidence of the functional role of TEs in the germline, suggests the existence of a broad range of host-TE dynamics going beyond the arms race model.

The Expansion and Diversification of Epigenetic Regulatory Networks Underpins Major Transitions in the Evolution of Land Plants

Epigenetic silencing is essential for regulating gene expression and cellular diversity in eukaryotes. While DNA and H3K9 methylation silence transposable elements (TEs), H3K27me3 marks deposited by the Polycomb repressive complex 2 (PRC2) silence varying proportions of TEs and genes across different lineages. Despite the major development role epigenetic silencing plays in multicellular eukaryotes, little is known about how epigenetic regulatory networks were shaped over evolutionary time. Here, we analyze epigenomes from diverse species across the green lineage to infer the chronological epigenetic recruitment of genes during land plant evolution. We first reveal the nature of plant heterochromatin in the unicellular chlorophyte microalga Chlorella sorokiniana and identify several genes marked with H3K27me3, highlighting the deep origin of PRC2-regulated genes in the green lineage. By incorporating genomic phylostratigraphy, we show how genes of differing evolutionary age occupy distinct epigenetic states in plants. While young genes tend to be silenced by H3K9 methylation, genes that emerged in land plants are preferentially marked with H3K27me3, some of which form part of a common network of PRC2-repressed genes across distantly related species. Finally, we analyze the potential recruitment of PRC2 to plant H3K27me3 domains and identify conserved DNA-binding sites of ancient transcription factor families known to interact with PRC2. Our findings shed light on the conservation and potential origin of epigenetic regulatory networks in the green lineage, while also providing insight into the evolutionary dynamics and molecular triggers that underlie the adaptation and elaboration of epigenetic regulation, laying the groundwork for its future consideration in other eukaryotic lineages.

Transposon insertion causes ctnnb2 transcript instability that results in the maternal effect zebrafish ichabod (ich) mutation

The maternal-effect mutation ichabod (ich) results in ventralized zebrafish embryos due to impaired induction of the dorsal canonical Wnt-signaling pathway. While previous studies linked the phenotype to reduced ctnnb2 transcript levels, the causative mutation remained unidentified. Using long-read sequencing, we discovered that the ich phenotype stems from the insertion of a non-autonomous CMC-Enhancer/Suppressor-mutator (CMC-EnSpm) transposon in the 3'UTR of the gene. Through reporter assays, we demonstrate that while wild type ctnnb2 mRNAs exhibit remarkably high stability throughout the early stages of development, the insertion of the transposon dramatically reduces transcript stability. Genome-wide mapping of the CMC-EnSpm transposons across multiple zebrafish strains also indicated ongoing transposition activity in the zebrafish genome. Our findings not only resolve the molecular basis of the ich mutation but also highlight the continuing mutagenic potential of endogenous transposons and reveal unexpected aspects of maternal transcript regulation during early zebrafish development.

Dynamic co-evolution of transposable elements and the piRNA pathway in African cichlid fishes

BACKGROUND

East African cichlid fishes have diversified in an explosive fashion, but the (epi)genetic basis of the phenotypic diversity of these fishes remains largely unknown. Although transposable elements (TEs) have been associated with phenotypic variation in cichlids, little is known about their transcriptional activity and epigenetic silencing. We set out to bridge this gap and to understand the interactions between TEs and their cichlid hosts.

RESULTS

Here, we describe dynamic patterns of TE expression in African cichlid gonads and during early development. Orthology inference revealed strong conservation of TE silencing factors in cichlids, and an expansion of piwil1 genes in Lake Malawi cichlids, likely driven by PiggyBac TEs. The expanded piwil1 copies have signatures of positive selection and retain amino acid residues essential for catalytic activity. Furthermore, the gonads of African cichlids express a Piwi-interacting RNA (piRNA) pathway that targets TEs. We define the genomic sites of piRNA production in African cichlids and find divergence in closely related species, in line with fast evolution of piRNA-producing loci.

CONCLUSIONS

Our findings suggest dynamic co-evolution of TEs and host silencing pathways in the African cichlid radiations. We propose that this co-evolution has contributed to cichlid genomic diversity.

Investigating the Impact of Whole-Genome Duplication on Transposable Element Evolution in Teleost Fishes

Transposable elements (TEs) can make up more than 50% of any given vertebrate's genome, with substantial variability in TE composition among lineages. TE variation is often linked to changes in gene regulation, genome size, and speciation. However, the role that genome duplication events have played in generating abrupt shifts in the composition of the mobilome over macroevolutionary timescales remains unclear. We investigated the degree to which the teleost genome duplication (TGD) shaped the diversification trajectory of the teleost mobilome. We integrate a new high coverage genome of Polypterus bichir with data from over 100 publicly available actinopterygian genomes to assess the macroevolutionary implications of genome duplication events on TE evolution in teleosts. Our results provide no evidence for a substantial shift in mobilome composition following the TGD event. Instead, the diversity of the teleost mobilome appears to have been shaped by a history of lineage-specific shifts in composition that are not correlated with commonly evoked drivers of diversification such as body size, water column usage, or latitude. Collectively, these results provide additional evidence for an emerging perspective that TGD did not catalyze bursts of diversification and innovation in the actinopterygian mobilome.

Genome sequence of the sugarcane aphid, Melanaphis sacchari (Hemiptera: Aphididae)

The sugarcane aphid, Melanaphis sacchari, is an agricultural pest that causes damage to plants in the Poaceae (the grasses) family, such as sorghum and sugarcane. In this study, we used nanopore long reads and a high-throughput chromosome conformation capture chromatin interaction maps to generate a chromosome-level assembly with a total length of 356.1 Mb, of which 85.5% (304.6 Mb) is contained within the 3 autosomes and the X chromosome. Repetitive sequences accounted for 16.29% of the chromosomes, and a total of 12,530 protein-coding genes were annotated, achieving 95.8% Benchmarking Universal Single-Copy Ortholog gene completeness. This offered a substantial improvement compared with previous low-quality genomic resources. A phylogenomic analysis by comparing M. sacchari with 24 published aphid genomes representing 3 aphid tribes revealed that M. sacchari belonged to the tribe Aphidini and maintained a conserved chromosome structure with other Aphidini species. The high-quality genomic resources reported in this study are useful for understanding the evolution of aphid genomes and studying pest management of M. sacchari.

Retrotransposon involves in photoperiodic spermatogenesis in Brandt's voles (Lasiopodomys brandtii) by co-transcription with flagellar genes

Photoperiod is a pivotal factor in affecting spermatogenesis in seasonal-breeding animals. Transposable elements have regulatory functions during spermatogenesis. However, whether it also functions in photoperiodic spermatogenesis in seasonal breeding animals is unknown. To explore this, we first annotated 5,501,822 transposons in the whole genome of Brandt's voles (Lasiopodomys brandtii), and revealed that LINEs were the most abundant, comprising 16.61 % of the genome. Following closely, SINEs accounted for 10.13 %, LTRs for 7.54 %, and DNA transposons for 0.70 %. Subsequently, we exposed male Brandt's voles to long-photoperiod (LP, 16 h/day) and short-photoperiod (SP, 8 h/day) from their embryonic stages, and obtained testes transcriptome at 4 and 10 weeks after birth. Differential expression and Pearson analysis indicated strongly positive correlations between the expression of differentially expressed retrotransposons and the adjacent genes. KO, KEGG and GSEA results showed that sperm flagellar genes were most enriched nearby the retrotransposons such as Dnah1, Dnah2, Dnah17, Dnali1. RT-PCR results showed that SINE/Alu_1213291 co-transcripted with Dnali1 gene. Our findings first reveal the regulatory function of transposons in photoperiodic spermatogenesis, providing insights into the role of photoperiod in seasonal reproduction in wild animals.

Transposable Element Diversity and Activity Patterns in Neotropical Salamanders

Transposable elements (TEs) compose a substantial proportion of the largest eukaryotic genomes. TE diversity has been hypothesized to be negatively correlated with genome size, yet empirical demonstrations of such a relationship in a phylogenetic context are largely lacking. Furthermore, the most abundant type of TEs in genomes varies across groups, and it is not clear if there are patterns of TE activity consistent with genome size among different taxa with large genome sizes. We use low-coverage sequencing of 16 species of Neotropical salamanders, which vary ∼7-fold in genome size, to estimate TE relative abundance and diversity for each species. We also compare the divergence of copies of each TE superfamily to estimate patterns of TE activity in each species. We find a negative relationship between TE diversity and genome size, which is consistent with the hypothesis that either competition among TEs or reduced selection against ectopic recombination may result in lower diversity in the largest genomes. We also find divergent activity patterns in the largest versus the smallest genomes, suggesting that the history of TE activity may explain differences in genome size. Our results suggest that both TE diversity and relative abundance may be predictable, at least within taxonomic groups.

Clam Genome and Transcriptomes Provide Insights into Molecular Basis of Morphological Novelties and Adaptations in Mollusks

Bivalve mollusks, comprising animals enclosed in two shell valves, are well-adapted to benthic life in many intertidal zones. Clams have evolved the buried lifestyle, which depends on their unique soft tissue structure and their wedge-shaped muscular foot and long extendible siphons. However, molecular mechanisms of adaptative phenotype evolution remain largely unknown. In the present study, we obtain the high-quality chromosome-level genome of Manila clam R. philippinarum, an economically important marine bivalve in many coastal areas. The genome is constructed by the Hi-C assisted assembly, which yields 19 chromosomes with a total of 1.17 Gb and BUSCO integrity of 92.23%. The de novo assembled genome has a contig N50 length of 307.7 kb and scaffold N50 of 59.5 Mb. Gene family expansion analysis reveals that a total of 24 single-copy gene families have undergone the significant expansion or contraction, including E3 ubiquitin ligase and dynein heavy chain. The significant expansion of transposable elements has been also identified, including long terminal repeats (LTR) and non-LTR retrotransposons. The comparative transcriptomics among different clam tissues reveals that extracellular matrix (ECM) receptors and neuroactive ligand receptors may play the important roles in tissue structural support and neurotransmission during their infaunal life. These findings of gene family expansion and tissue-specific expression may reflect the unique soft tissue structure of clams, suggesting the evolution of lineage-specific morphological novelties. The high-quality genome and transcriptome data of R. philippinarum will not only facilitate the genetic studies on clams but will also provide valuable information on morphological novelties in mollusks.

Characterization of hAT DNA transposon superfamily in the genome of Neotropical fish Apareiodon sp

DNA transposons are diverse in fish genomes and have been described to generate genomic evolutionary novelties. hAT transposable element data are scarce in Teleostei genomes, making it challenging to conduct comparative genomic studies to understand their neutrality or function. This study aimed to perform a genomic and molecular characterization of hAT copies to assess the diversity of these elements and associate changes in these sequences to genomic and karyotypic novelties in Apareiodon sp. The data revealed that hAT TEs are highly abundant in the Apareiodon sp. genome, with few possibly autonomous copies. Highly conserved sequences with likely functional transposases were observed in nine hAT elements. A great diversity of hAT subgroups was observed, especially from Ac, Charlie, Blackjack, Tip100, hAT6, and hAT5, and a similar wave of hAT genomic invasion was identified in the genome for these six groups of hAT sequences. The data also revealed a distinct number of microsatellites within degenerated hAT copies. hAT sites were demonstrated to be dispersed in the Apareiodon sp. chromosomes and not involved in W chromosome-specific region differentiation. In conclusion, the genomic analysis revealed a great diversity of hAT elements, possible autonomous copies, and differentiation of degenerated transposable elements into tandem sequences.

The emerging H3K9me3 chromatin landscape during zebrafish embryogenesis

The structural organization of eukaryotic genomes is contingent upon the fractionation of DNA into transcriptionally permissive euchromatin and repressive heterochromatin. However, we have a limited understanding of how these distinct states are first established during animal embryogenesis. Histone 3 lysine 9 trimethylation (H3K9me3) is critical to heterochromatin formation, and bulk establishment of this mark is thought to help drive large-scale remodeling of an initially naive chromatin state during animal embryogenesis. However, a detailed understanding of this process is lacking. Here, we leverage CUT&RUN to define the emerging H3K9me3 landscape of the zebrafish embryo with high sensitivity and temporal resolution. Despite the prevalence of DNA transposons in the zebrafish genome, we found that LTR transposons are preferentially targeted for embryonic H3K9me3 deposition, with different families exhibiting distinct establishment timelines. High signal-to-noise ratios afforded by CUT&RUN revealed new, emerging sites of low-amplitude H3K9me3 that initiated before the major wave of zygotic genome activation (ZGA). Early sites of establishment predominated at specific subsets of transposons and were particularly enriched for transposon sequences with maternal piRNAs and pericentromeric localization. Notably, the number of H3K9me3 enriched sites increased linearly across blastula development, while quantitative comparison revealed a >10-fold genome-wide increase in H3K9me3 signal at established sites over just 30 min at the onset of major ZGA. Continued maturation of the H3K9me3 landscape was observed beyond the initial wave of bulk establishment.

The genomes of all lungfish inform on genome expansion and tetrapod evolution

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri)1 and African2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers.

Comparative analysis of transposable elements dynamics in fish with different sex chromosome systems

Transposable elements (TEs) are widespread genomic components with substantial roles in genome evolution and sex chromosome differentiation. In this study, we compared the TE composition of three closely related fish with different sex chromosome systems: Megaleporinus elongatus (Z1Z1Z2Z2/Z1W1Z2W2), Megaleporinus macrocephalus (ZZ/ZW) (both with highly differentiated W sex chromosomes), and Leporinus friderici (without heteromorphic sex chromosomes). We created custom TE libraries for each species using clustering methods and manual annotation and prediction, and we predicted TE temporal dynamics through divergence-based analysis. The TE abundance ranged from 16% to 21% in the three mobilomes, with L. friderici having the lowest overall. Despite the recent amplification of TEs in all three species, we observed differing expansion activities, particularly between the two genera. Both Megaleporinus recently experienced high retrotransposon activity, with a reduction in DNA TEs, which could have implications in sex chromosome composition. In contrast, L. friderici showed the opposite pattern. Therefore, despite having similar TE compositions, Megaleporinus and Leporinus exhibit distinct TE histories that likely evolved after their separation, highlighting a rapid TE expansion over short evolutionary periods.

The first chromosome-level genome assembly and transcriptome sequencing provide insights into cantharidin production of the blister beetles

Blister beetles (Coleoptera: Meloidae) produce a natural defensive toxin cantharidin (CTD), which has been used for various cancer treatments and other diseases. Currently, the lack of chromosome-level reference genomes in Meloidae limits further understanding of the mechanism of CTD biosynthesis and environmental adaptation. In this study, the chromosome-level genome assembly of Mylabris phalerata was generated based on PacBio and Hi-C sequencing. This reference genome was about 136.68 Mb in size with contig N50 of 9.17 Mb and composed of 12 chromosomes. In comparison to six other Coleoptera insects, M. phalerata exhibited multiple expanded gene families enriched in juvenile hormone (JH) biosynthetic process pathway, farnesol dehydrogenase activity, and cytochrome P450, which may be related to CTD biosynthesis. Consistently, the transcriptomic analysis suggested the "terpenoid backbone biosynthesis" pathway and "the juvenile hormone" as putative core pathways of CTD biosynthesis and presented eight up-regulated differential expression genes in male adults as candidate genes. It is possible that the restricted feeding niche and lifestyle of M. phalerata were the cause of the gene family's contraction of odorant binding proteins. The ABC transporters (ABCs) related to exporting bound toxins out of the cell and the resistance to the self-secreted toxins (e.g. CTD) were also contracted, possibly due to other self-protection strategies in M. phalerata. A foundation of understanding CTD biosynthesis and environmental adaptation of blister beetles will be established by our reference genome and discoveries.

Differential Conservation and Loss of Chicken Repeat 1 (CR1) Retrotransposons in Squamates Reveal Lineage-Specific Genome Dynamics Across Reptiles

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniote vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; including ∼11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific genome dynamics have evolved over the course of squamate evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the chicken repeat 1 (CR1) retrotransposon, a TE family found in most tetrapod genomes which is the dominant TE in most reptiles. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

Expression of Transposable Elements throughout the Fasciola hepatica Trematode Life Cycle

BACKGROUND

Transposable elements (TEs) are major components of eukaryotic genomes. The extensive body of evidence suggests that although they were once considered "genomic parasites", transposons and their transcripts perform specific functions, such as regulation of early embryo development. Understanding the role of TEs in such parasites as trematodes is becoming critically important. Fasciola hepatica, a parasite affecting humans and livestock, undergoes a complex life cycle in diverse environments and hosts, and knowledge about its life cycle regulation is scarce so far.

METHODS

We summarized the data regarding the repetitive elements in F. hepatica and conducted bulk RNA-seq analysis across its life cycle stages. TE expression profiles were analyzed, focusing on differential expression and potential homology with previously described long non-coding RNAs (lncRNAs).

RESULTS

Differential expression analysis revealed stage-specific TE transcription patterns, notably peaking during egg and metacercariae stages. Some TEs showed homology with known lncRNAs and contained putative transcription factor binding sites. Interestingly, TE transcription levels were highest in eggs and metacercariae compared to adults, suggesting regulatory roles in trematode life cycle transitions.

CONCLUSIONS

These findings suggest that TEs may play roles in regulating trematode life cycle transitions. Moreover, TE homology with lncRNAs underscores their significance in gene regulation.

Genome Assembly of the Dyeing Poison Frog Provides Insights into the Dynamics of Transposable Element and Genome-Size Evolution

Genome size varies greatly across the tree of life and transposable elements are an important contributor to this variation. Among vertebrates, amphibians display the greatest variation in genome size, making them ideal models to explore the causes and consequences of genome size variation. However, high-quality genome assemblies for amphibians have, until recently, been rare. Here, we generate a high-quality genome assembly for the dyeing poison frog, Dendrobates tinctorius. We compare this assembly to publicly available frog genomes and find evidence for both large-scale conserved synteny and widespread rearrangements between frog lineages. Comparing conserved orthologs annotated in these genomes revealed a strong correlation between genome size and gene size. To explore the cause of gene-size variation, we quantified the location of transposable elements relative to gene features and find that the accumulation of transposable elements in introns has played an important role in the evolution of gene size in D. tinctorius, while estimates of insertion times suggest that many insertion events are recent and species-specific. Finally, we carry out population-scale mobile-element sequencing and show that the diversity and abundance of transposable elements in poison frog genomes can complicate genotyping from repetitive element sequence anchors. Our results show that transposable elements have clearly played an important role in the evolution of large genome size in D. tinctorius. Future studies are needed to fully understand the dynamics of transposable element evolution and to optimize primer or bait design for cost-effective population-level genotyping in species with large, repetitive genomes.

FishTEDB 2.0: an update fish transposable element (TE) database with new functions to facilitate TE research

We launched the initial version of FishTEDB in 2018, which aimed to establish an open-source, user-friendly, data-rich transposable element (TE) database. Over the past 5 years, FishTEDB 1.0 has gained approximately 10 000 users, accumulating more than 450 000 interactions. With the unveiling of extensive fish genome data and the increasing emphasis on TE research, FishTEDB needs to extend the richness of data and functions. To achieve the above goals, we introduced 33 new fish species to FishTEDB 2.0, encompassing a wide array of fish belonging to 48 orders. To make the updated database more functional, we added a genome browser to visualize the positional relationship between TEs and genes and the estimated TE insertion time in different species. In conclusion, we released a new version of the fish TE database, FishTEDB 2.0, designed to assist researchers in the future study of TE functions and promote the progress of biological theories related to TEs. Database URL: https://www.fishtedb.com/.

Transposable Element Expression Profiles in Premalignant Pigment Cell Lesions and Melanoma of Xiphophorus

Transposable elements (TEs) are characterized by their ability to change their genomic position. Through insertion or recombination leading to deletions and other chromosomal aberrations, they can cause genetic instability. The extent to which they thereby exert regulatory influence on cellular functions is unclear. To better characterize TEs in processes such as carcinogenesis, we used the well-established Xiphophorus melanoma model. By transcriptome sequencing, we show that an increasing total number in transposons correlates with progression of malignancy in melanoma samples from Xiphophorus interspecific hybrids. Further, by comparing the presence of TEs in the parental genomes of Xiphophorus maculatus and Xiphophorus hellerii, we could show that even in closely related species, genomic location and spectrum of TEs are considerably different.

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.

The emerging H3K9me3 chromatin landscape during zebrafish embryogenesis

The structural organization of eukaryotic genomes is contingent upon the fractionation of DNA into transcriptionally permissive euchromatin and repressive heterochromatin. However, we have a limited understanding of how these distinct states are first established during animal embryogenesis. Histone 3 lysine 9 trimethylation (H3K9me3) is critical to heterochromatin formation and bulk establishment of this mark is thought to help drive large-scale remodeling of an initially naive chromatin state during animal embryogenesis. However, a detailed understanding of this process is lacking. Here, we leverage CUT&RUN to define the emerging H3K9me3 landscape of the zebrafish embryo with high sensitivity and temporal resolution. Despite the prevalence of DNA transposons in the zebrafish genome, we found that LTR transposons are preferentially targeted for embryonic H3K9me3 deposition, with different families exhibiting distinct establishment timelines. High signal-to-noise ratios afforded by CUT&RUN revealed new, emerging sites of low-amplitude H3K9me3 that initiated before the major wave of zygotic genome activation (ZGA). Early sites of establishment predominated at specific subsets of transposons and were particularly enriched for transposon sequences with maternal piRNAs and pericentromeric localization. Notably, the number of H3K9me3 enriched sites increased linearly across blastula development, while quantitative comparison revealed a >10-fold genome-wide increase in H3K9me3 signal at established sites over just 30 minutes at the onset of ZGA. Continued maturation of the H3K9me3 landscape was observed beyond the initial wave of bulk establishment.

Transposable elements impact the population divergence of rice blast fungus Magnaporthe oryzae

Dynamic transposition of transposable elements (TEs) in fungal pathogens has significant impact on genome stability, gene expression, and virulence to the host. In Magnaporthe oryzae, genome plasticity resulting from TE insertion is a major driving force leading to the rapid evolution and diversification of this fungus. Despite their importance in M. oryzae population evolution and divergence, our understanding of TEs in this context remains limited. Here, we conducted a genome-wide analysis of TE transposition dynamics in the 11 most abundant TE families in M. oryzae populations. Our results show that these TEs have specifically expanded in recently isolated M. oryzae rice populations, with the presence/absence polymorphism of TE insertions highly concordant with population divergence on Geng/Japonica and Xian/Indica rice cultivars. Notably, the genes targeted by clade-specific TEs showed clade-specific expression patterns and are involved in the pathogenic process, suggesting a transcriptional regulation of TEs on targeted genes. Our study provides a comprehensive analysis of TEs in M. oryzae populations and demonstrates a crucial role of recent TE bursts in adaptive evolution and diversification of the M. oryzae rice-infecting lineage.

IMPORTANCE

Magnaporthe oryzae is the causal agent of the destructive blast disease, which caused massive loss of yield annually worldwide. The fungus diverged into distinct clades during adaptation toward the two rice subspecies, Xian/Indica and Geng/Japonica. Although the role of TEs in the adaptive evolution was well established, mechanisms underlying how TEs promote the population divergence of M. oryzae remain largely unknown. In this study, we reported that TEs shape the population divergence of M. oryzae by differentially regulating gene expression between Xian/Indica-infecting and Geng/Japonica-infecting populations. Our results revealed a TE insertion-mediated gene expression adaption that led to the divergence of M. oryzae population infecting different rice subspecies.

Multiple horizontal transfer events of a DNA transposon into turtles, fishes, and a frog

Horizontal transfer of transposable elements (HTT) has been reported across many species and the impact of such events on genome structure and function has been well described. However, few studies have focused on reptilian genomes, especially HTT events in Testudines (turtles). Here, as a consequence of investigating the repetitive content of Malaclemys terrapin terrapin (Diamondback turtle) we found a high similarity DNA transposon, annotated in RepBase as hAT-6_XT, shared between other turtle species, ray-finned fishes, and a frog. hAT-6_XT was notably absent in reptilian taxa closely related to turtles, such as crocodiles and birds. Successful invasion of DNA transposons into new genomes requires the conservation of specific residues in the encoded transposase, and through structural analysis, these residues were identified indicating some retention of functional transposition activity. We document six recent independent HTT events of a DNA transposon in turtles, which are known to have a low genomic evolutionary rate and ancient repeats.

Dynamic co-evolution of transposable elements and the piRNA pathway in African cichlid fishes

East African cichlid fishes have diversified in an explosive fashion, but the (epi)genetic basis of the phenotypic diversity of these fishes remains largely unknown. Although transposable elements (TEs) have been associated with phenotypic variation in cichlids, little is known about their transcriptional activity and epigenetic silencing. Here, we describe dynamic patterns of TE expression in African cichlid gonads and during early development. Orthology inference revealed an expansion of piwil1 genes in Lake Malawi cichlids, likely driven by PiggyBac TEs. The expanded piwil1 copies have signatures of positive selection and retain amino acid residues essential for catalytic activity. Furthermore, the gonads of African cichlids express a Piwi-interacting RNA (piRNA) pathway that target TEs. We define the genomic sites of piRNA production in African cichlids and find divergence in closely related species, in line with fast evolution of piRNA-producing loci. Our findings suggest dynamic co-evolution of TEs and host silencing pathways in the African cichlid radiations. We propose that this co-evolution has contributed to cichlid genomic diversity.

Genome sequencing and analysis of black flounder (Paralichthys orbignyanus) reveals new insights into Pleuronectiformes genomic size and structure

Black flounder (Paralichthys orbignyanus, Pleuronectiformes) is a commercially significant marine fish with promising aquaculture potential in Argentina. Despite extensive studies on Black flounder aquaculture, its limited genetic information available hampers the crucial role genetics plays in the development of this activity. In this study, we first employed Illumina sequencing technology to sequence the entire genome of Black flounder. Utilizing two independent libraries-one from a female and another from a male-with 150 bp paired-end reads, a mean insert length of 350 bp, and over 35 X-fold coverage, we achieved assemblies resulting in a genome size of ~ 538 Mbp. Analysis of the assemblies revealed that more than 98% of the core genes were present, with more than 78% of them having more than 50% coverage. This indicates a somehow complete and accurate genome at the coding sequence level. This genome contains 25,231 protein-coding genes, 445 tRNAs, 3 rRNAs, and more than 1,500 non-coding RNAs of other types. Black flounder, along with pufferfishes, seahorses, pipefishes, and anabantid fish, displays a smaller genome compared to most other teleost groups. In vertebrates, the number of transposable elements (TEs) is often correlated with genome size. However, it remains unclear whether the sizes of introns and exons also play a role in determining genome size. Hence, to elucidate the potential factors contributing to this reduced genome size, we conducted a comparative genomic analysis between Black flounder and other teleost orders to determine if the small genomic size could be explained by repetitive elements or gene features, including the whole genome genes and introns sizes. We show that the smaller genome size of flounders can be attributed to several factors, including changes in the number of repetitive elements, and decreased gene size, particularly due to lower amount of very large and small introns. Thus, these components appear to be involved in the genome reduction in Black flounder. Despite these insights, the full implications and potential benefits of genome reduction in Black flounder for reproduction and aquaculture remain incompletely understood, necessitating further research.

Evolution and expression patterns of the neo-sex chromosomes of the crested ibis

Bird sex chromosomes play a unique role in sex-determination, and affect the sexual morphology and behavior of bird species. Core waterbirds, a major clade of birds, share the common characteristics of being sexually monomorphic and having lower levels of inter-sexual conflict, yet their sex chromosome evolution remains poorly understood. Here, by we analyse of a chromosome-level assembly of a female crested ibis (Nipponia nippon), a typical core waterbird. We identify neo-sex chromosomes resulting from fusion of microchromosomes with ancient sex chromosomes. These fusion events likely occurred following the divergence of Threskiornithidae and Ardeidae. The neo-W chromosome of the crested ibis exhibits the characteristics of slow degradation, which is reflected in its retention of abundant gametologous genes. Neo-W chromosome genes display an apparent ovary-biased gene expression, which is largely driven by genes that are retained on the crested ibis W chromosome but lost in other bird species. These results provide new insights into the evolutionary history and expression patterns for the sex chromosomes of bird species.

De Novo Genome Assembly of the Whitespot Parrotfish (Scarus forsteni): A Valuable Scaridae Genomic Resource

Scarus forsteni, a whitespot parrotfish from the Scaridae family, is a herbivorous fish inhabiting coral reef ecosystems. The deterioration of coral reefs has highly affected the habitats of the parrotfish. The decline in genetic diversity of parrotfish emphasizes the critical importance of conserving their genetic variability to ensure the resilience and sustainability of marine ecosystems for future generations. In this study, a genome of S. forsteni was assembled de novo through using Illumina and Nanopore sequencing. The 1.71-Gb genome of S. forsteni, was assembled into 544 contigs (assembly level: contig). It exhibited an N50 length of 17.97 Mb and a GC content percentage of 39.32%. Our BUSCO analysis revealed that the complete protein of the S. forsteni genome had 98.10% integrity. Combined with structure annotation data, 34,140 (74.81%) genes were functionally annotated out of 45,638 predicted protein-coding genes. Upon comparing the genome size and TE content of teleost fishes, a roughly linear relationship was observed between these two parameters. However, TE content is not a decisive factor in determining the genome size of S. forsteni. Population history analysis results indicate that S. forsteni experienced two major population expansions, both of which occurred before the last interglacial period. In addition, through a comparative genomic analysis of the evolutionary relationship of other species, it was found that S. forsteni had the closest relationship with Cheilinus undulatus, another member of the Labridae family. Our expansion and contraction analysis of the gene family showed that the expansion genes were mainly associated with immune diseases, organismal systems, and cellular processes. At the same time, cell transcription and translation, sex hormone regulation, and other related pathways were also more prominent in the positive selection genes. The genomic sequence of S. forsteni offers valuable resources for future investigations on the conservation, evolution, and behavior of fish species.

Differential Conservation and Loss of CR1 Retrotransposons in Squamates Reveals Lineage-Specific Genome Dynamics across Reptiles

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniotic vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; ~11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific dynamics have evolved over the course of squamate evolution to constrain genome size across the order. Thus, squamates may represent a prime model for investigations into TE diversity and evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the CR1 retrotransposon, a TE family found in most tetrapod genomes. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds, and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

De novo assembly and characterization of a highly degenerated ZW sex chromosome in the fish Megaleporinus macrocephalus

BACKGROUND

Megaleporinus macrocephalus (piauçu) is a Neotropical fish within Characoidei that presents a well-established heteromorphic ZZ/ZW sex determination system and thus constitutes a good model for studying W and Z chromosomes in fishes. We used PacBio reads and Hi-C to assemble a chromosome-level reference genome for M. macrocephalus. We generated family segregation information to construct a genetic map, pool sequencing of males and females to characterize its sex system, and RNA sequencing to highlight candidate genes of M. macrocephalus sex determination.

RESULTS

The reference genome of M. macrocephalus is 1,282,030,339 bp in length and has a contig and scaffold N50 of 5.0 Mb and 45.03 Mb, respectively. In the sex chromosome, based on patterns of recombination suppression, coverage, FST, and sex-specific SNPs, we distinguished a putative W-specific region that is highly differentiated, a region where Z and W still share some similarities and is undergoing degeneration, and the PAR. The sex chromosome gene repertoire includes genes from the TGF-β family (amhr2, bmp7) and the Wnt/β-catenin pathway (wnt4, wnt7a), some of which are differentially expressed.

CONCLUSIONS

The chromosome-level genome of piauçu exhibits high quality, establishing a valuable resource for advancing research within the group. Our discoveries offer insights into the evolutionary dynamics of Z and W sex chromosomes in fish, emphasizing ongoing degenerative processes and indicating complex interactions between Z and W sequences in specific genomic regions. Notably, amhr2 and bmp7 are potential candidate genes for sex determination in M. macrocephalus.

Dispersal from the Qinghai-Tibet plateau by a high-altitude butterfly is associated with rapid expansion and reorganization of its genome

Parnassius glacialis is a typical "Out of the QTP" alpine butterfly that originated on the Qinghai-Tibet Plateau (QTP) and dispersed into relatively low-altitude mountainous. Here we assemble a chromosome-level genome of P. glacialis and resequence 9 populations in order to explore the genome evolution and local adaptation of this species. These results indicated that the rapid accumulation and slow unequal recombination of transposable elements (TEs) contributed to the formation of its large genome. Several ribosomal gene families showed extensive expansion and selective evolution through transposon-mediated processed pseudogenes. Additionally, massive structural variations (SVs) of TEs affected the genetic differentiation of low-altitude populations. These low-altitude populations might have experienced a genetic bottleneck in the past and harbor genes with selective signatures which may be responsible for the potential adaptation to low-altitude environments. These results provide a foundation for understanding genome evolution and local adaptation for "Out of the QTP" of P. glacialis.

Genome-wide SNP analysis provides insights into the XX/XY sex-determination system in silver barb (Barbonymus gonionotus)

Silver barb (Barbonymus gonionotus) is among the most economically important freshwater fish species in Thailand. It ranks fourth in economic value and third in production weight for fisheries and culture in Thailand. An XX/XY sex-determination system based on gynogenesis was previously reported for this fish. In this study, the molecular basis underlying the sex-determination system was further investigated. Genome-wide single-nucleotide polymorphism data were generated for 32 captive-bred silver barb individuals, previously scored by phenotypic sex, to identify sex-linked regions associated with sex determination. Sixty-three male-linked loci, indicating putative XY chromosomes, were identified. Male-specific loci were not observed, which indicates that the putative Y chromosome is young and the sex determination region is cryptic. A homology search revealed that most male-linked loci were homologous to the Mariner/Tc1 and Gypsy transposable elements and are probably the remnants of an initial accumulation of repeats on the Y chromosome from the early stages of sex chromosome differentiation. This research provides convincing insights into the mechanism of sex determination and reveals the potential sex determination regions in silver barb. The study provides the basic data necessary for increasing the commercial value of silver barbs through genetic improvements.

Comparative Analysis of Transposable Elements Reveals the Diversity of Transposable Elements in Decapoda and Their Effects on Genomic Evolution

Transposable elements (TEs) are mobile genetic elements that exist in the host genome and exert considerable influence on the evolution of the host genome. Since crustaceans, including decapoda, are considered ideal models for studying the relationship between adaptive evolution and TEs, TEs were identified and classified in the genomes of eight decapoda species and one diplostraca species (as the outgroup) using two strategies, namely homology-based annotation and de novo annotation. The statistics and classification of TEs showed that their proportion in the genome and their taxonomic composition in decapoda were different. Moreover, correlation analysis and transcriptome data demonstrated that there were more PIF-Harbinger TEs in the genomes of Eriocheir sinensis and Scylla paramamosain, and the expression patterns of PIF-Harbingers were significantly altered under air exposure stress conditions. These results signaled that PIF-Harbingers expanded in the genome of E. sinensis and S. paramamosain and might be related to their air exposure tolerance levels. Meanwhile, sequence alignment revealed that some Jockey-like sequences (JLSs) with high similarity to specific regions of the White spot syndrome virus (WSSV) genome existed in all eight decapod species. At the same time, phylogenetic comparison exposed that the phylogenetic tree constructed by JLSs was not in agreement with that of the species tree, and the distribution of each branch was significantly different. The abovementioned results signaled that these WSSV-specific JLSs might transfer horizontally and contribute to the emergence of WSSV. This study accumulated data for expanding research on TEs in decapod species and also provided new insights and future direction for the breeding of stress-resistant and disease-resistant crab breeds.

Widespread Deviant Patterns of Heterozygosity in Whole-Genome Sequencing Due to Autopolyploidy, Repeated Elements, and Duplication

Most population genomic tools rely on accurate single nucleotide polymorphism (SNP) calling and filtering to meet their underlying assumptions. However, genomic complexity, resulting from structural variants, paralogous sequences, and repetitive elements, presents significant challenges in assembling contiguous reference genomes. Consequently, short-read resequencing studies can encounter mismapping issues, leading to SNPs that deviate from Mendelian expected patterns of heterozygosity and allelic ratio. In this study, we employed the ngsParalog software to identify such deviant SNPs in whole-genome sequencing (WGS) data with low (1.5×) to intermediate (4.8×) coverage for four species: Arctic Char (Salvelinus alpinus), Lake Whitefish (Coregonus clupeaformis), Atlantic Salmon (Salmo salar), and the American Eel (Anguilla rostrata). The analyses revealed that deviant SNPs accounted for 22% to 62% of all SNPs in salmonid datasets and approximately 11% in the American Eel dataset. These deviant SNPs were particularly concentrated within repetitive elements and genomic regions that had recently undergone rediploidization in salmonids. Additionally, narrow peaks of elevated coverage were ubiquitous along all four reference genomes, encompassed most deviant SNPs, and could be partially associated with transposons and tandem repeats. Including these deviant SNPs in genomic analyses led to highly distorted site frequency spectra, underestimated pairwise FST values, and overestimated nucleotide diversity. Considering the widespread occurrence of deviant SNPs arising from a variety of sources, their important impact in estimating population parameters, and the availability of effective tools to identify them, we propose that excluding deviant SNPs from WGS datasets is required to improve genomic inferences for a wide range of taxa and sequencing depths.

Comparative Analysis of Transposable Elements in Strawberry Genomes of Different Ploidy Levels

Transposable elements (TEs) make up a large portion of plant genomes and play a vital role in genome structure, function, and evolution. Cultivated strawberry (Fragaria x ananassa) is one of the most important fruit crops, and its octoploid genome was formed through several rounds of genome duplications from diploid ancestors. Here, we built a pan-genome TE library for the Fragaria genus using ten published strawberry genomes at different ploidy levels, including seven diploids, one tetraploid, and two octoploids, and performed comparative analysis of TE content in these genomes. The TEs comprise 51.83% (F. viridis) to 60.07% (F. nilgerrensis) of the genomes. Long terminal repeat retrotransposons (LTR-RTs) are the predominant TE type in the Fragaria genomes (20.16% to 34.94%), particularly in F. iinumae (34.94%). Estimating TE content and LTR-RT insertion times revealed that species-specific TEs have shaped each strawberry genome. Additionally, the copy number of different LTR-RT families inserted in the last one million years reflects the genetic distance between Fragaria species. Comparing cultivated strawberry subgenomes to extant diploid ancestors showed that F. vesca and F. iinumae are likely the diploid ancestors of the cultivated strawberry, but not F. viridis. These findings provide new insights into the TE variations in the strawberry genomes and their roles in strawberry genome evolution.

Comprehensive detection of structural variation and transposable element differences between wild type laboratory lineages of C. elegans

Genomic structural variations (SVs) and transposable elements (TEs) can be significant contributors to genome evolution, altered gene expression, and risk of genetic diseases. Recent advancements in long-read sequencing have greatly improved the quality of de novo genome assemblies and enhanced the detection of sequence variants at the scale of hundreds or thousands of bases. Comparisons between two diverged wild isolates of Caenorhabditis elegans, the Bristol and Hawaiian strains, have been widely utilized in the analysis of small genetic variations. Genetic drift, including SVs and rearrangements of repeated sequences such as TEs, can occur over time from long-term maintenance of wild type isolates within the laboratory. To comprehensively detect both large and small structural variations as well as TEs due to genetic drift, we generated de novo genome assemblies and annotations for each strain from our lab collection using both long- and short-read sequencing and compared our assemblies and annotations with that of other lab wild type strains. Within our lab assemblies, we annotate over 3.1Mb of sequence divergence between the Bristol and Hawaiian isolates: 337,584 SNPs, 94,503 small insertion-deletions (<50bp), and 4,334 structural variations (>50bp). Further, we define the location and movement of specific DNA TEs between N2 Bristol and CB4856 Hawaiian wild type isolates. Specifically, we find the N2 Bristol genome has 20.6% more TEs from the Tc1/mariner family than the CB4856 Hawaiian genome. Moreover, we identified Zator elements as the most abundant and mobile TE family in the genome. Using specific TE sequences with unique SNPs, we also identify 38 TEs that moved intrachromosomally and 9 TEs that moved interchromosomally between the N2 Bristol and CB4856 Hawaiian genomes. By comparing the de novo genome assembly of our lab collection Bristol isolate to the VC2010 Bristol assembly, we also reveal that lab lineages display over 2 Mb of total variation: 1,162 SNPs, 1,528 indels, and 897 SVs with 95% of the variation due to SVs. Overall, our work demonstrates the unique contribution of SVs and TEs to variation and genetic drift between wild type laboratory strains assumed to be isogenic despite growing evidence of genetic drift and phenotypic variation.

Annotation of transposable elements in the transcriptome of the Neotropical brown stink bug Euschistus heros and its chromosomal distribution

Transposable elements (TEs) are DNA sequences capable of moving within the genome. Their distribution is very dynamic among organisms, and despite advances, there are still gaps in the understanding of the diversity and evolution of TEs in many insect species. In the case of Euschistus heros, considered the main stink bug in the soybean crop in Brazil, little is known about the participation of these elements. Therefore, the objective of the current work was to identify the different groups of transposable elements present in the E. heros transcriptome, evidencing their chromosomal distribution. Through RNA-Seq and de novo assembly, 60,009 transcripts were obtained, which were annotated locally via Blastn against specific databases. Of the 367 transcripts identified as TEs, 202 belong to Class II, with emphasis on the TIR order. Among Class I elements or retrotransposons, most were characterized as LINE. Phylogenetic analyses were performed with the protein domains, evidencing differences between Tc1-mariner sequences, which may be related to possible horizontal transfer events. The transposable elements that stood out in the transcriptome were selected for fluorescent in situ hybridization. DNA transposon probes hAT, Helitron, and Tc1-mariner showed mostly scattered signals, with the presence of some blocks. Retrotransposon probes Copia, Gypsy, Jockey, and RTE showed a more pulverized hybridization pattern, with the presence of small interstitial and/or terminal blocks. Studies like this one, integrating functional genomics and molecular cytogenetic tools, are essential to expanding knowledge about transcriptionally active mobile elements, and their behavior in the chromosomes.

Teleost genomic repeat landscapes in light of diversification rates and ecology

Repetitive DNA make up a considerable fraction of most eukaryotic genomes. In fish, transposable element (TE) activity has coincided with rapid species diversification. Here, we annotated the repetitive content in 100 genome assemblies, covering the major branches of the diverse lineage of teleost fish. We investigated if TE content correlates with family level net diversification rates and found support for a weak negative correlation. Further, we demonstrated that TE proportion correlates with genome size, but not to the proportion of short tandem repeats (STRs), which implies independent evolutionary paths. Marine and freshwater fish had large differences in STR content, with the most extreme propagation detected in the genomes of codfish species and Atlantic herring. Such a high density of STRs is likely to increase the mutational load, which we propose could be counterbalanced by high fecundity as seen in codfishes and herring.

Young SINEs in pig genomes impact gene regulation, genetic diversity, and complex traits

Transposable elements (TEs) are a major source of genetic polymorphisms and play a role in chromatin architecture, gene regulatory networks, and genomic evolution. However, their functional role in pigs and contributions to complex traits are largely unknown. We created a catalog of TEs (n = 3,087,929) in pigs and found that young SINEs were predominantly silenced by histone modifications, DNA methylation, and decreased accessibility. However, some transcripts from active young SINEs showed high tissue-specificity, as confirmed by analyzing 3570 RNA-seq samples. We also detected 211,067 dimorphic SINEs in 374 individuals, including 340 population-specific ones associated with local adaptation. Mapping these dimorphic SINEs to genome-wide associations of 97 complex traits in pigs, we found 54 candidate genes (e.g., ANK2 and VRTN) that might be mediated by TEs. Our findings highlight the important roles of young SINEs and provide a supplement for genotype-to-phenotype associations and modern breeding in pigs.

Transposable elements in early human embryo development and embryo models

Transposable elements (TEs), long discounted as 'selfish genomic elements,' are increasingly appreciated as the drivers of genomic evolution, genome organization, and gene regulation. TEs are particularly important in early embryo development, where advances in stem cell technologies, in tandem with improved computational and next-generation sequencing approaches, have provided an unprecedented opportunity to study the contribution of TEs to early mammalian development. Here, we summarize advances in our understanding of TEs in early human development and expand on how new stem cell-based embryo models can be leveraged to augment this understanding.

Consequences of polyploidy and divergence as revealed by cytogenetic mapping of tandem repeats in African clawed frogs (Xenopus, Pipidae)

Repetitive elements have been identified in several amphibian genomes using whole genome sequencing, but few studies have used cytogenetic mapping to visualize these elements in this vertebrate group. Here, we used fluorescence in situ hybridization and genomic data to map the U1 and U2 small nuclear RNAs and histone H3 in six species of African clawed frog (genus Xenopus), including, from subgenus Silurana, the diploid Xenopus tropicalis and its close allotetraploid relative X. calcaratus and, from subgenus Xenopus, the allotetraploid species X. pygmaeus, X. allofraseri, X. laevis, and X. muelleri. Results allowed us to qualitatively evaluate the relative roles of polyploidization and divergence in the evolution of repetitive elements because our focal species include allotetraploid species derived from two independent polyploidization events - one that is relatively young that gave rise to X. calcaratus and another that is older that gave rise to the other (older) allotetraploids. Our results demonstrated conserved loci number and position of signals in the species from subgenus Silurana; allotetraploid X. calcaratus has twice as many signals as diploid X. tropicalis. However, the content of repeats varied among the other allotetraploid species. We detected almost same number of signals in X. muelleri as in X. calcaratus and same number of signals in X. pygmaeus, X. allofraseri, X. laevis as in the diploid X. tropicalis. Overall, these results are consistent with the proposal that allopolyploidization duplicated these tandem repeats and that variation in their copy number was accumulated over time through reduction and expansion in a subset of the older allopolyploids.

Transposable Elements: Epigenetic Silencing Mechanisms or Modulating Tools for Vertebrate Adaptations? Two Sides of the Same Coin

Transposable elements constitute one of the main components of eukaryotic genomes. In vertebrates, they differ in content, typology, and family diversity and played a crucial role in the evolution of this taxon. However, due to their transposition ability, TEs can be responsible for genome instability, and thus silencing mechanisms were evolved to allow the coexistence between TEs and eukaryotic host-coding genes. Several papers are highlighting in TEs the presence of regulatory elements involved in regulating nearby genes in a tissue-specific fashion. This suggests that TEs are not sequences merely to silence; rather, they can be domesticated for the regulation of host-coding gene expression, permitting species adaptation and resilience as well as ensuring human health. This review presents the main silencing mechanisms acting in vertebrates and the importance of exploiting these mechanisms for TE control to rewire gene expression networks, challenging the general view of TEs as threatening elements.

The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order

Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.

Multiple and diversified transposon lineages contribute to early and recent bivalve genome evolution

BACKGROUND

Transposable elements (TEs) can represent one of the major sources of genomic variation across eukaryotes, providing novel raw materials for species diversification and innovation. While considerable effort has been made to study their evolutionary dynamics across multiple animal clades, molluscs represent a substantially understudied phylum. Here, we take advantage of the recent increase in mollusc genomic resources and adopt an automated TE annotation pipeline combined with a phylogenetic tree-based classification, as well as extensive manual curation efforts, to characterize TE repertories across 27 bivalve genomes with a particular emphasis on DDE/D class II elements, long interspersed nuclear elements (LINEs), and their evolutionary dynamics.

RESULTS

We found class I elements as highly dominant in bivalve genomes, with LINE elements, despite less represented in terms of copy number per genome, being the most common retroposon group covering up to 10% of their genome. We mined 86,488 reverse transcriptases (RVT) containing LINE coming from 12 clades distributed across all known superfamilies and 14,275 class II DDE/D-containing transposons coming from 16 distinct superfamilies. We uncovered a previously underestimated rich and diverse bivalve ancestral transposon complement that could be traced back to their most recent common ancestor that lived ~ 500 Mya. Moreover, we identified multiple instances of lineage-specific emergence and loss of different LINEs and DDE/D lineages with the interesting cases of CR1- Zenon, Proto2, RTE-X, and Academ elements that underwent a bivalve-specific amplification likely associated with their diversification. Finally, we found that this LINE diversity is maintained in extant species by an equally diverse set of long-living and potentially active elements, as suggested by their evolutionary history and transcription profiles in both male and female gonads.

CONCLUSIONS

We found that bivalves host an exceptional diversity of transposons compared to other molluscs. Their LINE complement could mainly follow a "stealth drivers" model of evolution where multiple and diversified families are able to survive and co-exist for a long period of time in the host genome, potentially shaping both recent and early phases of bivalve genome evolution and diversification. Overall, we provide not only the first comparative study of TE evolutionary dynamics in a large but understudied phylum such as Mollusca, but also a reference library for ORF-containing class II DDE/D and LINE elements, which represents an important genomic resource for their identification and characterization in novel genomes.

Chiropterans Are a Hotspot for Horizontal Transfer of DNA Transposons in Mammalia

Horizontal transfer of transposable elements (TEs) is an important mechanism contributing to genetic diversity and innovation. Bats (order Chiroptera) have repeatedly been shown to experience horizontal transfer of TEs at what appears to be a high rate compared with other mammals. We investigated the occurrence of horizontally transferred (HT) DNA transposons involving bats. We found over 200 putative HT elements within bats; 16 transposons were shared across distantly related mammalian clades, and 2 other elements were shared with a fish and two lizard species. Our results indicate that bats are a hotspot for horizontal transfer of DNA transposons. These events broadly coincide with the diversification of several bat clades, supporting the hypothesis that DNA transposon invasions have contributed to genetic diversification of bats.

Induction and application of human naive pluripotency

Over the past few decades, many attempts have been made to capture different states of pluripotency in vitro. Naive and primed pluripotent stem cells, corresponding to the pluripotency states of pre- and post-implantation epiblasts, respectively, have been well characterized in mice and can be interconverted in vitro. Here, we summarize the recently reported strategies to generate human naive pluripotent stem cells in vitro. We discuss their applications in studies of regulatory mechanisms involved in early developmental processes, including identification of molecular features, X chromosome inactivation modeling, transposable elements regulation, metabolic characteristics, and cell fate regulation, as well as potential for extraembryonic differentiation and blastoid construction for embryogenesis modeling. We further discuss the naive pluripotency-related research, including 8C-like cell establishment and disease modeling. We also highlight limitations of current naive pluripotency studies, such as imperfect culture conditions and inadequate responsiveness to differentiation signals.

High resolution genomes of multiple Xiphophorus species provide new insights into microevolution, hybrid incompatibility, and epistasis

Because of diverged adaptative phenotypes, fish species of the genus Xiphophorus have contributed to a wide range of research for a century. Existing Xiphophorus genome assemblies are not at the chromosomal level and are prone to sequence gaps, thus hindering advancement of the intra- and inter-species differences for evolutionary, comparative, and translational biomedical studies. Herein, we assembled high-quality chromosome-level genome assemblies for three distantly related Xiphophorus species, namely, X. maculatus, X. couchianus, and X. hellerii Our overall goal is to precisely assess microevolutionary processes in the clade to ascertain molecular events that led to the divergence of the Xiphophorus species and to progress understanding of genetic incompatibility to disease. In particular, we measured intra- and inter-species divergence and assessed gene expression dysregulation in reciprocal interspecies hybrids among the three species. We found expanded gene families and positively selected genes associated with live bearing, a special mode of reproduction. We also found positively selected gene families are significantly enriched in nonpolymorphic transposable elements, suggesting the dispersal of these nonpolymorphic transposable elements has accompanied the evolution of the genes, possibly by incorporating new regulatory elements in support of the Britten-Davidson hypothesis. We characterized inter-specific polymorphisms, structural variants, and polymorphic transposable element insertions and assessed their association to interspecies hybridization-induced gene expression dysregulation related to specific disease states in humans.

LTR Retroelements and Bird Adaptation to Arid Environments

TEs are known to be among the main drivers in genome evolution, leading to the generation of evolutionary advantages that favor the success of organisms. The aim of this work was to investigate the TE landscape in bird genomes to look for a possible relationship between the amount of specific TE types and environmental changes that characterized the Oligocene era in Australia. Therefore, the mobilome of 29 bird species, belonging to a total of 11 orders, was analyzed. Our results confirmed that LINE retroelements are not predominant in all species of this evolutionary lineage and highlighted an LTR retroelement dominance in species with an Australian-related evolutionary history. The bird LTR retroelement expansion might have happened in response to the Earth’s dramatic climate changes that occurred about 30 Mya, followed by a progressive aridification across most of Australian landmasses. Therefore, in birds, LTR retroelement burst might have represented an evolutionary advantage in the adaptation to arid/drought environments.

Pathways to polar adaptation in fishes revealed by long-read sequencing

Long-read sequencing is driving a new reality for genome science in which highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we utilized long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Ophthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance and vision have evidence of positive selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: haemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the haemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution.