Upgrading short read animal genome assemblies to chromosome level using comparative genomics and a universal probe set

Most recent initiatives to sequence and assemble new species’ genomes de novo fail to achieve the ultimate endpoint to produce contigs, each representing one whole chromosome. Even the best-assembled genomes (using contemporary technologies) consist of subchromosomal-sized scaffolds. To circumvent this problem, we developed a novel approach that combines computational algorithms to merge scaffolds into chromosomal fragments, PCR-based scaffold verification, and physical mapping to chromosomes. Multigenome-alignment-guided probe selection led to the development of a set of universal avian BAC clones that permit rapid anchoring of multiple scaffolds to chromosomes on all avian genomes. As proof of principle, we assembled genomes of the pigeon (Columbia livia) and peregrine falcon (Falco peregrinus) to chromosome levels comparable, in continuity, to avian reference genomes. Both species are of interest for breeding, cultural, food, and/or environmental reasons. Pigeon has a typical avian karyotype (2n = 80), while falcon (2n = 50) is highly rearranged compared to the avian ancestor. By using chromosome breakpoint data, we established that avian interchromosomal breakpoints appear in the regions of low density of conserved noncoding elements (CNEs) and that the chromosomal fission sites are further limited to long CNE “deserts.” This corresponds with fission being the rarest type of rearrangement in avian genome evolution. High-throughput multiple hybridization and rapid capture strategies using the current BAC set provide the basis for assembling numerous avian (and possibly other reptilian) species, while the overall strategy for scaffold assembly and mapping provides the basis for an approach that (provided metaphases can be generated) could be applied to any animal genome

Jurassic spark: Mapping the genomes of birds and other dinosaurs

The ultimate aim of a genome assembly is to create a contiguous length of sequence from the p- to q- terminus of each chromosome. Most assemblies are however highly fragmented, limiting their use in studies of gene mapping, phylogenomics and genomic organisation. To overcome these limitations, we developed a novel scaffold-to-chromosome anchoring method combining reference-assisted chromosome assembly (RACA) and fluorescence in situ hybridisation (FISH) to position scaffolds from de novo genomes onto chromosomes. Using RACA, scaffolds were ordered and orientated into ‘predicted chromosome fragments’ (PCFs) against a reference and outgroup genome. PCFs were verified using PCR prior to FISH mapping. A universal set of FISH probes developed through the selection of conserved regions were then used to map PCFs of peregrine falcon (Falco peregrinus Tunstall, 1771), pigeon (Columba livia Gmelin, 1789), ostrich (Struthio camelus Linnaeus, 1758), saker falcon (Falco cherrug Gray, 1834) the budgerigar (Melopsittacus undulatus Shaw, 1805). Using this approach, we were able to improve the N50 of genomes seven-fold. Results revealed that Interchromosomal breakpoint regions are limited to regions with low sequence conservation, shedding light on why most avian species have very stable karyotypes. Our combined FISH and bioinformatics approach represents a step-change in the mapping of genome assemblies, allowing comparative genomic research at a higher resolution than was previously possible. The universal probe set facilitates research into avian karyotype evolution and the role of chromosome rearrangements in adaptation and phenotypic diversity in birds. Indeed, they have been used on over 20 avian species plus non-avian reptiles (including turtles), shedding light into the evolution of dinosaur species. Non-avian dinosaurs remain subjects of intense biological enquiry while pervading popular culture and the creative arts. While organismal studies focus primarily on their morphology, relationships, likely behaviour, and ecology there have been few academic studies that have made extensive extrapolations about the nature of non-avian dinosaur genome structure prior to the emergence of modern birds. We have used multiple avian whole genome sequences assembled at a chromosomal level, to reconstruct the most likely gross genome organization of the overall genome structure of the diapsid ancestor and reconstruct the sequence of inter and intrachromosomal events that most likely occurred along the Archosauromorpha-Archosauria-Avemetatarsalia-Dinosauria-Theropoda-Maniraptora-Avialae lineage from the lepidosauromorph-archosauromorph divergence ~275 million years ago through to extant neornithine birds