Transcriptome of the dead: characterisation of immune genes and marker development from necropsy samples in a free-ranging marine mammal

Transcriptome analysis of cadmium exposure in kidney fibroblast cells of the North Atlantic Right Whale (Eubalaena glacialis)

An 8X15k oligonucleotide microarray was developed consisting of 2334 Eubalaena glacialis probes and 2166 Tursiops truncatus probes and used to measure the effects, at transcriptomic level, of cadmium exposure in right whale kidney fibroblast cells. Cells were exposed to three concentrations (1 μM, 0.1 μM, and 0.01 μM) of cadmium chloride (CdCl₂) for three exposure times (1, 4, and 24 h). Cells exposed to 1 μM CdCl₂ for 4 h and 24 h showed upregulated genes involved in protection from metal toxicity and oxidative stress, protein renaturation, apoptosis inhibition, as well as several regulators of cellular processes. Downregulated genes represented a suite of functions including cell proliferation, transcription regulation, actin polymerization, and stress fiber synthesis. The collection of differentially expressed genes in this study support proposed mechanisms of cadmium-induced apoptosis such as ubiquitin proteasome system disruption, Ca²⁺ homeostasis interference, mitochondrial membrane potential collapse, reactive oxygen species (ROS) production, and cell cycle arrest. The results also have confirmed the right whale microarray as a reproducible tool in measuring differentiated gene expression that could be a valuable asset for transcriptome analysis of other baleen whales and potential health assessment protocols.

An 85K SNP Array Uncovers Inbreeding and Cryptic Relatedness in an Antarctic Fur Seal Breeding Colony

High density single nucleotide polymorphism (SNP) arrays allow large numbers of individuals to be rapidly and cost-effectively genotyped at large numbers of genetic markers. However, despite being widely used in studies of humans and domesticated plants and animals, SNP arrays are lacking for most wild organisms. We developed a custom 85K Affymetrix Axiom array for an intensively studied pinniped, the Antarctic fur seal (Arctocephalus gazella). SNPs were discovered from a combination of genomic and transcriptomic resources and filtered according to strict criteria. Out of a total of 85,359 SNPs tiled on the array, 75,601 (88.6%) successfully converted and were polymorphic in 270 animals from a breeding colony at Bird Island in South Georgia. Evidence was found for inbreeding, with three genomic inbreeding coefficients being strongly intercorrelated and the proportion of the genome in runs of homozygosity being non-zero in all individuals. Furthermore, analysis of genomic relatedness coefficients identified previously unknown first-degree relatives and multiple second-degree relatives among a sample of ostensibly unrelated individuals. Such “cryptic relatedness” within fur seal breeding colonies may increase the likelihood of consanguineous matings and could therefore have implications for understanding fitness variation and mate choice. Finally, we demonstrate the cross-amplification potential of the array in three related pinniped species. Overall, our SNP array will facilitate future studies of Antarctic fur seals and has the potential to serve as a more general resource for the wider pinniped research community.

Transcriptome analysis revealed potential mechanisms of differences in physiological stress responses between caged male and female magpies

Background

Under caged conditions, birds are affected more severely by environmental stressors such as dietary structure, activity space, human disturbances, and pathogens, which may be reflected in the gene expression in peripheral blood or other tissues. Elucidating the molecular mechanism of these stress responses will help improve animal welfare.

Results

In the present study, the blood transcriptomes of six male and five female caged magpies (Pica pica) were sequenced, and a total of ~ 100 Gb in clean reads were generated using the Illumina HiSeq 2000 sequencer. A total of 420,291 unigenes were identified after assembly, of which 179,316 were annotated in five databases, 7471 were assigned to 269 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and 566 were assigned to the Clusters of Orthologous Groups (COG) functional classification “defense mechanisms”. Analysis of differentially expressed genes (DEGs) showed that 2657 unigenes were differentially expressed between males and females (q < 0.1), and these DEGs were assigned to 45 KEGG pathways involving stress resistance, immunity, energy metabolism, reproduction, lifespan regulation, and diseases. Further analysis revealed that females might be more sensitive to stress through upregulation of c-Jun N-terminal kinases (JNKs) and 5’AMP-activated protein kinase (AMPK), and were also possibly more sensitive to dynamic changes in energy. Females expressed higher major histocompatibility complex (MHC) class II levels than males, enhancing resistance to pathogens, and the DEGs related to reproduction included MAPK, CaMK, CPEB, and Cdc25. The genes related to stress, energy, and immunity were also likely related to the regulation of longevity. The upregulated JNKs in females might prolong lifespan and relieve antioxidant stress. Females may also activate the AMPK pathway and implement dietary restrictions to prolong lifespan, whereas males may upregulate SIRT1 and CRAB to increase lifespan.

Conclusions

Female magpies might be more sensitive to stress and dynamic changes in energy thus enhanced resistance to pathogens, and the genes related to stress, energy, and immunity were also possibly related to the regulation of longevity. Further confirmations with techniques such as RT-qPCR and western blot are necessary to validate the above arguments.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5804-0) contains supplementary material, which is available to authorized users.

Experimental validation of in silico predicted RAD locus frequencies using genomic resources and short read data from a model marine mammal

BACKGROUND

Restriction site-associated DNA sequencing (RADseq) has revolutionized the study of wild organisms by allowing cost-effective genotyping of thousands of loci. However, for species lacking reference genomes, it can be challenging to select the restriction enzyme that offers the best balance between the number of obtained RAD loci and depth of coverage, which is crucial for a successful outcome. To address this issue, PredRAD was recently developed, which uses probabilistic models to predict restriction site frequencies from a transcriptome assembly or other sequence resource based on either GC content or mono-, di- or trinucleotide composition. This program generates predictions that are broadly consistent with estimates of the true number of restriction sites obtained through in silico digestion of available reference genome assemblies. However, in practice the actual number of loci obtained could potentially differ as incomplete enzymatic digestion or patchy sequence coverage across the genome might lead to some loci not being represented in a RAD dataset, while erroneous assembly could potentially inflate the number of loci. To investigate this, we used genome and transcriptome assemblies together with RADseq data from the Antarctic fur seal (Arctocephalus gazella) to compare PredRAD predictions with empirical estimates of the number of loci obtained via in silico digestion and from de novo assemblies.

RESULTS

PredRAD yielded consistently higher predicted numbers of restriction sites for the transcriptome assembly relative to the genome assembly. The trinucleotide and dinucleotide models also predicted higher frequencies than the mononucleotide or GC content models. Overall, the dinucleotide and trinucleotide models applied to the transcriptome and the genome assemblies respectively generated predictions that were closest to the number of restriction sites estimated by in silico digestion. Furthermore, the number of de novo assembled RAD loci mapping to restriction sites was similar to the expectation based on in silico digestion.

CONCLUSIONS

Our study reveals generally high concordance between PredRAD predictions and empirical estimates of the number of RAD loci. This further supports the utility of PredRAD, while also suggesting that it may be feasible to sequence and assemble the majority of RAD loci present in an organism's genome.

De novo transcriptome assembly and RNA-Seq expression analysis in blood from beluga whales of Bristol Bay, AK

Assessing the health of marine mammal sentinel species is crucial to understanding the impacts of environmental perturbations on marine ecosystems and human health. In Arctic regions, beluga whales, Delphinapterus leucas, are upper level predators that may serve as a sentinel species, potentially forecasting impacts on human health. While gene expression profiling from blood transcriptomes has widely been used to assess health status and environmental exposures in human and veterinary medicine, its use in wildlife has been limited due to the lack of available genomes and baseline data. To this end we constructed the first beluga whale blood transcriptome de novo from samples collected during annual health assessments of the healthy Bristol Bay, AK stock during 2012-2014 to establish baseline information on the content and variation of the beluga whale blood transcriptome. The Trinity transcriptome assembly from beluga was comprised of 91,325 transcripts that represented a wide array of cellular functions and processes and was extremely similar in content to the blood transcriptome of another cetacean, the bottlenose dolphin. Expression of hemoglobin transcripts was much lower in beluga (25.6% of TPM, transcripts per million) than has been observed in many other mammals. A T12A amino acid substitution in the HBB sequence of beluga whales, but not bottlenose dolphins, was identified and may play a role in low temperature adaptation. The beluga blood transcriptome was extremely stable between sex and year, with no apparent clustering of samples by principle components analysis and <4% of genes differentially expressed (EBseq, FDR<0.05). While the impacts of season, sexual maturity, disease, and geography on the beluga blood transcriptome must be established, the presence of transcripts involved in stress, detoxification, and immune functions indicate that blood gene expression analyses may provide information on health status and exposure. This study provides a wealth of transcriptomic data on beluga whales and provides a sizeable pool of preliminary data for comparison with other studies in beluga whale.

De novo assembly of the ringed seal (Pusa hispida) blubber transcriptome: A tool that enables identification of molecular health indicators associated with PCB exposure

The ringed seal, Pusa hispida, is a keystone species in the Arctic marine ecosystem, and is proving a useful marine mammal for linking polychlorinated biphenyl (PCB) exposure to toxic injury. We report here the first de novo assembled transcriptome for the ringed seal (342,863 transcripts, of which 53% were annotated), which we then applied to a population of ringed seals exposed to a local PCB source in Arctic Labrador, Canada. We found an indication of energy metabolism imbalance in local ringed seals (n=4), and identified five significant gene transcript targets: plasminogen receptor (Plg-R(KT)), solute carrier family 25 member 43 receptor (Slc25a43), ankyrin repeat domain-containing protein 26-like receptor (Ankrd26), HIS30 (not yet annotated) and HIS16 (not yet annotated) that may represent indicators of PCB exposure and effects in marine mammals. The abundance profiles of these five gene targets were validated in blubber samples collected from 43 ringed seals using a qPCR assay. The mRNA transcript levels for all five gene targets, (Plg-R(KT), r²=0.43), (Slc25a43, r²=0.51), (Ankrd26, r²=0.43), (HIS30, r²=0.39) and (HIS16, r²=0.31) correlated with increasing levels of blubber PCBs. Results from the present study contribute to our understanding of PCB associated effects in marine mammals, and provide new tools for future molecular and toxicology work in pinnipeds.

Brain transcriptomes of harbor seals demonstrate gene expression patterns of animals undergoing a metabolic disease and a viral infection

Diseases of marine mammals can be difficult to diagnose because of their life history and protected status. Stranded marine mammals have been a particularly useful resource to discover and comprehend the diseases that plague these top predators. Additionally, advancements in high-throughput sequencing (HTS) has contributed to the discovery of novel pathogens in marine mammals. In this study, we use a combination of HTS and stranded harbor seals (Phoca vitulina) to better understand a known and unknown brain disease. To do this, we used transcriptomics to evaluate brain tissues from seven neonatal harbor seals that expired from an unknown cause of death (UCD) and compared them to four neonatal harbor seals that had confirmed phocine herpesvirus (PhV-1) infections in the brain. Comparing the two disease states we found that UCD animals showed a significant abundance of fatty acid metabolic transcripts in their brain tissue, thus we speculate that a fatty acid metabolic dysregulation contributed to the death of these animals. Furthermore, we were able to describe the response of four young harbor seals with PhV-1 infections in the brain. PhV-1 infected animals showed a significant ability to mount an innate and adaptive immune response, especially to combat viral infections. Our data also suggests that PhV-1 can hijack host pathways for DNA packaging and exocytosis. This is the first study to use transcriptomics in marine mammals to understand host and viral interactions and assess the death of stranded marine mammals with an unknown disease. Furthermore, we show the value of applying transcriptomics on stranded marine mammals for disease characterization.

Transcriptomic SNP discovery for custom genotyping arrays: impacts of sequence data, SNP calling method and genotyping technology on the probability of validation success

BACKGROUND

Single nucleotide polymorphism (SNP) discovery is an important goal of many studies. However, the number of 'putative' SNPs discovered from a sequence resource may not provide a reliable indication of the number that will successfully validate with a given genotyping technology. For this it may be necessary to account for factors such as the method used for SNP discovery and the type of sequence data from which it originates, suitability of the SNP flanking sequences for probe design, and genomic context. To explore the relative importance of these and other factors, we used Illumina sequencing to augment an existing Roche 454 transcriptome assembly for the Antarctic fur seal (Arctocephalus gazella). We then mapped the raw Illumina reads to the new hybrid transcriptome using BWA and BOWTIE2 before calling SNPs with GATK. The resulting markers were pooled with two existing sets of SNPs called from the original 454 assembly using NEWBLER and SWAP454. Finally, we explored the extent to which SNPs discovered using these four methods overlapped and predicted the corresponding validation outcomes for both Illumina Infinium iSelect HD and Affymetrix Axiom arrays.

RESULTS

Collating markers across all discovery methods resulted in a global list of 34,718 SNPs. However, concordance between the methods was surprisingly poor, with only 51.0 % of SNPs being discovered by more than one method and 13.5 % being called from both the 454 and Illumina datasets. Using a predictive modeling approach, we could also show that SNPs called from the Illumina data were on average more likely to successfully validate, as were SNPs called by more than one method. Above and beyond this pattern, predicted validation outcomes were also consistently better for Affymetrix Axiom arrays.

CONCLUSIONS

Our results suggest that focusing on SNPs called by more than one method could potentially improve validation outcomes. They also highlight possible differences between alternative genotyping technologies that could be explored in future studies of non-model organisms.

Genomic Methods Take the Plunge: Recent Advances in High-Throughput Sequencing of Marine Mammals

The dramatic increase in the application of genomic techniques to non-model organisms (NMOs) over the past decade has yielded numerous valuable contributions to evolutionary biology and ecology, many of which would not have been possible with traditional genetic markers. We review this recent progression with a particular focus on genomic studies of marine mammals, a group of taxa that represent key macroevolutionary transitions from terrestrial to marine environments and for which available genomic resources have recently undergone notable rapid growth. Genomic studies of NMOs utilize an expanding range of approaches, including whole genome sequencing, restriction site-associated DNA sequencing, array-based sequencing of single nucleotide polymorphisms and target sequence probes (e.g., exomes), and transcriptome sequencing. These approaches generate different types and quantities of data, and many can be applied with limited or no prior genomic resources, thus overcoming one traditional limitation of research on NMOs. Within marine mammals, such studies have thus far yielded significant contributions to the fields of phylogenomics and comparative genomics, as well as enabled investigations of fitness, demography, and population structure. Here we review the primary options for generating genomic data, introduce several emerging techniques, and discuss the suitability of each approach for different applications in the study of NMOs.

Sixteen kiwi (Apteryx spp) transcriptomes provide a wealth of genetic markers and insight into sex chromosome evolution in birds

Background

Kiwi represent the most basal extant avian lineage (paleognaths) and exhibit biological attributes that are unusual or extreme among living birds, such as large egg size, strong olfaction, nocturnality, flightlessness and long lifespan. Despite intense interest in their evolution and their threatened status, genomic resources for kiwi were virtually non-existent until the recent publication of a single genome. Here we present the most comprehensive kiwi transcriptomes to date, obtained via Illumina sequencing of whole blood and de novo assembly of mRNA sequences of eight individuals from each of the two rarest kiwi species, little spotted kiwi (LSK; Apteryx owenii) and rowi (A. rowi).

Results

Sequences obtained were orthologous with a wide diversity of functional genes despite the sequencing of a single tissue type. Individual and composite assemblies contain more than 7900 unique protein coding transcripts in each of LSK and rowi that show strong homology with chicken (Gallus gallus), including those associated with growth, development, disease resistance, reproduction and behavior. The assemblies also contain 66,909 SNPs that distinguish between LSK and rowi, 12,384 SNPs among LSK (associated with 3088 genes), and 29,313 SNPs among rowi (associated with 4953 genes). We found 3084 transcripts differentially expressed between LSK and rowi and 150 transcripts differentially expressed between the sexes. Of the latter, 83 could be mapped to chicken chromosomes with 95% syntenic with chromosome Z.

Conclusions

Our study has simultaneously sequenced multiple species, sexes, and individual kiwi at thousands of genes, and thus represents a significant leap forward in genomic resources available for kiwi. The expression pattern we observed among chromosome Z related genes in kiwi is similar to that observed in ostriches and emu, suggesting a common and ancestral pattern of sex chromosome homomorphy, recombination, and gene dosage among living paleognaths. The transcriptome assemblies described here will provide a rich resource for polymorphic marker development and studies of adaptation of these highly unusual and endangered birds.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2714-2) contains supplementary material, which is available to authorized users.

De novo assembly of the dual transcriptomes of a polymorphic raptor species and its malarial parasite

BACKGROUND

Studies of non-model species are important for understanding the molecular processes underpinning phenotypic variation under natural ecological conditions. The common buzzard (Buteo buteo; Aves: Accipitriformes) is a widespread and common Eurasian raptor with three distinct plumage morphs that differ in several fitness-related traits, including parasite infestation. To provide a genomic resource for plumage polymorphic birds in general and to search for candidate genes relating to fitness, we generated a transcriptome from a single dead buzzard specimen plus easily accessible, minimally invasive samples from live chicks.

RESULTS

We not only de novo assembled a near-complete buzzard transcriptome, but also obtained a significant fraction of the transcriptome of its malaria-like parasite, Leucocytozoon buteonis. By identifying melanogenesis-related transcripts that are differentially expressed in light ventral and dark dorsal feathers, but which are also expressed in other regions of the body, we also identified a suite of candidate genes that could be associated with fitness differences among the morphs. These include several immune-related genes, providing a plausible link between melanisation and parasite load. qPCR analysis of a subset of these genes revealed significant differences between ventral and dorsal feathers and an additional effect of morph.

CONCLUSION

This new resource provides preliminary insights into genes that could be involved in fitness differences between the buzzard colour morphs, and should facilitate future studies of raptors and their malaria-like parasites.

The degradation of proteins in pinniped skeletal muscle: viability of post-mortem tissue in physiological research

As marine divers, pinnipeds have a high capacity for exercise at depth while holding their breath. With finite access to oxygen, these species need to be capable of extended aerobic exercise and conservation of energy. Pinnipeds must deal with common physiological hurdles, such as hypoxia, exhaustion and acidosis, that are common to all exercising mammals. The physiological mechanisms in marine mammals used for managing oxygen and carbon dioxide have sparked much research, but access to animals and tissues is difficult and requires permits. Deceased animals that are either bycaught or stranded provide one potential source for tissues, but the validity of biochemical data from post-mortem samples has not been rigorously assessed. Tissues collected from stranded diving mammals may be a crucial source to add to our limited knowledge on the physiology of some of these animals and important to the conservation and management of these species. We aim to determine the reliability of biochemical assays derived from post-mortem tissue and to promote the immediate sampling of stranded animals for the purpose of physiological research. In this study, we mapped the temporal degradation of muscle enzymes from biopsied Northern elephant seals (Mirounga angustirostris) and highlight recommendations for storage protocols for the best preservation of tissue. We also compared the enzymatic activity of different muscle groups (pectoral and latissimus dorsi) in relation to locomotion and measured the effects of four freeze-thaw cycles on muscle tissue enzyme function. Results indicate that enzymatic activity fluctuates greatly, especially with varying storage temperature, storage time, species and muscle group being assayed. In contrast, proteins, such as myoglobin, remain relatively continuous in their increase at 4°C for 48 h. Stranded animals can be a valuable source of biochemical data, but enzyme assays should be used only with great caution in post-mortem tissues.

Transcriptome analysis of northern elephant seal (Mirounga angustirostris) muscle tissue provides a novel molecular resource and physiological insights

Background

The northern elephant seal, Mirounga angustirostris, is a valuable animal model of fasting adaptation and hypoxic stress tolerance. However, no reference sequence is currently available for this and many other marine mammal study systems, hindering molecular understanding of marine adaptations and unique physiology.

Results

We sequenced a transcriptome of M. angustirostris derived from muscle sampled during an acute stress challenge experiment to identify species-specific markers of stress axis activation and recovery. De novo assembly generated 164,966 contigs and a total of 522,699 transcripts, of which 68.70% were annotated using mouse, human, and domestic dog reference protein sequences. To reduce transcript redundancy, we removed highly similar isoforms in large gene families and produced a filtered assembly containing 336,657 transcripts. We found that a large number of annotated genes are associated with metabolic signaling, immune and stress responses, and muscle function. Preliminary differential expression analysis suggests a limited transcriptional response to acute stress involving alterations in metabolic and immune pathways and muscle tissue maintenance, potentially driven by early response transcription factors such as Cebpd.

Conclusions

We present the first reference sequence for Mirounga angustirostris produced by RNA sequencing of muscle tissue and cloud-based de novo transcriptome assembly. We annotated 395,102 transcripts, some of which may be novel isoforms, and have identified thousands of genes involved in key physiological processes. This resource provides elephant seal-specific gene sequences, complementing existing metabolite and protein expression studies and enabling future work on molecular pathways regulating adaptations such as fasting, hypoxia, and environmental stress responses in marine mammals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1253-6) contains supplementary material, which is available to authorized users.

Comparative transcriptomic analyses revealed divergences of two agriculturally important aphid species

BACKGROUND

Grain aphid (Sitobion avenae F) and pea aphid (Acyrthosiphon pisum) are two agriculturally important pest species, which cause significant yield losses to crop plants each year by inflicting damage both through the direct effects of feeding and by vectoring debilitating plant viruses. Although a close phylogenetic relationship between grain aphid and pea aphid was proposed, the biological variations between these two aphid species are obvious. While the host ranges of grain aphid is restricted to cereal crops and in particular wheat, that of pea aphid is wider, mainly colonizing leguminous plant species. Until now, the genetic factors underlying the divergence between grain aphid and pea aphid still remain unclear due to the limited genomic data of grain aphid available in public databases.

RESULTS

Based on a set of transcriptome data of grain aphid generated by using Roche 454 GS-FLX pyrosequencing, comparative analysis between this set of transcriptome data of grain aphid and mRNA sequences of pea aphid available in the public databases was performed. Compared with mRNA sequences of pea aphid, 4,857 unigenes were found to be specifically presented in the transcriptome of grain aphid under the rearing conditions described in this study. Furthermore, 3,368 orthologous pairs which could be calculated with both nonsynonymous (Ka) and synonymous (Ks) substitutions were used to infer their sequence divergences. The average differences in the coding, 5' and 3' untranslated regions of these orthologs were 10.53%, 21.29% and 18.96%, respectively. Moreover, of 340 orthologs which were identified to have evolved in response to positive selection based on the rates of Ka and Ks substitutions, 186 were predicted to be involved in secondary metabolism and xenobiotic metabolisms which might contribute to the divergence of these two aphid species.

CONCLUSIONS

The comprehensive transcriptome divergent sequence analysis between grain aphid and pea aphid provides an invaluable resource for the investigation of genes involved in host plant adaptation and evolution. Moreover, the demonstration of divergent transcriptome sequences between grain aphid and pea aphid pave the way for the investigation of the molecular mechanisms underpinning the biological variations of these two agriculturally important aphid species.

Brain transcriptome sequencing and assembly of three songbird model systems for the study of social behavior

Emberizid sparrows (emberizidae) have played a prominent role in the study of avian vocal communication and social behavior. We present here brain transcriptomes for three emberizid model systems, song sparrow Melospiza melodia, white-throated sparrow Zonotrichia albicollis, and Gambel's white-crowned sparrow Zonotrichia leucophrys gambelii. Each of the assemblies covered fully or in part, over 89% of the previously annotated protein coding genes in the zebra finch Taeniopygia guttata, with 16,846, 15,805, and 16,646 unique BLAST hits in song, white-throated and white-crowned sparrows, respectively. As in previous studies, we find tissue of origin (auditory forebrain versus hypothalamus and whole brain) as an important determinant of overall expression profile. We also demonstrate the successful isolation of RNA and RNA-sequencing from post-mortem samples from building strikes and suggest that such an approach could be useful when traditional sampling opportunities are limited. These transcriptomes will be an important resource for the study of social behavior in birds and for data driven annotation of forthcoming whole genome sequences for these and other bird species.

Transcriptomic analysis of endangered Chinese salamander: identification of immune, sex and reproduction-related genes and genetic markers

BACKGROUND

The Chinese salamander (Hynobius chinensis), an endangered amphibian species of salamander endemic to China, has attracted much attention because of its value of studying paleontology evolutionary history and decreasing population size. Despite increasing interest in the Hynobius chinensis genome, genomic resources for the species are still very limited. A comprehensive transcriptome of Hynobius chinensis, which will provide a resource for genome annotation, candidate genes identification and molecular marker development should be generated to supplement it.

PRINCIPAL FINDINGS

We performed a de novo assembly of Hynobius chinensis transcriptome by Illumina sequencing. A total of 148,510 nonredundant unigenes with an average length of approximately 580 bp were obtained. In all, 60,388 (40.66%) unigenes showed homologous matches in at least one database and 33,537 (22.58%) unigenes were annotated by all four databases. In total, 41,553 unigenes were categorized into 62 sub-categories by BLAST2GO search, and 19,468 transcripts were assigned to 140 KEGG pathways. A large number of unigenes involved in immune system, local adaptation, reproduction and sex determination were identified, as well as 31,982 simple sequence repeats (SSRs) and 460,923 putative single nucleotide polymorphisms (SNPs).

CONCLUSION

This dataset represents the first transcriptome analysis of the Chinese salamander (Hynobius chinensis), an endangered species, to be also the first time of hynobiidae. The transcriptome will provide valuable resource for further research in discovery of new genes, protection of population, adaptive evolution and survey of various pathways, as well as development of molecule markers in Chinese salamander; and reference information for closely related species.

Cross-amplification and validation of SNPs conserved over 44 million years between seals and dogs

High-density SNP arrays developed for humans and their companion species provide a rapid and convenient tool for generating SNP data in closely-related non-model organisms, but have not yet been widely applied to phylogenetically divergent taxa. Consequently, we used the CanineHD BeadChip to genotype 24 Antarctic fur seal (Arctocephalus gazella) individuals. Despite seals and dogs having diverged around 44 million years ago, 33,324 out of 173,662 loci (19.2%) could be genotyped, of which 173 were polymorphic and clearly interpretable. Two SNPs were validated using KASP genotyping assays, with the resulting genotypes being 100% concordant with those obtained from the high-density array. Two loci were also confirmed through in silico visualisation after mapping them to the fur seal transcriptome. Polymorphic SNPs were distributed broadly throughout the dog genome and did not differ significantly in proximity to genes from either monomorphic SNPs or those that failed to cross-amplify in seals. However, the nearest genes to polymorphic SNPs were significantly enriched for functional annotations relating to energy metabolism, suggesting a possible bias towards conserved regions of the genome.

Molecular characterisation and expression analysis of interferon gamma in response to natural Chlamydia infection in the koala, Phascolarctos cinereus

Interferon gamma (IFNγ) is a key Th1 cytokine, with a principal role in the immune response against intracellular organisms such as Chlamydia. Along with being responsible for significant morbidity in human populations, Chlamydia is also responsible for wide spread infection and disease in many animal hosts, with reports that many Australian koala subpopulations are endemically infected. An understanding of the role played by IFNγ in koala chlamydial diseases is important for the establishment of better prophylactic and therapeutic approaches against chlamydial infection in this host. A limited number of IFNγ sequences have been published from marsupials and no immune reagents to measure expression have been developed. Through preliminary analysis of the koala transcriptome, we have identified the full coding sequence of the koala IFNγ gene. Transcripts were identified in spleen and lymph node tissue samples. Phylogenetic analysis demonstrated that koala IFNγ is closely related to other marsupial IFNγ sequences and more distantly related to eutherian mammals. To begin to characterise the role of this important cytokine in the koala's response to chlamydial infection, we developed a quantitative real time PCR assay and applied it to a small cohort of koalas with and without active chlamydial disease, revealing significant differences in expression patterns between the groups. Description of the IFNγ sequence from the koala will not only assist in understanding this species' response to its most important pathogen but will also provide further insight into the evolution of the marsupial immune system.