Genomes of endangered great hammerhead and shortfin mako sharks reveal historic population declines and high levels of inbreeding in great hammerhead.

Despite increasing threats of extinction to Elasmobranchii (sharks and rays), whole genome-based conservation insights are lacking. Here, we present chromosome-level genome assemblies for the Critically Endangered great hammerhead (Sphyrna mokarran) and the Endangered shortfin mako (Isurus oxyrinchus) sharks, with genetic diversity and historical demographic comparisons to other shark species. The great hammerhead exhibited low genetic variation, with 8.7% of the 2.77 Gbp genome in runs of homozygosity (ROH) > 1 Mbp and 74.4% in ROH >100 kbp. The 4.98 Gbp shortfin mako genome had considerably greater diversity and <1% in ROH > 1 Mbp. Both these sharks experienced precipitous declines in effective population size (Ne) over the last 250 thousand years. While shortfin mako exhibited a large historical Ne that may have enabled the retention of higher genetic variation, the genomic data suggest a possibly more concerning picture for the great hammerhead, and a need for evaluation with additional individuals.

High-Quality Reference Genome for an Arid-Adapted Mammal, the Banner-Tailed Kangaroo Rat (Dipodomys spectabilis).

Kangaroo rats in the genus Dipodomys are found in a variety of habitat types in western North America, including deserts, arid and semiarid grasslands, and scrublands. Many Dipodomys species are experiencing strong population declines due to increasing habitat fragmentation, with two species listed as federally endangered in the United States. The precarious state of many Dipodomys populations, including those occupying extreme environments, make species of this genus valuable subjects for studying the impacts of habitat degradation and fragmentation on population genomic patterns and for characterizing the genomic bases of adaptation to harsh conditions. To facilitate exploration of such questions, we assembled and annotated a reference genome for the banner-tailed kangaroo rat (Dipodomys spectabilis) using PacBio HiFi sequencing reads, providing a more contiguous genomic resource than two previously assembled Dipodomys genomes. Using the HiFi data for D. spectabilis and publicly available sequencing data for two other Dipodomys species (Dipodomys ordii and Dipodomys stephensi), we demonstrate the utility of this new assembly for studies of congeners by conducting inference of historic effective population sizes (Ne) and linking these patterns to the species' current extinction risk statuses. The genome assembly presented here will serve as a valuable resource for population and conservation genomic studies of Dipodomys species, comparative genomic research within mammals and rodents, and investigations into genomic adaptation to extreme environments and changing landscapes.

Atypical Dermatophytosis in 12 North American Porcupines (Erethizon dorsatum) from the Northeastern United States 2010-2017.

Twelve wild North American porcupines (Erethizon dorsatum) out of a total of 44 of this species examined in an 8-year period were diagnosed with dermatopathies while being cared for at two wildlife rehabilitation clinics. Biopsy and necropsy were performed on seven and five animals, respectively. Atypical dermatophytosis was diagnosed in all cases. Lesions consisted of diffuse severe epidermal hyperkeratosis and mild hyperplasia with mild lymphoplasmacytic dermatitis and no folliculitis. Dermatophytes were noted histologically as hyphae and spores in hair shafts, and follicular and epidermal keratin. Trichophyton sp. was grown in 5/6 animals where culture was performed, with a molecular diagnosis of Arthroderma benhamiae/Trichophyton mentagrophytes in these five cases. Metagenomic analysis of formalin-fixed paraffin-embedded tissue samples from three cases identified fungi from 17 orders in phyla Basidiomycota and Ascomycota. Alteration of therapy from ketaconazole, which was unsuccessful in four out of five early cases, to terbinafine or nitraconazole led to the resolution of disease and recovery to release in four subsequent animals. In all, six animals were euthanized or died due to dermatopathy, no cases resolved spontaneously, and six cases were resolved with therapy. The work we present demonstrates an atypical lesion and anatomical distribution due to dermatophytosis in a series of free-ranging wild porcupines and the successful development of novel techniques for extracting and sequencing nucleic acids from fungus in archival formalin-fixed paraffin-embedded animal tissue.

White shark genome reveals ancient elasmobranch adaptations associated with wound healing and the maintenance of genome stability.

The white shark (Carcharodon carcharias; Chondrichthyes, Elasmobranchii) is one of the most publicly recognized marine animals. Here we report the genome sequence of the white shark and comparative evolutionary genomic analyses to the chondrichthyans, whale shark (Elasmobranchii) and elephant shark (Holocephali), as well as various vertebrates. The 4.63-Gbp white shark genome contains 24,520 predicted genes, and has a repeat content of 58.5%. We provide evidence for a history of positive selection and gene-content enrichments regarding important genome stability-related genes and functional categories, particularly so for the two elasmobranchs. We hypothesize that the molecular adaptive emphasis on genome stability in white and whale sharks may reflect the combined selective pressure of large genome sizes, high repeat content, high long-interspersed element retrotransposon representation, large body size, and long lifespans, represented across these two species. Molecular adaptation for wound healing was also evident, with positive selection in key genes involved in the wound-healing process, as well as Gene Ontology enrichments in fundamental wound-healing pathways. Sharks, particularly apex predators such as the white shark, are believed to have an acute sense of smell. However, we found very few olfactory receptor genes, very few trace amine-associated receptors, and extremely low numbers of G protein-coupled receptors. We did however, identify 13 copies of vomeronasal type 2 (V2R) genes in white shark and 10 in whale shark; this, combined with the over 30 V2Rs reported previously for elephant shark, suggests this gene family may underlie the keen odorant reception of chondrichthyans.

Comparative transcriptomics of elasmobranchs and teleosts highlight important processes in adaptive immunity and regional endothermy.

BACKGROUND: Comparative genomic and/or transcriptomic analyses involving elasmobranchs remain limited, with genome level comparisons of the elasmobranch immune system to that of higher vertebrates, non-existent. This paper reports a comparative RNA-seq analysis of heart tissue from seven species, including four elasmobranchs and three teleosts, focusing on immunity, but concomitantly seeking to identify genetic similarities shared by the two lamnid sharks and the single billfish in our study, which could be linked to convergent evolution of regional endothermy. RESULTS: Across seven species, we identified an average of 10,877 Swiss-Prot annotated genes from an average of 32,474 open reading frames within each species' heart transcriptome. About half of these genes were shared between all species while the remainder included functional differences between our groups of interest (elasmobranch vs. teleost and endotherms vs. ectotherms) as revealed by Gene Ontology (GO) and selection analyses. A repeatedly represented functional category, in both the uniquely expressed elasmobranch genes (total of 259) and the elasmobranch GO enrichment results, involved antibody-mediated immunity, either in the recruitment of immune cells (Fc receptors) or in antigen presentation, including such terms as "antigen processing and presentation of exogenous peptide antigen via MHC class II", and such genes as MHC class II, HLA-DPB1. Molecular adaptation analyses identified three genes in elasmobranchs with a history of positive selection, including legumain (LGMN), a gene with roles in both innate and adaptive immunity including producing antigens for presentation by MHC class II. Comparisons between the endothermic and ectothermic species revealed an enrichment of GO terms associated with cardiac muscle contraction in endotherms, with 19 genes expressed solely in endotherms, several of which have significant roles in lipid and fat metabolism. CONCLUSIONS: This collective comparative evidence provides the first multi-taxa transcriptomic-based perspective on differences between elasmobranchs and teleosts, and suggests various unique features associated with the adaptive immune system of elasmobranchs, pointing in particular to the potential importance of MHC Class II. This in turn suggests that expanded comparative work involving additional tissues, as well as genome sequencing of multiple elasmobranch species would be productive in elucidating the regulatory and genome architectural hallmarks of elasmobranchs.

Laser capture microdissection microscopy and genome sequencing of the avian malaria parasite, Plasmodium relictum.

Acquiring genomic material from avian malaria parasites for genome sequencing has proven problematic due to the nucleation of avian erythrocytes, which produces a large ratio of host to parasite DNA (∼1 million to 1 bp). We tested the ability of laser capture microdissection microscopy to isolate parasite cells from individual avian erythrocytes for four avian Plasmodium species, and subsequently applied whole genome amplification and Illumina sequencing methods to Plasmodium relictum (lineage pSGS1) to produce sequence reads of the P. relictum genome. We assembled ∼335 kbp of parasite DNA from this species, but were unable to completely avoid contamination by host DNA and other sources. However, it is clear that laser capture microdissection holds promise for the isolation of genomic material from haemosporidian parasites in intracellular life stages. In particular, laser capture microdissection may prove useful for isolating individual parasite species from co-infected hosts. Although not explicitly tested in this study, laser capture microdissection may also have important applications for isolation of rare parasite lineages and museum specimens for which no fresh material exists.

Mitochondrial genome of an Atlantic white shark (Carcharodon carcharias).

Here we report the first full mitochondrial genome sequence for a white shark caught in the Atlantic Ocean. The mitochondrial genome is 16,745 bp in length and contains 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a non-coding control region. The base composition of this mtDNA lineage is A: 30.7%, C: 26.9%, G: 13.8%, and T: 28.6%. In concordance with previous population genetic studies, the Atlantic caught individual forms a separate lineage from individuals caught on either side of the Pacific Ocean.

Comparative analysis of active retrotransposons in the transcriptomes of three species of heteromyid rodents.

Virtually every eukaryotic genome is replete with transposable elements (TEs). TE activity, or lack thereof, is of considerable evolutionary interest as TE insertions actively drive genome evolution by altering gene expression and/or function through different mechanisms. Herein, we take a comparative approach to better understand the variation in TE transcriptional activity. Our goals were to identify transcriptionally active TEs and to evaluate their relative expression levels in an effort to identify key determinants of TE activity. We do so in three related non-model rodent species (Dipodomys spectabilis, Chaetodipus baileyi, and Heteromys desmarestianus). We used these species (a) because rodents have long been productive models for the study of TEs; (b) because the known and dated relationships among these species permit strong phylogenetic inference (e.g., ancestral character states); and (c) because we have previously characterized key genes that underlie evolutionary adaptations in these species (e.g., osmoregulation). We used RNA-seq to characterize the transcriptomes of two different tissues--kidney and spleen--in each of these three species. Our data revealed a diversity of retrotransposons that were actively transcribed in these rodents (including LINEs, retroviruses/ERVs, and SINEs). We also identified tissue-specific differences in retrotransposon-activity in about half of the retrotransposons detected. We interpret these differences in TE activity in light of overall gene expression, with the transcribed SINEs enriched in differentially expressed sequences in all three species. These transcribed TEs may simply represent transpositionally active TEs. Alternatively, they may be TEs that have been co-opted to represent functionally important components of the host transcriptome (i.e., exaptations). Overall, our results contribute to a growing appreciation for the evolutionary interplay between TEs and their hosts.

Transcriptomic characterization of the immunogenetic repertoires of heteromyid rodents.

BACKGROUND: When populations evolve under disparate environmental conditions, they experience different selective pressures that shape patterns of sequence evolution and gene expression. These may be manifested in genetic and phenotypic differences such as a diverse immunogenetic repertoire in species from tropical latitudes that have greater and/or different parasite burdens than more temperate species. To test this idea, we compared the transcriptomes of one tropical species (Heteromys desmarestianus) and two species from temperate latitudes (Dipodomys spectabilis and Chaetodipus baileyi) from the Heteromyidae. We did so in a search for positive selection on sequences and/or differential expression, while controlling for phylogenetic history in our choice of species. RESULTS: We identified 127,812 contigs and annotated 34,878 of these, identifying immune genes associated with interleukins, cytokines, and the production of mast cells. We identified 632 genes that were upregulated in H. desmarestianus (8.7% of genes tested) and 492 (6.7%) that were downregulated. Gene ontology terms including "immune response" were associated with 31 (4.9%) of the 632 upregulated genes. We found preliminary evidence for positive selection on three genes (Palmitoyltransferase ZDHHC5 Ubiquitin-conjugating enzyme E2 N, Krueppel-like factor 10, and Spindle and kinetochore-associated protein 1) along the H. desmarestianus lineage. CONCLUSIONS: Overall our findings pinpoint genes in species from disparate environments that are on different evolutionary trajectories in terms of expression levels and/or nucleotide sequence. Our data indicate there are significant differences in the expression of genes among the spleen transcriptomes of these species and that a number of these differentially expressed genes do not show the same pattern of differential expression in another tissue type. This points to the possibility of expression differences between these species specific to the spleen transcriptome.

Natural selection and the genetic basis of osmoregulation in heteromyid rodents as revealed by RNA-seq.

One adaptation of ecological and evolutionary interest is the extraordinary ability of desert rodents to retain water during waste production. Much is known regarding the unique kidney physiology of kangaroo rats (Dipodomys spp.) and their ability to retain water during waste production, yet the genetic basis of these physiological adaptations is relatively unknown. Herein, we utilized RNA-seq data to conduct a comparative study to identify osmoregulatory genes expressed in heteromyid rodents. We sequenced kidney tissue from two temperate desert species (Dipodomys spectabilis and Chaetodipus baileyi) from two separate subfamilies of the Heteromyidae and compared these transcriptomes to a tropical mesic species (Heteromys desmarestianus) from a third subfamily. The evolutionary history of these subfamilies provided a robust phylogenetic control that allowed us to separate shared evolutionary history from convergence. Using two methods to detect differential expression (DE), we identified 1890 genes that showed consistent patterns of DE between the arid and mesic species. A three-species reciprocal BLAST analysis revealed 3511 sets of putative orthologues that, upon comparison to known Mus musculus sequences, revealed 323 annotated and full-length genic regions. Selection tests displayed evidence of positive selection (dn/ds > 1) on six genes in the two desert species and remained significant for one of these genes after correction for multiple testing. Thus, our data suggest that both the coding sequence and expression of genes have been shaped by natural selection to provide the genetic architecture for efficient osmoregulation in desert-adapted heteromyid rodents.

Microsatellite analyses across three diverse vertebrate transcriptomes (Acipenser fulvescens, Ambystoma tigrinum, and Dipodomys spectabilis).

Historically, many population genetics studies have utilized microsatellite markers sampled at random from the genome and presumed to be selectively neutral. Recent studies, however, have shown that microsatellites can occur in transcribed regions, where they are more likely to be under selection. In this study, we mined microsatellites from transcriptomes generated by 454-pyrosequencing for three vertebrate species: lake sturgeon (Acipenser fulvescens), tiger salamander (Ambystoma tigrinum), and kangaroo rat (Dipodomys spectabilis). We evaluated (i) the occurrence of microsatellites across species; (ii) whether particular gene ontology terms were over-represented in genes that contained microsatellites; (iii) whether repeat motifs were located in untranslated regions or coding sequences of genes; and (iv) in silico polymorphism. Microsatellites were less common in tiger salamanders than in either lake sturgeon or kangaroo rats. Across libraries, trinucleotides were found more frequently than any other motif type, presumably because they do not cause frameshift mutations. By evaluating variation across reads assembled to a given contig, we were able to identify repeat motifs likely to be polymorphic. Our study represents one of the first comparative data sets on the distribution of vertebrate microsatellites within expressed genes. Our results reinforce the idea that microsatellites do not always occur in noncoding DNA, but commonly occur in expressed genes.

Of contigs and quagmires: next-generation sequencing pitfalls associated with transcriptomic studies.

Molecular ecologists have good reasons to be excited about the newest DNA/RNA sequencing technologies. However, this exuberance should be tempered with a hefty dose of reality: new sequencing technologies come with significant new challenges. Herein, we offer a brief overview of some practical problems encountered during transcriptomics studies conducted in our laboratory, and of nontrivial issues that prospective practitioners should consider. These include template contamination (e.g. from xenobiotics) and the cutting-room floor problem, whereby most of the data are often unassembled, unannotated and unused. We also highlight computational requirements, including hardware, personnel time and associated skill sets. We are very optimistic about the future of molecular ecology, but we hope this cautionary overview will help neophytes better recognize some key challenges associated with new technologies.

A priori and a posteriori approaches for finding genes of evolutionary interest in non-model species: osmoregulatory genes in the kidney transcriptome of the desert rodent Dipodomys spectabilis (banner-tailed kangaroo rat).

One common goal in evolutionary biology is the identification of genes underlying adaptive traits of evolutionary interest. Recently next-generation sequencing techniques have greatly facilitated such evolutionary studies in species otherwise depauperate of genomic resources. Kangaroo rats (Dipodomys sp.) serve as exemplars of adaptation in that they inhabit extremely arid environments, yet require no drinking water because of ultra-efficient kidney function and osmoregulation. As a basis for identifying water conservation genes in kangaroo rats, we conducted a priori bioinformatics searches in model rodents (Mus musculus and Rattus norvegicus) to identify candidate genes with known or suspected osmoregulatory function. We then obtained 446,758 reads via 454 pyrosequencing to characterize genes expressed in the kidney of banner-tailed kangaroo rats (Dipodomys spectabilis). We also determined candidates a posteriori by identifying genes that were overexpressed in the kidney. The kangaroo rat sequences revealed nine different a priori candidate genes predicted from our Mus and Rattus searches, as well as 32 a posteriori candidate genes that were overexpressed in kidney. Mutations in two of these genes, Slc12a1 and Slc12a3, cause human renal diseases that result in the inability to concentrate urine. These genes are likely key determinants of physiological water conservation in desert rodents.

Comparative transcriptomics and gene expression in larval tiger salamander (Ambystoma tigrinum) gill and lung tissues as revealed by pyrosequencing.

Biologists are beginning to unravel the complexities of gene expression in model organisms by studying the transcriptome, the complement of genes that are transcribed in a given tissue. It is unclear, however, if findings from model systems apply to non-model organisms because of environmental effects on gene expression. Furthermore, there have been few efforts to quantify how transcriptome or gene expression varies across individuals and across tissues in natural environments. Herein, we describe transcriptomic profiling of gene expression in lung and gill tissue of three larval tiger salamanders. We do so with a hierarchical experimental design that captures variation in expression among genes, among tissues, and among individuals. Using 454 pyrosequencing, we produced high-quality sequence data of 59 megabases and assembled ~200,000 reads into 19,501 contigs. These contigs BLASTed to 3,599 transcripts, of which 721 were expressed in both tissues, 1,668 were unique to gill, and 1,210 unique to lung. Our data showed tissue-specific patterns in gene expression level with variation among transcripts and individuals. We identified genes and gene ontology terms related to respiration and compared their relative expression levels between gill and lung tissues. We also found evidence of exogenous genes associated with larval salamanders, and we identified ~1400 potential molecular markers (microsatellites and single nucleotide polymorphisms) that are associated with expressed genes. Given the tissue-specific differences we observed in transcriptomes, these data reinforce the idea that changes in gene expression serve as a primary mechanism underlying phenotypic plasticity.