Evolution and Diversity of Transposable Elements in Vertebrate Genomes

Sex-biased fertility degeneration induced by depletion of an auxiliary piRNA-pathway factor Qin in Bombyx mori

The PIWI-interacting RNA (piRNA) pathway is the major defense system for transposable elements (TEs) silencing in animal gonads, maintaining genomic integrity of germ cells and ensuring proper gametogenesis. An the piRNA-pathway factor, Qin, has been reported to participate in piRNA biogenesis in the lepidopteran model insect, Bombyx mori. Nevertheless, the physiological functions of Qin remain to be characterized. Here we demonstrated that Qin plays important roles in silkworm gonad development of both sexes. BmQin was predominantly expressed in gonads. Immunofluorescent staining revealed that BmQin is localized in the cytoplasm of both germ cells and somatic cells in gonads. Depletion of BmQin via CRISPR/Cas9 system induceed complete sterile in males, and partial sterile in females. Notably, mutants displayed severe defects in gonad development and gametogenesis. RNA-seq analysis revealed that the piRNA pathway was dysregulated in mutant gonads. In addition, apoptosis was significantly enhanced in mutant gonads. Our study revealed the physiological functions of BmQin in silkworm fertility and its auxiliary roles in the piRNA pathway in both male and female gonads.

Chromosome-level genome assembly reveals adaptive evolution of the invasive Amazon sailfin catfish (Pterygoplichthys pardalis)

Catfish represents a diverse lineage with variable number of chromosomes and complex relationships with humans. Although certain species pose significant invasive threats to native fish populations, comprehensive genomic investigations into the evolutionary adaptations that contribute to their invasion success are lacking. To address this gap, our study presents a high-quality genome assembly of the Amazon sailfin catfish (Pterygoplichthys pardalis), a member of the armored catfish family, along with a comprehensive comparative genomic analysis. By utilizing conserved genomic regions across different catfish species, we reconstructed the 29 ancestral chromosomes of catfish, including two microchromosomes (28 and 29) that show different fusion and breakage patterns across species. Our analysis shows that the Amazon sailfin catfish genome is notably larger (1.58 Gb) with more than 40,000 coding genes. The genome expansion was linked to early repetitive sequence expansions and recent gene duplications. Several expanded genes are associated with immune functions, including antigen recognition domains like the Ig-v-set domain and the tandem expansion of the CD300 gene family. We also identified specific insertions in CNEs (conserved non-coding elements) near genes involved in cellular processes and neural development. Additionally, rapidly evolving and positively selected genes in the Amazon sailfin catfish genome were found to be associated with collagen formation. Moreover, we identified multiple positively selected codons in hoxb9, which may lead to functional alterations. These findings provide insights into molecular adaptations in an invasive catfish that may underlie its widespread invasion success.

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.

Optimized toolkit for the manipulation of immortalized axolotl fibroblasts

The axolotl salamander model has broad utility for regeneration studies, but this model is limited by a lack of efficient cell-culture-based tools. The Axolotl Limb-1 (AL-1) fibroblast line, the only available immortalized axolotl cell line, was first published over 20 years ago, but many established molecular biology techniques, such as lipofectamine transfection, CRISPR-Cas9 mutagenesis, and antibiotic selection, work poorly or remain untested in AL-1 cells. Innovating technologies to manipulate AL-1 cells in culture and study their behavior following transplantation into the axolotl will complement in-vivo studies, decrease the number of animals used, and enable the faster, more streamlined investigation of regenerative biology questions. Here, we establish transfection, mutagenesis, antibiotic selection, and in-vivo transplantation techniques in axolotl AL-1 cells. These techniques will enable efficient culture with AL-1 cells and guide future tool development for the culture and manipulation of other salamander cell lines.

Multi-omics analysis reveals critical cis-regulatory roles of transposable elements in livestock genomes

Transposable elements (TEs) are important sources of genetic and regulatory variation, yet their functional roles in domesticated animals remain insufficiently explored. To address this gap, we comprehensively annotated TE types, ages, and distributions in the genomes of pig (Sus scrofa), cattle (Bos taurus), and chicken (Gallus gallus). Our analysis revealed species-specific patterns in TE abundance, amplification, and activity in modern genomes. By integrating transcriptomic and epigenomic data, we explored the impact of specific TE types on cis-regulatory elements (CREs) and constructed a TE expression atlas across five tissues in all three species. Our findings underscored the critical roles of tissue-specific TE expression and chromatin accessibility in regulating tissue-specific biological processes. Most notably, we developed a computational framework to uncover TE-mediated gene regulatory networks (TE-GRNs). Our findings provide valuable insights into the regulatory functions of TEs in livestock and offer a robust approach for studying TE-GRNs in diverse biological contexts.

Syntenic lncRNA locus exhibits DNA regulatory functions with sequence evolution

Syntenic long non-coding RNAs (lncRNAs) often show limited sequence conservation across species, prompting concern in the field. This study delves into functional signatures of syntenic lncRNAs between humans and zebrafish. Syntenic lncRNAs are highly expressed in zebrafish, with ∼90 % located near protein-coding genes, either in sense or antisense orientation. During early zebrafish development and in human embryonic stem cells (H1-hESC), syntenic lncRNA loci are enriched with cis-regulatory repressor signatures, influencing the expression of development-associated genes. In later zebrafish developmental stages and specific human cell lines, these syntenic lncRNA loci function as enhancers or transcription start sites (TSS) for protein-coding genes. Analysis of transposable elements (TEs) in syntenic lncRNA sequences revealed intriguing patterns: human lncRNAs are enriched in simple repeat elements, while their zebrafish counterparts show enrichment in LTR elements. This sequence evolution likely arises from post-rearrangement mutations that enhance DNA elements or cis-regulatory functions. It may also contribute to vertebrate innovation by creating novel transcription factor binding sites within the locus. This study highlights the conserved functionality of syntenic lncRNA loci through DNA elements, emphasizing their conserved roles across species despite sequence divergence.

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.

A chromosome-level, haplotype-resolved genome assembly and annotation for the Eurasian minnow (Leuciscidae: Phoxinus phoxinus) provide evidence of haplotype diversity

BACKGROUND

In this study, we present an in-depth analysis of the Eurasian minnow (Phoxinus phoxinus) genome, highlighting its genetic diversity, structural variations, and evolutionary adaptations. We generated an annotated haplotype-phased, chromosome-level genome assembly (2n = 50) by integrating high-fidelity (HiFi) long reads and chromosome conformation capture data (Hi-C).

RESULTS

We achieved a haploid size of 940 megabase pairs (Mbp) for haplome 1 and 929 Mbp for haplome 2 with high scaffold N50 values of 36.4 Mb and 36.6 Mb and BUSCO scores of 96.9% and 97.2%, respectively, indicating a highly complete genome assembly. We detected notable heterozygosity (1.43%) and a high repeat content (approximately 54%), primarily consisting of DNA transposons, which contribute to genome rearrangements and variations. We found substantial structural variations within the genome, including insertions, deletions, inversions, and translocations. These variations affect genes enriched in functions such as dephosphorylation, developmental pigmentation, phagocytosis, immunity, and stress response. In the annotation of protein-coding genes, 30,980 messenger RNAs and 23,497 protein-coding genes were identified with a high completeness score, which further underpins the high contiguity of our genome assemblies. We performed a gene family evolution analysis by comparing our proteome to 10 other teleost species, which identified immune system gene families that prioritize histone-based disease prevention over NB-LRR-related-based immune responses. Additionally, demographic analysis indicates historical fluctuations in the effective population size of P. phoxinus, likely correlating with past climatic changes.

CONCLUSIONS

This annotated, phased reference genome provides a crucial resource for resolving the taxonomic complexity within the genus Phoxinus and highlights the importance of haplotype-phased assemblies in understanding haplotype diversity in species characterized by high heterozygosity.

MCHelper automatically curates transposable element libraries across eukaryotic species

The number of species with high-quality genome sequences continues to increase, in part due to the scaling up of multiple large-scale biodiversity sequencing projects. While the need to annotate genic sequences in these genomes is widely acknowledged, the parallel need to annotate transposable element (TE) sequences that have been shown to alter genome architecture, rewire gene regulatory networks, and contribute to the evolution of host traits is becoming ever more evident. However, accurate genome-wide annotation of TE sequences is still technically challenging. Several de novo TE identification tools are now available, but manual curation of the libraries produced by these tools is needed to generate high-quality genome annotations. Manual curation is time-consuming, and thus impractical for large-scale genomic studies, and lacks reproducibility. In this work, we present the Manual Curator Helper tool MCHelper, which automates the TE library curation process. By leveraging MCHelper's fully automated mode with the outputs from three de novo TE identification tools, RepeatModeler2, EDTA, and REPET, in the fruit fly, rice, hooded crow, zebrafish, maize, and human, we show a substantial improvement in the quality of the TE libraries and genome annotations. MCHelper libraries are less redundant, with up to 65% reduction in the number of consensus sequences, have up to 11.4% fewer false positive sequences, and up to ∼48% fewer "unclassified/unknown" TE consensus sequences. Genome-wide TE annotations are also improved, including larger unfragmented insertions. Moreover, MCHelper is an easy-to-install and easy-to-use tool.

HROB Is Implicated in DNA Replication

DNA replication represents a series of precisely regulated events performed by a complex protein machinery that guarantees accurate duplication of the genetic information. Since DNA replication is permanently faced by a variety of exogenous and endogenous stressors, DNA damage response, repair and replication must be closely coordinated to maintain genomic integrity. HROB has been identified recently as a binding partner and activator of the Mcm8/9 helicase involved in DNA interstrand crosslink (ICL) repair. We identified HROB independently as a nuclear protein whose expression is co-regulated with various DNA replication factors. Accordingly, the HROB protein level showed a maximum in S phase and a downregulation in quiescence. Structural prediction and homology searches revealed that HROB is a largely intrinsically disordered protein bearing a helix-rich region and a canonical oligonucleotide/oligosaccharide-binding-fold motif that originated early in eukaryotic evolution. Employing a flow cytometry Förster resonance energy transfer (FRET) assay, we detected associations between HROB and proteins of the DNA replication machinery. Moreover, ectopic expression of HROB protein led to an almost complete shutdown of DNA replication. The available data imply a function for HROB during DNA replication across barriers such as ICLs.

Ecological diversification of sea catfishes is accompanied by genome-wide signatures of positive selection

Habitat transitions have shaped the evolutionary trajectory of many clades. Sea catfishes (Ariidae) have repeatedly undergone ecological transitions, including colonizing freshwaters from marine environments, leading to an adaptive radiation in Australia and New Guinea alongside non-radiating freshwater lineages elsewhere. Here, we generate and analyze one long-read reference genome and 66 short-read whole genome assemblies, in conjunction with genomic data for 54 additional species. We investigate how three major ecological transitions have shaped genomic variation among ariids over their ~ 50 million-year evolutionary history. Our results show that relatively younger freshwater lineages exhibit a higher incidence of positive selection than their more ancient marine counterparts. They also display a larger disparity in body shapes, a trend that correlates with a heightened occurrence of positive selection on genes associated with body size and elongation. Although positive selection in the Australia and New Guinea radiation does not stand out compared to non-radiating lineages overall, selection across the prolactin gene family during the marine-to-freshwater transition suggests that strong osmoregulatory adaptations may have facilitated their colonization and radiation. Our findings underscore the significant role of selection in shaping the genome and organismal traits in response to habitat shifts across macroevolutionary scales.

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.

Delineating the W Sex Chromosome in the Clam Shrimp, Eulimnadia texana

INTRODUCTION

Sex chromosomes have evolved independently across various lineages, often showing convergent degradation of the sex-limited chromosome. While extensively studied in model organisms with ancient sex chromosomal systems, the evolution of early-stage sex chromosomes remains poorly understood. Eulimnadia texana, a freshwater crustacean with a unique androdioecious breeding system (ZZ, ZW, and viable WW genotypes), provides a rare opportunity to study early sex chromosome evolution. This study examines E. texana's W chromosome for evidence of a small localized non-recombining region, characterized by a transposable element (TE) "hotspot," low gene density, and low GC content.

METHODS

Sex-linked markers were mapped onto the W chromosome (scaffold 1). TEs in the WW genome were identified using RepeatModeler and RepeatMasker. Statistical analyses compared TE distribution between the genome and scaffold 1, which was then divided into 20 equal-sized "bins" for finer-scale statistical analyses. Gene density and GC content were analyzed across these bins.

RESULTS

While no significant TE accumulation was found across the entire W chromosome compared to the remaining genome, a specific region (6.6-8.8 Mb, fourth bin) showed significantly higher TE accumulation. This region also exhibited low gene density and low GC content, indicative of reduced recombination.

CONCLUSION

Our findings suggest that E. texana's W chromosome contains a smaller region of crossover suppression, supporting the hypothesis that it is a proto-sex chromosome in early evolutionary development. This study provides valuable insights into early sex chromosome evolution and establishes E. texana as an ideal model for further investigation of evolutionary processes driving proto-sex chromosome differentiation.

IRescue: uncertainty-aware quantification of transposable elements expression at single cell level

Transposable elements (TEs) are mobile DNA repeats known to shape the evolution of eukaryotic genomes. In complex organisms, they exhibit tissue-specific transcription. However, understanding their role in cellular diversity across most tissues remains a challenge, when employing single-cell RNA sequencing (scRNA-seq), due to their widespread presence and genetic similarity. To address this, we present IRescue (Interspersed Repeats single-cell quantifier), a software capable of estimating the expression of TE subfamilies at the single-cell level. IRescue incorporates a unique UMI deduplication algorithm to rectify sequencing errors and employs an Expectation-Maximization procedure to effectively redistribute the counts of multi-mapping reads. Our study showcases the precision of IRescue through analysis of both simulated and real single cell and nuclei RNA-seq data from human colorectal cancer, brain, skin aging, and PBMCs during SARS-CoV-2 infection and recovery. By linking the expression patterns of TE signatures to specific conditions and biological contexts, we unveil insights into their potential roles in cellular heterogeneity and disease progression.

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.

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.

Discovery of numerous novel Helitron-like elements in eukaryote genomes using HELIANO

Helitron-like elements (HLEs) are widespread eukaryotic DNA transposons employing a rolling-circle transposition mechanism. Despite their prevalence in fungi, animals, and plant genomes, identifying Helitrons remains a formidable challenge. We introduce HELIANO, a software for annotating and classifying autonomous and non-autonomous HLE sequences from whole genomes. HELIANO overcomes several limitations of existing tools in speed and accuracy, demonstrated through benchmarking and its application to the complex genomes of frogs (Xenopus tropicalis and Xenopus laevis) and rice (Oryza sativa), where it uncovered numerous previously unidentified HLEs. In an extensive analysis of 404 eukaryote genomes, we found HLEs widely distributed across phyla, with exceptions in specific taxa. HELIANO's application led to the discovery of numerous new HLEs in land plants and identified 20 protein domains captured by certain autonomous HLE families. A comprehensive phylogenetic analysis further classified HLEs into two primary clades, HLE1 and HLE2, and revealed nine subgroups, some of which are enriched within specific taxa. The future use of HELIANO promises to improve the global analysis of HLEs across genomes, significantly advancing our understanding of this fascinating transposon superfamily.

A reference genome for the Harpy Eagle reveals steady demographic decline and chromosomal rearrangements in the origin of Accipitriformes

The Harpy Eagle (Harpia harpyja) is an iconic species that inhabits forested landscapes in Neotropical regions, with decreasing population trends mainly due to habitat loss, and currently classified as vulnerable. Here, we report on a chromosome-scale genome assembly for a female individual combining long reads, optical mapping, and chromatin conformation capture reads. The final assembly spans 1.35 Gb, with N50scaffold equal to 58.1 Mb and BUSCO completeness of 99.7%. We built the first extensive transposable element (TE) library for the Accipitridae to date and identified 7,228 intact TEs. We found a burst of an unknown TE ~ 13-22 million years ago (MYA), coincident with the split of the Harpy Eagle from other Harpiinae eagles. We also report a burst of solo-LTRs and CR1 retrotransposons ~ 31-33 MYA, overlapping with the split of the ancestor to all Harpiinae from other Accipitridae subfamilies. Comparative genomics with other Accipitridae, the closely related Cathartidae and Galloanserae revealed major chromosome-level rearrangements at the basal Accipitriformes genome, in contrast to a conserved ancient genome architecture for the latter two groups. A historical demography reconstruction showed a rapid decline in effective population size over the last 20,000 years. This reference genome serves as a crucial resource for future conservation efforts towards the Harpy Eagle.

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.

Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.

The Families of Non-LTR Transposable Elements within Neritimorpha and Other Gastropoda

Repeated sequences, especially transposable elements (TEs), are known to be abundant in some members of the important invertebrate class Gastropoda. TEs that do not have long terminal repeated sequences (non-LTR TEs) are frequently the most abundant type but have not been well characterised in any gastropod. Despite this, sequences in draft gastropod genomes are often described as non-LTR TEs, but without identification to family type. This study was conducted to characterise non-LTR TEs in neritimorph snails, using genomic skimming surveys of three species and the recently published draft genome of Theodoxus fluviatilis. Multiple families of non-LTR TEs from the I, Jockey, L1, R2 and RTE superfamilies were found, although there were notably few representatives of the first of these, which is nevertheless abundant in other Gastropoda. Phylogenetic analyses of amino acid sequences of the reverse transcriptase domain from the elements ORF2 regions found considerable interspersion of representatives of the four neritimorph taxa within non-LTR families and sub-families. In contrast, phylogenetic analyses of sequences from the elements' ORF1 region resolved the representatives from individual species as monophyletic. However, using either region, members of the two species of the Neritidae were closely related, suggesting their potential for investigation of phyletic evolution at the family level.

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.

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.

Gigantic animal cells suggest organellar scaling mechanisms across a 50-fold range in cell volume

The largest cells are orders of magnitude bigger than the smallest cells. Organelle content scales to maintain cell function, with different organelles increasing in volume, length, or number as cells increase in size. Scaling may also reflect functional demands placed on organelles by increased cell size. Amphibians exhibit exceptional diversity in cell size. Using transmission electron microscopy, we analyzed 3 species whose enterocyte cell volumes range from 228 to 10,593 μm3. We show that nuclear volume increases by an increase in radius while mitochondrial volume increases by an increase in total network length; the endoplasmic reticulum and Golgi apparatus, with their complex shapes, are intermediate. Notably, all 4 organelle types increase in total volume proportional to cell volume, despite variation in functional (i.e., metabolic, transport) demands. This pattern suggests that organellar building blocks are incorporated into more or larger organelles following the same rules across species that vary ~50-fold in cell sizes, consistent with a "limited precursor" model for organellar scaling that, in turn, assumes equivalent cytoplasmic concentrations of organellar building block proteins. Taken together, our results lead us to hypothesize that salamanders have evolved increased biosynthetic capacity to maintain functional protein concentrations despite huge cell volumes.

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.

Genome-scale phylogeny and comparative genomics of the fungal order Sordariales

The order Sordariales is taxonomically diverse, and harbours many species with different lifestyles and large economic importance. Despite its importance, a robust genome-scale phylogeny, and associated comparative genomic analysis of the order is lacking. In this study, we examined whole-genome data from 99 Sordariales, including 52 newly sequenced genomes, and seven outgroup taxa. We inferred a comprehensive phylogeny that resolved several contentious relationships amongst families in the order, and cleared-up intrafamily relationships within the Podosporaceae. Extensive comparative genomics showed that genomes from the three largest families in the dataset (Chaetomiaceae, Podosporaceae and Sordariaceae) differ greatly in GC content, genome size, gene number, repeat percentage, evolutionary rate, and genome content affected by repeat-induced point mutations (RIP). All genomic traits showed phylogenetic signal, and ancestral state reconstruction revealed that the variation of the properties stems primarily from within-family evolution. Together, the results provide a thorough framework for understanding genome evolution in this important group of fungi.

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.

Take a walk on the KRAB side

Canonical Krüppel-associated box (KRAB)-containing zinc finger proteins (KZFPs) act as major repressors of transposable elements (TEs) via the KRAB-mediated recruitment of the heterochromatin scaffold KRAB-associated protein (KAP)1. KZFP genes emerged some 420 million years ago in the last common ancestor of coelacanth, lungfish, and tetrapods, and dramatically expanded to give rise to lineage-specific repertoires in contemporary species paralleling their TE load and turnover. However, the KRAB domain displays sequence and function variations that reveal repeated diversions from a linear TE-KZFP trajectory. This Review summarizes current knowledge on the evolution of KZFPs and discusses how ancestral noncanonical KZFPs endowed with variant KRAB, SCAN or DUF3669 domains have been utilized to achieve KAP1-independent functions.

ChimeraTE: a pipeline to detect chimeric transcripts derived from genes and transposable elements

Transposable elements (TEs) produce structural variants and are considered an important source of genetic diversity. Notably, TE-gene fusion transcripts, i.e. chimeric transcripts, have been associated with adaptation in several species. However, the identification of these chimeras remains hindered due to the lack of detection tools at a transcriptome-wide scale, and to the reliance on a reference genome, even though different individuals/cells/strains have different TE insertions. Therefore, we developed ChimeraTE, a pipeline that uses paired-end RNA-seq reads to identify chimeric transcripts through two different modes. Mode 1 is the reference-guided approach that employs canonical genome alignment, and Mode 2 identifies chimeras derived from fixed or insertionally polymorphic TEs without any reference genome. We have validated both modes using RNA-seq data from four Drosophila melanogaster wild-type strains. We found ∼1.12% of all genes generating chimeric transcripts, most of them from TE-exonized sequences. Approximately ∼23% of all detected chimeras were absent from the reference genome, indicating that TEs belonging to chimeric transcripts may be recent, polymorphic insertions. ChimeraTE is the first pipeline able to automatically uncover chimeric transcripts without a reference genome, consisting of two running Modes that can be used as a tool to investigate the contribution of TEs to transcriptome plasticity.

Analyses of 600+ insect genomes reveal repetitive element dynamics and highlight biodiversity-scale repeat annotation challenges

Repetitive elements (REs) are integral to the composition, structure, and function of eukaryotic genomes, yet remain understudied in most taxonomic groups. We investigated REs across 601 insect species and report wide variation in RE dynamics across groups. Analysis of associations between REs and protein-coding genes revealed dynamic evolution at the interface between REs and coding regions across insects, including notably elevated RE-gene associations in lineages with abundant long interspersed nuclear elements (LINEs). We leveraged this large, empirical data set to quantify impacts of long-read technology on RE detection and investigate fundamental challenges to RE annotation in diverse groups. In long-read assemblies, we detected ∼36% more REs than short-read assemblies, with long terminal repeats (LTRs) showing 162% increased detection, whereas DNA transposons and LINEs showed less respective technology-related bias. In most insect lineages, 25%-85% of repetitive sequences were "unclassified" following automated annotation, compared with only ∼13% in Drosophila species. Although the diversity of available insect genomes has rapidly expanded, we show the rate of community contributions to RE databases has not kept pace, preventing efficient annotation and high-resolution study of REs in most groups. We highlight the tremendous opportunity and need for the biodiversity genomics field to embrace REs and suggest collective steps for making progress toward this goal.

The blackcap (Sylvia atricapilla) genome reveals a recent accumulation of LTR retrotransposons

Transposable elements (TEs) are mobile genetic elements that can move around the genome, and as such are a source of genomic variability. Based on their characteristics we can annotate TEs within the host genome and classify them into specific TE types and families. The increasing number of available high-quality genome references in recent years provides an excellent resource that will enhance the understanding of the role of recently active TEs on genetic variation and phenotypic evolution. Here we showcase the use of a high-quality TE annotation to understand the distinct effect of recent and ancient TE insertions on the evolution of genomic variation, within our study species the Eurasian blackcap (Sylvia atricapilla). We investigate how these distinct TE categories are distributed along the genome and evaluate how their coverage across the genome is correlated with four genomic features: recombination rate, gene coverage, CpG island coverage and GC content. We found within the recent TE insertions an accumulation of LTRs previously not seen in birds. While the coverage of recent TE insertions was negatively correlated with both GC content and recombination rate, the correlation with recombination rate disappeared and turned positive for GC content when considering ancient TE insertions.

Assessing the Expression of Long INterspersed Elements (LINEs) via Long-Read Sequencing in Diverse Human Tissues and Cell Lines

Transposable elements, such as Long INterspersed Elements (LINEs), are DNA sequences that can replicate within genomes. LINEs replicate using an RNA intermediate followed by reverse transcription and are typically a few kilobases in length. LINE activity creates genomic structural variants in human populations and leads to somatic alterations in cancer genomes. Long-read RNA sequencing technologies, including Oxford Nanopore and PacBio, can directly sequence relatively long transcripts, thus providing the opportunity to examine full-length LINE transcripts. This study focuses on the development of a new bioinformatics pipeline for the identification and quantification of active, full-length LINE transcripts in diverse human tissues and cell lines. In our pipeline, we utilized RepeatMasker to identify LINE-1 (L1) transcripts from long-read transcriptome data and incorporated several criteria, such as transcript start position, divergence, and length, to remove likely false positives. Comparisons between cancerous and normal cell lines, as well as human tissue samples, revealed elevated expression levels of young LINEs in cancer, particularly at intact L1 loci. By employing bioinformatics methodologies on long-read transcriptome data, this study demonstrates the landscape of L1 expression in tissues and cell lines.

Sequencing and Characterization of αs2-Casein Gene (CSN1S2) in the Old-World Camels Have Proven Genetic Variations Useful for the Understanding of Species Diversification

The CSN1S2 gene encodes αs2-casein, the third most abundant protein in camel milk. Despite its importance in foals, human nutrition, and dairy processing, the CSN1S2 gene in camels has received little attention. This study presents the first complete characterization of the CSN1S2 gene sequence in Old-World camels (Camelus bactrianus and Camelus dromedarius). Additionally, the gene promoter, consisting of 752 bp upstream of exon 1, was analyzed. The entire gene comprises 17 exons, ranging in length from 24 bp (exons 4, 8, 11, and 13) to 280 bp (exon 17). Interesting was the identification of the exon 12 in both species. The promoter analysis revealed 24 putative binding sites in the Bactrian camel and 22 in dromedary camel. Most of these sites were typical elements associated with milk protein, such as C/EBP-α, C/EBP-β, Oct-1, and AP1. The SNP discovery showed relatively high genetic diversity compared to other camel casein genes (CSN1S1, CSN2, and CSN3), with a total of 34 polymorphic sites across the two species. Particularly noteworthy is the transition g.311G>A in the CSN1S2 promoter, creating a new putative consensus binding site for a C/EBP-β in the Bactrian camel. At the exon level, two novel variants were found. One was detected in exon 6 of the Bactrian camel (g.3639C>G), resulting in an amino acid replacement, p.36Ile>Met. The second variant was found in noncoding exon 17 of dromedary CSN1S2 (g.1511G>T). Although this mutation occurs in the 3'-UnTranslated Region, it represents the first example of exonic polymorphism in the CSN1S2 for this species. This SNP also affects the binding sites of different microRNAs, including the seed sequence of the miRNA 4662a-3p, highlighting its role as a regulatory factor for CSN1S2 gene. A PCR-RFLP was set up for genotyping a dromedary Tunisian population (n = 157), and the minor allele frequency was found to be 0.27 for the G allele, indicating a potential yield improvement margin. The interspersed elements (INEs) analysis revealed 10 INEs covering 7.34% and 8.14% of the CSN1S2 sequence in the Bactrian and dromedary camels, respectively. Furthermore, six elements (A, B, F, H, I, and L) are shared among cattle and camels and are partially found in other ruminants, suggesting a common ancestral origin of these retrotransposons. Conversely, elements C, D, E, and G are specific to camels.

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.

Abandoning the Isochore Theory Can Help Explain Genome Compositional Organization in Fish

The organization of the genome nucleotide (AT/GC) composition in vertebrates remains poorly understood despite the numerous genome assemblies available. Particularly, the origin of the AT/GC heterogeneity in amniotes, in comparison to the homogeneity in anamniotes, is controversial. Recently, several exceptions to this dichotomy were confirmed in an ancient fish lineage with mammalian AT/GC heterogeneity. Hence, our current knowledge necessitates a reevaluation considering this fact and utilizing newly available data and tools. We analyzed fish genomes in silico with as low user input as possible to compare previous approaches to assessing genome composition. Our results revealed a disparity between previously used plots of GC% and histograms representing the authentic distribution of GC% values in genomes. Previous plots heavily reduced the range of GC% values in fish to comply with the alleged AT/GC homogeneity and AT-richness of their genomes. We illustrate how the selected sequence size influences the clustering of GC% values. Previous approaches that disregarded chromosome and genome sizes, which are about three times smaller in fish than in mammals, distorted their results and contributed to the persisting confusion about fish genome composition. Chromosome size and their transposons may drive the AT/GC heterogeneity apparent on mammalian chromosomes, whereas far less in fishes.

Abundance and Diversification of Repetitive Elements in Decapoda Genomes

Repetitive elements are a major component of DNA sequences due to their ability to propagate through the genome. Characterization of Metazoan repetitive profiles is improving; however, current pipelines fail to identify a significant proportion of divergent repeats in non-model organisms. The Decapoda order, for which repeat content analyses are largely lacking, is characterized by extremely variable genome sizes that suggest an important presence of repetitive elements. Here, we developed a new standardized pipeline to annotate repetitive elements in non-model organisms, which we applied to twenty Decapoda and six other Crustacea genomes. Using this new tool, we identified 10% more repetitive elements than standard pipelines. Repetitive elements were more abundant in Decapoda species than in other Crustacea, with a very large number of highly repeated satellite DNA families. Moreover, we demonstrated a high correlation between assembly size and transposable elements and different repeat dynamics between Dendrobranchiata and Reptantia. The patterns of repetitive elements largely reflect the phylogenetic relationships of Decapoda and the distinct evolutionary trajectories within Crustacea. In summary, our results highlight the impact of repetitive elements on genome evolution in Decapoda and the value of our novel annotation pipeline, which will provide a baseline for future comparative analyses.

maT and mosquito transposons in cnidarians: evolutionary history and intraspecific differences

Transposable elements exert a significant effect on the size and structure of eukaryotic genomes. Tc1/mariner superfamily elements represent the widely distributed and highly variable group of DNA transposons. Tc1/mariner elements include TLE/DD34-38E, MLE/DD34D, maT/DD37D, Visitor/DD41D, Guest/DD39D, mosquito/DD37E, and L18/DD37E families, all of which are well or less scarcely studied. However, more detailed research into the patterns of prevalence and diversity of Tc1/mariner transposons enables one to better understand the coevolution of the TEs and the eukaryotic genomes. We performed a detailed analysis of the maT/DD37D family in Cnidaria. The study of 77 genomic assemblies demonstrated that maT transposons are found in a limited number of cnidarian species belonging to classes Cubozoa (1 species), Hydrozoa (3 species) и Scyphozoa (5 species) only. The identified TEs were classified into 5 clades, with the representatives from Pelagiidae (class Scyphozoa) forming a separate clade of maT transposons, which has never been described previously. The potentially functional copies of maT transposons were identified in the hydrae. The phylogenetic analysis and the studies of distribution among the taxons and the evolutionary dynamics of the elements suggest that maT transposons of the cnidarians are the descendants of several independent invasion events occurring at different periods of time. We also established that the TEs of mosquito/DD37E family are found in Hydridae (class Hydrozoa) only. A comparison of maT and mosquito prevalence in two genomic assemblies of Hydra viridissima revealed obvious differences, thus demonstrating that each individual organism might carry a unique mobilome pattern. The results of the presented research make us better understand the diversity and evolution of Tc1/mariner transposons and their effect on the eukaryotic genomes.

Comparative genomics reveals insights into anuran genome size evolution

BACKGROUND

Amphibians, particularly anurans, display an enormous variation in genome size. Due to the unavailability of whole genome datasets in the past, the genomic elements and evolutionary causes of anuran genome size variation are poorly understood. To address this, we analyzed whole-genome sequences of 14 anuran species ranging in size from 1.1 to 6.8 Gb. By annotating multiple genomic elements, we investigated the genomic correlates of anuran genome size variation and further examined whether the genome size relates to habitat types.

RESULTS

Our results showed that intron expansions or contraction and Transposable Elements (TEs) diversity do not contribute significantly to genome size variation. However, the recent accumulation of transposable elements (TEs) and the lack of deletion of ancient TEs primarily accounted for the evolution of anuran genome sizes. Our study showed that the abundance and density of simple repeat sequences positively correlate with genome size. Ancestral state reconstruction revealed that genome size exhibits a taxon-specific pattern of evolution, with families Bufonidae and Pipidae experiencing extreme genome expansion and contraction events, respectively. Our result showed no relationship between genome size and habitat types, although large genome-sized species are predominantly found in humid habitats.

CONCLUSIONS

Overall, our study identified the genomic element and their evolutionary dynamics accounting for anuran genome size variation, thus paving a path to a greater understanding of the size evolution of the genome in amphibians.

Genome-wide characterization of long terminal repeat retrotransposons provides insights into trait evolution of four cucurbit species

Cucurbits are a diverse plant family that includes economically important crops, such as cucumber, watermelon, melon, and pumpkin. Knowledge of the roles that long terminal repeat retrotransposons (LTR-RTs) have played in diversification of cucurbit species is limited; to add to understanding of the roles of LTR-RTs, we assessed their distributions in four cucurbit species. We identified 381, 578, 1086, and 623 intact LTR-RTs in cucumber (Cucumis sativus L. var. sativus cv. Chinese Long), watermelon (Citrullus lanatus subsp. vulgaris cv. 97103), melon (Cucumis melo cv. DHL92), and Cucurbita (Cucurbita moschata var. Rifu), respectively. Among these LTR-RTs, the Ale clade of the Copia superfamily was the most abundant in all the four cucurbit species. Insertion time and copy number analysis revealed that an LTR-RT burst occurred approximately 2 million years ago in cucumber, watermelon, melon, and Cucurbita, and may have contributed to their genome size variation. Phylogenetic and nucleotide polymorphism analyses suggested that most LTR-RTs were formed after species diversification. Analysis of gene insertions by LTR-RTs revealed that the most frequent insertions were of Ale and Tekay and that genes related to dietary fiber synthesis were the most commonly affected by LTR-RTs in Cucurbita. These results increase our understanding of LTR-RTs and their roles in genome evolution and trait characterization in cucurbits.

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.

The genome assembly and annotation of the many-banded krait, Bungarus multicinctus

Snakes are a vital component of wildlife resources and are widely distributed across the globe. The many-banded krait Bungarus multicinctus is a highly venomous snake found across Southern Asia and central and southern China. Snakes are an ancient reptile group, and their genomes can provide important clues for understanding the evolutionary history of reptiles. Additionally, genomic resources play a crucial role in comprehending the evolution of all species. However, snake genomic resources are still scarce. Here, we present a highly contiguous genome of B. multicinctus with a size of 1.51 Gb. The genome contains a repeat content of 40.15%, with a total length exceeding 620 Mb. Additionally, we annotated a total of 24,869 functional genes. This research is of great significance for comprehending the evolution of B. multicinctus and provides genomic information on the genes involved in venom gland functions.

Genomics of cold adaptations in the Antarctic notothenioid fish radiation

Numerous novel adaptations characterise the radiation of notothenioids, the dominant fish group in the freezing seas of the Southern Ocean. To improve understanding of the evolution of this iconic fish group, here we generate and analyse new genome assemblies for 24 species covering all major subgroups of the radiation, including five long-read assemblies. We present a new estimate for the onset of the radiation at 10.7 million years ago, based on a time-calibrated phylogeny derived from genome-wide sequence data. We identify a two-fold variation in genome size, driven by expansion of multiple transposable element families, and use the long-read data to reconstruct two evolutionarily important, highly repetitive gene family loci. First, we present the most complete reconstruction to date of the antifreeze glycoprotein gene family, whose emergence enabled survival in sub-zero temperatures, showing the expansion of the antifreeze gene locus from the ancestral to the derived state. Second, we trace the loss of haemoglobin genes in icefishes, the only vertebrates lacking functional haemoglobins, through complete reconstruction of the two haemoglobin gene clusters across notothenioid families. Both the haemoglobin and antifreeze genomic loci are characterised by multiple transposon expansions that may have driven the evolutionary history of these genes.

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.

Genomic insights into metabolic flux in hummingbirds

Hummingbirds are very well adapted to sustain efficient and rapid metabolic shifts. They oxidize ingested nectar to directly fuel flight when foraging but have to switch to oxidizing stored lipids derived from ingested sugars during the night or long-distance migratory flights. Understanding how this organism moderates energy turnover is hampered by a lack of information regarding how relevant enzymes differ in sequence, expression, and regulation. To explore these questions, we generated a chromosome-scale genome assembly of the ruby-throated hummingbird (A. colubris) using a combination of long- and short-read sequencing, scaffolding it using existing assemblies. We then used hybrid long- and short-read RNA sequencing of liver and muscle tissue in fasted and fed metabolic states for a comprehensive transcriptome assembly and annotation. Our genomic and transcriptomic data found positive selection of key metabolic genes in nectivorous avian species and deletion of critical genes (SLC2A4, GCK) involved in glucostasis in other vertebrates. We found expression of a fructose-specific version of SLC2A5 putatively in place of insulin-sensitive SLC2A5, with predicted protein models suggesting affinity for both fructose and glucose. Alternative isoforms may even act to sequester fructose to preclude limitations from transport in metabolism. Finally, we identified differentially expressed genes from fasted and fed hummingbirds, suggesting key pathways for the rapid metabolic switch hummingbirds undergo.

Spindle-E is essential for gametogenesis in the silkworm, Bombyx mori

As a defense mechanism against transposable elements, the PIWI-interacting RNA (piRNA) pathway maintains genomic integrity and ensures proper gametogenesis in gonads. Numerous factors are orchestrated to ensure normal operation of the piRNA pathway. Spindle-E (Spn-E) gene was one of the first genes shown to participate in the piRNA pathway. In this study, we performed functional analysis of Spn-E in the model lepidopteran insect, Bombyx mori. Unlike the germline-specific expression pattern observed in Drosophila and mouse, BmSpn-E was ubiquitously expressed in all tissues tested, and it was highly expressed in gonads. Immunofluorescent staining showed that BmSpn-E was localized in both germ cells and somatic cells in ovary and was expressed in spermatocytes in testis. We used a binary transgenic CRISPR/Cas9 system to construct BmSpn-E mutants. Loss of BmSpn-E expression caused derepression of transposons in gonads. We also found that mutant gonads were much smaller than wild-type gonads and that the number of germ cells was considerably lower in mutant gonads. Quantitative real-time PCR analysis and TUNEL staining revealed that apoptosis was greatly enhanced in mutant gonads. Further, we found that the BmSpn-E mutation impacted gonadal development and gametogenesis at the early larval stage. In summary, our data provided the first evidence that BmSpn-E plays vital roles in gonadal development and gametogenesis in B. mori.

Whole-genome identification of transposable elements reveals the equine repetitive element insertion polymorphism in Chinese horses

Transposable elements (TEs) are diverse, abundant, and complicated in genomes. They not only can drive the genome evolution process but can also act as special resources for adaptation. However, little is known about the evolutionary processes that shaped horses. In this work, 126 horse assemblages involved in most horse breeds in China were used to investigate the patterns of TE variation for the first time. By using RepeatMasker and melt software, we found that the horse-specific short interspersed repetitive elements family, equine repetitive elements (ERE1), exhibited polymorphisms in horse genomes. Phylogenetic analysis based on these ERE1 loci (minor allele frequency ≥0.05) revealed three major horse groups, namely, those in northern China, southern China, and Qinghai-Tibetan, which mirrors the result determined by SNPs to some extent. The present ERE1 family emerged ~0.26 to 1.77 Mya ago, with an activity peak at ~0.49 Mya, which matches the early stage of the horse lineage and decreases after the divergence of Equus caballus and Equus ferus przewalskii. To detect the functional ERE1(s) associated with adaptation, locus-specific branch length, genome-wide association study, and absolute allele frequency difference analyses were conducted and resulted in two common protein-coding genes annotated by candidate ERE1s. They were clustered into the vascular smooth muscle contraction (p = 0.01, EDNRA) and apelin signalling pathways (p = 0.02, NRF1). Notably, ERE1 insertion into the EDNRA gene showed a higher association with adaptation among southern China horses and other horses in 15 populations and 451 individuals (p = 4.55 e-8). Our results provide a comprehensive understanding of TE variations to analyse the phylogenetic relationships and traits relevant to adaptive evolution in horses.

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.

Characterization of Two Transposable Elements and an Ultra-Conserved Element Isolated in the Genome of Zootoca vivipara (Squamata, Lacertidae)

Transposable elements (TEs) constitute a considerable fraction of eukaryote genomes representing a major source of genetic variability. We describe two DNA sequences isolated in the lizard Zootoca vivipara, here named Zv516 and Zv817. Both sequences are single-copy nuclear sequences, including a truncation of two transposable elements (TEs), SINE Squam1 in Zv516 and a Tc1/Mariner-like DNA transposon in Zv817. FISH analyses with Zv516 showed the occurrence of interspersed signals of the SINE Squam1 sequence on all chromosomes of Z. vivipara and quantitative dot blot indicated that this TE is present with about 4700 copies in the Z. vivipara genome. FISH and dot blot with Zv817 did not produce clear hybridization signals. Bioinformatic analysis showed the presence of active SINE Squam 1 copies in the genome of different lacertids, in different mRNAs, and intronic and coding regions of various genes. The Tc1/Mariner-like DNA transposon occurs in all reptiles, excluding Sphenodon and Archosauria. Zv817 includes a trait of 284 bp, representing an amniote ultra-conserved element (UCE). Using amniote UCE homologous sequences from available whole genome sequences of major amniote taxonomic groups, we performed a phylogenetic analysis which retrieved Prototheria as the sister group of Metatheria and Eutheria. Within diapsids, Testudines are the sister group to Aves + Crocodylia (Archosauria), and Sphenodon is the sister group to Squamata. Furthermore, large trait regions flanking the UCE are conserved at family level.