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.

Because of their near-global distribution, woodpeckers are often used as model systems in biogeographic and phylogeographic studies (Zink ; Fuchs , 2007; Klicka ; Navarro-Sigüenza ), as well as macroecological studies (Bock and Bock 1974; Blackburn ; Husak and Husak 2003; Ilsøe ). In North America, four woodpecker hybrid zones have also been studied for insights into avian speciation. These include flickers (Giller 1959; Smith 1987; Moore 1987; Moore and Price 1993; Seneviratne , 2016), sapsuckers (Sphyrapicus; Giller 1959; Smith 1987; Moore 1987; Moore and Price 1993; Seneviratne , 2016), Nuttall’s/Ladder-back woodpeckers (Manthey ) and Red-bellied/Golden-fronted woodpeckers (Giller 1959; Smith 1987; Moore 1987; Moore and Price 1993; Seneviratne , 2016). More recently, woodpeckers have gained attention for the high amount of repetitive DNA found in their genomes relative to other bird taxa (Sotero-Caio ; de Oliveira ; Bertocchi ), the result of high levels of genome-wide transposable elements (TEs), which are scarce in most bird genomes (Gao ). Manthey surveyed several woodpecker genomes, and found that TEs make up 17–31% of woodpecker genomes, compared to <10% for other bird species. Increasingly, researchers have suggested that TEs may play a critical role in driving avian evolution (Gao ; Suh ). However, only two woodpecker genomes (Picoides [Dryocopus] pubescens, Zhang ; and Dendrocopos noguchii, GCA_004320165.1) have been published, limiting our ability to undertake comparative genomics and to understand the architecture of TE evolution in woodpeckers.

The genus Melanperpes represents the largest radiation of New World woodpeckers (Aves: Picidae; Navarro-Sigüenza ). Melanerpes woodpeckers are found almost everywhere where forest occurs in the Americas. Various species range continuously from southern Canada to Argentina, with three species occurring in the West Indies. Members of the Red-bellied/Golden-fronted Woodpecker species complex (Melanerpes aurifrons/carolinus) are notable for the discord between plumage variation and phylogenetic structure, especially in Mexico and northern Mesoamerica, where various races show considerable plumage variation despite a lack of phylogeographic variation (García-Trejo ; Navarro-Sigüenza ; but see Barrowclough et al. (2018).

Here, we provide a high-quality, highly contiguous, reference genome for the Golden-fronted Woodpecker (Melanerpes aurifrons; Figure 1), which we expect to be foundational for research in avian hybrid zone dynamics, the evolution of reproductive isolation, and the role of TEs and other repetitive elements in driving molecular evolution and speciation in woodpeckers and other birds. Additionally, we believe our workflow will be of general interest to researchers looking to develop high quality reference genomes for non-model birds and other vertebrates. Our genome assembly was accomplished with only ∼46X coverage of long-read Oxford Nanopore and ∼52X coverage of short-read Illumina sequence. Furthermore, after troubleshooting the various steps in our workflow, we were able to go from raw data to completed assembly in less than two weeks, using a compute cluster with 72 CPU cores, 384 GB of RAM, and two NVIDIA Tesla P100 GPU accelerators. These sequencing and computational requirements should be fairly accessible to most research groups with modest budgets.

Figure 1

Male (left) and female (right) Golden-fronted Woodpecker (Melanerpes aurifrons). Photos by Bettina Arrigoni, cropped, and used under CC BY 2.0 license. Original photos available at: https://flickr.com/photos/69683857@N05/39849351035 and https://flickr.com/photos/69683857@N05/26752528708.

Methods

Specimen collection and DNA extraction

We collected an adult female Golden-fronted Woodpecker (Melanerpes aurifrons) in March 2019 in Dickens, Texas. This specimen and all associated genomic resources have been archived as a museum voucher at the Sam Noble Oklahoma Museum of Natural History (specimen number 24340, tissue number SMB682). This specimen is registered with NCBI as Biosample SAMN13719207. Scientific collecting was done under the following permits: Texas Scientific Collecting Permit SPR-0916-229, US Migratory Bird Treaty Act permit MB02276C-0, and University of Oklahoma IACUC permit R15-016A. Approximately ten milliliters of whole blood were transferred immediately to blood tubes coated with EDTA as an anticoagulant and kept cool on wet ice for 6 hr. Upon return to the lab, blood was aliquoted to several tubes and stored at -20°. We extracted high molecular weight DNA from a single aliquot using the Nanobind CBB Big DNA Kit (Circulomics, Inc). We followed the manufacturer’s high molecular weight (50–300+ kb) DNA protocol for cultured mammalian cells, with the exception that we diluted 10 uL of whole woodpecker blood into 190 uL of fresh PBS as sample input. This was done per recommendation of the manufacturer (personal communication) given that avian red blood cells are nucleated. We performed three rounds of HMW DNA extraction, with each extraction being used as input for an Oxford Nanopore sequencing library. For the initial MinION nanopore sequencing run, we treated the HMW DNA with Short Read Eliminator Kit (Circulomics, Inc.), following the manufacturer’s instructions. For the PromethION sequencing runs, this kit was not used in order to maximize DNA concentration for downstream library construction. For Illumina sequencing, we performed a standard genomic DNA extraction using the Quick-DNA miniprep kit (Zymo Research). Again, 10 uL of blood was diluted into 190 uL of PBS before following the manufacturer’s protocol.

Whole genome library preparation and sequencing

We first generated Oxford Nanopore long reads in-house on the MinION device using the SQK-LSK109 library preparation kit followed by sequencing on a R9.4.1 flow cell per the manufacturer’s instructions (Oxford Nanopore Technologies). After recovering less than 3 Gb of sequence data, we sent our remaining two HMW extracts to the UC Davis Core Lab to run on two PromethION flow cells. For the first PromethION run, a single library was prepared, loaded, and run for 48 hr. In an attempt to recover more sequence data, for the second PromethION flow cell, the HMW DNA was sheared using a Megaruptor (Diagenode, Inc.) set to a 50 kb target prior to a double library preparation allowing for a nuclease flush and a fresh library reload at the 24 hr mark of a 48 hr run. To generate the Illumina sequencing library, we used the NEBNext Ultra II FS DNA library kit (New England Biolabs, Inc); the initial enzymatic shearing step was accomplished via 10 min of incubation at 37°, after which we followed the manufacturer’s instructions. The library was indexed using NEBNext Multiplex Oligos for Illumina unique dual index kit (New England Biolabs, Inc). This library was run on one lane of a NovaSeq S-Prime flow cell at the Oklahoma Medical Research Foundation Clinical Genomics Center.

Read QC and trimming

Raw Nanopore reads (fast5) were converted into fastq format using Oxford Nanopore’s proprietary base-calling software guppy v.3.4 (https://community.nanoporetech.com). We evaluated Nanopore read quality in NanoPlot, and filtered reads (Q > 7, length > 10,000) using NanoFilt; both are part of the NanoPack distribution (De Coster ). We filtered and trimmed Illumina reads to remove adapters and low-quality bases using the standard settings in the bbduk program, which is part of BBMap v38.00 (http://sourceforge.net/projects/bbmap).

Genome size estimation

We used jellyfish v2.2.3 (Marçais and Kingsford 2011) to count the frequency of three distinct k-mers (17-mers, 26-mers, and 31-mers) in our trimmed Illumina sequencing reads. We plotted the resulting histogram file to establish the peak coverage at each k-mer. We then used the following formula from Liu to estimate genome size:

where G equals genome size, N equals the total number of base pairs in all reads, D equals the expected k-mer coverage depth, L equals the average read length, and K equals the k-mer size.

Genome assembly

We assembled an initial draft assembly using flye v2.6 (Lin ; Kolmogorov ). We corrected errors in that assembly using three iterations of racon v1.4.7 (Vaser ), which uses the original nanopore long reads to correct the assembly consensus output. Subsequently, we ran two iterations of short-read polishing in pilon v1.23 (Walker ). Pilon uses read alignment analysis to identify inconsistencies between a draft genome and an alignment of reads to that genome and attempts to correct base call errors, small and large indels and block substitutions, and to identify local mis-assemblies. We used minimap2 v2.17 (Li 2018) to map the Illumina short reads to our assembly to generate input alignments for pilon polishing. Finally, we removed redundant haploid and low coverage contigs from the polished assembly using purge_haplotigs v1.1 (Roach ). That workflow begins by generating a read depth histogram to establish coverage cutoff levels for low coverage, high coverage, and the midpoint in coverage between haploid and diploid coverage. Subsequent steps remove contigs if they fall outside this range. To prevent interference from repeat elements, a bed file containing repetitive element coordinates was provided to purge_haplotigs (see below). At each step, we evaluated the completeness of the genome assembly via BUSCO v3 benchmarking (Simão ) using the aves_odb9 dataset, which contains 4915 near-universal single-copy orthologs.

Identifying repetitive elements and annotating the mitochondrial genome

We estimated the portion of our genome comprised of repetitive elements using RepeatMasker v4.0.9p2. (http://www.repeatmasker.org/). To ensure avian-specific annotations, we utilized the ‘chicken’ species entry in the RepeatMasker database. The output of this was converted to a bed file for use in purge_haplotigs (see above). Using the Zebra Finch (Taeniopygia guttata), rather than chicken, resulted in less than 0.1% difference in the estimated repetitive elements frequencies reported by RepeatMasker.

To identify contigs corresponding to the mitochondrial genome, we performed BLASTn using an existing sequence of M. aurifrons NADH-dehydrogenase subunit II (MF766655) against our draft assembly, which recovered a single contig. Using that contig as a draft mitochondrial genome, we annotated all mitochondrial genes using the MitoAnnotator web server (http://mitofish.aori.u-tokyo.ac.jp/; Iwasaki ). As an additional metric to evaluate the robustness of our assembly, we assessed the quality of the mtDNA assembly by comparing all protein coding gene sequences to an annotated mitochondrial genome of the Great Spotted Woodpecker (Dendrocopos major, NC028174.1) and we confirmed proper secondary structure for all tRNAs using the tRNAscan-SE 2.0 web server (http://lowelab.ucsc.edu/tRNAscan-SE/index.html; Lowe and Chan 2016). Finally, we performed a Google Scholar search to look for other avian genomes published in 2019 to compare our results to other high-quality wild bird reference genomes.

Reference assisted pseudomolecule scaffolding

We generated pseudomolecule chromosomal scaffolds using a reference genome assisted approach in RaGOO v1.1 (Alonge ). RaGOO attempts to cluster, order, and orient assembly contigs based on a Minimap2 (Li 2018) alignment of those contigs to a reference genome. There is no chromosomal-level genome available for woodpeckers and allies (Order: Piciformes) and all bird species with chromosomal-level genomes are roughly similar in phylogenetic distance (Zhang ). Therefore, we chose to use the Anna’s Hummingbird chromosome-assembled genome (Calypte anna bCalAnn1_v1.p; GCA_003957555.2; Korlach ) as our reference. To evaluate the validity of our reference-guided scaffold we measured genome synteny by aligning the RaGOO-scaffolded assembly to the Anna’s Hummingbird reference genome as well as to two additional chromosomal-scale bird genomes: Zebra Finch (Taeniopygia guttata, GCA_008822105.1) and Kakapo (Strigops habroptila, GCA_004027225.1). We used the nucmer module in MUMMER v4.0.0b2 (Kurtz ) to perform the alignments which were subsequently filtered using MUMMER’s delta_filter module with many-to-many alignments allowed, minimum alignment length equal to 500, and minimum alignment identity equal to 80%. We filtered this output to only select mapped clusters equaling at least 3000 base pairs in the RaGOO (query) assembly. To visualize genome synteny we generated circos plots in Circa (http://omgenomics.com/circa); to improve visual clarity, only chromosomes 1 through 14 (including 4A, 4B and 5A) and the Z and W chromosomes were used.

Data availability

The genome assembly and raw reads have been deposited to NCBI (project PRJNA598863; WWNC00000000). Output from RaGOO-scaffolded assembly is available from figshare. Code to replicate our analyses is available at: https://github.com/wileygenomics/Nanopore_Avian_Genome_Assembly. Supplemental material available at figshare: https://doi.org/10.25387/g3.11952324.

Results and Discussion

Sequencing run statistics and genome assembly results

We generated a total of 3.9 M Oxford Nanopore reads for a total of 63.99 Gb of long read sequence data (Table 1). The initial MinION sequencing run produced only limited data but had longer reads than either the unsheared or sheared HMW runs on the PromethION (Table 1). Megaruptor shearing resulted in slightly lower median read length but increased read length N50. Comparing our two PromethION runs, HMW shearing, along with reloading a fresh library after the initial 24-hour runtime, increased the number of reads recovered by 18% and total base pairs yielded by 40%. These results, along with the slightly higher read length N50, suggest that users should consider HMW DNA shearing and a nuclease flush prior to Oxford Nanopore library preparation for metazoan whole genome assemblies.

Table 1

Run statistics for three Oxford Nanopore sequencing runs

Sequencing Run	Reads	Median Read Length	Read Length N50	Median Read Qual	Total Bases
MinION Run	101,989	15,105	39,203	10.5	2.29 X 109
PromethION Run 1	1.77 X 106	9,478	30.742	10.6	27.43 X 109
PromethION Run 2	2.09 X 106	9,178	34,314	10.5	34.27 X 109

Our initial flye assembly resulted in 1519 contigs which were organized into 1456 scaffolds with 41x mean coverage. The total length of the assembly was 1.35 Gb with a contig L50 of 14.5 Mb and a scaffold L50 of 15.9 Mb. The largest contig was 4.28 Mb. Fifty percent of the assembly was recovered in just 33 contigs (L50) – those contigs were at least 14.5 Mb in size (N50). Ninety percent of the assembly was recovered in 136 contigs (L90) all of which were at least 1.4 Mb in size (N90). In contrast, our final genome assembly after polishing and haplotig purging (see below) resulted in 441 unscaffolded contigs with an estimated genome size of 1.309 Gb, plus a complete circularized mitochondrial genome (N = 16,844). The largest contig was 46.8 Mb. Contig L50 was 28 contigs and contig N50 was 16 Mb.

Our Illumina sequencing generated a total of 490.7 M paired end reads for a total of 74.1 Gb of short read sequence data. After bbduk trimming, we retained 71.2 Gb of Illumina data from 487.8 M reads. Histograms for k-mer coverage (Supplementary Figure S1) suggest an average k-mer depth ranging from 46 (k-mer = 17 bp) to 41 (k-mer = 31 bp). This results in an estimated genome size of 1.38 – 1.40 Gb, depending on the k-mer value (Table 2). For resulting analyses, we consider 1.38 Gb to be the best estimate of genome size. This suggests that our assembly recovered 94.2% of the true Golden-fronted Woodpecker genome, whereas a recent review suggests that on-average, bird genomes generated with short reads alone often recover only 71–89% of the full genome (Peona ). Given our genome size estimate, our Oxford Nanopore reads were at an average depth of 43.4, while our trimmed Illumina reads were at an average sequencing depth of 51.6.

Table 2

Results of genome size estimate using k-mer analysis

k-mer size	k-mer coverage	Estimated genome size
17	46	1.377 Gb
26	42	1.404 Gb
31	41	1.379 Gb

Racon correction eliminates gaps and merges small contigs

Consensus correction using racon eliminated gaps in scaffolds, meaning that all recovered contigs were ungapped. Racon greatly reduced the number of contigs, but this was apparently mostly achieved by merging small contigs with overlapping regions from larger ones. Thus, racon had the effect of slightly reducing assembly size with only slight increases in contig L50, L90, N50, and N90 (Table 3). The largest effect of racon correction was observed in the first iteration, by the third round, little change was observed.

Table 3

Assembly contiguity statistics for the initial assembly and after multiple iterations of consensus correction

Stage	# of contigs	Assembly size	Max contig size	contig L50	contig L90	contig N50	contig N90
flye (no correction)	1519	1,353 Mb	42.8 Mb	33 fragments	136 contigs	14.5 Mb	1.4 Mb
racon 1	847	1344 Mb	46.8 Mb	29 contigs	117 contigs	15.9 Mb	1.5 Mb
racon 2	823	1343 Mb	46.8 Mb	29 contigs	117 contigs	15.9 Mb	1.5 Mb
racon 3	808	1343 Mb	46.8 Mb	29 contigs	117 contigs	15.9 Mb	1.5 Mb
pilon 2	808	1346 Mb	46.8 Mb	29 contigs	117 contigs	15.9 Mb	1.5 Mb
purge_haplotigs	441	1309 Mb	46.8 Mb	28 contigs	100 contigs	16 Mb	2.3 Mb

Polishing and haplotig purging reduce genome assembly errors

Whereas racon improved the contiguity of the assembly, it achieved only modest improvement in BUSCO benchmarking (Table 4). On the other hand, BUSCO scores greatly improved (76.5–92.6%) after pilon polishing. A second iteration of pilon resulted in only negligible improvement in BUSCO benchmarking. The improvement in BUSCO scores after polishing is likely due to correction of basecalling, indel, and local assembly errors that obscured identification of orthologs during AUGUSTUS annotation in the BUSCO pipeline. Haplotig purging had negligible impact on BUSCO scores, with a 0.2% reduction in duplications and a 0.1% increase in missing orthologs. Instead, purge_haplotigs identified 366 contigs as redundant, resulting in a more accurate genome representation while maintaining the overall genome assembly. Our final assembly recovered 92.6% of the BUSCO orthologs as complete.

Table 4

BUSCO summarized benchmarking at various stages of genome assembly

Stage	Complete (Single, Duplicate)	Fragmented	Missing
flye (no correction)	76.2% (74.7%, 1.5%)	10.6%	13.2%
racon 3	76.5% (75.1%, 1.4%)	11.0%	12.5%
pilon 1	92.6% (90.8%, 1.8%)	4.5%	2.9%
pilon 2	92.7% (90.8%, 1.9%)	4.5%	2.8%
purge_haplotigs	92.6% (90.9%, 1.7%)	4.5%	2.9%

Repetitive elements in the assembly

RepeatMasker estimated that 305 Mb (25.8%) of our assembly is comprised of repetitive elements. Nearly all of these elements (93.9%) – a total of 287 Mb (21.9% of the genome) – are part of the CR1 transposable element (TE) superfamily.

Birds have smaller genomes than most other tetrapods (Gregory ), likely due to the metabolic cost of flight (Wright ). The streamlining of avian genomes appears to be driven by a substantial reduction in the frequency of TEs in the genome (Gao ). For most tetrapods, TEs represent a sizeable portion of the genome. For example, in mammals TEs typically account for 1/3 to 1/2 of the entire genome (Platt ), while in birds’ closest living relatives, the crocodiles, TEs represent about 1/3 of the genome (Green ). In contrast, birds have substantially lower frequency of TEs – typically ranging from 4–10% of the total genome size (Sotero-Caio ). Woodpeckers and allies are the exception. Manthey demonstrated that TEs occurred at higher frequency across several genera of woodpeckers (17–30%), due principally to expansion of CR1 TEs (15–17% of total TEs). Here, CR1 TEs represent nearly 22% of the Golden-fronted Woodpecker genome. Manthey also observed that their TE identification protocol may underestimate the true TE frequency by 2–5%, so a direct comparison should be considered tentative. However, our result does represent the highest CR1 proportion reported for any avian genome to date.

Mitogenome assembly and annotation results

The mitochondrial genome assembly (16,844 base pairs; Figure 2) recovered all 13 protein-coding genes, 12S and 16S ribosomal RNAs, and 22 tRNAs, all in the standard order observed in other woodpeckers and allies (Mindell ; Eo 2017). Pairwise alignment of each of these features shows all features starting and ending coincident with homologous features in the Great Spotted Woodpecker mitochondrial genome, except in a small number of instances where our initial annotation failed to properly assess protein-coding genes with incomplete (-T) stop codons, which are common in vertebrate mitochondrial genomes (Ojala ). In nearly all cases where this occurred, the adjacent 3` base would result in a valid stop codon, so this “error” seems to be a limitation of the MitoAnnotator algorithm rather than a measure of error in our mitogenome assembly. Our initial flye assembly recovered the mitochondrion as a single contig, with a sequence length of 16,858 bases, and 99.6% pairwise identity to the final polished mitochondrion sequence. Differences were due to homopolymer indel errors that resulted in frameshift errors. This comparison provides independent evidence to suggest that racon consensus correction plus pilon short read polishing is able to correct most errors in long read genome assemblies.

Figure 2

Circularized annotated mitochondrial genome assembly for the Golden-fronted Woodpecker (Melanerpes aurifrons). Figure generated by MitoAnnotator [36].

Synteny with other avian genomes and the usefulness of reference-based scaffolding

In general, synteny plots show that most of our Anna’s Hummingbird guided reference-based scaffolded pseudo-chromosomes had largely one-to-one synteny plots with chromosomes of all three bird genomes evaluated (Figure 3). In some instances, the chromosomal numbers do not agree. For example, our pseudo-chromosome 5A maps almost exclusively to the Kakapo chromosome 4. This is likely to be simply the result of differences between chromosome numbering conventions. Both the Zebra Finch and Anna’s Hummingbird synteny plot show that a portion of those birds’ W chromosomes map to our Golden-fronted Woodpecker pseudo-chromosome 3, suggesting that this pseudo-molecule may be misassembled. We agree with Alonge that reference-guided assemblies will be generally useful for comparative genomics studies of animal biodiversity, such as to map genome-wide phylogenetic markers, or to map the locations of structural variants and/or F outliers, etc. However, given the relatively low number of contigs in our assembly, researchers may prefer to do those analyses on the 441 unscaffolded contigs rather than on the pseudo-chromosomes.

Figure 3

Ideograms showing synteny between Golden-fronted Woodpecker pseudo-chromosomes and chromosomes from three different bird species.

Comparisons to recently published avian genomes

There is increasing interest in obtaining reference genomes for wild bird research. While a full chromosome level assembly would be ideal, many research questions would benefit from a highly contiguous but less-than-chromosome assembly. How should research groups with limited resources go about generating such genomes? Comparing assembly statistics for some recently published avian genomes (Table 5), we find that hybrid (PacBio or Oxford Nanopore long reads + Illumina short reads) assemblies are generally more contiguous than assemblies generated with only Illumina sequencing, including those with third generation Illumina libraries (e.g., 10X and Hi-C) with our genome possibly being more complete than many other published genomes (Peona ). As de novo genome assembly algorithms such as flye (Kolmogorov ), RedBean (Ruan and Li 2020), Shasta (https://chanzuckerberg.github.io/shasta/), and Canu (Koren ) that are capable with dealing with noisy Oxford Nanopore data continue to be developed, we expect to see more researchers pursuing similar, low overhead, genome assemblies for their own studies of animal biodiversity genomics.

Table 5

Comparison of various genome assembly statistics for recently published wild bird genomes. * N50 scaffold not reported for assemblies comprised only of ungapped contigs. ** Statistics only reported for scaffolded contigs

Species	Reference	Method	Assembly size	# of contigs/# of scaffolds	N50 contig	N50 scaffold	BUSCO complete	% Repeats
Pavo cristatus (Galliformes)	(Dhar et al. 2019)	Nanopore (low coverage) + short read	932 Mb	685,241/ 15,025	14.7 kb	0.23 Mb	Not reported	7.3%
Melospiza melodia (Passeriformes)	(Louha et al. 2019)	Paired end + HiC reads	978 Mb	Not reported	31.7 kb	5.6 Mb	87.5%	7.4%
Rhegmatorhina melanosticta (Passeriformes)	(Coelho et al. 2019)	Paired end + 10× linked reads	1.03 Gb	Not reported / 715	137 kb	3.3 Mb	89.2%	Not reported
Numida meleagris (Galliformes)	(Vignal et al. 2019)	Paired end + Mate pair	1.04Gb	Not reported / 2,739	234 kb	7.8 Mb	90.7%	19.5%
Colinus virginanus (Galliformes)	(Salter et al. 2019)	Paired end + HiC reads	866 Mb	Not reported **/1,512	Not reported **	66.8 Mb	90.8%	Not reported
Eremophila alpestris (Passeriformes)	(Mason et al. 2019)	Paired end + mate pair	1.04 Gb	Not reported ** / 2,708	Not reported**	10.6 Mb	94.5%	Not reported
Grus nigricollis (Gruiformes)	(Zhou et al. 2019)	Nanopore + paired end	1.33 Gb	1,837/ NA*	17.9 Mb	NA*	97.7%	8.1%
Melanerpes aurifrons (Piciformes)	This study	Nanopore + paired end	1.30 Gb/	441 / NA*	16.0 Mb	NA*	92.6%	25.8%