De novo assembly of the Indo-Pacific humpback dolphin leucocyte transcriptome to identify putative genes involved in the aquatic adaptation and immune response

Chromosome-level genome provides insights into environmental adaptability and innate immunity in the common dolphin (delphinus delphis)

The common dolphin (Delphinus delphis) is widely distributed worldwide and well adapted to various habitats. Animal genomes store clues about their pasts, and can reveal the genes underlying their evolutionary success. Here, we report the first high-quality chromosome-level genome of D. delphis. The assembled genome size was 2.56 Gb with a contig N50 of 63.85 Mb. Phylogenetically, D. delphis was close to Tursiops truncatus and T. aduncus. The genome of D. delphis exhibited 428 expanded and 1,885 contracted gene families, and 120 genes were identified as positively selected. The expansion of the HSP70 gene family suggested that D. delphis has a powerful system for buffering stress, which might be associated with its broad adaptability, longevity, and detoxification capacity. The expanded IFN-α and IFN-ω gene families, as well as the positively selected genes encoding tripartite motif-containing protein 25, peptidyl-prolyl cis-trans isomerase NIMA-interacting 1, and p38 MAP kinase, were all involved in pathways for antiviral, anti-inflammatory, and antineoplastic mechanisms. The genome data also revealed dramatic fluctuations in the effective population size during the Pleistocene. Overall, the high-quality genome assembly and annotation represent significant molecular resources for ecological and evolutionary studies of Delphinus and help support their sustainable treatment and conservation.

Transcriptome profiling of blood from common bottlenose dolphins (Tursiops truncatus) in the northern Gulf of Mexico to enhance health assessment capabilities

Following the 2010 Deepwater Horizon disaster and subsequent unusual mortality event, adverse health impacts have been reported in bottlenose dolphins in Barataria Bay, LA including impaired stress response and reproductive, pulmonary, cardiac, and immune function. These conditions were primarily diagnosed through hands-on veterinary examinations and analysis of standard diagnostic panels. In human and veterinary medicine, gene expression profiling has been used to identify molecular mechanisms underlying toxic responses and disease states. Identification of molecular markers of exposure or disease may enable earlier detection of health effects or allow for health evaluation when the use of specialized methodologies is not feasible. To date this powerful tool has not been applied to augment the veterinary data collected concurrently during dolphin health assessments. This study examined transcriptomic profiles of blood from 76 dolphins sampled in health assessments during 2013–2018 in the waters near Barataria Bay, LA and Sarasota Bay, FL. Gene expression was analyzed in conjunction with the substantial suite of health data collected using principal component analysis, differential expression testing, over-representation analysis, and weighted gene co-expression network analysis. Broadly, transcript profiles of Barataria Bay dolphins indicated a shift in immune response, cytoskeletal alterations, and mitochondrial dysfunction, most pronounced in dolphins likely exposed to Deepwater Horizon oiling. While gene expression profiles in Barataria Bay dolphins were altered compared to Sarasota Bay for all years, profiles from 2013 exhibited the greatest alteration in gene expression. Differentially expressed transcripts included genes involved in immunity, inflammation, reproductive failure, and lung or cardiac dysfunction, all of which have been documented in dolphins from Barataria Bay following the Deepwater Horizon oil spill. The genes and pathways identified in this study may, with additional research and validation, prove useful as molecular markers of exposure or disease to assist wildlife veterinarians in evaluating the health of dolphins and other cetaceans.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

The Landscape of Genome-Wide and Gender-Specific Microsatellites in Indo-Pacific Humpback Dolphin and Potential Applications in Cetacean Resource Investigation

Microsatellites are one of the important genome characterizations that can be a valuable resource for variety identification, genetic diversity, phylogenetic analysis, as well as comparative and conservation genomics research. Here, we developed comprehensive microsatellites through genome-wide mining for the threatened cetacean Indo-Pacific humpback dolphin (Sousa chinensis). We found 87,757 microsatellites with 2–6 bp nucleotide motifs, showing that about 32.5 microsatellites per megabase comprises microsatellites sequences. Approximately 97.8% of the markers developed in this study were consistent with the published identified markers. About 75.3% microsatellites were with dinucleotide motifs, followed by tetranucleotide motifs (17.4%), sharing the same composition pattern as other cetaceans. The microsatellites were not evenly distributed in the S. chinensis genome, mainly in non-coding regions, with only about 0.5% of the markers located in coding regions. The microsatellite-containing genes were mainly functionally enriched in the methylation process, probably demonstrating the potential impacts of microsatellites on biological functions. Polymorphic microsatellites were developed between different genders of S. chinensis, which was expected to lay the foundation for genetic diversity investigation in cetaceans. The specific markers for a male Indo-Pacific humpback dolphin will provide comprehensive and representative male candidate markers for sex identification, providing a potential biomolecular tool for further analysis of population structure and social behavior of wild populations, population trend evaluation, and species conservation management.

Transcriptomic Profiling and Microsatellite Identification in Cobia (Rachycentron canadum), Using High-Throughput RNA Sequencing

Cobia (Rachycentron canadum) is a marine teleost species with great productive potential worldwide. However, the genomic information currently available for this species in public databases is limited. Such lack of information hinders gene expression assessments that might bring forward novel insights into the physiology, ecology, evolution, and genetics of this potential aquaculture species. In this study, we report the first de novo transcriptome assembly of R. canadum liver, improving the availability of novel gene sequences for this species. Illumina sequencing of liver transcripts generated 1,761,965,794 raw reads, which were filtered into 1,652,319,304 high-quality reads. De novo assembly resulted in 101,789 unigenes and 163,096 isoforms, with an average length of 950.61 and 1,617.34 nt, respectively. Moreover, we found that 126,013 of these transcripts bear potentially coding sequences, and 125,993 of these elements (77.3%) correspond to functionally annotated genes found in six different databases. We also identified 701 putative ncRNA and 35,414 putative lncRNA. Interestingly, homologues for 410 of these putative lncRNAs have already been observed in previous analyses with Danio rerio, Lates calcarifer, Seriola lalandi dorsalis, Seriola dumerili, or Echeneis naucrates. Finally, we identified 7894 microsatellites related to cobia's putative lncRNAs. Thus, the information derived from the transcriptome assembly described herein will likely assist future nutrigenomics and breeding programs involving this important fish farming species.

Whole genomes reveal multiple candidate genes and pathways involved in the immune response of dolphins to a highly infectious virus

Wildlife species are challenged by various infectious diseases that act as important demographic drivers of populations and have become a great conservation concern particularly under growing environmental changes. The new era of whole genome sequencing provides new opportunities and avenues to explore the role of genetic variants in the plasticity of immune responses, particularly in non-model systems. Cetacean morbillivirus (CeMV) has emerged as a major viral threat to cetacean populations worldwide, contributing to the death of thousands of individuals of multiple dolphin and whale species. To understand the genomic basis of immune responses to CeMV, we generated and analysed whole genomes of 53 Indo-Pacific bottlenose dolphins (Tursiops aduncus) exposed to Australia's largest known CeMV-related mortality event that killed at least 50 dolphins from three different species. The genomic data set consisted of 10,168,981 SNPs anchored onto 23 chromosome-length scaffolds and 77 short scaffolds. Whole genome analysis indicated that levels of inbreeding in the dolphin population did not influence the outcome of an individual. Allele frequency estimates between survivors and nonsurvivors of the outbreak revealed 15,769 candidate SNPs, of which 689 were annotated to 295 protein coding genes. These included 50 genes with functions related to innate and adaptive immune responses, and cytokine signalling pathways and genes thought to be involved in immune responses to other morbilliviruses. Our study characterised genomic regions and pathways that may contribute to CeMV immune responses in dolphins. This represents a stride towards clarifying the complex interactions of the cetacean immune system and emphasises the value of whole genome data sets in understanding genetic elements that are essential for species conservation, including disease susceptibility and adaptation.

Integrated Full-Length Transcriptome and RNA-Seq to Identify Immune System Genes from the Skin of Sperm Whale (Physeter macrocephalus)

Cetaceans are a group of secondary aquatic mammals whose ancestors returned to the ocean from land, and during evolution, their immune systems adapted to the aquatic environment. Their skin, as the primary barrier to environmental pathogens, supposedly evolved to adapt to a new living environment. However, the immune system in the skin of cetaceans and the associated molecular mechanisms are still largely unknown. To better understand the immune system, we extracted RNA from the sperm whale's (Physeter macrocephalus) skin and performed PacBio full-length sequencing and RNA-seq sequencing. We obtained a total of 96,350 full-length transcripts with an average length of 1705 bp and detected 5150 genes that were associated with 21 immune-related pathways by gene annotation enrichment analysis. Moreover, we found 89 encoding genes corresponding to 33 proteins were annotated in the NOD-like receptor (NLR)-signaling pathway, including NOD1, NOD2, RIP2, and NF-kB genes, which were discussed in detail and predicted to play essential roles in the immune system of the sperm whale. Furthermore, NOD1 was highly conservative during evolution by the sequence comparison and phylogenetic tree. These results provide new information about the immune system in the skin of cetaceans, as well as the evolution of immune-related genes.

Phylogenomic Resolution of the Cetacean Tree of Life Using Target Sequence Capture

The evolution of cetaceans, from their early transition to an aquatic lifestyle to their subsequent diversification, has been the subject of numerous studies. However, although the higher-level relationships among cetacean families have been largely settled, several aspects of the systematics within these groups remain unresolved. Problematic clades include the oceanic dolphins (37 spp.), which have experienced a recent rapid radiation, and the beaked whales (22 spp.), which have not been investigated in detail using nuclear loci. The combined application of high-throughput sequencing with techniques that target specific genomic sequences provide a powerful means of rapidly generating large volumes of orthologous sequence data for use in phylogenomic studies. To elucidate the phylogenetic relationships within the Cetacea, we combined sequence capture with Illumina sequencing to generate data for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}3200 protein-coding genes for 68 cetacean species and their close relatives including the pygmy hippopotamus. By combining data from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$>$\end{document}38,000 exons with existing sequences from 11 cetaceans and seven outgroup taxa, we produced the first comprehensive comparative genomic data set for cetaceans, spanning 6,527,596 aligned base pairs (bp) and 89 taxa. Phylogenetic trees reconstructed with maximum likelihood and Bayesian inference of concatenated loci, as well as with coalescence analyses of individual gene trees, produced mostly concordant and well-supported trees. Our results completely resolve the relationships among beaked whales as well as the contentious relationships among oceanic dolphins, especially the problematic subfamily Delphinidae. We carried out Bayesian estimation of species divergence times using MCMCTree and compared our complete data set to a subset of clocklike genes. Analyses using the complete data set consistently showed less variance in divergence times than the reduced data set. In addition, integration of new fossils (e.g., Mystacodon selenensis) indicates that the diversification of Crown Cetacea began before the Late Eocene and the divergence of Crown Delphinidae as early as the Middle Miocene. [Cetaceans; phylogenomics; Delphinidae; Ziphiidae; dolphins; whales.]

Distinct evolution of toll-like receptor signaling pathway genes in cetaceans

BACKGROUND

The relatively rapid spread and diversity of marine pathogens posed an initial and ongoing challenge for cetaceans (whales, dolphins, and porpoises), descendants of terrestrial mammals that transitioned from land to sea approximately 56 million years ago. Toll-like receptors (TLRs) play important roles in regulating immunity against pathogen infections by detecting specific molecular patterns and activating a wide range of downstream signaling pathways. The ever-increasing catalogue of mammalian genomes offers unprecedented opportunities to reveal genetic changes associated with evolutionary and ecological processes.

OBJECTIVE

This study aimed to explore the molecular evolution of TLR signaling pathway genes in cetaceans.

METHODS

Genes involved in the TLR signaling pathway were retrieved by BLAST searches using human coding sequences as queries. We tested each gene for positive selection along the cetacean branches using PAML and Hyphy. Physicochemical property changes of amino acids at all positively selected residues were assessed by TreeSAAP and visualized with WebLogo. Bovine and dolphin TLR4 was assessed using human embryonic kidney cell line HEK293, which lacks TLR4 and its co-receptor MD-2.

RESULTS

We demonstrate that eight TLR signaling pathway genes are under positive selection in cetaceans. These include key genes in the response to Gram-negative bacteria: TLR4, CD14, and LY96 (MD-2). Moreover, 41 out of 65 positively selected sites were inferred to harbor substitution that dramatically changes the physicochemical properties of amino acids, with most of them situated in or adjacent to functional regions. We also found strong evidence that positive selection occurred in the lineage of the Yangtze finless porpoise, likely reflecting relatively recent adaptions to a freshwater milieu. Species-specific differences in TLR4 response were observed between cetacean and terrestrial species. Cetacean TLR4 was significantly less responsive to lipopolysaccharides from a terrestrial E. coli strain, possibly a reflection of the arms race of host-pathogen co-evolution faced by cetaceans in an aquatic environment.

CONCLUSION

This study provides further impetus for studies on the evolution and function of the cetacean immune system.

Whole Genome Sequencing of Chinese White Dolphin (Sousa chinensis) for High-Throughput Screening of Antihypertensive Peptides

Chinese white dolphin (Sousa chinensis), also known as the Indo-Pacific humpback dolphin, has been classified as “Vulnerable” on the IUCN Red List of Threatened Species. It is a special cetacean species that lives in tropical and subtropical nearshore waters, with significant differences from other cetaceans. Here, we sequenced and assembled a draft genome of the Chinese white dolphin with a total length of 2.3 Gb and annotation of 18,387 protein-coding genes. Genes from certain expanded families are potentially involved in DNA replication and repairing, suggesting that they may be related to adaptation of this marine mammal to nearshore environments. We also discovered that its historical population had undergone a remarkable bottleneck incident before the Mindel glaciation. In addition, a comparative genomic survey on antihypertensive peptides (AHTPs) among five representative mammals with various residential habitats (such as remarkable differences in exogenous ion concentrations and sea depth) revealed that these small bioactive peptides were highly conserved among these examined mammals, and they had the most abundant hits in collagen subunit proteins, especially for two putative AHTP peptides Gly-Leu-Pro (GLP) and Leu-Gly-Pro (LGP). Our genome assembly will be a valuable resource for further genetic researches on adaptive ecology and conservation biology of cetaceans, and for in-depth investigations into bioactive peptides in aquatic and terrestrial mammals for development of peptide-based drugs to treat various human cardiovascular diseases.

The genome resources for conservation of Indo-Pacific humpback dolphin, Sousa chinensis

The Indo-Pacific humpback dolphin (Sousa chinensis), is a threatened marine mammal and belongs to the First Order of the National Key Protected Wild Aquatic Animals List in China. However, limited genomic information is available for studies of its population genetics and biological conservation. Here, we have assembled a genomic sequence of this species using a whole genome shotgun (WGS) sequencing strategy after a pilot low coverage genome survey. The total assembled genome size was 2.34 Gb: with a contig N50 of 67 kb and a scaffold N50 of 9 Mb (107.6-fold sequencing coverage). The S. chinensis genome contained 24,640 predicted protein-coding genes and had approximately 37% repeated sequences. The completeness of the genome assembly was evaluated by benchmarking universal single copy orthologous genes (BUSCOs): 94.3% of a total 4,104 expected mammalian genes were identified as complete, and 2.3% were identified as fragmented. This newly produced high-quality assembly and annotation of the genome will greatly promote the future studies of the genetic diversity, conservation and evolution.

Comparative transcriptome analysis reveals the role of p53 signalling pathway during red-spotted grouper nervous necrosis virus infection in Lateolabrax japonicus brain cells

Nervous necrosis virus (NNV) is one of the fish pathogens that have caused mass mortalities of many marine and freshwater fishes in the world. To better comprehend the molecular immune mechanism of sea perch (Lateolabrax japonicus) against NNV infection, the comparative transcriptome analysis of red-spotted grouper nervous necrosis virus (RGNNV)-infected or mock-infected L. japonicus brain (LJB) cells was performed via RNA sequencing technology. Here, 1,969 up-regulated genes and 9,858 down-regulated genes, which were widely implicated in immune response pathways, were identified. Furthermore, we confirmed that p53 signalling pathway was repressed at 48 hr post-RGNNV infection, as indicated by up-regulation of Mdm2 and down-regulation of p53 and its downstream target genes, including Bax, Casp8 and CytC. Overexpression of L. japonicus p53 (Ljp53) significantly inhibited RGNNV replication and up-regulated the expression of apoptosis-related genes, whereas the down-regulation caused by pifithrin-α led to the opposite effect, suggesting Ljp53 might promote cell apoptosis to repress virus replication. Luciferase assay indicated that Ljp53 could enhance the promoter activities of zebrafish interferon (IFN)1, indicating that Ljp53 could exert its anti-RGNNV activities by enforcing the type I IFN response. This study revealed the potential antiviral role of p53 during NNV infection.

High-quality whole-genome sequence of an abundant Holarctic odontocete, the harbour porpoise (Phocoena phocoena)

The harbour porpoise (Phocoena phocoena) is a highly mobile cetacean found across the Northern hemisphere. It occurs in coastal waters and inhabits basins that vary broadly in salinity, temperature and food availability. These diverse habitats could drive subtle differentiation among populations, but examination of this would be best conducted with a robust reference genome. Here, we report the first harbour porpoise genome, assembled de novo from an individual originating in the Kattegat Sea (Sweden). The genome is one of the most complete cetacean genomes currently available, with a total size of 2.39 Gb and 50% of the total length found in just 34 scaffolds. Using 122 of the longest scaffolds, we were able to show high levels of synteny with the genome of the domestic cattle (Bos taurus). Our draft annotation comprises 22,154 predicted genes, which we further annotated through matches to the NCBI nucleotide database, GO categorization and motif prediction. Within the predicted genes, we have confirmed the presence of >20 genes or gene families that have been associated with adaptive evolution in other cetaceans. Overall, this genome assembly and draft annotation represent a crucial addition to the genomic resources currently available for the study of porpoises and Phocoenidae evolution, phylogeny and conservation.

Bioaccumulation of organic pollutants in Indo-Pacific humpback dolphin: A review on current knowledge and future prospects

Indo-Pacific humpback dolphin (Sousa chinensis) are chronically exposed to organic pollutants since they inhabit shallow coastal waters that are often impacted by anthropogenic activities. The aim of this review was to evaluate existing knowledge on the occurrence of organic pollutants in Indo-Pacific humpback dolphins, identify knowledge gaps, and offer recommendations for future research directions. We discussed the trends in the bioaccumulation of organic pollutants in Indo-Pacific humpback dolphins focusing on sources, physicochemical properties, and usage patterns. Furthermore, we examined factors that influence bioaccumulation such as gender, age, dietary intake and tissue-specific distribution. Studies on bioaccumulation in Indo-Pacific humpback dolphin remain scarce, despite high concentrations above 13,000 ng/g lw we previously detected for PFOS, ∑PBDE and chlorinated paraffins. The maximum concentration of organochlorines detected was 157,000 ng/g wt. Furthermore, variations in bioaccumulation were shown to be caused by factors such as usage patterns and physicochemical properties of the pollutant. However, restrictions in sampling inhibit investigations on exposure pathway and toxicity of organic pollutants in Indo-Pacific humpback dolphin. We proposed the use of biopsy sampling, predictive bioaccumulation and toxicity modeling, and monitoring other emerging contaminants such as microplastics and pharmaceuticals for future health risk assessment on this critically endangered marine mammal species.

New Technologies for Monitoring Marine Mammal Health

Omics technologies have been developed in recent decades and applied to different subjects, although the greatest advancements have been achieved in human biology and disease. Genome sequencing and the exploration of its coding and noncoding regions are rapidly yielding meaningful answers to diverse questions, relating genome information to protein activity to environmental changes. In the past, marine mammal genetic and transcriptional studies have been restricted due to the lack of reference genomes. But the advance of high-throughput sequencing is revolutionizing the life sciences technologies. As long-lived organisms, at the top of the food chain, marine mammals play an important role in marine ecosystems and while their protected status is in favor of conservation of the species, it also complicates the researcher's approach to traditional measurements of health. Omics data generated by high-throughput technologies will represent an important key for improving the scientific basis for understanding both marine mammal and environment health.

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 sequencing, assembly and analysis of eight different transcriptomes from the Malayan pangolin

Pangolins are scale-covered mammals, containing eight endangered species. Maintaining pangolins in captivity is a significant challenge, in part because little is known about their genetics. Here we provide the first large-scale sequencing of the critically endangered Manis javanica transcriptomes from eight different organs using Illumina HiSeq technology, yielding ~75 Giga bases and 89,754 unigenes. We found some unigenes involved in the insect hormone biosynthesis pathway and also 747 lipids metabolism-related unigenes that may be insightful to understand the lipid metabolism system in pangolins. Comparative analysis between M. javanica and other mammals revealed many pangolin-specific genes significantly over-represented in stress-related processes, cell proliferation and external stimulus, probably reflecting the traits and adaptations of the analyzed pregnant female M. javanica. Our study provides an invaluable resource for future functional works that may be highly relevant for the conservation of pangolins.

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.

Integrative analysis of transcriptomics and proteomics of skeletal muscles of the Chinese indigenous Shaziling pig compared with the Yorkshire breed

Background

The Shaziling pig (Sus scrofa) is a well-known indigenous breed in China. One of its main advantages over European breeds is its high meat quality. However, little genetic information is available for the Shaziling pig. To screen for differentially expressed genes and proteins that might be responsible for the meat quality, the longissimus dorsi muscles from Shaziling and Yorkshire pig breeds were investigated using an integrative analysis of transcriptomics and proteomics, involving high-throughput sequencing, the two-dimensional gel electrophoresis, and mass spectrometry.

Results

Sequencing produced 79,320 unigenes by de novo assembly, and 488 differentially expressed genes in the longissimus dorsi muscle of Shaziling pig compared with the Yorkshire breed were identified. Gene Ontology term enrichment of biological functions and Kyoto Encyclopedia of Genes and Genomes analysis showed that the gene products were mainly involved in metabolism, protein binding, and regulation of skeletal muscle development. At the protein level, 23 differentially expressed proteins were identified, which were potentially associated with fatty acid metabolism, the glycolytic pathway, and skeletal muscle growth. Eight differentially expressed genes were confirmed by real-time PCR. These results give an insight into the mechanisms underlying the formation of skeletal muscle in the Shaziling pig.

Conclusions

Certain differentially expressed genes and proteins are involved in fatty acid metabolism, intramuscular fat deposition, and skeletal muscle growth in the Shaziling pig. These results provide candidate genes for improving meat quality and will promote further transcriptomic research in Shaziling pigs.

Electronic supplementary material

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

De novo characterization of the liver transcriptome of javelin goby Synechogobius hasta and analysis of its transcriptomic profile following waterborne copper exposure

Previous studies have investigated the physiological responses to chronic copper (Cu) exposure in the liver of Synechogobius hasta; however, little information is available on the underlying molecular mechanisms. In an effort to better understand the mechanisms of Cu toxicity and to illuminate global gene expression patterns modulated by Cu exposure, we obtained the liver transcriptome information of S. hasta by RNA sequencing (RNA-seq) technology and also investigated the differential expression of genes following waterborne Cu exposure. Using the Illumina sequencing platform, as many as 60,217 unigenes were generated, with 815 bp of average length and 1298 bp of unigene N50 after filtering and assembly. For functional annotation analysis, 34,860, 31,526, 31,576, 25,808, 11,542, and 21,721 unigenes were annotated to the NR, NT, Swiss-Prot, KEGG, COG, and GO databases, respectively, and total annotation unigenes were 37,764. After 30 days of exposure to 55 μg Cu/l, a total of 292 and 1076 genes were significantly up- and down-regulated, respectively. By KEGG analysis, 660 had a specific pathway annotation. Subsequent bioinformatics analysis revealed that the differentially expressed genes were mainly related to lipid metabolism, immune system, apoptosis, and signal transduction, suggesting that these signaling pathways may be regulated by Cu exposure. The present study provides comprehensive sequence information for subsequent gene expression studies regarding S. hasta, and the transcriptome profiling after Cu exposure is also expected to improve our understanding of the molecular toxicology of Cu.

A phylogenomic analysis of the role and timing of molecular adaptation in the aquatic transition of cetartiodactyl mammals

Recent studies have reported multiple cases of molecular adaptation in cetaceans related to their aquatic abilities. However, none of these has included the hippopotamus, precluding an understanding of whether molecular adaptations in cetaceans occurred before or after they split from their semi-aquatic sister taxa. Here, we obtained new transcriptomes from the hippopotamus and humpback whale, and analysed these together with available data from eight other cetaceans. We identified more than 11 000 orthologous genes and compiled a genome-wide dataset of 6845 coding DNA sequences among 23 mammals, to our knowledge the largest phylogenomic dataset to date for cetaceans. We found positive selection in nine genes on the branch leading to the common ancestor of hippopotamus and whales, and 461 genes in cetaceans compared to 64 in hippopotamus. Functional annotation revealed adaptations in diverse processes, including lipid metabolism, hypoxia, muscle and brain function. By combining these findings with data on protein-protein interactions, we found evidence suggesting clustering among gene products relating to nervous and muscular systems in cetaceans. We found little support for shared ancestral adaptations in the two taxa; most molecular adaptations in extant cetaceans occurred after their split with hippopotamids.

De Novo assembly of the Japanese flounder (Paralichthys olivaceus) spleen transcriptome to identify putative genes involved in immunity

Background

Japanese flounder (Paralichthys olivaceus) is an economically important marine fish in Asia and has suffered from disease outbreaks caused by various pathogens, which requires more information for immune relevant genes on genome background. However, genomic and transcriptomic data for Japanese flounder remain scarce, which limits studies on the immune system of this species. In this study, we characterized the Japanese flounder spleen transcriptome using an Illumina paired-end sequencing platform to identify putative genes involved in immunity.

Methodology/Principal Findings

A cDNA library from the spleen of P. olivaceus was constructed and randomly sequenced using an Illumina technique. The removal of low quality reads generated 12,196,968 trimmed reads, which assembled into 96,627 unigenes. A total of 21,391 unigenes (22.14%) were annotated in the NCBI Nr database, and only 1.1% of the BLASTx top-hits matched P. olivaceus protein sequences. Approximately 12,503 (58.45%) unigenes were categorized into three Gene Ontology groups, 19,547 (91.38%) were classified into 26 Cluster of Orthologous Groups, and 10,649 (49.78%) were assigned to six Kyoto Encyclopedia of Genes and Genomes pathways. Furthermore, 40,928 putative simple sequence repeats and 47, 362 putative single nucleotide polymorphisms were identified. Importantly, we identified 1,563 putative immune-associated unigenes that mapped to 15 immune signaling pathways.

Conclusions/Significance

The P. olivaceus transciptome data provides a rich source to discover and identify new genes, and the immune-relevant sequences identified here will facilitate our understanding of the mechanisms involved in the immune response. Furthermore, the plentiful potential SSRs and SNPs found in this study are important resources with respect to future development of a linkage map or marker assisted breeding programs for the flounder.

Whole transcriptome analysis using next-generation sequencing of sterile-cultured Eisenia andrei for immune system research

Recently, earthworms have become a useful model for research into the immune system, and it is expected that results obtained using this model will shed light on the sophisticated vertebrate immune system and the evolution of the immune response, and additionally help identify new biomolecules with therapeutic applications. However, for earthworms to be used as a genetic model of the invertebrate immune system, basic molecular and genetic resources, such as an expressed sequence tag (EST) database, must be developed for this organism. Next-generation sequencing technologies have generated EST libraries by RNA-seq in many model species. In this study, we used Illumina RNA-sequence technology to perform a comprehensive transcriptome analysis using an RNA sample pooled from sterile-cultured Eisenia andrei. All clean reads were assembled de novo into 41,423 unigenes using the Trinity program. Using this transcriptome data, we performed BLAST analysis against the GenBank non-redundant (NR) database and obtained a total of 12,285 significant BLAST hits. Furthermore, gene ontology (GO) analysis assigned 78 unigenes to 24 immune class GO terms. In addition, we detected a unigene with high similarity to beta-1,3-glucuronyltransferase 1 (GlcAT-P), which mediates a glucuronyl transfer reaction during the biosynthesis of the carbohydrate epitope HNK-1 (human natural killer-1, also known as CD57), a marker of NK cells. The identified transcripts will be used to facilitate future research into the immune system using E. andrei.

De novo assembly and characterization of narrow-ridged finless porpoise renal transcriptome and identification of candidate genes involved in osmoregulation

During the evolutionary transition from land to water, cetaceans have undergone numerous critical challenges, with osmoregulation being the major one. Two subspecies of the narrow-ridged finless porpoise (Neophocaena asiaeorientalis), the freshwater Yangtze finless porpoise (N. a. asiaeorientalis, NAA) and the marine East Asian finless porpoise (N. a. sunameri, NAS), provide excellent subjects to understand the genetic basis of osmoregulatory divergence between freshwater and marine mammals. The kidney plays an important and well-established role in osmoregulation in marine mammals and thus, herein, we utilized RNA-seq to characterize the renal transcriptome and preliminarily analyze the divergence between the NAA and the NAS. Approximately 48.98 million clean reads from NAS and 49.40 million clean reads from NAA were obtained by RNA-Seq. And 73,449 (NAS) and 68,073 (NAA) unigenes were assembled. Among these annotations, 22,231 (NAS) and 21,849 (NAA) unigenes were annotated against the NCBI nr protein database. The ion channel complex GO term and four pathways were detected as relevant to osmoregulation by GO and KEGG pathway classification of these annotated unigenes. Although the endangered status of the study species prevented analysis of biological replicates, we identified nine differentially expressed genes (DEGs) that may be vital in the osmoregulation of the narrow-ridged finless porpoise and worthwhile for future studies. Of these DEGs, the differential expression and distribution of the aquaporin-2 (AQP2) in the collecting duct were verified using immunohistochemical experiments. Together, this work is the first report of renal transcriptome sequencing in cetaceans, and it will provide a valuable resource for future molecular genetics studies on cetacean osmoregulation.