Species and population specific gene expression in blood transcriptomes of marine turtles

Shreya M. Banjeree; Jamie Adkins Stoll; Camryn D. Allen; Jennifer M. Lynch; Heather S. Harris; Lauren Kenyon; Richard E. Connon; Eleanor J. Sterling; Eugenia Naro-Maciel; Kathryn McFadden; Margaret Lamont; James Benge; Nadia B. Fernandez; Jeffrey A. Seminoff; Scott R. Benson; Rebecca L. Lewison; Tomoharu Eguchi; Tammy M. Summers; Jessy R. Hapdei; Marc R. Rice; Summer Martin; T. Todd Jones; Peter H. Dutton; George H. Balazs; Lisa M. Komoroske

doi:10.1186/s12864-021-07656-5

Home/Archive/USGS

USGS · 70229231

Species and population specific gene expression in blood transcriptomes of marine turtles

View original source ↗Export metadata ↓

Abstract

Background Transcriptomic data has demonstrated utility to advance the study of physiological diversity and organisms’ responses to environmental stressors. However, a lack of genomic resources and challenges associated with collecting high-quality RNA can limit its application for many wild populations. Minimally invasive blood sampling combined with de novo transcriptomic approaches has great potential to alleviate these barriers. Here, we advance these goals for marine turtles by generating high quality de novo blood transcriptome assemblies to characterize functional diversity and compare global transcriptional profiles between tissues, species, and foraging aggregations. Results We generated high quality blood transcriptome assemblies for hawksbill ( Eretmochelys imbricata ) , loggerhead ( Caretta caretta ), green ( Chelonia mydas ), and leatherback ( Dermochelys coriacea ) turtles. The functional diversity in assembled blood transcriptomes was comparable to those from more traditionally sampled tissues. A total of 31.3% of orthogroups identified were present in all four species, representing a core set of conserved genes expressed in blood and shared across marine turtle species. We observed strong species-specific expression of these genes, as well as distinct transcriptomic profiles between green turtle foraging aggregations that inhabit areas of greater or lesser anthropogenic disturbance. Conclusions Obtaining global gene expression data through non-lethal, minimally invasive sampling can greatly expand the applications of RNA-sequencing in protected long-lived species such as marine turtles. The distinct differences in gene expression signatures between species and foraging aggregations provide insight into the functional genomics underlying the diversity in this ancient vertebrate lineage. The transcriptomic resources generated here can be used in further studies examining the evolutionary ecology and anthropogenic impacts on marine turtles.

Explore related subjects

BMC Genomics Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Shreya M. Banjeree, Jamie Adkins Stoll, Camryn D. Allen, Jennifer M. Lynch, Heather S. Harris, Lauren Kenyon, Richard E. Connon, Eleanor J. Sterling, Eugenia Naro-Maciel, Kathryn McFadden, Margaret Lamont, James Benge, Nadia B. Fernandez, Jeffrey A. Seminoff, Scott R. Benson, Rebecca L. Lewison, Tomoharu Eguchi, Tammy M. Summers, Jessy R. Hapdei, Marc R. Rice, Summer Martin, T. Todd Jones, Peter H. Dutton, George H. Balazs, Lisa M. Komoroske. 2021-05-13. Species and population specific gene expression in blood transcriptomes of marine turtles. https://doi.org/10.1186/s12864-021-07656-5

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2020

An empirical comparison of population genetic analyses using microsatellite and SNP data for a species of conservation concern

Background Use of genomic tools to characterize wildlife populations has increased in recent years. In the past, genetic characterization has been accomplished with more traditional genetic tools (e.g., microsatellites). The explosion of genomic methods and the subsequent creation of large SNP datasets has led to the promise of increased precision in population genetic parameter estimates and identification of demographically and evolutionarily independent groups, as well as questions about the future usefulness of the more traditional genetic tools. At present, few empirical comparisons of population genetic parameters and clustering analyses performed with microsatellites and SNPs have been conducted. Results Here we used microsatellite and SNP data generated from Gunnison sage-grouse ( Centrocercus minimus ) samples to evaluate concordance of the results obtained from each dataset for common metrics of genetic diversity ( H O , H E , F IS , A R ) and differentiation ( F ST , G ST , D Jost ). Additionally, we evaluated clustering of individuals using putatively neutral (SNPs and microsatellites), putatively adaptive, and a combined dataset of putatively neutral and adaptive loci. We took particular interest in the conservation implications of any differences. Generally, we found high concordance between microsatellites and SNPs for H E , F IS , A R , and all differentiation estimates. Although there was strong correlation between metrics from SNPs and microsatellites, the magnitude of the diversity and differentiation metrics were quite different in some cases. Clustering analyses also showed similar patterns, though SNP data was able to cluster individuals into more distinct groups. Importantly, clustering analyses with SNP data suggest strong demographic independence among the six distinct populations of Gunnison sage-grouse with some indication of evolutionary independence in two or three populations; a finding that was not revealed by microsatellite data. Conclusion We demonstrate that SNPs have three main advantages over microsatellites: more precise estimates of population-level diversity, higher power to identify groups in clustering methods, and the ability to consider local adaptation. This study adds to a growing body of work comparing the use of SNPs and microsatellites to evaluate genetic diversity and differentiation for a species of conservation concern with relatively high population structure and using the most common method of obtaining SNP genotypes for non-model organisms.

Arizona, Colorado, New Mexico, Utah↗

USGS2019

Effect of growth rate on transcriptomic responses to immune stimulation in wild-type, domesticated, and GH-transgenic coho salmon

Background Transcriptomic responses to immune stimulation were investigated in coho salmon (Oncorhynchus kisutch) with distinct growth phenotypes. Wild-type fish were contrasted to strains with accelerated growth arising either from selective breeding (i.e. domestication) or genetic modification. Such distinct routes to accelerated growth may have unique implications for relationships and/or trade-offs between growth and immune function. Results RNA-Seq was performed on liver and head kidney in four ‘growth response groups’ injected with polyinosinic-polycytidylic acid (Poly I:C; viral mimic), peptidoglycan (PGN; bacterial mimic) or PBS (control). These groups were: 1) ‘W’: wild-type, 2) ‘TF’: growth hormone (GH) transgenic salmon with ~ 3-fold higher growth-rate than W, 3) ‘TR’: GH transgenic fish ration restricted to possess a growth-rate equal to W, and 4) ‘D’: domesticated non-transgenic fish showing growth-rate intermediate to W and TF. D and TF showed a higher similarity in transcriptomic response compared to W and TR. Several immune genes showed constitutive expression differences among growth response groups, including perforin 1 and C-C motif chemokine 19-like. Among the affected immune pathways, most were up-regulated by Poly I:C and PGN. In response to PGN, the c-type lectin receptor signalling pathway responded uniquely in TF and TR. In response to stimulation with both immune mimics, TR responded more strongly than other groups. Further, group-specific pathway responses to PGN stimulation included NOD-like receptor signalling in W and platelet activation in TR. TF consistently showed the most attenuated immune response relative to W, and more DEGs were apparent in TR than TF and D relative to W, suggesting that a non-satiating ration coupled with elevated circulating GH levels may cause TR to possess enhanced immune capabilities. Alternatively, TF and D salmon are prevented from acquiring the same level of immune response as TR due to direction of energy to high overall somatic growth. Further study of the effects of ration restriction in growth-modified fishes is warranted. Conclusions These findings improve our understanding of the pleiotropic effects of growth modification on the immunological responses of fish, revealing unique immune pathway responses depending on the mechanism of growth acceleration and nutritional availability.

BMC Genomics↗

USGS2018

New insights into the phylogenetics and population structure of the prairie falcon (Falco mexicanus)

Background Management requires a robust understanding of between- and within-species genetic variability, however such data are still lacking in many species. For example, although multiple population genetics studies of the peregrine falcon ( Falco peregrinus ) have been conducted, no similar studies have been done of the closely-related prairie falcon ( F. mexicanus ) and it is unclear how much genetic variation and population structure exists across the species’ range. Furthermore, the phylogenetic relationship of F. mexicanus relative to other falcon species is contested. We utilized a genomics approach (i.e., genome sequencing and assembly followed by single nucleotide polymorphism genotyping) to rapidly address these gaps in knowledge. Results We sequenced the genome of a single female prairie falcon and generated a 1.17 Gb (gigabases) draft genome assembly. We generated maximum likelihood phylogenetic trees using complete mitochondrial genomes as well as nuclear protein-coding genes. This process provided evidence that F. mexicanus is an outgroup to the clade that includes the peregrine falcon and members of the subgenus Hierofalco. We annotated > 16,000 genes and almost 600,000 high-quality single nucleotide polymorphisms (SNPs) in the nuclear genome, providing the raw material for a SNP assay design featuring > 140 gene-associated markers and a molecular-sexing marker. We subsequently genotyped ~ 100 individuals from California (including the San Francisco East Bay Area, Pinnacles National Park and the Mojave Desert) and Idaho (Snake River Birds of Prey National Conservation Area). We tested for population structure and found evidence that individuals sampled in California and Idaho represent a single panmictic population. Conclusions Our study illustrates how genomic resources can rapidly shed light on genetic variability in understudied species and resolve phylogenetic relationships. Furthermore, we found evidence of a single, randomly mating population of prairie falcons across our sampling locations. Prairie falcons are highly mobile and relatively rare long-distance dispersal events may promote gene flow throughout the range. As such, California’s prairie falcons might be managed as a single population, indicating that management actions undertaken to benefit the species at the local level have the potential to influence the species as a whole.

BMC Genomics↗