F1000Research F1000Research 2046-1402 F1000 Research Limited London, UK 10.12688/f1000research.12636.1 Data Note Articles Bioinformatics Genomics Draft genome of tule elk Cervus elaphus nannodes

[version 1; peer review: 1 approved, 1 approved with reservations]

 Mizzi Jessica E. Data Curation Formal Analysis Software Writing – Original Draft Preparation Writing – Review & Editing https://orcid.org/0000-0002-6604-4618 1 Lounsberry Zachary T. Resources Writing – Review & Editing 2 Brown C. Titus Resources Supervision Writing – Review & Editing https://orcid.org/0000-0001-6001-2677 3 Sacks Benjamin N. Conceptualization Funding Acquisition Resources Supervision Writing – Review & Editing https://orcid.org/0000-0003-0143-6589 a 4 Microbiology Graduate Group, University of California, Davis, CA, 95616, USA Veterinary Genetics Laboratory, University of California, Davis, CA, 95616, USA Department of Population Health and Reproduction, School of Veterinary Medicine, University of California, Davis, CA, 95616, USA Department of Population Health and Reproduction and Mammalian Ecology and Conservation Unit, Veterinary Genetics Laboratory, School of Veterinary Medicine, University of California, Davis, CA, 95616, USA bnsacks@ucdavis.edu

Competing interests: No competing interests were disclosed.

 15 9 2017 2017 6 1691 19 6 2026 Copyright: © 2017 Mizzi JE et al. 2017 This is an open access article distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper presents the first draft genome of the tule elk ( Cervus elaphus nannodes), a sub species native to California that underwent an extreme genetic bottleneck in the late 1800s. The genome was generated from Illumina HiSeq 3000 whole genome sequencing of four individuals, resulting in the assembly of 2.395 billion base pairs (Gbp) over 602,862 contigs over 500 bp and N50 = 6,885 bp. This genome provides a resource to facilitate future genomic research on elk and other cervids.

 Cervus elaphus nannodes genome draft mammalian genome assembly tule elk California Department of Fish and Wildlife, Big Game Management Program P1580009 Support for this project was provided by a grant to BNS from the California Department of Fish and Wildlife, FY1516 Big Game Management Program (Grant ID P1580009). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Introduction

To date, the closest genomic resource for elk ( Cervus elaphus) is a full mitochondrial assembly of white-tailed deer ( Odocoileus virginianus), a distantly related cervid 1 . The present paper presents the first de novo genomic draft of the tule elk ( C. elaphus nannodes). This California-endemic elk subspecies underwent a major genetic bottleneck when its numbers were reduced to as few as three individuals in the 1870s 2, 3 . Although their numbers have increased to >5,000 today 4 , the historical bottleneck nevertheless left its mark on the elk’s genome, rendering it more homozygous than other elk subspecies.

Our motivation for generating a genomic resource for the tule elk was to create a reference for identifying single nucleotide polymorphisms (SNPs) to develop assays to monitor elk population abundance and for related population genetic applications. Due to the relatively low coverage generated in this work (40X overall with an average of 10X coverage from each individual), we used the MEGAHIT metagenome assembler, which has been found to perform well on low-quality or low-coverage DNA sequencing in bacteria 5 .

 Methods Sample collection and library prep

Elk were selected from four geographically distinct populations across northern California to maximize genomic diversity (San Luis Reservoir, California Valley, American Canyon, and the San Luis National Wildlife Refuge 4 ). Genomic DNA was extracted from skin biopsies, which were obtained by the California Department of Fish and Wildlife as part of their elk management activities 4 . We extracted DNA from skin using Qiagen DNeasy blood & tissue kits (QIAGEN Inc., Valencia, CA), according to the manufacturer’s instructions. The DNA was then fragmented via sonication using a Bioruptor (Diagenode, Denville, NJ) to 300 to 400 base pairs (bp) prior to adapter ligation. After verification of fragment size range using agarose gel electrophoresis, NEBNext® Ultra™ DNA Library Prep Kit for Illumina® (New England Biolabs, Inc., Ipswich, MA) was used to ligate Illumina adapters. Multiplexed libraries were prepared using NEBNext Multiplex Oligos for Illumina (New England Biolabs) to individually barcode each of four individual elk. Barcodes were annealed using low-cycle polymerase chain reactions during library preparation. To assess library quality, trace analysis was performed using a Bioanalyzer 2100 (Agilent, Santa Clara, CA) and fluorometric DNA quantitation of libraries was performed using a Qubit fluorometer (Invitrogen, Carlsbad, CA) prior to equilibrating sample concentrations and pooling for sequencing. After library quality control, four samples (one from each population) were pooled in equimolar concentrations and submitted for paired-end sequencing. Samples were sequenced on an Illumina HiSeq 3000 at the DNA Technologies and Expression Analysis Core of the UC Davis Genome Center.

 Bioinformatics processing

Sequencing quality on demultiplexed reads was evaluated using FastQC v0.11.3 (RRID:SCR_014583) 6 . The Illumina TruSeq3-PE sequencing adapters were removed using Trimmomatic v0.30 (RRID:SCR_011848) 7 with the ILLUMINACLIP parameter set to TruSeq3-PE.fa:2:40:15. The TruSeq3-PE.fa sequence was downloaded from https://anonscm.debian.org/cgit/debian-med/trimmomatic.git/plain/adapters/TruSeq3-PE.fa. LEADING, TRAILING, and SLIDING parameters were set to 2, resulting in the removal of bases with a quality score of 2 or less according to a phred33 quality scoring matrix. The SLIDINGWINDOW parameter of 4:2 was used to clip reads once the quality score fell below 2 within the window. The MINLENGTH parameter set to 25 dropped any reads that fell below that length due to quality trimming. The demultiplexed, quality-filtered reads were interleaved using the interleave-reads.py script in khmer v2.0 (RRID:SCR_001156) 8 . The assembly was performed using MEGAHIT v1.0.5 9 on interleaved quality filtered reads. Genome statistical analysis was done using QUAST v3.0 (RRID:SCR_001228) 10 . All code used is publicly available at https://github.com/dib-lab/2017-tule-elk/.

 Results

We obtained 377,980,276 demultiplexed 150 bp paired-end raw reads, containing a total of 113.394 Gbp of sequence, or approximately 40X coverage of the approximately 3 Gbp tule elk genome. Sequence assembly resulted in the generation of a total genome sequence size of 2.395 Gbp. Reads were assembled into 602,862 contiguous sequences ("contigs") averaging 3,973 bp in length with a minimum contig length of 201 bp. The G+C content of the genome was 41.55%. The N50 was 6,885 bp and maximum contig length was 72,391 bp. Additional assembly statistics are available in Table 1. No contigs (e.g. under a certain size or likely to reflect repeats) were removed from the assembly.

Quality metrics on tule elk ( <italic toggle="yes">Cervus elaphus nannodes</italic>) assembly, as generated with QUAST v3.0.
Metric Tule elk assembly
# contigs (≥ 200 bp) 1,367,218
# contigs ≥ 500 bp 602,862
# contigs (≥ 1000 bp) 460,702
# contigs (≥ 5000 bp) 160,229
# contigs (≥ 10000 bp) 51,790
# contigs (≥ 25000 bp) 2,606
# contigs (≥ 50000 bp) 36
Total length (≥ 200 bp) 2,607,088,486
Total length (≥ 1000 bp) 2,295,163,580
Total length (≥ 5000 bp) 1,531,314,985
Total length (≥ 10000 bp) 771,863,493
Total length (≥ 25000 bp) 80,157,993
Total length (≥ 50000 bp) 2,056,962
Largest contig 72,391
Total length 2,395,105,945
GC 41.55%
N50 6,885
N75 3,646
L50 103,346
L75 222,107
# N's per 100 kbp 0

This genome can serve as the basis for further genomic work on tule elk and other cervids, such as the development of a SNP assay to track elk population movement across increasingly developed northern Californian terrain. Furthermore, it is the first whole genome assembly available from the family Cervidae, providing a useful interim reference genome for bioinformatic analyses on other deer and elk species.

 Data availability

Raw reads are available in the SRA under the BioProject ID PRJNA345218. The genome draft is available at https://doi.org/10.6084/m9.figshare.5382565.v1 11 .

Code used in this study have been archived at http://doi.org/10.5281/zenodo.887935 12

Acknowledgements

JM would like to thank Luiz Irber, Camille Scott, and Lisa Johnson of the DIB lab at UC Davis for assistance with bioinformatics processing. We also thank C. Langner and J. Hobbs of the California Department of Fish and Wildlife for providing samples.

 Seabury CM Bhattarai EK Taylor JF : Correction: Genome-Wide Polymorphism and Comparative Analyses in the White-Tailed Deer ( Odocoileus Virginianus): A Model for Conservation Genomics. PLoS One. Public Library of Science.2011;6(2): e15811. 21283515 10.1371/annotation/80c6965b-ffa0-4100-8a20-3f41a08b2894 3034728 McCullough DR : The Tule Elk: Its History, Behavior, and Ecology. elibrary.ru.1969. Reference Source Sacks BN Lounsberry ZT Kalani T : Development and Characterization of 15 Polymorphic Dinucleotide Microsatellite Markers for Tule Elk Using HiSeq3000. J Hered. 2016;107(7):666–669. 27638816 10.1093/jhered/esw069 Hobbs JH : Draft Conservation and Management Plan for Elk. Report to the California Department of Fish and Wildlife.2014;41. Seitz A Nieselt K : Improving ancient DNA genome assembly. PeerJ. 2017;5:e3126. 28392981 10.7717/peerj.3126 5384568 Andrews S : FastQC: a quality control tool for high throughput sequence data.2010. Reference Source Bolger AM Lohse M Usadel B : Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. 24695404 10.1093/bioinformatics/btu170 4103590 Crusoe MR Alameldin HF Awad S : The khmer software package: enabling efficient nucleotide sequence analysis [version 1; referees: 2 approved, 1 approved with reservations]. F1000Res. 2015;4:900. 26535114 10.12688/f1000research.6924.1 4608353 Li D Luo R Liu CM : MEGAHIT v1.0: A fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods. 2016;102:3–11. 27012178 10.1016/j.ymeth.2016.02.020 Gurevich A Saveliev V Vyahhi N : QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29(8):1072–75. 23422339 10.1093/bioinformatics/btt086 3624806 Mizzi J Lounsberry ZT Brown CT : Tule Elk Draft Genome. figshare. 2017. Data Source Mizzi J : dib-lab/2017-tule-elk: Ready for publication. Zenodo. 2017. Data Source 10.5256/f1000research.13682.r27606 Reviewer response for version 1 Brauning Rudiger 1 Referee https://orcid.org/0000-0001-5068-9166 Invermay Agricultural Centre, AgResearch, Mosgiel, New Zealand

Competing interests: No competing interests were disclosed.

 16 11 2017 Copyright: © 2017 Brauning R 2017 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve-with-reservations

The authors describe the generation of a draft assembly for tule elk in the style of a brief genome announcement. For SNP detection and primer design this assembly is fine. It could e.g. be used in combination with Genotyping by Sequencing on additional individuals.

Materials and methods are sound and provided in full.

However a quick search of NCBI's taxonomy resource reveals that since June 2017 there is a genome assembly for red deer available https://www.ncbi.nlm.nih.gov/genome/10790. The authors therefore cannot claim to present the first whole genome assembly from the family Cervidae. Please change that statement.

Suggested further improvements:

Results

I would have liked to see a figure for the total amount of sequence after filtering as a simple way of showing how good or bad the sequence run was.

Table 1's readability would be improved by getting all figures to align right.

I'd also recommend to add another assembly metric to look at the gene content; either using something like BUSCO or by mapping the refseq sequences of a related, well annotated species (e.g. cattle) against the draft genome.

Methods

Sample collection and library prep

I see that each individual has two tissue samples. The authors entered a sample ID into the 'tissue' field of NCBI's BioSample database. I'd recommend removing this and adding the animal ID in the 'isolate' field.

Please expand the entries in the 'isolation source' field. It says e.g. "Am. Cyn" which probably means American Canyon.

Bioinformatics processing

Checking the code I believe the statement "LEADING, TRAILING, and SLIDING parameters were set to 2" should read "LEADING and TRAILING parameters were set to 2".

Are sufficient details of methods and materials provided to allow replication by others?

Yes

Is the rationale for creating the dataset(s) clearly described?

Yes

Are the datasets clearly presented in a useable and accessible format?

Yes

Are the protocols appropriate and is the work technically sound?

Yes

Reviewer Expertise:

NA

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.

Mizzi Jessica UC Davis, USA

Competing interests: I declare no competing interests.

 5 12 2017

Thank you for your review of this paper. Version 2 has been edited to reflect the presence of the red deer genome and a citation to that genome has been made. Table 1 has been reformatted for readability. The changes you’ve requested to the NCBI BioSample entry have been made. The trimmomatic code in the Bioinformatics Processing section has been edited to remove the erroneous “SLIDING” parameter. We’ve added text to the first sentence of the results section that describes the quality of sequence data in terms of standard quality scores. We opted not to provide details on the gene content relative to a related genome as we felt this could be done more comprehensively in the future once the red deer genome has been published and peer-reviewed.

 10.5256/f1000research.13682.r26607 Reviewer response for version 1 Olsen Steve 1 Referee Infectious Bacterial Diseases Research Unit, National Animal Disease Center, ARS-USDA, Ames, IA, USA

Competing interests: No competing interests were disclosed.

 25 10 2017 Copyright: © 2017 Olsen S 2017 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve

This article describes the generation of a draft genome (40X coverage from 4 animals) of the tule elk (Cervus elaphus nannodes). The research methods are fairly standard for the Illumina sequencing used.  At 602,862 contigs, the genome is very prelminary and will require quite a bit of additional work in order for it to be applicable to a wide range of applications.  The report basically falls into a category of a genome announcement.

Are sufficient details of methods and materials provided to allow replication by others?

Yes

Is the rationale for creating the dataset(s) clearly described?

Yes

Are the datasets clearly presented in a useable and accessible format?

Yes

Are the protocols appropriate and is the work technically sound?

Yes

Reviewer Expertise:

NA

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

Mizzi Jessica UC Davis, USA

Competing interests: I declare no competing interests.

 5 12 2017

Thank you for your review of this paper.