Retrieving the near-complete genome of a threatened bird from wild frozen samples

Conservation genomics often relies on non-invasive methods to obtain DNA fragments which limit the power of multi-omic analyses for threatened species. Collecting samples from frozen dead animals in the wild provides an alternative approach to obtaining high-quality nucleic acids. Here, we report multi-omic analyses based on a well-preserved great bustard individual (Otis tarda, Otidiformes) of a recent death found in the mountainous region in Gansu, China. We generated a near-complete genome assembly (OTswu) having only 18 gaps scattering in 8 out of the 40 assembled chromosomes. Unlikely most other bird genome assemblies, OTswu contains complete chromosome models (2n = 80). We demonstrated that the great bustard genome likely retained the ancestral avian karyotype. We also characterized the DNA methylation landscapes of OTswu which are strongly correlated with GC content and gene expression. Our phylogenomic analysis suggested Otidiformes and Musophagiformes were sister groups that diverged from each other 46.3 million years ago. The genetic diversity of great bustard was found the lowest among the four available Otidiformes genomes, possibly due to population declines during past glacial periods. As one of the heaviest migratory birds, great bustard possesses several expanded gene families related to cardiac contraction, actin contraction, calcium ion signaling transduction, as well as positively selected genes enriching for metabolism. Finally, we identified an extremely young evolutionary stratum on the sex chromosome, a rare case among Neoaves. Together, our study combining long-read sequencing and RNA-seq technology provides a working strategy for conducting multi-omic analyses for threatened species by retrieving high-quality nucleic acids from dead animals frozen in the wild.