Conservation and divergence in the frog immunome: pyrosequencing and de novo assembly of immune tissue transcriptomes

Molecular diversity and functional implication of amphibian interferon complex: Remarking immune adaptation in vertebrate evolution

Cross-species comparison of vertebrate genomes has unraveled previously unknown complexities of interferon (IFN) systems in amphibian species. Recent genomic curation revealed that amphibian species have evolved expanded repertoires of four types of intron-containing IFN genes akin to those seen in jawed fish, intronless type I IFNs and intron-containing type III IFNs akin to those seen in amniotes, as well as uniquely intronless type III IFNs. This appears to be the case with at least ten analyzed amphibian species; with distinct species encoding diverse repertoires of these respective IFN gene subsets. Amphibians represent a key stage in vertebrate evolution, and in this context offer a unique perspective into the divergent and converged pathways leading to the emergence of distinct IFN families and groups. Recent studies have begun to unravel the roles of amphibian IFNs during these animals' immune responses in general and during their antiviral responses, in particular. However, the pleiotropic potentials of these highly expanded amphibian IFN repertoires warrant further studies. Based on recent reports and our omics analyses using Xenopus models, we posit that amphibian IFN complex may have evolved novel functions, as indicated by their extensive molecular diversity. Here, we provide an overview and an update of the present understanding of the amphibian IFN complex in the context of the evolution of vertebrate immune systems. A greater understanding of the amphibian IFN complex will grant new perspectives on the evolution of vertebrate immunity and may yield new measures by which to counteract the global amphibian declines.

Immune priming prior to pathogen exposure sheds light on the relationship between host, microbiome and pathogen in disease

Dynamic interactions between host, pathogen and host-associated microbiome dictate infection outcomes. Pathogens including Batrachochytrium dendrobatidis (Bd) threaten global biodiversity, but conservation efforts are hindered by limited understanding of amphibian host, Bd and microbiome interactions. We conducted a vaccination and infection experiment using Eastern hellbender salamanders (Cryptobranchus alleganiensis alleganiensis) challenged with Bd to observe infection, skin microbial communities and gene expression of host skin, pathogen and microbiome throughout the experiment. Most animals survived high Bd loads regardless of their vaccination status and vaccination did not affect pathogen load, but host gene expression differed based on vaccination. Oral vaccination (exposure to killed Bd) stimulated immune gene upregulation while topically and sham-vaccinated animals did not significantly upregulate immune genes. In early infection, topically vaccinated animals upregulated immune genes but orally and sham-vaccinated animals downregulated immune genes. Bd increased pathogenicity-associated gene expression in late infection when Bd loads were highest. The microbiome was altered by Bd, but there was no correlation between anti-Bd microbe abundance or richness and pathogen burden. Our observations suggest that hellbenders initially generate a vigorous immune response to Bd, which is ineffective at controlling disease and is subsequently modulated. Interactions with antifungal skin microbiota did not influence disease progression.

Transcriptome characterization and SSR discovery in the giant spiny frog Quasipaa spinosa

The giant spiny frog Quasipaa spinosa (Amphibia: Ranidae) is a large unique frog species found mainly in southern China with a low amount of fat and high protein, and it has become one of the most important aquaculture animal species in China. To better understand its genetic background and screen potential molecular markers for artificial breeding and species conservation, we constructed an expression profile of Q. spinosa with high-throughput RNA sequencing and acquired potential SSR markers. Approximately 81.7 Gb of data and 93,887 unigenes were generated. The transcriptome contains 2085 (80.7 %) complete BUSCOs, suggesting that our assembly methods were effective and accurate.These unigenes were functionally classified using 7 functional databases, yielding 17,482 Pfam-, 12,752 Sting-, 17,526 KEGG-, 24,341 Swiss-Prot-, 28,604 Nr-, 16,287 GO- and 12,752 COG-annotated unigenes. Among several amphibian species, Q. spinosa unigenes had the highest number of hits to Xenopus tropicalis (35.25 %), followed by Xenopus laevis (12.68 %). 1417 unigenes were assigned to the immune system. In addition, a total of 33,019 candidate SSR markers were identified from the constructed library. Further tests with 20 loci and 118 large-scale breeding specimens gathered from four culture farms in China showed that 15 (75 %) loci were polymorphic, with the number of alleles per locus varying from 3 to 9 (mean of 4.3). The PIC values for the SSR markers ranged from 0.19 to 0.82, with an average value of 0.43, indicating moderate polymorphism in Q. spinosa. The transcriptomic profile and SSR repertoire obtained in the present study will facilitate population genetic studies and the selective breeding of amphibian species.

Xenopus Interferon Complex: Inscribing the Amphibiotic Adaption and Species-Specific Pathogenic Pressure in Vertebrate Evolution?

Several recent studies have revealed previously unknown complexity of the amphibian interferon (IFN) system. Being unique in vertebrate animals, amphibians not only conserve and multiply the fish-like intron-containing IFN genes, but also rapidly evolve amniote-like intronless IFN genes in each tested species. We postulate that the amphibian IFN system confers an essential model to study vertebrate immune evolution in molecular and functional diversity to cope with unprecedented pathophysiological requirement during terrestrial adaption. Studies so far have ascribed a potential role of these IFNs in immune regulation against intracellular pathogens, particularly viruses; however, many knowledge gaps remain elusive. Based on recent reports about IFN's multifunctional properties in regulation of animal physiological and defense responses, we interpret that amphibian IFNs may evolve novel function pertinent to their superior molecular diversity. Such new function revealed by the emerging studies about antifungal and developmental regulation of amphibian IFNs will certainly promote our understanding of immune evolution in vertebrates to address current pathogenic threats causing amphibian decline.

Frog Skin Innate Immune Defences: Sensing and Surviving Pathogens

Amphibian skin is a mucosal surface in direct and continuous contact with a microbially diverse and laden aquatic and/or terrestrial environment. As such, frog skin is an important innate immune organ and first line of defence against pathogens in the environment. Critical to the innate immune functions of frog skin are the maintenance of physical, chemical, cellular, and microbiological barriers and the complex network of interactions that occur across all the barriers. Despite the global decline in amphibian populations, largely as a result of emerging infectious diseases, we understand little regarding the cellular and molecular mechanisms that underlie the innate immune function of amphibian skin and defence against pathogens. In this review, we discuss the structure, cell composition and cellular junctions that contribute to the skin physical barrier, the antimicrobial peptide arsenal that, in part, comprises the chemical barrier, the pattern recognition receptors involved in recognizing pathogens and initiating innate immune responses in the skin, and the contribution of commensal microbes on the skin to pathogen defence. We briefly discuss the influence of environmental abiotic factors (natural and anthropogenic) and pathogens on the immunocompetency of frog skin defences. Although some aspects of frog innate immunity, such as antimicrobial peptides are well-studied; other components and how they contribute to the skin innate immune barrier, are lacking. Elucidating the complex network of interactions occurring at the interface of the frog's external and internal environments will yield insight into the crucial role amphibian skin plays in host defence and the environmental factors leading to compromised barrier integrity, disease, and host mortality.

Identification of critical sex-biased genes in Andrias davidianus by de novo transcriptome

The Chinese giant salamander Andrias davidianus is a protected amphibian with high nutritional and economic value. Understanding its sex determination mechanism is important for improving culture techniques and sex control in breeding. However, little information on the characterization of critical genes involved in sex is available. Herein, sequencing of ovary and test produced 40,783,222 and 46,128,902 raw reads, respectively, which were jointly assembled into 80,497 unigenes. Of these, 36,609 unigenes were annotated, of which 8907 were female-biased and 10,385 were male-biased. Several sex-related pathways were observed, including the Wnt signaling pathway. After elevated temperature and estrogen exposure, neomale and neofemale specimens were identified by a female-specific marker for the first time. RT-qPCR analysis showed the expression profile of ten selected sex-biased genes to be exhibited consistently in male and neomale and in female and neofemale, with the exception of the Amh and TfIIIa genes. Results suggested that these genes may play important roles in A. davidianus sex determination and gonad development. This provides a basis for further investigation of the molecular mechanisms of sex determination in amphibians.

Transcriptome analyses of immune tissues from three Japanese frogs (genus Rana ) reveals their utility in characterizing major histocompatibility complex class II

Background

In Japan and East Asia, endemic frogs appear to be tolerant or not susceptible to chytridiomycosis, a deadly amphibian disease caused by the chytrid fungus Batrachochytridium dendrobatidis (Bd). Japanese frogs may have evolved mechanisms of immune resistance to pathogens such as Bd. This study characterizes immune genes expressed in various tissues of healthy Japanese Rana frogs.

Results

We generated transcriptome data sets of skin, spleen and blood from three adult Japanese Ranidae frogs (Japanese brown frog Rana japonica, the montane brown frog Rana ornativentris, and Tago’s brown frog Rana tagoi tagoi) as well as whole body of R. japonica and R. ornativentris tadpoles. From this, we identified tissue- and stage-specific differentially expressed genes; in particular, the spleen was most enriched for immune-related genes. A specific immune gene, major histocompatibility complex class IIB (MHC-IIB), was further characterized due to its role in pathogen recognition. We identified a total of 33 MHC-IIB variants from the three focal species (n = 7 individuals each), which displayed evolutionary signatures related to increased MHC variation, including balancing selection. Our supertyping analyses of MHC-IIB variants from Japanese frogs and previously studied frog species identified potential physiochemical properties of MHC-II that may be important for recognizing and binding chytrid-related antigens.

Conclusions

This is one of the first studies to generate transcriptomic resources for Japanese frogs, and contributes to further understanding the immunogenetic factors associated with resistance to infectious diseases in amphibians such as chytridiomycosis. Notably, MHC-IIB supertyping analyses identified unique functional properties of specific MHC-IIB alleles that may partially contribute to Bd resistance, and such properties provide a springboard for future experimental validation.

Electronic supplementary material

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

Mitigating amphibian chytridiomycoses in nature

Amphibians across the planet face the threat of population decline and extirpation caused by the disease chytridiomycosis. Despite consensus that the fungal pathogens responsible for the disease are conservation issues, strategies to mitigate their impacts in the natural world are, at best, nascent. Reducing risk associated with the movement of amphibians, non-amphibian vectors and other sources of infection remains the first line of defence and a primary objective when mitigating the threat of disease in wildlife. Amphibian-associated chytridiomycete fungi and chytridiomycosis are already widespread, though, and we therefore focus on discussing options for mitigating the threats once disease emergence has occurred in wild amphibian populations. All strategies have shortcomings that need to be overcome before implementation, including stronger efforts towards understanding and addressing ethical and legal considerations. Even if these issues can be dealt with, all currently available approaches, or those under discussion, are unlikely to yield the desired conservation outcome of disease mitigation. The decision process for establishing mitigation strategies requires integrated thinking that assesses disease mitigation options critically and embeds them within more comprehensive strategies for the conservation of amphibian populations, communities and ecosystems.This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

Transcriptome analysis illuminates the nature of the intracellular interaction in a vertebrate-algal symbiosis

During embryonic development, cells of the green alga Oophila amblystomatis enter cells of the salamander Ambystoma maculatum forming an endosymbiosis. Here, using de novo dual-RNA seq, we compared the host salamander cells that harbored intracellular algae to those without algae and the algae inside the animal cells to those in the egg capsule. This two-by-two-way analysis revealed that intracellular algae exhibit hallmarks of cellular stress and undergo a striking metabolic shift from oxidative metabolism to fermentation. Culturing experiments with the alga showed that host glutamine may be utilized by the algal endosymbiont as a primary nitrogen source. Transcriptional changes in salamander cells suggest an innate immune response to the alga, with potential attenuation of NF-κB, and metabolic alterations indicative of modulation of insulin sensitivity. In stark contrast to its algal endosymbiont, the salamander cells did not exhibit major stress responses, suggesting that the host cell experience is neutral or beneficial.

Selective Landscapes in newt Immune Genes Inferred from Patterns of Nucleotide Variation

Host–pathogen interactions may result in either directional selection or in pressure for the maintenance of polymorphism at the molecular level. Hence signatures of both positive and balancing selection are expected in immune genes. Because both overall selective pressure and specific targets may differ between species, large-scale population genomic studies are useful in detecting functionally important immune genes and comparing selective landscapes between taxa. Such studies are of particular interest in amphibians, a group threatened worldwide by emerging infectious diseases. Here, we present an analysis of polymorphism and divergence of 634 immune genes in two lineages of Lissotriton newts: L. montandoni and L. vulgaris graecus. Variation in newt immune genes has been shaped predominantly by widespread purifying selection and strong evolutionary constraint, implying long-term importance of these genes for functioning of the immune system. The two evolutionary lineages differ in the overall strength of purifying selection which can partially be explained by demographic history but may also signal differences in long-term pathogen pressure. The prevalent constraint notwithstanding, 23 putative targets of positive selection and 11 putative targets of balancing selection were identified. The latter were detected by composite tests involving the demographic model and further validated in independent population samples. Putative targets of balancing selection encode proteins which may interact closely with pathogens but include also regulators of immune response. The identified candidates will be useful for testing whether genes affected by balancing selection are more prone to interspecific introgression than other genes in the genome.

Research proceedings on amphibian model organisms

Model organisms have long been important in biology and medicine due to their specific characteristics. Amphibians, especially Xenopus, play key roles in answering fundamental questions on developmental biology, regeneration, genetics, and toxicology due to their large and abundant eggs, as well as their versatile embryos, which can be readily manipulated and developed in vivo. Furthermore, amphibians have also proven to be of considerable benefit in human disease research due to their conserved cellular developmental and genomic organization. This review gives a brief introduction on the progress and limitations of these animal models in biology and human disease research, and discusses the potential and challenge of Microhyla fissipes as a new model organism.

Sorting duplicated loci disentangles complexities of polyploid genomes masked by genotyping by sequencing

Many plants and animals of polyploid origin are currently enjoying a genomics explosion enabled by modern sequencing and genotyping technologies. However, routine filtering of duplicated loci in most studies using genotyping by sequencing introduces an unacceptable, but often overlooked, bias when detecting selection. Retained duplicates from ancient whole-genome duplications (WGDs) may be found throughout genomes, whereas retained duplicates from recent WGDs are concentrated at distal ends of some chromosome arms. Additionally, segmental duplicates can be found at distal ends or nearly anywhere in a genome. Evidence shows that these duplications facilitate adaptation through one of two pathways: neo-functionalization or increased gene expression. Filtering duplicates removes distal ends of some chromosomes, and distal ends are especially known to harbour adaptively important genes. Thus, filtering of duplicated loci impoverishes the interpretation of genomic data as signals from contiguous duplicated genes are ignored. We review existing strategies to genotype and map duplicated loci; we focus in detail on an overlooked strategy of using gynogenetic haploids (1N) as a part of new genotyping by sequencing studies. We provide guidelines on how to use this haploid strategy for studies on polyploid-origin vertebrates including how it can be used to screen duplicated loci in natural populations. We conclude by discussing areas of research that will benefit from better inclusion of polyploid loci; we particularly stress the sometimes overlooked fact that basing genomic studies on dense maps provides value added in the form of locating and annotating outlier loci or colocating outliers into islands of divergence.

Comparative transcriptome analyses of seven anurans reveal functions and adaptations of amphibian skin

Animal skin, which is the tissue that directly contacts the external surroundings, has evolved diverse functions to adapt to various environments. Amphibians represent the transitional taxon from aquatic to terrestrial life. Exploring the molecular basis of their skin function and adaptation is important to understand the survival and evolutionary mechanisms of vertebrates. However, comprehensive studies on the molecular mechanisms of skin functions in amphibians are scarce. In this study, we sequenced the skin transcriptomes of seven anurans belonging to three families and compared the similarities and differences in expressed genes and proteins. Unigenes and pathways related to basic biological processes and special functions, such as defense, immunity, and respiration, were enriched in functional annotations. A total of 108 antimicrobial peptides were identified. The highly expressed genes were similar in species of the same family but were different among families. Additionally, the positively selected orthologous groups were involved in biosynthesis, metabolism, immunity, and defense processes. This study is the first to generate extensive transcriptome data for the skin of seven anurans and provides unigenes and pathway candidates for further studies on amphibian skin function and adaptation.

Using "Omics" and Integrated Multi-Omics Approaches to Guide Probiotic Selection to Mitigate Chytridiomycosis and Other Emerging Infectious Diseases

Emerging infectious diseases in wildlife are responsible for massive population declines. In amphibians, chytridiomycosis caused by Batrachochytrium dendrobatidis, Bd, has severely affected many amphibian populations and species around the world. One promising management strategy is probiotic bioaugmentation of antifungal bacteria on amphibian skin. In vivo experimental trials using bioaugmentation strategies have had mixed results, and therefore a more informed strategy is needed to select successful probiotic candidates. Metagenomic, transcriptomic, and metabolomic methods, colloquially called "omics," are approaches that can better inform probiotic selection and optimize selection protocols. The integration of multiple omic data using bioinformatic and statistical tools and in silico models that link bacterial community structure with bacterial defensive function can allow the identification of species involved in pathogen inhibition. We recommend using 16S rRNA gene amplicon sequencing and methods such as indicator species analysis, the Kolmogorov-Smirnov Measure, and co-occurrence networks to identify bacteria that are associated with pathogen resistance in field surveys and experimental trials. In addition to 16S amplicon sequencing, we recommend approaches that give insight into symbiont function such as shotgun metagenomics, metatranscriptomics, or metabolomics to maximize the probability of finding effective probiotic candidates, which can then be isolated in culture and tested in persistence and clinical trials. An effective mitigation strategy to ameliorate chytridiomycosis and other emerging infectious diseases is necessary; the advancement of omic methods and the integration of multiple omic data provide a promising avenue toward conservation of imperiled species.

4Pipe4--A 454 data analysis pipeline for SNP detection in datasets with no reference sequence or strain information

Background

Next-generation sequencing datasets are becoming more frequent, and their use in population studies is becoming widespread. For non-model species, without a reference genome, it is possible from a panel of individuals to identify a set of SNPs that can be used for further population genotyping. However the lack of a reference genome to which the sequenced data could be compared makes the finding of SNPs more troublesome. Additionally when the data sources (strains) are not identified (e.g. in datasets of pooled individuals), the problem of finding reliable variation in these datasets can become much more difficult due to the lack of specialized software for this specific task.

Results

Here we describe 4Pipe4, a 454 data analysis pipeline particularly focused on SNP detection when no reference or strain information is available. It uses a command line interface to automatically call other programs, parse their outputs and summarize the results. The variation detection routine is built-in in the program itself. Despite being optimized for SNP mining in 454 EST data, it is flexible enough to automate the analysis of genomic data or even data from other NGS technologies. 4Pipe4 will output several HTML formatted reports with metrics on many of the most common assembly values, as well as on all the variation found. There is also a module available for finding putative SSRs in the analysed datasets.

Conclusions

This program can be especially useful for researchers that have 454 datasets of a panel of pooled individuals and want to discover and characterize SNPs for subsequent individual genotyping with customized genotyping arrays. In comparison with other SNP detection approaches, 4Pipe4 showed the best validation ratio, retrieving a smaller number of SNPs but with a considerably lower false positive rate than other methods.

4Pipe4’s source code is available at https://github.com/StuntsPT/4Pipe4.

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems.

De novo Transcriptome Assemblies of Rana (Lithobates) catesbeiana and Xenopus laevis Tadpole Livers for Comparative Genomics without Reference Genomes

In this work we studied the liver transcriptomes of two frog species, the American bullfrog (Rana (Lithobates) catesbeiana) and the African clawed frog (Xenopus laevis). We used high throughput RNA sequencing (RNA-seq) data to assemble and annotate these transcriptomes, and compared how their baseline expression profiles change when tadpoles of the two species are exposed to thyroid hormone. We generated more than 1.5 billion RNA-seq reads in total for the two species under two conditions as treatment/control pairs. We de novo assembled these reads using Trans-ABySS to reconstruct reference transcriptomes, obtaining over 350,000 and 130,000 putative transcripts for R. catesbeiana and X. laevis, respectively. Using available genomics resources for X. laevis, we annotated over 97% of our X. laevis transcriptome contigs, demonstrating the utility and efficacy of our methodology. Leveraging this validated analysis pipeline, we also annotated the assembled R. catesbeiana transcriptome. We used the expression profiles of the annotated genes of the two species to examine the similarities and differences between the tadpole liver transcriptomes. We also compared the gene ontology terms of expressed genes to measure how the animals react to a challenge by thyroid hormone. Our study reports three main conclusions. First, de novo assembly of RNA-seq data is a powerful method for annotating and establishing transcriptomes of non-model organisms. Second, the liver transcriptomes of the two frog species, R. catesbeiana and X. laevis, show many common features, and the distribution of their gene ontology profiles are statistically indistinguishable. Third, although they broadly respond the same way to the presence of thyroid hormone in their environment, their receptor/signal transduction pathways display marked differences.

A de novo Assembly of the Common Frog (Rana temporaria) Transcriptome and Comparison of Transcription Following Exposure to Ranavirus and Batrachochytrium dendrobatidis

Amphibians are experiencing global declines and extinctions, with infectious diseases representing a major factor. In this study we examined the transcriptional response of metamorphic hosts (common frog, Rana temporaria) to the two most important amphibian pathogens: Batrachochytrium dendrobatidis (Bd) and Ranavirus. We found strong up-regulation of a gene involved in the adaptive immune response (AP4S1) at four days post-exposure to both pathogens. We detected a significant transcriptional response to Bd, covering the immune response (innate and adaptive immunity, complement activation, and general inflammatory responses), but relatively little transcriptional response to Ranavirus. This may reflect the higher mortality rates found in wild common frogs infected with Ranavirus as opposed to Bd. These data provide a valuable genomic resource for the amphibians, contribute insight into gene expression changes after pathogen exposure, and suggest potential candidate genes for future host-pathogen research.

Linking genetic and environmental factors in amphibian disease risk

A central question in evolutionary biology is how interactions between organisms and the environment shape genetic differentiation. The pathogen Batrachochytrium dendrobatidis (Bd) has caused variable population declines in the lowland leopard frog (Lithobates yavapaiensis); thus, disease has potentially shaped, or been shaped by, host genetic diversity. Environmental factors can also influence both amphibian immunity and Bd virulence, confounding our ability to assess the genetic effects on disease dynamics. Here, we used genetics, pathogen dynamics, and environmental data to characterize L. yavapaiensis populations, estimate migration, and determine relative contributions of genetic and environmental factors in predicting Bd dynamics. We found that the two uninfected populations belonged to a single genetic deme, whereas each infected population was genetically unique. We detected an outlier locus that deviated from neutral expectations and was significantly correlated with mortality within populations. Across populations, only environmental variables predicted infection intensity, whereas environment and genetics predicted infection prevalence, and genetic diversity alone predicted mortality. At one locality with geothermally elevated water temperatures, migration estimates revealed source-sink dynamics that have likely prevented local adaptation. We conclude that integrating genetic and environmental variation among populations provides a better understanding of Bd spatial epidemiology, generating more effective conservation management strategies for mitigating amphibian declines.

TRUFA: A User-Friendly Web Server for de novo RNA-seq Analysis Using Cluster Computing

Application of next-generation sequencing (NGS) methods for transcriptome analysis (RNA-seq) has become increasingly accessible in recent years and are of great interest to many biological disciplines including, eg, evolutionary biology, ecology, biomedicine, and computational biology. Although virtually any research group can now obtain RNA-seq data, only a few have the bioinformatics knowledge and computation facilities required for transcriptome analysis. Here, we present TRUFA (TRanscriptome User-Friendly Analysis), an open informatics platform offering a web-based interface that generates the outputs commonly used in de novo RNA-seq analysis and comparative transcriptomics. TRUFA provides a comprehensive service that allows performing dynamically raw read cleaning, transcript assembly, annotation, and expression quantification. Due to the computationally intensive nature of such analyses, TRUFA is highly parallelized and benefits from accessing high-performance computing resources. The complete TRUFA pipeline was validated using four previously published transcriptomic data sets. TRUFA’s results for the example datasets showed globally similar results when comparing with the original studies, and performed particularly better when analyzing the green tea dataset. The platform permits analyzing RNA-seq data in a fast, robust, and user-friendly manner. Accounts on TRUFA are provided freely upon request at https://trufa.ifca.es.

Transcriptome analysis revealed positive selection of immune-related genes in tilapia

High-throughput sequencing of transcriptome promises a new approach for detecting evolutionary divergence among species. Up to now, the information about evolution of immune genes in cultured fish, especially in tilapias which would aid to understand the molecular basis of immune phenotypic differentiation is still lack. Thus, in the present study, we used high-throughput sequencing to obtain large amount of gene sequences in blue tilapia and characterized the diversity of orthologs among Nile tilapia, blue tilapia and zebrafish. A total of 52,424,506 raw reads, representing 31,404 unigenes were obtained from blue tilapia cDNA library of mixed tissues, including brain, pituitary, gill, heart, liver, spleen, kidney, intestine, muscle, testis and ovary. Based on Ks value, we calculated that the divergence time between Nile tilapia and blue tilapia is 2.93 million years ago. And the tilapias are both apart from zebrafish in 197 million years ago. Furthermore, the positive selected genes were identified by calculating of Ka/Ks ratio. Several immune-related genes were identified as positively selected genes, such as Notch2 and nfatc3b. Considering that these genes play crucial role in immune regulating function, the immune system genes met a great variation under environment selection in tilapias which suggests fast evolution in immune system of cultured tilapias.

Conservation genetics and genomics of amphibians and reptiles

Amphibians and reptiles as a group are often secretive, reach their greatest diversity often in remote tropical regions, and contain some of the most endangered groups of organisms on earth. Particularly in the past decade, genetics and genomics have been instrumental in the conservation biology of these cryptic vertebrates, enabling work ranging from the identification of populations subject to trade and exploitation, to the identification of cryptic lineages harboring critical genetic variation, to the analysis of genes controlling key life history traits. In this review, we highlight some of the most important ways that genetic analyses have brought new insights to the conservation of amphibians and reptiles. Although genomics has only recently emerged as part of this conservation tool kit, several large-scale data sources, including full genomes, expressed sequence tags, and transcriptomes, are providing new opportunities to identify key genes, quantify landscape effects, and manage captive breeding stocks of at-risk species.