The genome sequence of the eastern grey squirrel, Sciurus carolinensis Gmelin, 1788

We present a genome assembly from an individual male Sciurus carolinensis (the eastern grey squirrel; Vertebrata; Mammalia; Eutheria; Rodentia; Sciuridae). The genome sequence is 2.82 gigabases in span. The majority of the assembly (92.3%) is scaffolded into 21 chromosomal-level scaffolds, with both X and Y sex chromosomes assembled.


Background
The eastern grey squirrel, Sciurus carolinensis, is native to eastern North America, where it plays important roles in forest regeneration through its habit of caching food nuts and seeds (Corbet & Hill, 1991) 1 . In North America, S. carolinensis has been introduced outside its native range such that it is now found from the Canadian Pacific northwest to Florida. S. carolinensis was introduced to Britain (in 1876), Ireland (in 1911), Italy (in 1948, South Africa (before 1900), Australia (in 1880s, extirpated in 1973 and Pitcairn island (in 1987) (see https://www. cabi.org/isc/datasheet/49075). S. carolinensis, which thrives in urban parklands and gardens, is classed as invasive in Europe and on Pitcairn island. In Britain and Ireland the expansion of S. carolinensis populations has driven decline in populations of the native red squirrel, Sciurus vulgaris, which we have also assembled (Mead et al., 2020). The negative impact of S. carolinensis is through interspecific competition, leading to competitive exclusion of S. vulgaris, and by their carriage of squirrelpox virus, to which they are resistant but S. vulgaris are not (Chantrey et al., 2014) (Darby et al., 2014. The S. carolinensis genome will aid analyses of resistance and susceptibility to squirrelpox, as well as to the genomics of invasiveness.

Genome sequence report
The genome was sequenced from DNA extracted from a naturally deceased male S. carolinensis collected as part of a squirrel monitoring project run by the Wildlife Trust for Lancashire, Manchester and North Merseyside. A total of 74fold coverage in Pacific Biosciences single-molecule long reads (N50 28 kb) and 40-fold coverage in 10X Genomics read clouds (from molecules with an estimated N50 of 19 kb) were generated. Primary assembly contigs were scaffolded with chromosome conformation HiC data (42-fold coverage). A contamination check identified a small number of lowcoverage contigs that were likely to have derived from an apicomplexan parasite infecting the squirrel (Léveillé et al., 2020); these were removed. Subsequent manual assembly curation corrected 272 missing/misjoins and removed three haplotypic duplications, reducing the scaffold number by 19% and increasing the scaffold N50 by 242% The final assembly has a total length of 2.82 Gb in 752 sequence scaffolds with a scaffold N50 of 148.2 Mb (Table 1). The majority, 92.3%, of the assembly sequence was assigned to 21 chromosomal-level scaffolds representing 19 autosomes (numbered by sequence length), and the X and Y sex chromosomes (Figure 1- Figure 5; Table 2) plus 13 unlocalised scaffolds (assigned to chromosomes but with ambiguous placement). The assembly has a BUSCO (Simão et al., 2015) completeness of 93.7% using the mammalia_odb9 reference set. The primary assembly is a large-scale mosaic of both haplotypes (i.e. is not 1 See https://animaldiversity.org/accounts/Sciurus_carolinensis/: fully phased) and we have therefore also deposited the contigs corresponding to the alternate haplotype. The S. carolinensis mSciCar1 genome sequence is largely collinear with that of S. vulgaris mSciVul1 (Figure 4).

Methods
The eastern grey squirrel specimen was collected by the Wildlife Trust for Lancashire, Manchester and North Merseyside as part of an ongoing programme of recovery of dead squirrels. A full tissue dissection and preservation in 80% ethanol was undertaken and the specimen accessioned by the Natural History Museum, London.
DNA was extracted using an agarose plug extraction from spleen tissue following the Bionano Prep Animal Tissue DNA Isolation Soft Tissue Protocol 2 . Pacific Biosciences CLR long read and 10X Genomics read cloud sequencing libraries were constructed according to the manufacturers' instructions. Sequencing was performed by the Scientific Operations core at the Wellcome Sanger Institute on Pacific Biosciences SEQUEL I (single molecule long read) and Illumina HiSeq X (10X Genomics Chromium). HiC data were generated using the Dovetail v1.0 kit and sequenced on HiSeq X.
See Table 3 for software versions and sources. Assembly was carried out using Falcon-unzip (Chin et al., 2016), haplotypic duplication was identified and removed with purge_dups (Guan et al., 2020) and a first round of scaffolding carried out with 10X Genomics read clouds using scaff10x. Scaffolding with Hi-C data was carried out using SALSA2. The Hi-C scaffolded assembly was polished with arrow using the PacBio data, then polished with the 10X Genomics Illumina data by aligning to the assembly with longranger align, calling variants with     The genome sequence is released openly for reuse. The S. carolinensis genome sequencing initiative is part of the Wellcome Sanger Institute's "25 genomes for 25 years" project 3 . It is also part of the Vertebrate Genome Project (VGP) 4 and the Darwin Tree of Life (DToL) project 5 . The specimen has been preserved in ethanol and deposited with the Natural History Museum, London under registration number NHMUK ZD 2019.214, where it will remain accessible to the research community for posterity. All raw data and the assembly have been deposited in the ENA. The genome will be annotated and presented through the Ensembl pipeline at the European Bioinformatics Institute. Raw data and assembly accession identifiers are reported in Table 1.