Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the ‘venom-ome’ and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 ‘venom-ome-specific toxins’ (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

The snow leopard, Panthera uncia, is an elusive high-altitude specialist that inhabits vast, inaccessible habitat across Asia. We conducted the first range-wide genetic assessment of snow leopards based on noninvasive scat surveys. Thirty-three microsatellites were genotyped and a total of 683 bp of mitochondrial DNA sequenced in 70 individuals. Snow leopards exhibited low genetic diversity at microsatellites (A_N = 5.8, H_O = 0.433, H_E = 0.568), virtually no mtDNA variation, and underwent a bottleneck in the Holocene (∼8000 years ago) coinciding with increased temperatures, precipitation, and upward treeline shift in the Tibetan Plateau. Multiple analyses supported 3 primary genetic clusters: (1) Northern (the Altai region), (2) Central (core Himalaya and Tibetan Plateau), and (3) Western (Tian Shan, Pamir, trans-Himalaya regions). Accordingly, we recognize 3 subspecies, Panthera uncia irbis (Northern group), Panthera uncia uncia (Western group), and Panthera uncia uncioides (Central group) based upon genetic distinctness, low levels of admixture, unambiguous population assignment, and geographic separation. The patterns of variation were consistent with desert-basin “barrier effects” of the Gobi isolating the northern subspecies (Mongolia), and the trans-Himalaya dividing the central (Qinghai, Tibet, Bhutan, and Nepal) and western subspecies (India, Pakistan, Tajikistan, and Kyrgyzstan). Hierarchical Bayesian clustering analysis revealed additional subdivision into a minimum of 6 proposed management units: western Mongolia, southern Mongolia, Tian Shan, Pamir-Himalaya, Tibet-Himalaya, and Qinghai, with spatial autocorrelation suggesting potential connectivity by dispersing individuals up to ∼400 km. We provide a foundation for global conservation of snow leopard subspecies, and set the stage for in-depth landscape genetics and genomic studies.

The development of high-quality chromosomally assigned reference genomes constitutes a key feature for understanding genome architecture of a species and is critical for the discovery of the genetic blueprints of traits of biological significance. South American camelids serve people in extreme environments and are important fiber and companion animals worldwide. Despite this, the alpaca reference genome lags far behind those available for other domestic species. Here we produced a chromosome-level improved reference assembly for the alpaca genome using the DNA of the same female Huacaya alpaca as in previous assemblies. We generated 190X Illumina short-read, 8X Pacific Biosciences long-read and 60X Dovetail Chicago® chromatin interaction scaffolding data for the assembly, used testis and skin RNAseq data for annotation, and cytogenetic map data for chromosomal assignments. The new assembly VicPac3.1 contains 90% of the alpaca genome in just 103 scaffolds and 76% of all scaffolds are mapped to the 36 pairs of the alpaca autosomes and the X chromosome. Preliminary annotation of the assembly predicted 22,462 coding genes and 29,337 isoforms. Comparative analysis of selected regions of the alpaca genome, such as the major histocompatibility complex, the region involved in the Minute Chromosome Syndrome and candidate genes for high-altitude adaptations, reveal unique features of the alpaca genome. The alpaca reference genome VicPac3.1 presents a significant improvement in completeness, contiguity and accuracy over VicPac2 and is an important tool for the advancement of genomics research in all New World camelids.

Disorders of sex development (DSD) and reproduction are not uncommon among horses, though knowledge about their molecular causes is sparse. Here we characterized a ~200 kb homozygous deletion in chromosome 29 at 29.7–29.9 Mb. The region contains AKR1C genes which function as ketosteroid reductases in steroid hormone biosynthesis, including androgens and estrogens. Mutations in AKR1C genes are associated with human DSDs. Deletion boundaries, sequence properties and gene content were studied by PCR and whole genome sequencing of select deletion homozygotes and control animals. Deletion analysis by PCR in 940 horses, including 622 with DSDs and reproductive problems and 318 phenotypically normal controls, detected 67 deletion homozygotes of which 79% were developmentally or reproductively abnormal. Altogether, 8–9% of all abnormal horses were homozygous for the deletion, with the highest incidence (9.4%) among cryptorchids. The deletion was found in ~4% of our phenotypically normal cohort, ~1% of global warmblood horses and ponies, and ~7% of draught breeds of general horse population as retrieved from published data. Based on the abnormal phenotype of the carriers, the functionally relevant gene content, and the low incidence in general population, we consider the deletion in chromosome 29 as a risk factor for equine DSDs and reproductive disorders.