Construction of Red Fox Chromosomal Fragments from the Short-Read Genome Assembly

Whole Genomes Inform Genetic Rescue Strategy for Montane Red Foxes in North America

A few iconic examples have proven the value of facilitated gene flow for counteracting inbreeding depression and staving off extinction; yet, the practice is often not implemented for fear of causing outbreeding depression. Using genomic sequencing, climatic niche modeling, and demographic reconstruction, we sought to assess the risks and benefits of using translocations as a tool for recovery of endangered montane red fox (Vulpes vulpes) populations in the western United States. We demonstrated elevated inbreeding and homozygosity of deleterious alleles across all populations, but especially those isolated in the Cascade and Sierra Nevada ranges. Consequently, translocations would be expected to increase population growth by masking deleterious recessive alleles. Demographic reconstructions further indicated shallow divergences of less than a few thousand years among montane populations, suggesting low risk of outbreeding depression. These genomic-guided findings set the stage for future management, the documentation of which will provide a roadmap for recovery of other data-deficient taxa.

Maintenance of a narrow hybrid zone between native and introduced red foxes (Vulpes vulpes) despite conspecificity and high dispersal capabilities

Human-facilitated introductions of nonnative populations can lead to secondary contact between allopatric lineages, resulting in lineage homogenisation or the formation of stable hybrid zones maintained by reproductive barriers. We investigated patterns of gene flow between the native Sacramento Valley red fox (Vulpes vulpes patwin) and introduced conspecifics of captive-bred origin in California's Central Valley. Considering their recent divergence (20-70 kya), we hypothesised that any observed barriers to gene flow were primarily driven by pre-zygotic (e.g. behavioural differences) rather than post-zygotic (e.g. reduced hybrid fitness) barriers. We also explored whether nonnative genes could confer higher fitness in the human-dominated landscape resulting in selective introgression into the native population. Genetic analysis of red foxes (n = 682) at both mitochondrial (cytochrome b + D-loop) and nuclear (19,051 SNPs) loci revealed narrower cline widths than expected under a simulated model of unrestricted gene flow, consistent with the existence of reproductive barriers. We identified several loci with reduced introgression that were previously linked to behavioural divergence in captive-bred and domestic canids, supporting pre-zygotic, yet possibly hereditary, barriers as a mechanism driving the narrowness and stability of the hybrid zone. Several loci with elevated gene flow from the nonnative into the native population were linked to genes associated with domestication and adaptation to human-dominated landscapes. This study contributes to our understanding of hybridisation dynamics in vertebrates, particularly in the context of species introductions and landscape changes, underscoring the importance of considering how multiple mechanisms may be maintaining lineages at the species and subspecies level.

Chromosome-level genome assembly of the Arctic fox (Vulpes lagopus) using PacBio sequencing and Hi-C technology

The Arctic fox (Vulpes lagopus) is the only fox species occurring in the Arctic and has adapted to its extreme climatic conditions. Currently, the molecular basis of its adaptation to the extreme climate has not been characterized. Here, we applied PacBio sequencing and chromosome structure capture technique to assemble the first V. lagopus genome assembly, which is assembled into chromosome fragments. The genome assembly has a total length of 2.345 Gb with a contig N50 of 31.848 Mb and a scaffold N50 of 131.537 Mb, consisting of 25 pseudochromosomal scaffolds. The V. lagopus genome had approximately 32.33% repeat sequences. In total, 21,278 protein-coding genes were predicted, of which 99.14% were functionally annotated. Compared with 12 other mammals, V. lagopus was most closely related to V. Vulpes with an estimated divergence time of ~7.1 Ma. The expanded gene families and positively selected genes potentially play roles in the adaptation of V. lagopus to Arctic extreme environment. This high-quality assembled genome will not only promote future studies of genetic diversity and evolution in foxes and other canids but also provide important resources for conservation of Arctic species.

Using Ancestry-Informative SNPs to Quantify Introgression of European Alleles into North American Red Foxes

A recent study demonstrated that British red foxes introduced to the mid-Atlantic coastal plain (ACP) of the eastern United States during the late 18th century successfully interbred with indigenous American red foxes despite half a million year's divergence. However, a large disparity in frequency of European mitochondria (27%) versus Y chromosomes (1%) left unclear the magnitude of genetic exchange. We sought to quantify genomic introgression using 35 autosomal and 5 X-chromosome ancestry-informative markers (AIMs) in conjunction with diagnostic Y chromosome single nucleotide polymorphism (Y-SNP) markers to characterize the modern state of red foxes in the eastern United States and to gain insight into the potential role of reproductive barriers. European admixture was highest in the ACP and apparently restricted to the central eastern United States. We estimated only slightly (and nonsignificantly) European ancestry in autosomal than X-chromosome markers. European ancestry from autosomal and X-chromosome markers (36.4%) was higher than the corresponding mitochondrial (mt) DNA estimate (26.4%) in the ACP. Only 1 of 124 males (<1%) in the ACP had European Y chromosomes, which was similar to the neighboring regions, in which 2 of 99 (2%) males carried a European Y chromosome (the same haplotype). Although we could not rule out drift as the cause of low European Y-chromosome frequency, results were also consistent with F1 male infertility. In the future, more extensive genomic sequencing will enable a more thorough investigation of possible barrier genes on the X chromosome as well as throughout the genome.

Hypothalamic transcriptome of tame and aggressive silver foxes (Vulpes vulpes) identifies gene expression differences shared across brain regions

The underlying neurological events accompanying dog domestication remain elusive. To reconstruct the domestication process in an experimental setting, silver foxes (Vulpes vulpes) have been deliberately bred for tame vs aggressive behaviors for more than 50 generations at the Institute for Cytology and Genetics in Novosibirsk, Russia. The hypothalamus is an essential part of the hypothalamic-pituitary-adrenal axis and regulates the fight-or-flight response, and thus, we hypothesized that selective breeding for tameness/aggressiveness has shaped the hypothalamic transcriptomic profile. RNA-seq analysis identified 70 differentially expressed genes (DEGs). Seven of these genes, DKKL1, FBLN7, NPL, PRIMPOL, PTGRN, SHCBP1L and SKIV2L, showed the same direction expression differences in the hypothalamus, basal forebrain and prefrontal cortex. The genes differentially expressed across the three tissues are involved in cell division, differentiation, adhesion and carbohydrate processing, suggesting an association of these processes with selective breeding. Additionally, 159 transcripts from the hypothalamus demonstrated differences in the abundance of alternative spliced forms between the tame and aggressive foxes. Weighted gene coexpression network analyses also suggested that gene modules in hypothalamus were significantly associated with tame vs aggressive behavior. Pathways associated with these modules include signal transduction, interleukin signaling, cytokine-cytokine receptor interaction and peptide ligand-binding receptors (eg, G-protein coupled receptor [GPCR] ligand binding). Current studies show the selection for tameness vs aggressiveness in foxes is associated with unique hypothalamic gene profiles partly shared with other brain regions and highlight DEGs involved in biological processes such as development, differentiation and immunological responses. The role of these processes in fox and dog domestication remains to be determined.

The Red Fox Y-Chromosome in Comparative Context

While the number of mammalian genome assemblies has proliferated, Y-chromosome assemblies have lagged behind. This discrepancy is caused by biological features of the Y-chromosome, such as its high repeat content, that present challenges to assembly with short-read, next-generation sequencing technologies. Partial Y-chromosome assemblies have been developed for the cat (Feliscatus), dog (Canislupusfamiliaris), and grey wolf (Canislupuslupus), providing the opportunity to examine the red fox (Vulpesvulpes) Y-chromosome in the context of closely related species. Here we present a data-driven approach to identifying Y-chromosome sequence among the scaffolds that comprise the short-read assembled red fox genome. First, scaffolds containing genes found on the Y-chromosomes of cats, dogs, and wolves were identified. Next, analysis of the resequenced genomes of 15 male and 15 female foxes revealed scaffolds containing male-specific k-mers and patterns of inter-sex copy number variation consistent with the heterogametic chromosome. Analyzing variation across these two metrics revealed 171 scaffolds containing 3.37 Mbp of putative Y-chromosome sequence. The gene content of these scaffolds is consistent overall with that of the Y-chromosome in other carnivore species, though the red fox Y-chromosome carries more copies of BCORY2 and UBE1Y than has been reported in related species and fewer copies of SRY than in other canids. The assignment of these scaffolds to the Y-chromosome serves to further characterize the content of the red fox draft genome while providing resources for future analyses of canid Y-chromosome evolution.

A Near Chromosome Assembly of the Dromedary Camel Genome

The dromedary camel is an economically and socially important species of livestock in many parts of the world, being used for transport and the production of milk and meat. Much like cattle and horses, the camel may be found in industrial farming conditions as well as used in sporting. Camel racing is a multi-million dollar industry, with some specimens being valued at upward of 9.5 million USD. Despite its apparent value to humans, the dromedary camel is a neglected species in genomics. While cattle and other domesticated species have had much attention in terms of genome assembly, the camel has only been assembled to scaffold level, which does not give a clear indication of the order or chromosomal location of sequenced fragments. In this study, the Reference Assistant Chromosome Assembly (RACA) algorithm was implemented to use read-pair information of camel scaffolds, aligned with the cattle and human genomes in order to organize and orient these scaffolds in a near-chromosome level assembly. This method generated 72 large size fragments (N50 54.36 Mb). These predicted chromosome fragments (PCFs) were then compared with comparative maps of camel and cytogenetic map of alpaca chromosomes, allowing us to further upgrade the assembly. This dromedary camel assembly will be an invaluable tool to verify future camel assemblies generated with chromatin conformation or/and long read technologies. This study provides the first near-chromosome assembly of the dromedary camel, thus adding this economically important species to a growing pool of knowledge regarding the genome structure of domesticated livestock.

A near-chromosome-scale genome assembly of the gemsbok (Oryx gazella): an iconic antelope of the Kalahari desert

Background

The gemsbok (Oryx gazella) is one of the largest antelopes in Africa. Gemsbok are heterothermic and thus highly adapted to live in the desert, changing their feeding behavior when faced with extreme drought and heat. A high-quality genome sequence of this species will assist efforts to elucidate these and other important traits of gemsbok and facilitate research on conservation efforts.

Findings

Using 180 Gbp of Illumina paired-end and mate-pair reads, a 2.9 Gbp assembly with scaffold N50 of 1.48 Mbp was generated using SOAPdenovo. Scaffolds were extended using Chicago library sequencing, which yielded an additional 114.7 Gbp of DNA sequence. The HiRise assembly using SOAPdenovo + Chicago library sequencing produced a scaffold N50 of 47 Mbp and a final genome size of 2.9 Gbp, representing 90.6% of the estimated genome size and including 93.2% of expected genes according to Benchmarking Universal Single-Copy Orthologs analysis. The Reference-Assisted Chromosome Assembly tool was used to generate a final set of 47 predicted chromosome fragments with N50 of 86.25 Mbp and containing 93.8% of expected genes. A total of 23,125 protein-coding genes and 1.14 Gbp of repetitive sequences were annotated using de novo and homology-based predictions.

Conclusions

Our results provide the first high-quality, chromosome-scale genome sequence assembly for gemsbok, which will be a valuable resource for studying adaptive evolution of this species and other ruminants.

Insights into Emydid Turtle Cytogenetics: The European Pond Turtle as a Model Species

Our knowledge of Testudines evolution is limited by the lack of modern cytogenetic data. Compared to other reptiles, there is little information even on chromosome banding, let alone molecular cytogenetic data. Here, we provide detailed information on the karyotype of the European pond turtle Emys orbicularis, a model Emydidae, employing both chromosome banding and molecular cytogenetics. We provide a high-resolution G-banded karyotype and a map of rDNA genes and telomeric sequences using fluorescence in situ hybridization. We test hypotheses of sex-determining mechanisms in Emys by comparative genomic hybridization to determine if Emys has a cryptic sex-specific region. Our results provide valuable data to guide future efforts on genome sequencing and anchoring in Emydidae and for understanding karyotype evolution in Testudines.