Teaching transposon classification as a means to crowd source the curation of repeat annotation - a tardigrade perspective.

BACKGROUND: The advancement of sequencing technologies results in the rapid release of hundreds of new genome assemblies a year providing unprecedented resources for the study of genome evolution. Within this context, the significance of in-depth analyses of repetitive elements, transposable elements (TEs) in particular, is increasingly recognized in understanding genome evolution. Despite the plethora of available bioinformatic tools for identifying and annotating TEs, the phylogenetic distance of the target species from a curated and classified database of repetitive element sequences constrains any automated annotation effort. Moreover, manual curation of raw repeat libraries is deemed essential due to the frequent incompleteness of automatically generated consensus sequences. RESULTS: Here, we present an example of a crowd-sourcing effort aimed at curating and annotating TE libraries of two non-model species built around a collaborative, peer-reviewed teaching process. Manual curation and classification are time-consuming processes that offer limited short-term academic rewards and are typically confined to a few research groups where methods are taught through hands-on experience. Crowd-sourcing efforts could therefore offer a significant opportunity to bridge the gap between learning the methods of curation effectively and empowering the scientific community with high-quality, reusable repeat libraries. CONCLUSIONS: The collaborative manual curation of TEs from two tardigrade species, for which there were no TE libraries available, resulted in the successful characterization of hundreds of new and diverse TEs in a reasonable time frame. Our crowd-sourcing setting can be used as a teaching reference guide for similar projects: A hidden treasure awaits discovery within non-model organisms.

Comparative analysis of bats and rodents' genomes suggests a relation between non-LTR retrotransposons, cancer incidence, and ageing.

The presence in nature of species showing drastic differences in lifespan and cancer incidence has recently increased the interest of the scientific community. In particular, the adaptations and the genomic features underlying the evolution of cancer-resistant and long-lived organisms have recently focused on transposable elements (TEs). In this study, we compared the content and dynamics of TE activity in the genomes of four rodent and six bat species exhibiting different lifespans and cancer susceptibility. Mouse, rat, and guinea pig genomes (short-lived and cancer-prone organisms) were compared with that of naked mole rat (Heterocephalus glaber) which is a cancer-resistant organism and the rodent with the longest lifespan. The long-lived bats of the genera Myotis, Rhinolophus, Pteropus and Rousettus were instead compared with Molossus molossus, which is one of the organisms with the shortest lifespan among the order Chiroptera. Despite previous hypotheses stating a substantial tolerance of TEs in bats, we found that long-lived bats and the naked mole rat share a marked decrease of non-LTR retrotransposons (LINEs and SINEs) accumulation in recent evolutionary times.

An annotated chromosome-scale reference genome for Eastern black-eared wheatear (Oenanthe melanoleuca).

Pervasive convergent evolution and in part high incidences of hybridization distinguish wheatears (songbirds of the genus Oenanthe) as a versatile system to address questions at the forefront of research on the molecular bases of phenotypic and species diversification. To prepare the genomic resources for this venture, we here generated and annotated a chromosome-scale assembly of the Eastern black-eared wheatear (Oenanthe melanoleuca). This species is part of the Oenanthe hispanica complex that is characterized by convergent evolution of plumage coloration and high rates of hybridization. The long-read-based male nuclear genome assembly comprises 1.04 Gb in 32 autosomes, the Z chromosome, and the mitogenome. The assembly is highly contiguous (contig N50, 12.6 Mb; scaffold N50, 70 Mb), with 96% of the genome assembled at the chromosome level and 95.5% benchmarking universal single-copy orthologs (BUSCO) completeness. The nuclear genome was annotated with 18,143 protein-coding genes and 31,333 mRNAs (annotation BUSCO completeness, 98.0%), and about 10% of the genome consists of repetitive DNA. The annotated chromosome-scale reference genome of Eastern black-eared wheatear provides a crucial resource for research into the genomics of adaptation and speciation in an intriguing group of passerines.

Satellite DNA evolution in Corvoidea inferred from short and long reads.

Satellite DNA (satDNA) is a fast-evolving portion of eukaryotic genomes. The homogeneous and repetitive nature of such satDNA causes problems during the assembly of genomes, and therefore it is still difficult to study it in detail in nonmodel organisms as well as across broad evolutionary timescales. Here, we combined the use of short- and long-read data to explore the diversity and evolution of satDNA between individuals of the same species and between genera of birds spanning ~40 millions of years of bird evolution using birds-of-paradise (Paradisaeidae) and crow (Corvus) species. These avian species highlighted the presence of a GC-rich Corvoidea satellitome composed of 61 satellite families and provided a set of candidate satDNA monomers for being centromeric on the basis of length, abundance, homogeneity and transcription. Surprisingly, we found that the satDNA of crow species rapidly diverged between closely related species while the satDNA appeared more similar between birds-of-paradise species belonging to different genera.

A beginner's guide to manual curation of transposable elements.

BACKGROUND: In the study of transposable elements (TEs), the generation of a high confidence set of consensus sequences that represent the diversity of TEs found in a given genome is a key step in the path to investigate these fascinating genomic elements. Many algorithms and pipelines are available to automatically identify putative TE families present in a genome. Despite the availability of these valuable resources, producing a library of high-quality full-length TE consensus sequences largely remains a process of manual curation. This know-how is often passed on from mentor-to-mentee within research groups, making it difficult for those outside the field to access this highly specialised skill. RESULTS: Our manuscript attempts to fill this gap by providing a set of detailed computer protocols, software recommendations and video tutorials for those aiming to manually curate TEs. Detailed step-by-step protocols, aimed at the complete beginner, are presented in the Supplementary Methods. CONCLUSIONS: The proposed set of programs and tools presented here will make the process of manual curation achievable and amenable to all researchers and in special to those new to the field of TEs.

Recurrent chromosome reshuffling and the evolution of neo-sex chromosomes in parrots.

The karyotype of most birds has remained considerably stable during more than 100 million years' evolution, except for some groups, such as parrots. The evolutionary processes and underlying genetic mechanism of chromosomal rearrangements in parrots, however, are poorly understood. Here, using chromosome-level assemblies of four parrot genomes, we uncover frequent chromosome fusions and fissions, with most of them occurring independently among lineages. The increased activities of chromosomal rearrangements in parrots are likely associated with parrot-specific loss of two genes, ALC1 and PARP3, that have known functions in the repair of double-strand breaks and maintenance of genome stability. We further find that the fusion of the ZW sex chromosomes and chromosome 11 has created a pair of neo-sex chromosomes in the ancestor of parrots, and the chromosome 25 has been further added to the sex chromosomes in monk parakeet. Together, the combination of our genomic and cytogenetic analyses characterizes the complex evolutionary history of chromosomal rearrangements and sex chromosomes in parrots.

The avian W chromosome is a refugium for endogenous retroviruses with likely effects on female-biased mutational load and genetic incompatibilities.

It is a broadly observed pattern that the non-recombining regions of sex-limited chromosomes (Y and W) accumulate more repeats than the rest of the genome, even in species like birds with a low genome-wide repeat content. Here, we show that in birds with highly heteromorphic sex chromosomes, the W chromosome has a transposable element (TE) density of greater than 55% compared to the genome-wide density of less than 10%, and contains over half of all full-length (thus potentially active) endogenous retroviruses (ERVs) of the entire genome. Using RNA-seq and protein mass spectrometry data, we were able to detect signatures of female-specific ERV expression. We hypothesize that the avian W chromosome acts as a refugium for active ERVs, probably leading to female-biased mutational load that may influence female physiology similar to the 'toxic-Y' effect in Drosophila males. Furthermore, Haldane's rule predicts that the heterogametic sex has reduced fertility in hybrids. We propose that the excess of W-linked active ERVs over the rest of the genome may be an additional explanatory variable for Haldane's rule, with consequences for genetic incompatibilities between species through TE/repressor mismatches in hybrids. Together, our results suggest that the sequence content of female-specific W chromosomes can have effects far beyond sex determination and gene dosage. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Digest: Population genomics reveals convergence toward melanism in different island populations.

Distinct traits between mainland and island populations provide an excellent opportunity to study the evolution and genetic basis of these traits. In this issue, Walsh et al. unraveled the evolution of black plumage color that arose independently in two island populations of the white-winged fairywren. They also described the first steps in understanding the genetic underpinnings of this trait.

Genetic Barriers to Historical Gene Flow between Cryptic Species of Alpine Bumblebees Revealed by Comparative Population Genomics.

Evidence is accumulating that gene flow commonly occurs between recently diverged species, despite the existence of barriers to gene flow in their genomes. However, we still know little about what regions of the genome become barriers to gene flow and how such barriers form. Here, we compare genetic differentiation across the genomes of bumblebee species living in sympatry and allopatry to reveal the potential impact of gene flow during species divergence and uncover genetic barrier loci. We first compared the genomes of the alpine bumblebee Bombus sylvicola and a previously unidentified sister species living in sympatry in the Rocky Mountains, revealing prominent islands of elevated genetic divergence in the genome that colocalize with centromeres and regions of low recombination. This same pattern is observed between the genomes of another pair of closely related species living in allopatry (B. bifarius and B. vancouverensis). Strikingly however, the genomic islands exhibit significantly elevated absolute divergence (dXY) in the sympatric, but not the allopatric, comparison indicating that they contain loci that have acted as barriers to historical gene flow in sympatry. Our results suggest that intrinsic barriers to gene flow between species may often accumulate in regions of low recombination and near centromeres through processes such as genetic hitchhiking, and that divergence in these regions is accentuated in the presence of gene flow.

Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird-of-paradise.

Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat-rich and GC-rich regions (genomic "dark matter") limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long-read, linked-read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC-rich microchromosomes and the repeat-rich W chromosome. Telomere-to-telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes.

The tuatara genome reveals ancient features of amniote evolution.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

Publisher Correction: The tuatara genome reveals ancient features of amniote evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Discovery and population genomics of structural variation in a songbird genus.

Structural variation (SV) constitutes an important type of genetic mutations providing the raw material for evolution. Here, we uncover the genome-wide spectrum of intra- and interspecific SV segregating in natural populations of seven songbird species in the genus Corvus. Combining short-read (N = 127) and long-read re-sequencing (N = 31), as well as optical mapping (N = 16), we apply both assembly- and read mapping approaches to detect SV and characterize a total of 220,452 insertions, deletions and inversions. We exploit sampling across wide phylogenetic timescales to validate SV genotypes and assess the contribution of SV to evolutionary processes in an avian model of incipient speciation. We reveal an evolutionary young (~530,000 years) cis-acting 2.25-kb LTR retrotransposon insertion reducing expression of the NDP gene with consequences for premating isolation. Our results attest to the wealth and evolutionary significance of SV segregating in natural populations and highlight the need for reliable SV genotyping.

Dynamic evolutionary history and gene content of sex chromosomes across diverse songbirds.

Songbirds have a species number close to that of mammals and are classic models for studying speciation and sexual selection. Sex chromosomes are hotspots of both processes, yet their evolutionary history in songbirds remains unclear. We characterized genomes of 11 songbird species, with 5 genomes of bird-of-paradise species. We conclude that songbird sex chromosomes have undergone four periods of recombination suppression before species radiation, producing a gradient of pairwise sequence divergence termed 'evolutionary strata'. The latest stratum was probably due to a songbird-specific burst of retrotransposon CR1-E1 elements at its boundary, instead of the chromosome inversion generally assumed for suppressing sex-linked recombination. The formation of evolutionary strata has reshaped the genomic architecture of both sex chromosomes. We find stepwise variations of Z-linked inversions, repeat and guanine-cytosine (GC) contents, as well as W-linked gene loss rate associated with the age of strata. A few W-linked genes have been preserved for their essential functions, indicated by higher and broader expression of lizard orthologues compared with those of other sex-linked genes. We also find a different degree of accelerated evolution of Z-linked genes versus autosomal genes among species, potentially reflecting diversified intensity of sexual selection. Our results uncover the dynamic evolutionary history of songbird sex chromosomes and provide insights into the mechanisms of recombination suppression.

Impact of non-LTR retrotransposons in the differentiation and evolution of anatomically modern humans.

BACKGROUND: Transposable elements are biologically important components of eukaryote genomes. In particular, non-LTR retrotransposons (N-LTRrs) played a key role in shaping the human genome throughout evolution. In this study, we compared retrotransposon insertions differentially present in the genomes of Anatomically Modern Humans, Neanderthals, Denisovans and Chimpanzees, in order to assess the possible impact of retrotransposition in the differentiation of the human lineage. RESULTS: We first identified species-specific N-LTRrs and established their distribution in present day human populations. These analyses shortlisted a group of N-LTRr insertions that were found exclusively in Anatomically Modern Humans. These insertions are associated with an increase in the number of transcriptional/splicing variants of those genes they inserted in. The analysis of the functionality of genes containing human-specific N-LTRr insertions reflects changes that occurred during human evolution. In particular, the expression of genes containing the most recent N-LTRr insertions is enriched in the brain, especially in undifferentiated neurons, and these genes associate in networks related to neuron maturation and migration. Additionally, we identified candidate N-LTRr insertions that have likely produced new functional variants exclusive to modern humans, whose genomic loci show traces of positive selection. CONCLUSIONS: Our results strongly suggest that N-LTRr impacted our differentiation as a species, most likely inducing an increase in neural complexity, and have been a constant source of genomic variability all throughout the evolution of the human lineage.

How complete are "complete" genome assemblies?-An avian perspective.

The genomics revolution has led to the sequencing of a large variety of nonmodel organisms often referred to as "whole" or "complete" genome assemblies. But how complete are these, really? Here, we use birds as an example for nonmodel vertebrates and find that, although suitable in principle for genomic studies, the current standard of short-read assemblies misses a significant proportion of the expected genome size (7% to 42%; mean 20 ± 9%). In particular, regions with strongly deviating nucleotide composition (e.g., guanine-cytosine-[GC]-rich) and regions highly enriched in repetitive DNA (e.g., transposable elements and satellite DNA) are usually underrepresented in assemblies. However, long-read sequencing technologies successfully characterize many of these underrepresented GC-rich or repeat-rich regions in several bird genomes. For instance, only ~2% of the expected total base pairs are missing in the last chicken reference (galGal5). These assemblies still contain thousands of gaps (i.e., fragmented sequences) because some chromosomal structures (e.g., centromeres) likely contain arrays of repetitive DNA that are too long to bridge with currently available technologies. We discuss how to minimize the number of assembly gaps by combining the latest available technologies with complementary strengths. At last, we emphasize the importance of knowing the location, size and potential content of assembly gaps when making population genetic inferences about adjacent genomic regions.

Correction to: Transposable Elements Activity is Positively Related to Rate of Speciation in Mammals.

The original version of the article unfortunately contained tagging error in Given and Surname of all the authors.

Transposable Elements Activity is Positively Related to Rate of Speciation in Mammals.

Transposable elements (TEs) play an essential role in shaping eukaryotic genomes and generating variability. Speciation and TE activity bursts could be strongly related in mammals, in which simple gradualistic models of differentiation do not account for the currently observed species variability. In order to test this hypothesis, we designed two parameters: the Density of insertion (DI) and the Relative rate of speciation (RRS). DI is the ratio between the number of TE insertions in a genome and its size, whereas the RRS is a conditional parameter designed to identify potential speciation bursts. Thus, by analyzing TE insertions in mammals, we defined the genomes as "hot" (high DI) and "cold" (low DI). Then, comparing TE activity among 29 taxonomical families of the whole Mammalia class, 16 intra-order pairs of mammalian species, and four superorders of Eutheria, we showed that taxa with high rates of speciation are associated with "hot" genomes, whereas taxa with low ones are associated with "cold" genomes. These results suggest a remarkable correlation between TE activity and speciation, also being consistent with patterns describing variable rates of differentiation and accounting for the different time frames of the speciation bursts.