Abundance Of Transposable Elements Research Articles

Varying recombination landscapes between individuals are driven by polymorphic transposable elements.

Meiotic recombination is a prominent force shaping genome evolution, and understanding the causes for varying recombination landscapes within and between species has remained a central, though challenging, question. Recombination rates are widely observed to negatively associate with the abundance of transposable elements (TEs), selfish genetic elements that move between genomic locations. While such associations are usually interpreted as recombination influencing the efficacy of selection at removing TEs, accumulating findings suggest that TEs could instead be the cause rather than the consequence. To test this prediction, we formally investigated the influence of polymorphic, putatively active TEs on recombination rates. We developed and benchmarked a novel approach that uses PacBio long-read sequencing to efficiently, accurately, and cost-effectively identify crossovers (COs), a key recombination product, among large numbers of pooled recombinant individuals. By applying this approach to Drosophila strains with distinct TE insertion profiles, we found that polymorphic TEs, especially RNA-based TEs and TEs with local enrichment of repressive marks, reduce the occurrence of COs. Such an effect leads to different CO frequencies between homologous sequences with and without TEs, contributing to varying CO maps between individuals. The suppressive effect of TEs on CO is further supported by two orthogonal approaches-analyzing the distributions of COs in panels of recombinant inbred lines in relation to TE polymorphism and applying marker-assisted estimations of CO frequencies to isogenic strains with and without transgenically inserted TEs. Our investigations reveal how the constantly changing mobilome can actively modify recombination landscapes, shaping genome evolution within and between species.

Prevalent Fast Evolution of Genes Involved in Heterochromatin Functions.

Heterochromatin is a gene-poor and repeat-rich genomic compartment universally found in eukaryotes. Despite its low transcriptional activity, heterochromatin plays important roles in maintaining genome stability, organizing chromosomes, and suppressing transposable elements. Given the importance of these functions, it is expected that genes involved in heterochromatin regulation would be highly conserved. Yet, a handful of these genes were found to evolve rapidly. To investigate whether these previous findings are anecdotal or general to genes modulating heterochromatin, we compile an exhaustive list of 106 candidate genes involved in heterochromatin functions and investigate their evolution over short and long evolutionary time scales in Drosophila. Our analyses find that these genes exhibit significantly more frequent evolutionary changes, both in the forms of amino acid substitutions and gene copy number change, when compared to genes involved in Polycomb-based repressive chromatin. While positive selection drives amino acid changes within both structured domains with diverse functions and intrinsically disordered regions, purifying selection may have maintained the proportions of intrinsically disordered regions of these proteins. Together with the observed negative associations between the evolutionary rate of these genes and the genomic abundance of transposable elements, we propose an evolutionary model where the fast evolution of genes involved in heterochromatin functions is an inevitable outcome of the unique functional roles of heterochromatin, while the rapid evolution of transposable elements may be an effect rather than cause. Our study provides an important global view of the evolution of genes involved in this critical cellular domain and provides insights into the factors driving the distinctive evolution of heterochromatin.

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.

Prevalent fast evolution of genes involved in heterochromatin functions.

Heterochromatin is a gene-poor and repeat-rich genomic compartment universally found in eukaryotes. Despite its low transcriptional activity, heterochromatin plays important roles in maintaining genome stability, organizing chromosomes, and suppressing transposable elements (TEs). Given the importance of these functions, it is expected that the genes involved in heterochromatin regulation would be highly conserved. Yet, a handful of these genes were found to evolve rapidly. To investigate whether these previous findings are anecdotal or general to genes modulating heterochromatin, we compile an exhaustive list of 106 candidate genes involved in heterochromatin functions and investigate their evolution over short and long evolutionary time scales in Drosophila. Our analyses find that these genes exhibit significantly more frequent evolutionary changes, both in the forms of amino acid substitutions and gene copy number change, when compared to genes involved in Polycomb-based repressive chromatin. While positive selection drives amino acid changes within both structured domains with diverse functions and intrinsically disordered regions (IDRs), purifying selection may have maintained the proportions of IDRs of these proteins. Together with the observed negative associations between evolutionary rates of these genes and genomic TE abundance, we propose an evolutionary model where the fast evolution of genes involved in heterochromatin functions is an inevitable outcome of the unique functional roles of heterochromatin, while the rapid evolution of TEs may be an effect rather than cause. Our study provides an important global view of the evolution of genes involved in this critical cellular domain and provides insights into the factors driving the distinctive evolution of heterochromatin.

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.

Consistent accumulation of transposable elements in species of the Hawaiian Tetragnatha spiny-leg adaptive radiation across the archipelago chronosequence

Abstract The ecological and phenotypic diversity observed in oceanic island radiations presents an evolutionary paradox: a high level of genetic variation is typically required for diversification, but species colonizing a new island commonly suffer from founder effects. This reduction in population size leads to lower genetic diversity, which ultimately results in a reduction in the efficiency of natural selection. Then, what is the source of genetic variation which acts as the raw material for ecological and phenotypic diversification in oceanic archipelagos? Transposable elements (TEs) are mobile genetic elements that have been linked to the generation of genetic diversity, and evidence suggests that TE activity and accumulation along the genome can result from reductions in population size. Here, we use the Hawaiian spiny-leg spider radiation (Tetragnatha) to test whether TE accumulation increases due to demographic processes associated with island colonization. We sequenced and quantified TEs in 23 individuals representing 16 species from the spiny-leg radiation and 4 individuals from its sister radiation, the Hawaiian web-building Tetragnatha. Our results show that founder effects resulting from colonization of new islands have not resulted in TE accumulation over evolutionary time. Specifically, we found no evidence for increase in abundance of specific TE superfamilies, nor an accumulation of ‘young TEs’ in lineages which have recently colonized a new island or are present in islands with active volcanoes. We also found that the DNA/hAT transposon superfamily is by far the most abundant TE superfamily in the Tetragnatha radiation. This work shows that there is no clear trend of increasing TE abundance for the spiny-leg radiation across the archipelago chronosequence, and TE accumulation is not affected by population oscillations associated with island colonization events. Therefore, despite their known role in the generation of genetic diversity, TE activity does not appear to be the mechanism explaining the evolutionary paradox of insular diversification in the Tetragnatha spiny-leg radiation.

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.

Technology to the rescue: how to uncover the role of transposable elements in preimplantation development.

Transposable elements (TEs) are highly expressed in preimplantation development. Preimplantation development is the phase when the cells of the early embryo undergo the first cell fate choice and change from being totipotent to pluripotent. A range of studies have advanced our understanding of TEs in preimplantation, as well as their epigenetic regulation and functional roles. However, many questions remain about the implications of TE expression during early development. Challenges originate first due to the abundance of TEs in the genome, and second because of the limited cell numbers in preimplantation. Here we review the most recent technological advancements promising to shed light onto the role of TEs in preimplantation development. We explore novel avenues to identify genomic TE insertions and improve our understanding of the regulatory mechanisms and roles of TEs and their RNA and protein products during early development.

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.

Phased assembly of neo-sex chromosomes reveals extensive Y degeneration and rapid genome evolution in Rumex hastatulus.

Y chromosomes are thought to undergo progressive degeneration due to stepwise loss of recombination and subsequent reduction in selection efficiency. However, the timescales and evolutionary forces driving degeneration remain unclear. To investigate the evolution of sex chromosomes on multiple timescales, we generated a high-quality phased genome assembly of the massive older (<10MYA) and neo (<200,000 years) sex chromosomes in the XYY cytotype of the dioecious plant Rumex hastatulus and a hermaphroditic outgroup R. salicifolius. Our assemblies, supported by fluorescence in situ hybridization, confirmed that the neo-sex chromosomes were formed by two key events: an X-autosome fusion and a reciprocal translocation between the homologous autosome and the Y chromosome. The enormous sex-linked regions of the X (296 MB) and two Y chromosomes (503 MB) both evolved from large repeat-rich genomic regions with low recombination; however, the complete loss of recombination on the Y still led to over 30% gene loss and major rearrangements. In the older sex-linked region, there has been a significant increase in transposable element abundance, even into and near genes. In the neo sex-linked regions, we observed evidence of extensive rearrangements without gene degeneration and loss. Overall, we inferred significant degeneration during the first 10 million years of Y chromosome evolution but not on very short timescales. Our results indicate that even when sex chromosomes emerge from repetitive regions of already-low recombination, the complete loss of recombination on the Y chromosome still leads to a substantial increase in repetitive element content and gene degeneration.

Annotation and Comparative Genomics of Prokaryotic Transposable Elements.

The data generated in nearly 30years of bacterial genome sequencing has revealed the abundance of transposable elements (TE) and their importance in genome and transcript remodeling through the mediation of DNA insertions and deletions, structural rearrangements, and regulation of gene expression. Furthermore, what we have learned from studying transposition mechanisms and their regulation in bacterial TE is fundamental to our current understanding of TE in other organisms because much of what has been observed in bacteria is conserved in all domains of life. However, unlike eukaryotic TE, prokaryotic TE sequester and transmit important classes of genes that impact host fitness, such as resistance to antibiotics and heavy metals and virulence factors affecting animals and plants, among other acquired traits. This provides dynamism and plasticity to bacteria, which would otherwise be propagated clonally. The insertion sequences (IS), the simplest form of prokaryotic TE, are autonomous and compact mobile genetic elements. These can be organized into compound transposons, in which two similar IS can flank any DNA segment and render it transposable. Other more complex structures, called unit transposons, can be grouped into four major families (Tn3, Tn7, Tn402, Tn554) with specific genetic characteristics. This chapter will revisit the prominent structural features of these elements, focusing on a genomic annotation framework and comparative analysis. Relevant aspects of TE will also be presented, stressing their key position in genome impact and evolution, especially in the emergence of antimicrobial resistance and other adaptive traits.

Quantifying the role of genome size and repeat content in adaptive variation and the architecture of flowering time in Amaranthus tuberculatus.

Genome size variation, largely driven by repeat content, is poorly understood within and among populations, limiting our understanding of its significance for adaptation. Here we characterize intraspecific variation in genome size and repeat content across 186 individuals of Amaranthus tuberculatus, a ubiquitous native weed that shows flowering time adaptation to climate across its range and in response to agriculture. Sequence-based genome size estimates vary by up to 20% across individuals, consistent with the considerable variability in the abundance of transposable elements, unknown repeats, and rDNAs across individuals. The additive effect of this variation has important phenotypic consequences-individuals with more repeats, and thus larger genomes, show slower flowering times and growth rates. However, compared to newly-characterized gene copy number and polygenic nucleotide changes underlying variation in flowering time, we show that genome size is a marginal contributor. Differences in flowering time are reflected by genome size variation across sexes and marginally, habitats, while polygenic variation and a gene copy number variant within the ATP synthesis pathway show consistently stronger environmental clines than genome size. Repeat content nonetheless shows non-neutral distributions across the genome, and across latitudinal and environmental gradients, demonstrating the numerous governing processes that in turn influence quantitative genetic variation for phenotypes key to plant adaptation.

Insights into the Superrosids phylogeny and flavonoid synthesis from the telomere-to-telomere gap-free genome assembly of Penthorum chinense Pursh.

The completion of the first telomere-to-telomere (T2T) genome assembly of Penthorum chinense Pursh (PC), a prominent medicinal plant in China, represents a significant achievement. This assembly spans a length of 257.5Mb and consists of nine chromosomes. PC's notably smaller genome size in Saxifragales, compared to that of Paeonia ostii, can be attributed to the low abundance of transposable elements. By utilizing single-copy genes from 30 species, including 28 other Superrosids species, we successfully resolved a previously debated Superrosids phylogeny. Our findings unveiled Saxifragales as the sister group to the core rosids, with both being the sister group to Vitales. Utilizing previously characterized cytochrome P450 (CYP) genes, we predicted the compound classes that most CYP genes of PC are involved in synthesizing, providing insight into PC's potential metabolic diversity. Metabolomic and transcriptomic data revealed that the richest sources of the three most noteworthy medicinal components in PC are young leaves and flowers. We also observed higher activity of upstream genes in the flavonoid synthesis pathway in these plant parts. Additionally, through weighted gene co-expression network analysis, we identified gene regulatory networks associated with the three medicinal components. Overall, these findings deepen our understanding of PC, opening new avenues for further research and exploration.

In Silico Chromosome Mapping of the Male-Specific/Linked Loci in the Jade Perch (Scortum barcoo) Suggests Chromosome 19 as the Putative Y Sex Chromosome

Jade perch (Scortum barcoo) has an XX/XY sex-determination system (SDS); however, its sex chromosomes and sex-determining region remain unknown. The recent availability of the jade perch chromosome-level genomic data provides a valuable resource for pinpointing the location of functional genes and the whole genomic structure. In this study, we conducted. In silico chromosome mapping of male-specific/linked loci of jade perch and identified a potential 11.18 Mb male-linked region localized on chromosome 19 (SBA19). Repeat annotations of the male-linked region revealed an abundance of transposable elements, particularly Ty3/Gypsy and novel repeats. Sequence analysis of this region identified a remnant of amh gene, which is considered a potential candidate for SDS in many teleosts. A duplicate copy of amh remnant was located at SBA6. These duplicated amh copies were highly similar to those of XX/XY SDS in teleosts, in which one copy of amh was identified on the Y sex chromosome. Taken all together, we hypothesize SBA19 as the putative sex chromosome and the 11.18 Mb male-linked region to be a potential male-determining region.

Welcome to the big leaves: Best practices for improving genome annotation in non-model plant genomes.

Robust standards to evaluate quality and completeness are lacking in eukaryotic structural genome annotation, as genome annotation software is developed using model organisms and typically lacks benchmarking to comprehensively evaluate the quality and accuracy of the final predictions. The annotation of plant genomes is particularly challenging due to their large sizes, abundant transposable elements, and variable ploidies. This study investigates the impact of genome quality, complexity, sequence read input, and method on protein-coding gene predictions. The impact of repeat masking, long-read and short-read inputs, and de novo and genome-guided protein evidence was examined in the context of the popular BRAKER and MAKER workflows for five plant genomes. The annotations were benchmarked for structural traits and sequence similarity. Benchmarks that reflect gene structures, reciprocal similarity search alignments, and mono-exonic/multi-exonic gene counts provide a more complete view of annotation accuracy. Transcripts derived from RNA-read alignments alone are not sufficient for genome annotation. Gene prediction workflows that combine evidence-based and ab initio approaches are recommended, and a combination of short and long reads can improve genome annotation. Adding protein evidence from de novo assemblies, genome-guided transcriptome assemblies, or full-length proteins from OrthoDB generates more putative false positives as implemented in the current workflows. Post-processing with functional and structural filters is highly recommended. While the annotation of non-model plant genomes remains complex, this study provides recommendations for inputs and methodological approaches. We discuss a set of best practices to generate an optimal plant genome annotation and present a more robust set of metrics to evaluate the resulting predictions.

H4K16ac activates the transcription of transposable elements and contributes to their cis-regulatory function

Mammalian genomes harbor abundant transposable elements (TEs) and their remnants, with numerous epigenetic repression mechanisms enacted to silence TE transcription. However, TEs are upregulated during early development, neuronal lineage, and cancers, although the epigenetic factors contributing to the transcription of TEs have yet to be fully elucidated. Here, we demonstrate that the male-specific lethal (MSL)-complex-mediated histone H4 acetylation at lysine 16 (H4K16ac) is enriched at TEs in human embryonic stem cells (hESCs) and cancer cells. This in turn activates transcription of subsets of full-length long interspersed nuclear elements (LINE1s, L1s) and endogenous retrovirus (ERV) long terminal repeats (LTRs). Furthermore, we show that the H4K16ac-marked L1 and LTR subfamilies display enhancer-like functions and are enriched in genomic locations with chromatin features associated with active enhancers. Importantly, such regions often reside at boundaries of topologically associated domains and loop with genes. CRISPR-based epigenetic perturbation and genetic deletion of L1s reveal that H4K16ac-marked L1s and LTRs regulate the expression of genes in cis. Overall, TEs enriched with H4K16ac contribute to the cis-regulatory landscape at specific genomic locations by maintaining an active chromatin landscape at TEs.

Chromosome-level genomes of two armyworms, Mythimna separata and Mythimna loreyi, provide insights into the biosynthesis and reception of sex pheromones.

Mythimna separata and Mythimna loreyi are global pests of gramineous cereals, heavily controlled with synthetic insecticides. Here, we generated two high-quality chromosome-level genome assemblies for M. separata (688 Mb) and M. loreyi (683 Mb). Our analysis identified Z and W chromosomes, with few genes and abundant transposable elements (TEs) found on the W chromosome. We also observed a recent explosion of long interspersed nuclear elements (LINEs), which contributed to the larger genomes of Mythimna. The two armyworms diverged ~10.5 MYA, with only three chromosomes have intrachromosomal rearrangements. Additionally, we observed a tandem repeat expansion of α-amylase genes in Mythimna, which may promote the digestion of carbohydrates and exacerbate their damage to crops. Furthermore, we inferred the sex pheromone biosynthesis pathway for M. separata, M. loreyi and Spodoptera frugiperda. We discovered that M. loreyi and S. frugiperda synthesized the same major constituents of sex pheromones through different pathways. Specifically, the double bonds in the dominant sex pheromone components of S. frugiperda were generated by Δ9- and Δ11-desaturase, while they were generated by Δ11-desaturase and chain-shortening reactions in M. loreyi. We also identified pheromone receptor (PR) genes and inferred their corresponding components. These findings provide a better understanding of sex pheromone communication and promote the development of a new pest control strategy involving pheromone traps, which are more effective and environmentally friendly than current strategies.

RBBP4 is an epigenetic barrier for the induced transition of pluripotent stem cells into totipotent 2C-like cells.

Cellular totipotency is critical for whole-organism generation, yet how totipotency is established remains poorly illustrated. Abundant transposable elements (TEs) are activated in totipotent cells, which is critical for embryonic totipotency. Here, we show that the histone chaperone RBBP4, but not its homolog RBBP7, is indispensable for maintaining the identity of mouse embryonic stem cells (mESCs). Auxin-induced degradation of RBBP4, but not RBBP7, reprograms mESCs to the totipotent 2C-like cells. Also, loss of RBBP4 enhances transition from mESCs to trophoblast cells. Mechanistically, RBBP4 binds to the endogenous retroviruses (ERVs) and functions as an upstream regulator by recruiting G9a to deposit H3K9me2 on ERVL elements, and recruiting KAP1 to deposit H3K9me3 on ERV1/ERVK elements, respectively. Moreover, RBBP4 facilitates the maintenance of nucleosome occupancy at the ERVK and ERVL sites within heterochromatin regions through the chromatin remodeler CHD4. RBBP4 depletion leads to the loss of the heterochromatin marks and activation of TEs and 2C genes. Together, our findings illustrate that RBBP4 is required for heterochromatin assembly and is a critical barrier for inducing cell fate transition from pluripotency to totipotency.

The gyrfalcon (Falco rusticolus) genome.

High-quality genome assemblies are characterized by high-sequence contiguity, completeness, and a low error rate, thus providing the basis for a wide array of studies focusing on natural species ecology, conservation, evolution, and population genomics. To provide this valuable resource for conservation projects and comparative genomics studies on gyrfalcon (Falco rusticolus), we sequenced and assembled the genome of this species using third-generation sequencing strategies and optical maps. Here, we describe a highly contiguous and complete genome assembly comprising 20 scaffolds and 13 contigs with a total size of 1.193 Gbp, including 8,064 complete Benchmarking Universal Single-Copy Orthologs (BUSCOs) of the total 8,338 BUSCO groups present in the library aves_odb10. Of these BUSCO genes, 96.7% were complete, 96.1% were present as a single copy, and 0.6% were duplicated. Furthermore, 0.8% of BUSCO genes were fragmented and 2.5% (210) were missing. A de novo search for transposable elements (TEs) identified 5,716 TEs that masked 7.61% of the F. rusticolus genome assembly when combined with publicly available TE collections. Long interspersed nuclear elements, in particular, the element Chicken-repeat 1 (CR1), were the most abundant TEs in the F. rusticolus genome. A de novo first-pass gene annotation was performed using 293,349 PacBio Iso-Seq transcripts and 496,195 transcripts derived from the assembly of 42,429,525 Illumina PE RNA-seq reads. In all, 19,602 putative genes, of which 59.31% were functionally characterized and associated with Gene Ontology terms, were annotated. A comparison of the gyrfalcon genome assembly with the publicly available assemblies of the domestic chicken (Gallus gallus), zebra finch (Taeniopygia guttata), and hummingbird (Calypte anna) revealed several genome rearrangements. In particular, nine putative chromosome fusions were identified in the gyrfalcon genome assembly compared with those in the G. gallus genome assembly. This genome assembly, its annotation for TEs and genes, and the comparative analyses presented, complement and strength the base of high-quality genome assemblies and associated resources available for comparative studies focusing on the evolution, ecology, and conservation of Aves.

Intracellular common gardens reveal niche differentiation in transposable element community during bacterial adaptive evolution

The distribution and abundance of transposable elements across the tree of life have significantly shaped the evolution of cellular organisms, but the underlying mechanisms shaping these ecological patterns remain elusive. Here we establish a “common garden” approach to study causal ecological interactions between a xenogeneic conditional lethal sacB gene and the community of transposable insertion sequences (ISs) in a multipartite prokaryote genome. Xenogeneic sacB of low, medium, or high GC content was individually inserted into three replicons of a model bacterium Sinorhizobium fredii, and exhibited replicon- and GC-dependent variation in genetic stability. This variation was largely attributable to multidimensional niche differentiation for IS community members. The transposition efficiency of major active ISs depended on the nucleoid-associated xenogeneic silencer MucR. Experimentally eliminating insertion activity of specific ISs by deleting MucR strongly demonstrated a dominant role of niche differentiation among ISs. This intracellular common garden approach in the experimental evolution context allows not only for evaluating genetic stability of natural and synthetic xenogeneic genes of different sequence signatures in host cells but also for tracking and testing causal relationships in unifying ecological principles in genome ecology.