Low diversity, activity, and density of transposable elements in five avian genomes

The genome of a globally invasive passerine, the common myna, Acridotheres tristis

In an era of global climate change, biodiversity conservation is receiving increased attention. Conservation efforts are greatly aided by genetic tools and approaches, which seek to understand patterns of genetic diversity and how they impact species health and their ability to persist under future climate regimes. Invasive species offer vital model systems in which to investigate questions regarding adaptive potential, with a particular focus on how changes in genetic diversity and effective population size interact with novel selection regimes. The common myna (Acridotheres tristis) is a globally invasive passerine and is an excellent model species for research both into the persistence of low-diversity populations and the mechanisms of biological invasion. To underpin research on the invasion genetics of this species, we present the genome assembly of the common myna. We describe the genomic landscape of this species, including genome wide allelic diversity, methylation, repeats, and recombination rate, as well as an examination of gene family evolution. Finally, we use demographic analysis to identify that some native regions underwent a dramatic population increase between the two most recent periods of glaciation, and reveal artefactual impacts of genetic bottlenecks on demographic analysis.

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.

Taming transposable elements in livestock and poultry: a review of their roles and applications

Livestock and poultry play a significant role in human nutrition by converting agricultural by-products into high-quality proteins. To meet the growing demand for safe animal protein, genetic improvement of livestock must be done sustainably while minimizing negative environmental impacts. Transposable elements (TE) are important components of livestock and poultry genomes, contributing to their genetic diversity, chromatin states, gene regulatory networks, and complex traits of economic value. However, compared to other species, research on TE in livestock and poultry is still in its early stages. In this review, we analyze 72 studies published in the past 20 years, summarize the TE composition in livestock and poultry genomes, and focus on their potential roles in functional genomics. We also discuss bioinformatic tools and strategies for integrating multi-omics data with TE, and explore future directions, feasibility, and challenges of TE research in livestock and poultry. In addition, we suggest strategies to apply TE in basic biological research and animal breeding. Our goal is to provide a new perspective on the importance of TE in livestock and poultry genomes.

Transposable elements in the transcriptome of the velvetbean caterpillar Anticarsia gemmatalis Hübner, 1818 (Lepidoptera: Erebidae)

Transposable elements (TEs) are DNA sequences that possess the ability to move from one genomic location to another. These sequences contribute to a significant fraction of the genomes of most eukaryotes and can impact their architecture and regulation. In this paper, we present the first data related to the identification and characterization of TEs present in the transcriptome of Anticarsia gemmatalis. Approximately, 835 transcripts showed significant similarity to TEs and (or) characteristic domains. Retrotransposons accounted for 71.2% (595 sequences) of the identified elements, while DNA transposons were less abundant, with 240 annotations (28.8%). TEs were classified into 30 superfamilies, with SINE3/5S and Gypsy being the most abundant. Based on the sequences of TEs found in the transcriptome, we were able to locate conserved regions in the chromosomes of this species. The analysis of differential expression of TEs in susceptible and resistant strains, challenged and not challenged with Bacillus thuringiensis ( Bt) from in silico analysis, indicated that exposure to Bt can regulate the transcription of mobile genetic elements in the velvetbean caterpillar. Thus, these data contribute significantly to the knowledge of the structure and composition of these elements in the genome of this species, and suggest the role of stress on their expression.

Dinosaurs: Comparative Cytogenomics of Their Reptile Cousins and Avian Descendants

Reptiles known as dinosaurs pervade scientific and popular culture, while interest in their genomics has increased since the 1990s. Birds (part of the crown group Reptilia) are living theropod dinosaurs. Chromosome-level genome assemblies cannot be made from long-extinct biological material, but dinosaur genome organization can be inferred through comparative genomics of related extant species. Most reptiles apart from crocodilians have both macro- and microchromosomes; comparative genomics involving molecular cytogenetics and bioinformatics has established chromosomal relationships between many species. The capacity of dinosaurs to survive multiple extinction events is now well established, and birds now have more species in comparison with any other terrestrial vertebrate. This may be due, in part, to their karyotypic features, including a distinctive karyotype of around n = 40 (~10 macro and 30 microchromosomes). Similarity in genome organization in distantly related species suggests that the common avian ancestor had a similar karyotype to e.g., the chicken/emu/zebra finch. The close karyotypic similarity to the soft-shelled turtle (n = 33) suggests that this basic pattern was mostly established before the Testudine-Archosaur divergence, ~255 MYA. That is, dinosaurs most likely had similar karyotypes and their extensive phenotypic variation may have been mediated by increased random chromosome segregation and genetic recombination, which is inherently higher in karyotypes with more and smaller chromosomes.

Ineffective silencing of transposable elements on an avian W Chromosome

One of the defining features of transposable elements (TEs) is their ability to move to new locations in the host genome. To minimise the potentially deleterious effects of de novo TE insertions, hosts have evolved several mechanisms to control TE activity, including recombination-mediated removal and epigenetic silencing; however, increasing evidence suggests that silencing of TEs is often incomplete. The crow family experienced a recent radiation of LTR retrotransposons (LTRs), offering an opportunity to gain insight into the regulatory control of young, potentially still active TEs. We quantified the abundance of TE-derived transcripts across several tissues in 15 Eurasian crows (Corvus (corone) spp.) raised under common garden conditions and find evidence for ineffective TE suppression on the female-specific W Chromosome. Using RNA-seq data, we show that ~ 9.5% of all transcribed TEs had considerably greater (average: 16-fold) transcript abundance in female crows, and that more than 85% of these female-biased TEs originated on the W Chromosome. After accounting for differences in TE density among chromosomal classes, W-linked TEs were significantly more highly expressed than TEs residing on other chromosomes, consistent with ineffective silencing on the former. Together, our results suggest that the crow W Chromosome acts as a source of transcriptionally active TEs, with possible negative fitness consequences for female birds analogous to Drosophila (an X/Y system), where overexpression of Y-linked TEs is associated with male-specific aging and fitness loss ('toxic Y').

Comparative Population Genomics of Cryptic Speciation and Adaptive Divergence in Bicknell's and Gray-Cheeked Thrushes (Aves: Catharus bicknelli and Catharus minimus)

Cryptic speciation may occur when reproductive isolation is recent or the accumulation of morphological differences between sister lineages is slowed by stabilizing selection preventing phenotypic differentiation. In North America, Bicknell's Thrush (Catharus bicknelli) and its sister species, the Gray-cheeked Thrush (Catharus minimus), are parapatrically breeding migratory songbirds, distinguishable in nature only by subtle differences in song and coloration, and were recognized as distinct species only in the 1990s. Previous molecular studies have estimated that the species diverged approximately 120,000-420,000 YBP and found very low levels of introgression despite their similarity and sympatry in the spring (prebreeding) migration. To further clarify the history, genetic divergence, genomic structure, and adaptive processes in C. bicknelli and C. minimus, we sequenced and assembled high-coverage reference genomes of both species and resequenced genomes from population samples of C. bicknelli, C. minimus, and two individuals of the Swainson's Thrush (Catharus ustulatus). The genome of C. bicknelli exhibits markedly higher abundances of transposable elements compared with other Catharus and chicken. Demographic and admixture analyses confirm moderate genome-wide differentiation (Fst ≈ 0.10) and limited gene flow between C. bicknelli and C. minimus, but suggest a more recent divergence than estimates based on mtDNA. We find evidence of rapid evolution of the Z-chromosome and elevated divergence consistent with natural selection on genomic regions near genes involved with neuronal processes in C. bicknelli. These genomes are a useful resource for future investigations of speciation, migration, and adaptation in Catharus thrushes.

Characterization and expression pattern of ZB and PS transposons in zebrafish

Despite comprising much of the genome, transposons were once thought of as junk. However, transposons play many roles in the eukaryotic genome, such as providing new proteins as domesticated genes, expressing during germline-soma differentiation, function in DNA rearrangement in the offspring, and so on. We sought to describe the distribution and structural organization of the two autonomous transposons (ZB and PS) in the zebrafish genome and examine their expression patterns in embryos and adult tissues. The intact copy of ZB and PS was queried by BLAST on NCBI and ENSEMBL using default parameters. Of the copies with coverage and identity, more than 90 % were downloaded to do structural analysis. Spatial and temporal expression patterns were detected by qRT-PCR and Whole-mount in situ hybridization (WISH). There are 19 intact copies of ZB, encoding 341 amino acid residues with DD34E catalytic domain and flanked by 201bp TIRs, and seven intact PS copies, containing 425 amino acid residues with DD35D catalytic domain flanked by 28bp TIRs, were detected in the genome of zebrafish respectively. Analysis of genomic insertions indicated that both ZB and PS transposons are prone to be retained in the intron and intergenic regions of the zebrafish genome. The sense and antisense transcripts of ZB and PS were detected during embryonic development stages and exhibited similar expression patterns. The difference is that the sense strand transcript of ZB was explicitly expressed in midbrain-hindbrain boundary (MHB) and otic vesicle (OV), and pharyngeal arches and pharyngeal pouches (PA&PP) at 48 hpf. In adult zebrafish, the expressions of ZB and PS in muscle and brain are much higher than in other tissues. Our study results indicate that ZB and PS transposons may be involved in the embryonic development and regulation of somatic cells of certain adult tissues, such as the brain and muscle.

Multiple Invasions of Visitor, a DD41D Family of Tc1/mariner Transposons, throughout the Evolution of Vertebrates

Although the DD41D (named as Visitor, VS) family of Tc1/mariner transposons was discovered in Arthropods and Mollusca, the evolution profile of this family is still largely unknown. We found that VS is widespread in the animal kingdom, including 140 species of 18 orders in invertebrates and 30 species of 12 orders in vertebrates, and one land plant species. Our data revealed multiple horizontal transfer events in both invertebrates and vertebrates and invasion into multiple lineages of mammals, including Chiroptera (seven species), Dasyuromorphia/Marsupialia (one species), Didelphimorphia/Marsupialia (one species), Diprotodontia/Marsupialia (two species), and Primates (one species). Phylogenetic analysis revealed a close relationship of VSs to DD37D/maT and DD34D/mariner and confirmed that VSs with the DD40D signature identified previously are not a distinct family but originated from DD41D/VS. Age analysis revealed that the most recent invasion of VSs was found in ray-finned fishes and a toad, followed by relatively young invasions in bats and marsupials, whereas VSs in mammals, jawless fishes, and lizards were mainly represented by ancient copies, suggesting old age. Phylogenetic analyses and comparison of pairwise distances between VSs and recombination-activating gene 1 (RAG1) support horizontal transfer events of VSs in vertebrates. The intact VSs from bats were nonfunctional as determined by the transposition activity assay. Some vertebrate lineages and species were identified as the hot hosts of Tc1/mariner transposons. Overall, our study presents the evolution profile of VSs and suggests that VSs play roles in diversifying and shaping the genomes of diverse animal lineages.

Intruder (DD38E), a recently evolved sibling family of DD34E/Tc1 transposons in animals

BACKGROUND

A family of Tc1/mariner transposons with a characteristic DD38E triad of catalytic amino acid residues, named Intruder (IT), was previously discovered in sturgeon genomes, but their evolutionary landscapes remain largely unknown.

RESULTS

Here, we comprehensively investigated the evolutionary profiles of ITs, and evaluated their cut-and-paste activities in cells. ITs exhibited a narrow taxonomic distribution pattern in the animal kingdom, with invasions into two invertebrate phyla (Arthropoda and Cnidaria) and three vertebrate lineages (Actinopterygii, Agnatha, and Anura): very similar to that of the DD36E/IC family. Some animal orders and species seem to be more hospitable to Tc1/mariner transposons, one order of Amphibia and seven Actinopterygian orders are the most common orders with horizontal transfer events and have been invaded by all four families (DD38E/IT, DD35E/TR, DD36E/IC and DD37E/TRT) of Tc1/mariner transposons, and eight Actinopterygii species were identified as the major hosts of these families. Intact ITs have a total length of 1.5-1.7 kb containing a transposase gene flanked by terminal inverted repeats (TIRs). The phylogenetic tree and sequence identity showed that IT transposases were most closely related to DD34E/Tc1. ITs have been involved in multiple events of horizontal transfer in vertebrates and have invaded most lineages recently (< 5 million years ago) based on insertion age analysis. Accordingly, ITs presented high average sequence identity (86-95%) across most vertebrate species, suggesting that some are putatively active. ITs can transpose in human HeLa cells, and the transposition efficiency of consensus TIRs was higher than that of the TIRs of natural isolates.

CONCLUSIONS

We conclude that DD38E/IT originated from DD34E/Tc1 and can be detected in two invertebrate phyla (Arthropoda and Cnidaria), and in three vertebrate lineages (Actinopterygii, Agnatha and Anura). IT has experienced multiple HT events in animals, dominated by recent amplifications in most species and has high identity among vertebrate taxa. Our reconstructed IT transposon vector designed according to the sequence from the "cat" genome showed high cut-and-paste activity. The data suggest that IT has been acquired recently and is active in many species. This study is meaningful for understanding the evolution of the Tc1/mariner superfamily members and their hosts.

The chicken model organism for epigenomic research

The chicken model organism has advanced the areas of developmental biology, virology, immunology, oncology, epigenetic regulation of gene expression, conservation biology, and genomics of domestication. Further, the chicken model organism has aided in our understanding of human disease. Through the recent advances in high-throughput sequencing and bioinformatic tools, researchers have successfully identified sequences in the chicken genome that have human orthologs, improving mammalian genome annotation. In this review, we highlight the importance of chicken as an animal model in basic and pre-clinical research. We will present the importance of chicken in poultry epigenetics and in genomic studies that trace back to their ancestor, the last link between human and chicken in the tree of life. There are still many genes of unknown function in the chicken genome yet to be characterized. By taking advantage of recent sequencing technologies, it is possible to gain further insight into the chicken epigenome.

A Highly Contiguous Genome for the Golden-Fronted Woodpecker (Melanerpes aurifrons) via Hybrid Oxford Nanopore and Short Read Assembly

Woodpeckers are found in nearly every part of the world and have been important for studies of biogeography, phylogeography, and macroecology. Woodpecker hybrid zones are often studied to understand the dynamics of introgression between bird species. Notably, woodpeckers are gaining attention for their enriched levels of transposable elements (TEs) relative to most other birds. This enrichment of TEs may have substantial effects on molecular evolution. However, comparative studies of woodpecker genomes are hindered by the fact that no high-contiguity genome exists for any woodpecker species. Using hybrid assembly methods combining long-read Oxford Nanopore and short-read Illumina sequencing data, we generated a highly contiguous genome assembly for the Golden-fronted Woodpecker (Melanerpes aurifrons). The final assembly is 1.31 Gb and comprises 441 contigs plus a full mitochondrial genome. Half of the assembly is represented by 28 contigs (contig L50), each of these contigs is at least 16 Mb in size (contig N50). High recovery (92.6%) of bird-specific BUSCO genes suggests our assembly is both relatively complete and relatively accurate. Over a quarter (25.8%) of the genome consists of repetitive elements, with 287 Mb (21.9%) of those elements assignable to the CR1 superfamily of transposable elements, the highest proportion of CR1 repeats reported for any bird genome to date. Our assembly should improve comparative studies of molecular evolution and genomics in woodpeckers and allies. Additionally, the sequencing and bioinformatic resources used to generate this assembly were relatively low-cost and should provide a direction for development of high-quality genomes for studies of animal biodiversity.

Alu insertion polymorphisms shared by Papio baboons and Theropithecus gelada reveal an intertwined common ancestry

Background

Baboons (genus Papio) and geladas (Theropithecus gelada) are now generally recognized as close phylogenetic relatives, though morphologically quite distinct and generally classified in separate genera. Primate specific Alu retrotransposons are well-established genomic markers for the study of phylogenetic and population genetic relationships. We previously reported a computational reconstruction of Papio phylogeny using large-scale whole genome sequence (WGS) analysis of Alu insertion polymorphisms. Recently, high coverage WGS was generated for Theropithecus gelada. The objective of this study was to apply the high-throughput "poly-Detect" method to computationally determine the number of Alu insertion polymorphisms shared by T. gelada and Papio, and vice versa, by each individual Papio species and T. gelada. Secondly, we performed locus-specific polymerase chain reaction (PCR) assays on a diverse DNA panel to complement the computational data.

Results

We identified 27,700 Alu insertions from T. gelada WGS that were also present among six Papio species, with nearly half (12,956) remaining unfixed among 12 Papio individuals. Similarly, each of the six Papio species had species-indicative Alu insertions that were also present in T. gelada. In general, P. kindae shared more insertion polymorphisms with T. gelada than did any of the other five Papio species. PCR-based genotype data provided additional support for the computational findings.

Conclusions

Our discovery that several thousand Alu insertion polymorphisms are shared by T. gelada and Papio baboons suggests a much more permeable reproductive barrier between the two genera then previously suspected. Their intertwined evolution likely involves a long history of admixture, gene flow and incomplete lineage sorting.

Germline-restricted chromosome (GRC) is widespread among songbirds

An unusual supernumerary chromosome has been reported for two related avian species, the zebra and Bengalese finches. This large, germline-restricted chromosome (GRC) is eliminated from somatic cells and spermatids and transmitted via oocytes only. Its origin, distribution among avian lineages, and function were mostly unknown so far. Using immunolocalization of key meiotic proteins, we found that GRCs of varying size and genetic content are present in all 16 songbird species investigated and absent from germline genomes of all eight examined bird species from other avian orders. Results of fluorescent in situ hybridization of microdissected GRC probes and their sequencing indicate that GRCs show little homology between songbird species and contain a variety of repetitive elements and unique sequences with paralogs in the somatic genome. Our data suggest that the GRC evolved in the common ancestor of all songbirds and underwent significant changes in the extant descendant lineages.

Time lapse: A glimpse into prehistoric genomics

For the purpose of this review, ‘time-lapse’ refers to the reconstruction of ancestral (in this case dinosaur) karyotypes using genome assemblies of extant species. Such reconstructions are only usually possible when genomes are assembled to ‘chromosome level’ i.e. a complete representation of all the sequences, correctly ordered contiguously on each of the chromosomes. Recent paleontological evidence is very clear that birds are living dinosaurs, the latest example of dinosaurs emerging from a catastrophic extinction event. Non-avian dinosaurs (ever present in the public imagination through art, and broadcast media) emerged some 240 million years ago and have displayed incredible phenotypic diversity. Here we report on our recent studies to infer the overall karyotype of the Theropod dinosaur lineage from extant avian chromosome level genome assemblies. Our work first focused on determining the likely karyotype of the avian ancestor (most likely a chicken-sized, two-legged, feathered, land dinosaur from the Jurassic period) finding karyotypic similarity to the chicken. We then took the work further to determine the likely karyotype of the bird-lizard ancestor and the chromosomal changes (chiefly translocations and inversions) that occurred between then and modern birds. A combination of bioinformatics and cross-species fluorescence in situ hybridization (zoo-FISH) uncovered a considerable number of translocations and fissions from a ‘lizard-like’ genome structure of 2n = 36–46 to one similar to that of soft-shelled turtles (2n = 66) from 275 to 255 million years ago (mya). Remarkable karyotypic similarities between some soft-shelled turtles and chicken suggests that there were few translocations from the bird-turtle ancestor (plus ∼7 fissions) through the dawn of the dinosaurs and pterosaurs, through the theropod linage and on to most to modern birds. In other words, an avian-like karyotype was in place about 240mya when the dinosaurs and pterosaurs first emerged. We mapped 49 chromosome inversions from then to the present day, uncovering some gene ontology enrichment in evolutionary breakpoint regions. This avian-like karyotype with its many (micro)chromosomes provides the basis for variation (the driver of natural selection) through increased random segregation and recombination. It may therefore contribute to the ability of dinosaurs to survive multiple extinction events, emerging each time as speciose and diverse.

Reconstruction of the diapsid ancestral genome permits chromosome evolution tracing in avian and non-avian dinosaurs

Genomic organisation of extinct lineages can be inferred from extant chromosome-level genome assemblies. Here, we apply bioinformatic and molecular cytogenetic approaches to determine the genomic structure of the diapsid common ancestor. We then infer the events that likely occurred along this lineage from theropod dinosaurs through to modern birds. Our results suggest that most elements of a typical ‘avian-like’ karyotype (40 chromosome pairs, including 30 microchromosomes) were in place before the divergence of turtles from birds ~255 mya. This genome organisation therefore predates the emergence of early dinosaurs and pterosaurs and the evolution of flight. Remaining largely unchanged interchromosomally through the dinosaur–theropod route that led to modern birds, intrachromosomal changes nonetheless reveal evolutionary breakpoint regions enriched for genes with ontology terms related to chromatin organisation and transcription. This genomic structure therefore appears highly stable yet contributes to a large degree of phenotypic diversity, as well as underpinning adaptive responses to major environmental disruptions via intrachromosomal repatterning.

Ancient diapsids diverged into the lineages leading to turtles and birds over 250 million years ago. Here, the authors use genomic and molecular cytogenetic analyses of modern species to infer the genome structure of the diapsid common ancestor (DCA) and the changes occurring along the lineage to birds through theropod dinosaurs.