RepeatFiller newly identifies megabases of aligning repetitive sequences and improves annotations of conserved non-exonic elements

Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample

Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.

The evolution history of an allotetraploid mangrove tree analysed with a new tool Allo4D

Mangrove species are broadly classified as true mangroves and mangrove associates. The latter are amphibious plants that can survive in the intertidal zone and reproduce naturally in terrestrial environments. Their widespread distribution and extensive adaptability make them ideal research materials for exploring adaptive evolution. In this study, we de novo assembled two genomes of mangrove associates (the allotetraploid Barringtonia racemosa (2n = 4x = 52) and diploid Barringtonia asiatica (2n = 2x = 26)) to investigate the role of allopolyploidy in the evolutionary history of mangrove species. We developed a new allotetraploid-dividing tool Allo4D to distinguish between allotetraploid scaffold-scale subgenomes and verified its accuracy and reliability using real and simulated data. According to the two subgenomes of allotetraploid B. racemosa divided using Allo4D, the allopolyploidization event was estimated to have occurred approximately one million years ago (Mya). We found that B. racemosa, B. asiatica, and Diospyros lotus shared a whole genome duplication (WGD) event during the K-Pg (Cretaceous-Paleozoic) period. K-Pg WGD and recent allopolyploidization events contributed to the speciation of B. racemosa and its adaptation to coastal habitats. We found that genes in the glucosinolates (GSLs) pathway, an essential pathway in response to various biotic and abiotic stresses, expanded rapidly in B. racemosa during polyploidization. In summary, this study provides a typical example of the adaptation of allopolyploid plants to extreme environmental conditions. The newly developed tool, Allo4D, can effectively divide allotetraploid subgenomes and explore the evolutionary history of polyploid plants, especially for species whose ancestors are unknown or extinct.

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Segregation between an ornamental and a disease driver gene provides insights into pigment cell regulation

Genetic interactions are adaptive within a species. Hybridization can disrupt such species-specific genetic interactions and creates novel interactions that alter the hybrid progeny overall fitness. Hybrid incompatibility, which refers to degenerative genetic interactions that decrease the overall hybrid survival, is one of the results from combining two diverged genomes in hybrids. The discovery of spontaneous lethal tumorigenesis and underlying genetic interactions in select hybrids between diverged Xiphophorus species showed that lethal pathological process can result from degenerative genetic interactions. Such genetic interactions leading to lethal phenotype are thought to shield gene flow between diverged species. However, hybrids between certain Xiphophorus species do not develop such tumors. Here we report the identification of a locus residing in the genome of one Xiphophorus species that represses an oncogene from a different species. Our finding provides insights into normal and pathological pigment cell development, regulation and molecular mechanism in hybrid incompatibility.

Significance

The Dobzhansky-Muller model states epistatic interactions occurred between genes in diverged species underlies hybrid incompatibility. There are a few vertebrate interspecies hybrid cases that support the Dobzhansky-Muller model. This study reports a fish hybrid system where incompatible genetic interactions are involved in neuronal regulation of pigment cell biology, and also identified a novel point of regulation for pigment cells.

PlantFUNCO: Integrative Functional Genomics Database Reveals Clues into Duplicates Divergence Evolution

Evolutionary epigenomics and, more generally, evolutionary functional genomics, are emerging fields that study how non-DNA-encoded alterations in gene expression regulation are an important form of plasticity and adaptation. Previous evidence analyzing plants' comparative functional genomics has mostly focused on comparing same assay-matched experiments, missing the power of heterogeneous datasets for conservation inference. To fill this gap, we developed PlantFUN(ctional)CO(nservation) database, which is constituted by several tools and two main resources: interspecies chromatin states and functional genomics conservation scores, presented and analyzed in this work for three well-established plant models (Arabidopsis thaliana, Oryza sativa, and Zea mays). Overall, PlantFUNCO elucidated evolutionary information in terms of cross-species functional agreement. Therefore, providing a new complementary comparative-genomics source for assessing evolutionary studies. To illustrate the potential applications of this database, we replicated two previously published models predicting genetic redundancy in A. thaliana and found that chromatin states are a determinant of paralogs degree of functional divergence. These predictions were validated based on the phenotypes of mitochondrial alternative oxidase knockout mutants under two different stressors. Taking all the above into account, PlantFUNCO aim to leverage data diversity and extrapolate molecular mechanisms findings from different model organisms to determine the extent of functional conservation, thus, deepening our understanding of how plants epigenome and functional noncoding genome have evolved. PlantFUNCO is available at https://rocesv.github.io/PlantFUNCO.

Chromosome level genome assembly of the Etruscan shrew Suncus etruscus

Suncus etruscus is one of the world's smallest mammals, with an average body mass of about 2 grams. The Etruscan shrew's small body is accompanied by a very high energy demand and numerous metabolic adaptations. Here we report a chromosome-level genome assembly using PacBio long read sequencing, 10X Genomics linked short reads, optical mapping, and Hi-C linked reads. The assembly is partially phased, with the 2.472 Gbp primary pseudohaplotype and 1.515 Gbp alternate. We manually curated the primary assembly and identified 22 chromosomes, including X and Y sex chromosomes. The NCBI genome annotation pipeline identified 39,091 genes, 19,819 of them protein-coding. We also identified segmental duplications, inferred GO term annotations, and computed orthologs of human and mouse genes. This reference-quality genome will be an important resource for research on mammalian development, metabolism, and body size control.

Integrating gene annotation with orthology inference at scale

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

The genome of the pygmy right whale illuminates the evolution of rorquals

BACKGROUND

Baleen whales are a clade of gigantic and highly specialized marine mammals. Their genomes have been used to investigate their complex evolutionary history and to decipher the molecular mechanisms that allowed them to reach these dimensions. However, many unanswered questions remain, especially about the early radiation of rorquals and how cancer resistance interplays with their huge number of cells. The pygmy right whale is the smallest and most elusive among the baleen whales. It reaches only a fraction of the body length compared to its relatives and it is the only living member of an otherwise extinct family. This placement makes the pygmy right whale genome an interesting target to update the complex phylogenetic past of baleen whales, because it splits up an otherwise long branch that leads to the radiation of rorquals. Apart from that, genomic data of this species might help to investigate cancer resistance in large whales, since these mechanisms are not as important for the pygmy right whale as in other giant rorquals and right whales.

RESULTS

Here, we present a first de novo genome of the species and test its potential in phylogenomics and cancer research. To do so, we constructed a multi-species coalescent tree from fragments of a whole-genome alignment and quantified the amount of introgression in the early evolution of rorquals. Furthermore, a genome-wide comparison of selection rates between large and small-bodied baleen whales revealed a small set of conserved candidate genes with potential connections to cancer resistance.

CONCLUSIONS

Our results suggest that the evolution of rorquals is best described as a hard polytomy with a rapid radiation and high levels of introgression. The lack of shared positive selected genes between different large-bodied whale species supports a previously proposed convergent evolution of gigantism and hence cancer resistance in baleen whales.

Prevalent Introgression Underlies Convergent Evolution in the Diversification of Pungitius Sticklebacks

New mutations and standing genetic variations contribute significantly to repeated phenotypic evolution in sticklebacks. However, less is known about the role of introgression in this process. We analyzed taxonomically and geographically comprehensive genomic data from Pungitius sticklebacks to decipher the extent of introgression and its consequences for the diversification of this genus. Our results demonstrate that introgression is more prevalent than suggested by earlier studies. Although gene flow was generally bidirectional, it was often asymmetric and left unequal genomic signatures in hybridizing species, which might, at least partly, be due to biased hybridization and/or population size differences. In several cases, introgression of variants from one species to another was accompanied by transitions of pelvic and/or lateral plate structures-important diagnostic traits in Pungitius systematics-and frequently left signatures of adaptation in the core gene regulatory networks of armor trait development. This finding suggests that introgression has been an important source of genetic variation and enabled phenotypic convergence among Pungitius sticklebacks. The results highlight the importance of introgression of genetic variation as a source of adaptive variation underlying key ecological and taxonomic traits. Taken together, our study indicates that introgression-driven convergence likely explains the long-standing challenges in resolving the taxonomy and systematics of this small but phenotypically highly diverse group of fish.

Loss of a gluconeogenic muscle enzyme contributed to adaptive metabolic traits in hummingbirds

Hummingbirds possess distinct metabolic adaptations to fuel their energy-demanding hovering flight, but the underlying genomic changes are largely unknown. Here, we generated a chromosome-level genome assembly of the long-tailed hermit and screened for genes that have been specifically inactivated in the ancestral hummingbird lineage. We discovered that FBP2 (fructose-bisphosphatase 2), which encodes a gluconeogenic muscle enzyme, was lost during a time period when hovering flight evolved. We show that FBP2 knockdown in an avian muscle cell line up-regulates glycolysis and enhances mitochondrial respiration, coincident with an increased mitochondria number. Furthermore, genes involved in mitochondrial respiration and organization have up-regulated expression in hummingbird flight muscle. Together, these results suggest that FBP2 loss was likely a key step in the evolution of metabolic muscle adaptations required for true hovering flight.

The Landscape of Liver Chromatin Accessibility and Conserved Non-coding Elements in Larimichthys crocea, Nibea albiflora, and Lateolabrax maculatus

Large yellow croaker (Larimichthys crocea), yellow drum (Nibea albiflora), and Chinese seabass (Lateolabrax maculatus) are important economic marine fishes in China. The conserved non-coding elements (CNEs) in the liver tissues of the three kinds of fish are directly or indirectly involved in the regulation of gene expression and affect liver functions. However, the fishes' CNEs and even chromatin accessibility landscape have not been effectively investigated. Hence, this study established the landscapes of the fishes' genome-wide chromatin accessibility and CNEs by detecting regions of the open chromatin in their livers using an assay for transposase-accessible chromatin by high-throughput sequencing (ATAC-seq) and comparative genomics approach. The results showed that Smad1, Sp1, and Foxl1 transcription factor binding motifs were considerably enriched in the chromatin accessibility landscape in the liver of the three species, and the three transcription factors (TFs) had a wide range of common targets. The hypothetical gene set was targeted by one, two, or all three TFs, which was much higher than would be expected for an accidental outcome. The gene sets near the CNEs were mainly enriched through processes such as a macromolecule metabolic process and ribonucleoprotein complex biogenesis. The active CNEs were found in the promoter regions of genes such as ap1g1, hax1, and ndufs2. And 5 CNEs were predicted to be highly conserved active enhancers. These results demonstrated that Smad1, Sp1, and Foxl1 might be related to the liver function in the three fishes. In addition, we found a series of ATAC-seq-labeled CNEs located in the gene promoter regions, and highly conserved H3k27ac + -labeled CNEs located in the liver function genes. The highly conserved nature of these regulatory elements suggests that they play important roles in the liver in fish. This study mined the landscape of chromatin accessibility and CNEs of three important economic fishes to fill the knowledge gaps in this field. Moreover, the work provides useful data for the industrial application and theoretical research of these three fish species.

Contradictory Phylogenetic Signals in the Laurasiatheria Anomaly Zone

Relationships among laurasiatherian clades represent one of the most highly disputed topics in mammalian phylogeny. In this study, we attempt to disentangle laurasiatherian interordinal relationships using two independent genome-level approaches: (1) quantifying retrotransposon presence/absence patterns, and (2) comparisons of exon datasets at the levels of nucleotides and amino acids. The two approaches revealed contradictory phylogenetic signals, possibly due to a high level of ancestral incomplete lineage sorting. The positions of Eulipotyphla and Chiroptera as the first and second earliest divergences were consistent across the approaches. However, the phylogenetic relationships of Perissodactyla, Cetartiodactyla, and Ferae, were contradictory. While retrotransposon insertion analyses suggest a clade with Cetartiodactyla and Ferae, the exon dataset favoured Cetartiodactyla and Perissodactyla. Future analyses of hitherto unsampled laurasiatherian lineages and synergistic analyses of retrotransposon insertions, exon and conserved intron/intergenic sequences might unravel the conflicting patterns of relationships in this major mammalian clade.

Chromosome-length genome assembly and linkage map of a critically endangered Australian bird: the helmeted honeyeater

Background

The helmeted honeyeater (Lichenostomus melanops cassidix) is a Critically Endangered bird endemic to Victoria, Australia. To aid its conservation, the population is the subject of genetic rescue. To understand, monitor, and modulate the effects of genetic rescue on the helmeted honeyeater genome, a chromosome-length genome and a high-density linkage map are required.

Results

We used a combination of Illumina, Oxford Nanopore, and Hi-C sequencing technologies to assemble a chromosome-length genome of the helmeted honeyeater, comprising 906 scaffolds, with length of 1.1 Gb and scaffold N50 of 63.8 Mb. Annotation comprised 57,181 gene models. Using a pedigree of 257 birds and 53,111 single-nucleotide polymorphisms, we obtained high-density linkage and recombination maps for 25 autosomes and Z chromosome. The total sex-averaged linkage map was 1,347 cM long, with the male map being 6.7% longer than the female map. Recombination maps revealed sexually dimorphic recombination rates (overall higher in males), with average recombination rate of 1.8 cM/Mb. Comparative analyses revealed high synteny of the helmeted honeyeater genome with that of 3 passerine species (e.g., 32 Hi-C scaffolds mapped to 30 zebra finch autosomes and Z chromosome). The genome assembly and linkage map suggest that the helmeted honeyeater exhibits a fission of chromosome 1A into 2 chromosomes relative to zebra finch. PSMC analysis showed a ∼15-fold decline in effective population size to ∼60,000 from mid- to late Pleistocene.

Conclusions

The annotated chromosome-length genome and high-density linkage map provide rich resources for evolutionary studies and will be fundamental in guiding conservation efforts for the helmeted honeyeater.

Gene losses in the common vampire bat illuminate molecular adaptations to blood feeding

Vampire bats are the only mammals that feed exclusively on blood. To uncover genomic changes associated with this dietary adaptation, we generated a haplotype-resolved genome of the common vampire bat and screened 27 bat species for genes that were specifically lost in the vampire bat lineage. We found previously unknown gene losses that relate to reduced insulin secretion (FFAR1 and SLC30A8), limited glycogen stores (PPP1R3E), and a unique gastric physiology (CTSE). Other gene losses likely reflect the biased nutrient composition (ERN2 and CTRL) and distinct pathogen diversity of blood (RNASE7) and predict the complete lack of cone-based vision in these strictly nocturnal bats (PDE6H and PDE6C). Notably, REP15 loss likely helped vampire bats adapt to high dietary iron levels by enhancing iron excretion, and the loss of CYP39A1 could have contributed to their exceptional cognitive abilities. These findings enhance our understanding of vampire bat biology and the genomic underpinnings of adaptations to blood feeding.

Interspecies transcriptomics identify genes that underlie disproportionate foot growth in jerboas

Despite the great diversity of vertebrate limb proportion and our deep understanding of the genetic mechanisms that drive skeletal elongation, little is known about how individual bones reach different lengths in any species. Here, we directly compare the transcriptomes of homologous growth cartilages of the mouse (Mus musculus) and bipedal jerboa (Jaculus jaculus), the latter of which has "mouse-like" arms but extremely long metatarsals of the feet. Intersecting gene-expression differences in metatarsals and forearms of the two species revealed that about 10% of orthologous genes are associated with the disproportionately rapid elongation of neonatal jerboa feet. These include genes and enriched pathways not previously associated with endochondral elongation as well as those that might diversify skeletal proportion in addition to their known requirements for bone growth throughout the skeleton. We also identified transcription regulators that might act as "nodes" for sweeping differences in genome expression between species. Among these, Shox2, which is necessary for proximal limb elongation, has gained expression in jerboa metatarsals where it has not been detected in other vertebrates. We show that Shox2 is sufficient to increase mouse distal limb length, and a nearby putative cis-regulatory region is preferentially accessible in jerboa metatarsals. In addition to mechanisms that might directly promote growth, we found evidence that jerboa foot elongation may occur in part by de-repressing latent growth potential. The genes and pathways that we identified here provide a framework to understand the modular genetic control of skeletal growth and the remarkable malleability of vertebrate limb proportion.

Finding and Characterizing Repeats in Plant Genomes

Plant genomes contain a particularly high proportion of repeated structures of various types. This chapter proposes a guided tour of the available software that can help biologists to scan automatically for these repeats in sequence data or check hypothetical models intended to characterize their structures. Since transposable elements (TEs) are a major source of repeats in plants, many methods have been used or developed for this broad class of sequences. They are representative of the range of tools available for other classes of repeats and we have provided two sections on this topic (for the analysis of genomes or directly of sequenced reads), as well as a selection of the main existing software. It may be hard to keep up with the profusion of proposals in this dynamic field and the rest of the chapter is devoted to the foundations of an efficient search for repeats and more complex patterns. We first introduce the key concepts of the art of indexing and mapping or querying sequences. We end the chapter with the more prospective issue of building models of repeat families. We present the Machine Learning approach first, seeking to build predictors automatically for some families of ET, from a set of sequences known to belong to this family. A second approach, the linguistic (or syntactic) approach, allows biologists to describe themselves and check the validity of models of their favorite repeat family.

Evolutionary conservation in noncoding genomic regions

Humans may share more genomic commonalities with other species than previously thought. According to current estimates, ~5% of the human genome is functionally constrained, which is a much larger fraction than the ~1.5% occupied by annotated protein-coding genes. Hence, ~3.5% of the human genome comprises likely functional conserved noncoding elements (CNEs) preserved among organisms, whose common ancestors existed throughout hundreds of millions of years of evolution. As whole-genome sequencing emerges as a standard procedure in genetic analyses, interpretation of variations in CNEs, including the elucidation of mechanistic and functional roles, becomes a necessity. Here, we discuss the phenomenon of noncoding conservation via four dimensions (sequence, regulatory conservation, spatiotemporal expression, and structure) and the potential significance of CNEs in phenotype variation and disease.

Convergent Losses of TLR5 Suggest Altered Extracellular Flagellin Detection in Four Mammalian Lineages

Toll-like receptors (TLRs) play an important role for the innate immune system by detecting pathogen-associated molecular patterns. TLR5 encodes the major extracellular receptor for bacterial flagellin and frequently evolves under positive selection, consistent with coevolutionary arms races between the host and pathogens. Furthermore, TLR5 is inactivated in several vertebrates and a TLR5 stop codon polymorphism is widespread in human populations. Here, we analyzed the genomes of 120 mammals and discovered that TLR5 is convergently lost in four independent lineages, comprising guinea pigs, Yangtze river dolphin, pinnipeds, and pangolins. Validated inactivating mutations, absence of protein-coding transcript expression, and relaxed selection on the TLR5 remnants confirm these losses. PCR analysis further confirmed the loss of TLR5 in the pinniped stem lineage. Finally, we show that TLR11, encoding a second extracellular flagellin receptor, is also absent in these four lineages. Independent losses of TLR5 and TLR11 suggest that a major pathway for detecting flagellated bacteria is not essential for different mammals and predicts an impaired capacity to sense extracellular flagellin.

The mole genome reveals regulatory rearrangements associated with adaptive intersexuality

Linking genomic variation to phenotypical traits remains a major challenge in evolutionary genetics. In this study, we use phylogenomic strategies to investigate a distinctive trait among mammals: the development of masculinizing ovotestes in female moles. By combining a chromosome-scale genome assembly of the Iberian mole, Talpa occidentalis, with transcriptomic, epigenetic, and chromatin interaction datasets, we identify rearrangements altering the regulatory landscape of genes with distinct gonadal expression patterns. These include a tandem triplication involving CYP17A1, a gene controlling androgen synthesis, and an intrachromosomal inversion involving the pro-testicular growth factor gene FGF9, which is heterochronically expressed in mole ovotestes. Transgenic mice with a knock-in mole CYP17A1 enhancer or overexpressing FGF9 showed phenotypes recapitulating mole sexual features. Our results highlight how integrative genomic approaches can reveal the phenotypic impact of noncoding sequence changes.

A highly contiguous genome assembly of the bat hawkmoth Hyles vespertilio (Lepidoptera: Sphingidae)

BACKGROUND

Adapted to different ecological niches, moth species belonging to the Hyles genus exhibit a spectacular diversity of larval color patterns. These species diverged ∼7.5 million years ago, making this rather young genus an interesting system to study a wide range of questions including the process of speciation, ecological adaptation, and adaptive radiation.

RESULTS

Here we present a high-quality genome assembly of the bat hawkmoth Hyles vespertilio, the first reference genome of a member of the Hyles genus. We generated 51× Pacific Biosciences long reads with an average read length of 8.9 kb. Pacific Biosciences reads longer than 4 kb were assembled into contigs, resulting in a 651.4-Mb assembly consisting of 530 contigs with an N50 value of 7.5 Mb. The circular mitochondrial contig has a length of 15,303 bp. The H. vespertilio genome is very repeat-rich and exhibits a higher repeat content (50.3%) than other Bombycoidea species such as Bombyx mori (45.7%) and Manduca sexta (27.5%). We developed a comprehensive gene annotation workflow to obtain consensus gene models from different evidence including gene projections, protein homology, transcriptome data, and ab initio predictions. The resulting gene annotation is highly complete with 94.5% of BUSCO genes being completely present, which is higher than the BUSCO completeness of the B. mori (92.2%) and M. sexta (90%) annotations.

CONCLUSIONS

Our gene annotation strategy has general applicability to other genomes, and the H. vespertilio genome provides a valuable molecular resource to study a range of questions in this genus, including phylogeny, incomplete lineage sorting, speciation, and hybridization. A genome browser displaying the genome, alignments, and annotations is available at https://genome-public.pks.mpg.de/cgi-bin/hgTracks?db=HLhylVes1.

RepeatFiller newly identifies megabases of aligning repetitive sequences and improves annotations of conserved non-exonic elements

BACKGROUND

Transposons and other repetitive sequences make up a large part of complex genomes. Repetitive sequences can be co-opted into a variety of functions and thus provide a source for evolutionary novelty. However, comprehensively detecting ancestral repeats that align between species is difficult because considering all repeat-overlapping seeds in alignment methods that rely on the seed-and-extend heuristic results in prohibitively high runtimes.

RESULTS

Here, we show that ignoring repeat-overlapping alignment seeds when aligning entire genomes misses numerous alignments between repetitive elements. We present a tool, RepeatFiller, that improves genome alignments by incorporating previously undetected local alignments between repetitive sequences. By applying RepeatFiller to genome alignments between human and 20 other representative mammals, we uncover between 22 and 84 Mb of previously undetected alignments that mostly overlap transposable elements. We further show that the increased alignment coverage improves the annotation of conserved non-exonic elements, both by discovering numerous novel transposon-derived elements that evolve under constraint and by removing thousands of elements that are not under constraint in placental mammals.

CONCLUSIONS

RepeatFiller contributes to comprehensively aligning repetitive genomic regions, which facilitates studying transposon co-option and genome evolution. Source code: https://github.com/hillerlab/GenomeAlignmentTools.