Transcriptome-wide identification, molecular evolution and expression analysis of Toll-like receptor family in a Tibet fish, Gymnocypris przewalskii

Comprehensive transcriptomic analysis of Tibetan Schizothoracinae fish Gymnocypris przewalskii reveals how it adapts to a high altitude aquatic life

Background

Understanding the genetic basis of adaptation to high altitude life is of paramount importance for preserving and managing genetic diversity in highland animals. This objective has been addressed mainly in terrestrial fauna but rarely in aquatic animals. Tibetan Schizothoracinae fish is the ideal model system in evolutionary biology, carrying key insights into evolutionary genetics of speciation and adaptation at high altitude. Gymnocypris przewalskii is the newly formed Schizothoracinae fish species in the Tibetan Plateau, inhabits chronic cold, extreme saline and alkaline aquatic environment in Lake Qinghai, thus evolving the unique genomic signatures to adapt extremely severe environments.

Results

To characterize its genomic features, we assembled de novo transcriptome of G. przewalskii from Lake Qinghai. Intriguingly, by comparative genomic analyses of G. przewalskii and 8 other fish species, we identified potential expansions in gene families related to energy metabolism, transport and developmental functions, possibly underlying the adaptation to these environmental stresses. Through comprehensive molecular evolution analyses, we found that sets of genes controlling mitochondrion, ion homoeostasis, acid-base balance and innate immunity show significant signals of positive selection. Compared to previous studies on highland fishes, we failed to identify any positively selected genes related to hypoxia response.

Conclusions

Our findings provide comprehensive insights into the genetic basis of teleost fish that underlie their adaptation to extreme high altitude aquatic life on the Tibetan Plateau.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0925-z) contains supplementary material, which is available to authorized users.

Molecular cloning and characterization of LrTLR4, analysis of its inductive expression and associated down-stream signaling molecules following lipopolysaccharide stimulation and Gram-negative bacterial infection

Toll-like receptors (TLRs) play key roles in innate immunity from lower to higher vertebrates. Among various TLR types, TLR4 was reported to recognize LPS in higher vertebrates resulting in the activation of down-stream signaling pathway. Except in some teleosts, function of TLR4 in most fish species including rohu (Labeo rohita) a commercially important fish species in the South-East Asian countries remained unknown. To investigate it, full-length cDNA of Labeo rohita TLR4 (LrTLR4) was cloned, and it consisted of 2729 bp, with a single ORF of 2469 bp encoding a polypeptide of 822 aa with a predicted molecular mass of 94.753 kDa. Structurally, LrTLR4 consisted of 25 LRRs (leucine rich repeat regions), one TM (trans-membrane) domain and one TIR (Toll/interleukin-1 receptor) domain, and was similar to higher vertebrate's TLR4. Phylogenetically, LrTLR4 exhibited highest (85%) identity with the common carp TLR4b amino acids sequence, and formed a separate subgroup in the phylogenetic tree. LrTLR4 was widely expressed in all tested organs/tissues, and amidst the tissues highest expression was detected in blood and the lowest in eye. In response to LPS-stimulation, LrTLR4 was induced with the activation of MyD88-dependent and TRIF-dependent signaling pathway resulting in pro-inflammatory cytokines (interleukin 6 and 8) and type I IFN gene expression. Infection of rohu with a Gram-negative fish pathogen (Aeromonas hydrophila), also activated LrTLR4. Together, these findings suggest the important role of TLR4 in LPS sensing and augmentation of innate immunity against Gram-negative bacterial infection in fish.

Identification and functional analysis of the toll-like receptor 20.2 gene in grass carp, Ctenopharyngodon idella

We characterized and identified the cDNA sequence of Toll-like receptor 20.2 in Ctenopharyngodon idella (gctlr20.2); it consisted of 3197 bp, with an open reading frame of 2835 bp that encoded a 944 amino acid polypeptide. Relatively, high expression levels of gctlr20.2 were observed in the spleen, head kidney, liver and brain tissues, with lower expression levels in the trunk kidney, intestine and heart tissues. In vivo and in vitro, after being challenged with Aeromonas hydrophila or grass carp reovirus (GCRV), gctlr20.2 expression was induced in C. idella kidney cells stimulated with lipopolysaccharide, flagellin or polyinosinic-polycytidylic acid. Overexpression of gctlr20.2 increased the expression of il1β, il8 and tnf-α, but not ifn, and also increased the activity of the nf-κB signal pathway. Silencing, via siRNA-tlr20.2, inhibited gctlr20.2 transcription by 65.7% and down-regulated the expression of inflammatory cytokine genes, but not tnf-α. This study increases understanding of the immune system in C. idella.

Identification and characterization of a novel Toll-like receptor 4 homologue in blunt snout bream, Megalobrama amblycephala

Toll-like receptors (TLRs) are central players in the innate immune system in response to a wide range of pathogen infection. Among various TLRs, TLR4 plays a key role in recognition of bacterial lipopolysaccharides (LPS). In the present study, we identified and characterized a novel TLR4 homologue (maTLR4b) in blunt snout bream (Megalobrama amblycephala) which was significantly distinct from previously reported M. amblycephala TLR4 (tentatively named maTLR4a). The results showed that the complete cDNA sequence of maTLR4b was 3261 bp with an open reading frame encoding a polypeptide of 820 amino acids, and that its genomic sequence was 3793 bp, which had 3 exons. Structurally, the deduced maTLR4b protein showed a typical TLR domain architecture, including a signal peptide, eight leucine-rich repeats (LRRs) in the extracellular region, a transmembrane domain, and a Toll-Interleukin 1 receptor (TIR) domain in the cytoplasmic region. Phylogenetic analysis revealed that all TLR4s from teleost fish formed a monophyletic clade. Both maTLR4a and maTLR4b were divided into two distinct branches, and showed the highest level of similarity with the grass carp TLR4.2 and TLR4.4 homologue, respectively. MaTLR4b was constitutively expressed in all healthy tissues tested although at different levels. After LPS stimulation, the expression levels were significantly up-regulated in spleen, and peaked at 4 h between maTLR4a and maTLR4b, but with a distinct and complementary expression patterns. Taken together, these results suggested that maTLR4b is indeed a functional homologue of TLR4 in other species, which may play vital role in innate immune.

Identification and characterization of three TLR1 subfamily members from the orange-spotted grouper, Epinephelus coioides

Toll-like receptors (TLRs), which play important roles in host defense against pathogen infection, are the most intensively studied pattern recognition receptors (PRRs). In this study, we identified three novel TLR1 subfamily members, including TLR1 (EcTLR1b), TLR2 (EcTLR2b) and TLR14 (EcTLR14), from the orange-spotted grouper (Epinephelus coioides). EcTLR1b and EcTLR2b displayed low sequence identity with the previously reported grouper TLR1 (EcTLR1a) and TLR2 (EcTLR2a), respectively. The open reading frames (ORFs) of EcTLR1b, EcTLR2b and EcTLR14 contain 2484 bp, 2394 bp and 2640 bp, which encode the corresponding 827 amino acids (aa), 797 aa and 879 aa, respectively. All three TLRs have leucine-rich repeat (LRR) domains (including an LRR-NT (except for EcTLR1b), several LRR motifs and an LRR-CT), a trans-membrane region and a Toll/interleukin-1 receptor (TIR) domain. The TIR domains of the three TLRs exhibited conserved boxes, namely box1, box2 and box3, and their 3D models were similar to those of human TLR1 or TLR2. Sequence alignment demonstrated that the TIR domains of the three TLRs shared higher sequence identity with those of other species than the full-length receptors. Phylogenetic analysis indicated that EcTLR1s and EcTLR2s are characterized by their differing evolutionary status, whereas EcTLR14 was found to be in the same group as other piscine TLR14/18s. The three TLRs were ubiquitously expressed in seven tested tissues of healthy groupers, although their expression profiles were different. Post Cryptocaryon irritans infection, TLR1s expression was up-regulated in the gills. The expression of TLR2b was mainly increased in the spleen, but decreased in the gills, which was similar to the expression pattern of TLR2a post C. irritans infection. Unlike EcTLR1b and EcTLR2b, however, the grouper TLR14 transcript was substantially induced in both tissues post challenge. These findings may be helpful in understanding the innate immune mechanism of host anti-parasite infection.

Molecular cloning and expression of toll-like receptor 4 (tlr4) in the blunt snout bream (Megalobrama amblycephala)

Toll-like receptors (TLRs) play a pivotal role in teleost innate immune system. In this study, Megalobrama amblycephala (ma) tlr4 gene was cloned, its putative polypeptide product characterized, and expression analysed. Matlr4 cDNA is 2862 bp long, with an open reading frame of 2364 bp encoding 787 amino acids. MaTlr4 is a typical TLR protein, including the extracellular part with nine leucine-rich repeat motifs, a transmembrane region and a cytoplasmic Toll/interleukin-1 receptor domain. MaTlr4 has the highest level of identity (94%) and similarity (97%) with the grass carp Tlr4.2 homolog. This was also corroborated by the phylogenetic analysis, which placed MaTlr4 in a cluster with other cyprinid homologs. Matlr4 mRNA was ubiquitously expressed in all examined tissues and during all sampled developmental stages. The observed peak in matlr4 mRNA expression during gastrula and somite stages is in good agreement with its proposed role in the development of the neural system. Temporal expression patterns of matlr4 and maMyD88 mRNAs and proteins were analyzed in liver, spleen, head kidney, trunk kidney and intestine after Aeromonas hydrophila infection. And mRNA expression varied between different time-points. Both MaTlr4 and MaMyD88 protein expressions at 12 hpi were significantly enhanced in head kidney and intestine. These results indicate that matlr4 is involved in the immune response in M. amblycephala, and that it is indeed a functional homologue of tlr4s described in other animal species.

Positive Darwinian selection within interferon regulatory factor genes of Gymnocypris przewalskii (Cyprinidae) on the Tibetan Plateau

Tibetan Plateau (TP) had experienced phased uplift, resulting in inhospitable environment of low temperature, hypoxia and high ultraviolet radiation for Tibetan wildlife. Many organisms can well adapt to TP, it is of ecological and evolutionary interest to untangle how organisms adapt to extreme environment on TP through evolution. Previous studies mainly focused on hypoxia and metabolism related genes, but we know little about the evolutionary history of immune genes in Tibetan wildlife. In this study, we first identified 10 interferon regulatory factor (IRF) genes from Tibetan naked carp Gymnocypris przewalskii. Within this gene family, IRF3, IRF5, IRF7 and IRF8 contained positive selection sites. Evidences indicated that positive selection may lead to IRF genes functional alternations, presumably driving genes towards adaptation to the environmental changes. Taken together, our results suggested 4 candidate genes as interesting targets for further experimental confirmation of their functional variations and contributions to high altitude adaptation in Tibet fish.

Identification and characterization of a novel Toll-like receptor 2 homologue in the large yellow croaker Larimichthys crocea

Toll-like receptors (TLRs) are key components of innate immunity that play significant roles in immune defence against pathogen invasion. In the present study, we identified a novel TLR2 homologue (LycTLR2b) in large yellow croaker (Larimichthys crocea) that shared low sequence identity with the previously reported large yellow croaker TLR2 (tentatively named LycTLR2a). The full-length cDNA of LycTLR2b was 2926 nucleotides (nt) long and encoded a protein consisting of 797 amino acids (aa). The deduced LycTLR2b protein exhibited a typical TLR domain architecture including a signal peptide, seven leucine-rich repeats (LRRs) in the extracellular region, a transmembrane domain, and a Toll-Interleukin 1 receptor (TIR) domain in the cytoplasmic region. Phylogenetic analysis showed that both LycTLR2a and LycTLR2b fall into a major clade formed by all TLR2 sequences, and are divided into two distinct branches. Genomic organization revealed that the LycTLR2b gene lacks intron, which is similar to zebrafish and human TLR2 genes, whereas the LycTLR2a gene contains multiple introns, as found in damselfish TLR2a and Fugu TLR2 genes. Syntenic analysis suggested that the occurrence of LycTLR2a and LycTLR2b may result from a relatively recent genome duplication event. LycTLR2b mRNA was constitutively expressed in all tissues examined although at different levels. Following bacterial vaccine challenge, LycTLR2b expression levels were significantly up-regulated in both spleen and head kidney tissues. Taken together, these results indicated that two different TLR2 homologues, which may play roles in antibacterial immunity, exist in large yellow croaker.