Mitochondrial Genome Sequence of the Land Snail Oreohelix idahoensis

We announce the complete mitochondrial genome sequence of Oreohelix idahoensis, a threatened land snail endemic to the Pacific Northwest of the United States. The circular genome is 14.2 kb and contains 13 protein-coding genes, 2 rRNA genes, and 21 tRNA genes.

M ountain snails (Oreohelix) are the most diverse group of land snails in North America (79 species) (1), and Oreohelix spp. are the dominant malacological fauna of the Rocky Mountains (2). While the genus is widely distributed, most species in this group are endemic to single mountains or individual limestone outcrops within a mountain. The ecological specialization of these snails to small islands of limestone in montane ecosystems makes them vulnerable to extirpation, and many species are considered endangered or highly threatened by the International Union for Conservation of Nature and Natural Resources (IUCN; https://www.iucnredlist.org/) and local state governments. However, conservation officials are hesitant to develop conservation plans for many Oreohelix species, as recent genetic studies indicate that there is a lack of phylogenetic evidence for many existing taxonomic units (3) and the possible existence of cryptic species (4). These findings suggest that our understanding of oreohelicid diversity is limited and highlight the need for genomic resources to unravel species relationships in this group so that conservation plans can be developed based on robust species relationships.
To address this need, we sent a tissue sample of a single Oreohelix idahoensis snail collected from Lucile, ID, to the HudsonAlpha Institute for Biotechnology (Birmingham, AL) for total genomic DNA isolation and library preparation. DNA was extracted using a MagAttract high-molecular-weight DNA kit from Qiagen (product number [PN] 67563), and a 10ϫ linked read library was prepared using a Chromium v.2 genome reagent kit (10xGenomics; PN 120258), with 2.5 ng as the starting material. The linkedread library was then sequenced on three HiSeq X 150-bp paired-end lanes, resulting in 50.2ϫ coverage of the genome.
After assembly, we oriented and linearized the O. idahoensis mitochondrial genome with all other available stylommatophoran and systellommatophoran mitochondrial genomes (37 species) using the cyclic gap-free alignment tool MARS v.1.0 (7) and then aligned sequences using MAFFT (v.7.307) (8). The genome was then annotated with the MITOS Web server (v.1.0) (9), and annotations were confirmed after visually comparing BLAST alignments of the genome to all existing stylommatophoran and systellommato-phoran genomes. We then constructed a maximum likelihood phylogeny using the aligned sequences with RAxML (v.8.2.9) (10), specifying a GTRGAMMA model of nucleotide evolution and evaluating nodal support by performing 100 bootstrap replicates, to determine the phylogenetic placement of O. idahoensis.
Data availability. The genome sequence has been deposited in GenBank under the accession number MK290736, and the reads used for assembly are available in the SRA under the accession number PRJNA553530.

ACKNOWLEDGMENTS
We are grateful to Joel Sauder of the Idaho Department of Fish and Game for facilitating permits and providing logistic support while conducting fieldwork (permit number 160302). We also thank Andrew Rankin for helpful comments on the manuscript and providing the COI gene sequence. This work was funded by an NSF Graduate Research Fellowship program award to T. Mason Linscott (NSF 1144152), an initial startup funding award from IBEST to Christine E. Parent, and an Institutional Development Award (IDeA) from the National