Agassiz’s desert tortoise (Gopherus agassizii) is a long-lived species native to the Mojave Desert and is listed as threatened under the US Endangered Species Act. To aid conservation efforts for preserving the genetic diversity of this species, we generated a whole genome reference sequence with an annotation based on deep transcriptome sequences of adult skeletal muscle, lung, brain, and blood. The draft genome assembly for G. agassizii has a scaffold N50 length of 252 kbp and a total length of 2.4 Gbp. Genome annotation reveals 20,172 protein-coding genes in the G. agassizii assembly, and that gene structure is more similar to chicken than other turtles. We provide a series of comparative analyses demonstrating (1) that turtles are among the slowest-evolving genome-enabled reptiles, (2) amino acid changes in genes controlling desert tortoise traits such as shell development, longevity and osmoregulation, and (3) fixed variants across the Gopherus species complex in genes related to desert adaptations, including circadian rhythm and innate immune response. This G. agassizii genome reference and annotation is the first such resource for any tortoise, and will serve as a foundation for future analysis of the genetic basis of adaptations to the desert environment, allow for investigation into genomic factors affecting tortoise health, disease and longevity, and serve as a valuable resource for additional studies in this species complex.

Data Availability: All genomic and transcriptomic sequence files are available from the NIH-NCBI BioProject database (accession numbers PRJNA352725, PRJNA352726, and PRJNA281763). All genome assembly, transcriptome assembly, predicted protein, transcript, genome annotation, repeatmasker, phylogenetic trees, .vcf and GO enrichment files are available on Harvard Dataverse (doi:10.7910/DVN/EH2S9K).

Background: The green anole lizard, Anolis carolinensis, is a key species for both laboratory and field-based studies of evolutionary genetics, development, neurobiology, physiology, behavior, and ecology. As the first non-avian reptilian genome sequenced, A. carolinesis is also a prime reptilian model for comparison with other vertebrate genomes. The public databases of Ensembl and NCBI have provided a first generation gene annotation of the anole genome that relies primarily on sequence conservation with related species. A second generation annotation based on tissue-specific transcriptomes would provide a valuable resource for molecular studies.

Results: Here we provide an annotation of the A. carolinensis genome based on de novo assembly of deep transcriptomes of 14 adult and embryonic tissues. This revised annotation describes 59,373 transcripts, compared to 16,533 and 18,939 currently for Ensembl and NCBI, and 22,962 predicted protein-coding genes. A key improvement in this revised annotation is coverage of untranslated region (UTR) sequences, with 79% and 59% of transcripts containing 5’ and 3’ UTRs, respectively. Gaps in genome sequence from the current A. carolinensis build (Anocar2.0) are highlighted by our identification of 16,542 unmapped transcripts, representing 6,695 orthologues, with less than 70% genomic coverage.

Conclusions: Incorporation of tissue-specific transcriptome sequence into the A. carolinensis genome annotation has markedly improved its utility for comparative and functional studies. Increased UTR coverage allows for more accurate predicted protein sequence and regulatory analysis. This revised annotation also provides an atlas of gene expression specific to adult and embryonic tissues.

In species with highly heteromorphic sex chromosomes, the degradation of one of the sex chromosomes will result in unequal gene expression between the sexes (e.g. between XX females and XY males) and between the sex chromosomes and the autosomes. Dosage compensation is a process whereby genes on the sex chromosomes achieve equal gene expression. We compared genome-wide levels of transcription between males and females, and between the X chromosome and the autosomes in the green anole, Anolis carolinensis. We present evidence for dosage compensation between the sexes, and between the sex chromosomes and the autosomes. When dividing the X chromosome into regions based on linkage groups, we discovered that genes in the first reported X-linked region, anole linkage group b (LGb), exhibit complete dosage compensation, although the rest of the X-linked genes exhibit incomplete dosage compensation. Our data further suggest that the mechanism of this dosage compensation is upregulation of the X chromosome in males. We report that approximately 10% of coding genes, most of which are on the autosomes, are differentially expressed between males and females. In addition, genes on the X chromosome exhibited higher ratios of nonsynonymous to synonymous substitution than autosomal genes, consistent with the fast-X effect. Our results from the green anole add an additional observation of dosage compensation in a species with XX/XY sex determination.

Bacterial lipopolysaccharides (LPS) are structural components of the outer membranes of Gram-negative bacteria and also are potent inducers of inflammation in mammals. Higher vertebrates are extremely sensitive to LPS, but lower vertebrates, like fish, are resistant to their systemic toxic effects. However, the effects of LPS on the fish intestinal mucosa remain unknown. Edwardsiella ictaluri is a primitive member of the Enterobacteriaceae family that causes enteric septicemia in channel catfish (Ictalurus punctatus). E. ictaluri infects and colonizes deep lymphoid tissues upon oral or immersion infection. Both gut and olfactory organs are the primary sites of invasion. At the systemic level, E. ictaluri pathogenesis is relatively well characterized, but our knowledge about E. ictaluri intestinal interaction is limited. Recently, we observed that E. ictaluri oligo-polysaccharide (O-PS) LPS mutants have differential effects on the intestinal epithelia of orally inoculated catfish. Here we evaluate the effects of E. ictaluri O-PS LPS mutants by using a novel catfish intestinal loop model and compare it to the rabbit ileal loop model inoculated with Salmonella enterica serovar Typhimurium LPS. We found evident differences in rabbit ileal loop and catfish ileal loop responses to E. ictaluri and S. Typhimurium LPS. We determined that catfish respond to E. ictaluri LPS but not to S. Typhimurium LPS. We also determined that E. ictaluri inhibits cytokine production and induces disruption of the intestinal fish epithelia in an O-PS-dependent fashion. The E. ictaluri wild type and ΔwibT LPS mutant caused intestinal tissue damage and inhibited proinflammatory cytokine synthesis, in contrast to E. ictaluri Δgne and Δugd LPS mutants. We concluded that the E. ictaluri O-PS subunits play a major role during pathogenesis, since they influence the recognition of the LPS by the intestinal mucosal immune system of the catfish. The LPS structure of E. ictaluri mutants is needed to understand the mechanism of interaction.

Background: Lizards are evolutionarily the most closely related vertebrates to humans that can lose and regrow an entire appendage. Regeneration in lizards involves differential expression of hundreds of genes that regulate wound healing, musculoskeletal development, hormonal response, and embryonic morphogenesis. While microRNAs are able to regulate large groups of genes, their role in lizard regeneration has not been investigated.

Results: MicroRNA sequencing of green anole lizard (Anolis carolinensis) regenerating tail and associated tissues revealed 350 putative novel and 196 known microRNA precursors. Eleven microRNAs were differentially expressed between the regenerating tail tip and base during maximum outgrowth (25 days post autotomy), including miR-133a, miR-133b, and miR-206, which have been reported to regulate regeneration and stem cell proliferation in other model systems. Three putative novel differentially expressed microRNAs were identified in the regenerating tail tip.

Conclusions: Differentially expressed microRNAs were identified in the regenerating lizard tail, including known regulators of stem cell proliferation. The identification of 3 putative novel microRNAs suggests that regulatory networks, either conserved in vertebrates and previously uncharacterized or specific to lizards, are involved in regeneration. These findings suggest that differential regulation of microRNAs may play a role in coordinating the timing and expression of hundreds of genes involved in regeneration.

Introduction: Nutrient availability, assimilation, and allocation can have important and lasting effects on the immune system development of growing animals. Though carotenoid pigments have immunostimulatory properties in many animals, relatively little is known regarding how they influence the immune system during development. Moreover, studies linking carotenoids to health at any life stage have largely been restricted to birds and mammals. We investigated the effects of carotenoid supplementation on multiple aspects of immunity in juvenile veiled chameleons (Chamaeleo calyptratus). We supplemented half of the chameleons with lutein (a xanthophyll carotenoid) for 14 weeks during development and serially measured multiple aspects of immune function, including: agglutination and lysis performance of plasma, wound healing, and plasma nitric oxide concentrations before and after wounding.

Results: Though lutein supplementation effectively elevated circulating carotenoid concentrations throughout the developmental period, we found no evidence that carotenoid repletion enhanced immune function at any point. However, agglutination and lysis scores increased, while baseline nitric oxide levels decreased, as chameleons aged.

Conclusions: Taken together, our results indicate that body mass and age, but not carotenoid access, may play an important role in immune performance of growing chameleons. Hence, studying well-understood physiological processes in novel taxa can provide new perspectives on alternative physiological processes and nutrient function.

Introduction: Immune function is a vital physiological process that is often suppressed during times of resource scarcity due to investments in other physiological systems. While energy is the typical currency that has been examined in such trade-offs, limitations of other resources may similarly lead to trade-offs that affect immune function. Specifically, water is a critical resource with profound implications for organismal ecology, yet its availability can fluctuate at local, regional, and even global levels. Despite this, the effect of osmotic state on immune function has received little attention.

Results: Using agglutination and lysis assays as measures of an organism's plasma concentration of natural antibodies and capacity for foreign cell destruction, respectively, we tested the independent effects of osmotic state, digestive state, and energy balance on innate immune function in free-ranging and laboratory populations of the Gila monster, Heloderma suspectum. This desert-dwelling lizard experiences dehydration and energy resource fluctuations on a seasonal basis. Dehydration was expected to decrease innate immune function, yet we found that dehydration increased lysis and agglutination abilities in both lab and field studies, a relationship that was not simply an effect of an increased concentration of immune molecules. Laboratory-based differences in digestive state were not associated with lysis or agglutination metrics, although in our field population, a loss of fat stores was correlated with an increase in lysis.

Conclusions: Depending on the life history of an organism, osmotic state may have a greater influence on immune function than energy availability. Thus, consideration of osmotic state as a factor influencing immune function will likely improve our understanding of ecoimmunology and the disease dynamics of a wide range of species.

Lizards, which are amniote vertebrates like humans, are able to lose and regenerate a functional tail. Understanding the molecular basis of this process would advance regenerative approaches in amniotes, including humans. We have carried out the first transcriptomic analysis of tail regeneration in a lizard, the green anole Anolis carolinensis, which revealed 326 differentially expressed genes activating multiple developmental and repair mechanisms. Specifically, genes involved in wound response, hormonal regulation, musculoskeletal development, and the Wnt and MAPK/FGF pathways were differentially expressed along the regenerating tail axis. Furthermore, we identified 2 microRNA precursor families, 22 unclassified non-coding RNAs, and 3 novel protein-coding genes significantly enriched in the regenerating tail. However, high levels of progenitor/stem cell markers were not observed in any region of the regenerating tail. Furthermore, we observed multiple tissue-type specific clusters of proliferating cells along the regenerating tail, not localized to the tail tip. These findings predict a different mechanism of regeneration in the lizard than the blastema model described in the salamander and the zebrafish, which are anamniote vertebrates. Thus, lizard tail regrowth involves the activation of conserved developmental and wound response pathways, which are potential targets for regenerative medical therapies.

Agassiz’s desert tortoise (Gopherus agassizii) is a long-lived species native to the Mojave Desert and is listed as threatened under the US Endangered Species Act. To aid conservation efforts for preserving the genetic diversity of this species, we generated a whole genome reference sequence with an annotation based on deep transcriptome sequences of adult skeletal muscle, lung, brain, and blood. The draft genome assembly for G. agassizii has a scaffold N50 length of 252 kbp and a total length of 2.4 Gbp. Genome annotation reveals 20,172 protein-coding genes in the G. agassizii assembly, and that gene structure is more similar to chicken than other turtles. We provide a series of comparative analyses demonstrating (1) that turtles are among the slowest-evolving genome-enabled reptiles, (2) amino acid changes in genes controlling desert tortoise traits such as shell development, longevity and osmoregulation, and (3) fixed variants across the Gopherus species complex in genes related to desert adaptations, including circadian rhythm and innate immune response. This G. agassizii genome reference and annotation is the first such resource for any tortoise, and will serve as a foundation for future analysis of the genetic basis of adaptations to the desert environment, allow for investigation into genomic factors affecting tortoise health, disease and longevity, and serve as a valuable resource for additional studies in this species complex.