Conserved regulation of the lymphocyte-specific expression of lck in the Fugu and mammals

Environmental salinity influences the branchial expression of TCR pathway related genes based on transcriptome of a catadromous fish

Environmental salinity not only affects the physiological processes such as osmoregulation and hormonal control, but also changes the immune system in fishes. Studies are limited in fish on the roles of the T cell receptor (TCR)-related genes in relation to changes in environmental salinity. A large group of salinity-challenged transcripts was obtained in gills of marbled eel (Anguilla marmorata). Moreover, bioinformatic ways were used to identify the enriched TCR pathway related genes which were significantly different expressed in fresh water (FW), brackish water (BW) and seawater (SW). Meanwhile, the RT-qPCR results were validated and consistent with the RNA-seq results. TCR a, TCR b, CD45, CD28, PI3K, LCK and LAT were up-regulated when the salinity increases in BW and SW, which connected with the related signaling pathways (Ras-MAPK and PKC pathway). CD4 and Zap70 were down-regulated when the salinity increases in BW and SW, which connected with the PLC pathway. The research offers a novel viewpoint to explore the immune pathways including the TCR pathway in fish based on transcriptome.

Lyn is involved in host defense against S. agalactiae infection and BCR signaling in Nile tilapia (Oreochromis niloticus)

Lyn, a member of Src protein kinase family, plays a crucial role in immune reactions against pathogenic infection. In this study, Lyn from Nile tilapia (Oreochromis niloticus) (OnLyn) was identified and characterized at expression pattern against bacterial infection, and regulation function in BCR signaling. The open reading frame of OnLyn contained 1536 bp of nucleotide sequence encoded a protein of 511 amino acids. The OnLyn protein was highly conversed to other species Lyn, including SH3, SH2 and a catalytic Tyr kinase (TyrKc) domain. Transcriptional expression analysis revealed that OnLyn was detected in all examined tissues and was highly expressed in the head kidney. The up-regulation OnLyn expression was observed in the head kidney and spleen following challenge with Streptococcus agalactiae (S. agalactiae) in vivo, and was also displayed in head kidney leukocytes challenge with S. agalactiae and LPS in vitro. In addition, after induction with mouse anti-OnIgM mAb in vitro, the OnLyn expression and phosphorylation of OnLyn (Y⁵⁰⁷) were significantly up-regulated in the head kidney leukocytes. Moreover, after treatment with AZD0530 and mouse anti-OnIgM monoclonal antibody, the down-regulation of cytoplasmic free-Ca²⁺ concentration was detected in the head kidney leukocytes in vitro. Taken together, the findings of this study revealed that OnLyn might play potential roles in BCR signaling and get involved in host defense against bacterial infection in Nile tilapia.

Transcriptionally promiscuous "blurry" promoters in Tc1/mariner transposons allow transcription in distantly related genomes

Background

We have recently described a peculiar feature of the promoters in two Drosophila Tc1-like elements, Bari1 and Bari3. The AT-richness and the presence of weak core-promoter motifs make these promoters, that we have defined “blurry”, able to activate transcription of a reporter gene in cellular systems as diverse as fly, human, yeast and bacteria. In order to clarify whether the blurry promoter is a specific feature of the Bari transposon family, we have extended this study to promoters isolated from three additional DNA transposon and from two additional LTR retrotransposons.

Results

Here we show that the blurry promoter is also a feature of two vertebrate transposable elements, Sleeping Beauty and Hsmar1, belonging to the Tc1/mariner superfamily. In contrast, this feature is not shared by the promoter of the hobo transposon, which belongs to the hAT superfamily, nor by LTR retrotransposon-derived promoters, which, in general, do not activate transcription when introduced into non-related genomes.

Conclusions

Our results suggest that the blurry promoter could be a shared feature of the members of the Tc1/mariner superfamily with possible evolutionary and biotechnological implications.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0155-6) contains supplementary material, which is available to authorized users.

Grouper (Epinephelus coioides) TCR signaling pathway was involved in response against Cryptocaryon irritans infection

T cell activation is a complicated process accompanying with the activation of T cell receptor (TCR) signaling pathway, which is not well described in teleost fish. The initiation of this pathway depends on the interaction of membrane TCR co-receptors (e.g. CD4/8, CD3 and CD45) and a series of cytoplasmic protein tyrosine kinases (e.g. Lck, Fyn and ZAP70). Cyptocaryon irritans is a ciliate pathogen of marine fish white spot disease causing huge economic lost in marine aquaculture. This parasite can infect fish gill and skin and is considered to be a good pathogen model for fish gill and skin mucosal immunity. Our previous studies showed the locally mucosal antibody response was important for fish defense against this parasite. While how TCR signaling pathway involved in T cell activation to help B cell activation in C. irritans infected fish is still not known. In the present study, we cloned a grouper TCR co-receptor gene EcCD3ε (537 bp) and its three kinase genes, including EcLck (1512 bp), EcFyn (1605 bp) and EcZAP70 (1893 bp). Homology analysis showed that they all shared the highest identity with corresponding genes from Takifugu rubripes (EcCD3ε 41%, EcLck 88%, EcFyn 98% and EcZAP70 93%), and their conserved motifs involved in the signaling transduction were analyzed. The tissue distribution analysis showed these four genes were high expressed in thymus, and it is interesting to find their comparative high expression in skin, gill and midgut mucosal immune tissues. In C. irritans infected grouper, the expression of three TCR co-receptors (EcCD4-1, EcCD3ε and EcCD45) and three kinases (EcLck, EcFyn and EcZAP70) was tested in skin, gill, head kidney and spleen at 0, 12 h, 24 h, 2 d, 3 d, 5 d and 7 d. All six genes were significantly up-regulated in skin at most tested time points, which indicate the possibility of skin local T cell activation to support the local antibody response. Compared to three TCR co-receptors, significantly up-regulation of three kinases were seen in the spleen, and the spleen fold changes of these three kinases were much higher than head kidney, which indicates spleen maybe the major systematic immune organs for T cell activation in C. irritans infected fish.

Molecular cloning and characterization of lymphocyte cell kinase from humphead snapper (Lutjanus sanguineus)

Lymphocyte cell kinase (LCK) belongs to the Src family of tyrosine kinases, which involves in the proliferation control of lymphocytes. In this study, we cloned the LCK gene of humphead snapper (Lutjanus sanguineus) (designed as LsLCK). Sequence analysis showed that the full-length cDNA of LsLCK was 2279 bp, contained a 1506-bp open reading frame (ORF), encoding a polypeptide of 501 amino acids. The deduced amino acid possessed the typical structural features of known LCK proteins, including four Src homology (SH) domains arranged as the SH1 domain followed by a regulatory C-terminal tail (COOH-domain), SH2 and SH3 adapter domains and SH4 domain which required for membrane attachment and CD4/CD8 binding. Fluorescent quantitative real-time PCR analysis indicated that LsLCK transcripts were expressed mainly in thymus, spleen and head kidney in healthy fish. Moreover, the mRNA expressions in these tissues were significantly up-regulated after challenge with Vibrio harveyi. The results of immunohistochemistry showed that LsLCK protein localized distinctly in cytoplasm of cell in thymus, spleen and head kidney. Taken together, these findings indicated that LsLCK may play an important role in the immune response of humphead snapper against bacterial infection.

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.

A linkage map of transcribed single nucleotide polymorphisms in rohu (Labeo rohita) and QTL associated with resistance to Aeromonas hydrophila

BACKGROUND

Production of carp dominates world aquaculture. More than 1.1 million tonnes of rohu carp, Labeo rohita (Hamilton), were produced in 2010. Aeromonas hydrophila is a bacterial pathogen causing aeromoniasis in rohu, and is a major problem for carp production worldwide. There is a need to better understand the genetic mechanisms affecting resistance to this disease, and to develop tools that can be used with selective breeding to improve resistance. Here we use a 6 K SNP array to genotype 21 full-sibling families of L. rohita that were experimentally challenged intra-peritoneally with a virulent strain of A. hydrophila to scan the genome for quantitative trait loci associated with disease resistance.

RESULTS

In all, 3193 SNPs were found to be informative and were used to create a linkage map and to scan for QTL affecting resistance to A. hydrophila. The linkage map consisted of 25 linkage groups, corresponding to the number of haploid chromosomes in L. rohita. Male and female linkage maps were similar in terms of order, coverage (1384 and 1393 cM, respectively) and average interval distances (1.32 and 1.35 cM, respectively). Forty-one percent of the SNPs were annotated with gene identity using BLAST (cut off E-score of 0.001). Twenty-one SNPs mapping to ten linkage groups showed significant associations with the traits hours of survival and dead or alive (P <0.05 after Bonferroni correction). Of the SNPs showing significant or suggestive associations with the traits, several were homologous to genes of known immune function or were in close linkage to such genes. Genes of interest included heat shock proteins (70, 60, 105 and "small heat shock proteins"), mucin (5b precursor and 2), lectin (receptor and CD22), tributyltin-binding protein, major histocompatibility loci (I and II), complement protein component c7-1, perforin 1, ubiquitin (ligase, factor e4b isoform 2 and conjugation enzyme e2 c), proteasome subunit, T-cell antigen receptor and lymphocyte specific protein tyrosine kinase.

CONCLUSIONS

A panel of markers has been identified that will be validated for use with both genomic and marker-assisted selection to improve resistance of L. rohita to A. hydrophila.

Fish T cells: recent advances through genomics

This brief review is intended to provide a concise overview of the current literature concerning T cells, advances in identifying distinct T cell functional subsets, and in distinguishing effector cells from memory cells. We compare and contrast a wealth of recent progress made in T cell immunology of teleost, elasmobranch, and agnathan fish, to knowledge derived from mammalian T cell studies. From genome studies, fish clearly have most components associated with T cell function and we can speculate on the presence of putative T cell subsets, and the ability to detect their differentiation to form memory cells. Some recombinant proteins for T cell associated cytokines and antibodies for T cell surface receptors have been generated that will facilitate studying the functional roles of teleost T cells during immune responses. Although there is still a long way to go, major advances have occurred in recent years for investigating T cell responses, thus phenotypic and functional characterization is on the near horizon.

T cell diversity and TcR repertoires in teleost fish

In vertebrates, the diverse and extended range of antigenic motifs is matched to large populations of lymphocytes. The concept of immune repertoire was proposed to describe this diversity of lymphocyte receptors--IG and TR--required for the recognition specificity. Immune repertoires have become useful tools to describe lymphocyte and receptor populations during the immune system development and in pathological situations. In teleosts, the presence of conventional T cells was first proposed to explain graft rejection and optimized specific antibody production. The discovery of TR genes definitely established the reality of conventional T cells in fish. The development of genomic and EST databases recently led to the description of several key T cell markers including CD4, CD8, CD3, CD28, CTLA4, as well as important cytokines, suggesting the existence of different T helper (Th) subtypes, similar to the mammalian Th1, Th2 and Th17. Over the last decade, repertoire studies have demonstrated that both public and private responses occur in fish as they do in mammals, and in vitro specific cytotoxicity assays have been established. While such typical features of T cells are similar in both fish and mammals, the structure of particular repertoires such as the one of gut intra-epithelial lymphocytes seems to be very different. Future studies will further reveal the particular characteristics of teleost T cell repertoires and adaptive responses.

Description of an elasmobranch TCR coreceptor: CD8α from Rhinobatos productus

Cell-mediated immunity plays an essential role for the control and eradication of intracellular pathogens. To learn more about the evolutionary origins of the first signal (Signal 1) for T-cell activation, we cloned CD8α from an elasmobranch, Rhinobatos productus. Similar to full-length CD8α cDNAs from other vertebrates, Rhpr-CD8α (1800bp) encodes a 219 amino acid open reading frame composed of a signal peptide, an extracellular IgSF V domain and a stalk/hinge region followed by a well-conserved transmembrane domain and cytoplasmic tail. Overall, the mature Rhpr-CD8α protein (201 aa) displays ∼ 30% amino acid identity with mammalian CD8α including absolute conservation of cysteine residues involved in the IgSf V domain fold and dimerization of CD8αα and CD8αβ. One prominent feature is the absence of the LCK association motif (CXC) that is needed for achieving signal 1 in tetrapods. Both elasmobranch and teleost CD8α protein sequences possess a similar but distinctly different motif (CXH) in the cytoplasmic tail. The overall genomic structure of CD8α has been conserved during the course of vertebrate evolution both for the number of exons and phase of splicing. Finally, quantitative RTPCR demonstrated that elasmobranch CD8α is expressed in lymphoid-rich tissues similar to CD8 in other vertebrates. The results from this study indicate the existence of CD8 prior to the emergence of the gnathostomes (>450 MYA) while providing evidence that the canonical LCK association motif in mammals is likely a derived characteristic of tetrapod CD8α, suggesting potential differences for T-cell education and activation in the various gnathostomes.

Enforced expression of Simian virus 40 large T-antigen leads to testicular germ cell tumors in zebrafish

Testicular germ cell tumors (TGCTs) are the most common malignancy in young men. However, there are few in vivo animal models that have been developed to study this disease. We have used the pufferfish (fugu) lymphocyte-specific protein tyrosine kinase (flck) promoter, which has been shown to enforce high-level expression in the testes of transgenic mice, to express Simian virus 40 large T-antigen in zebrafish testes. Zebrafish that express T-antigen develop TGCTs after a long latency of >1 year. Although overt TGCTs are only evident in 20% of the fish, occult TGCTs can be detected in 90% of the transgenic fish by 36 month of age. The TGCTs resemble the human disease in terms of morphology and gene expression pattern, and can be transplanted to healthy wild-type recipient fish. In addition, enforced expression of the zebrafish stem cell leukemia (scl) gene in the zebrafish testes also generated TGCTs in transgenic fish. These results demonstrate the feasibility of studying TGCTs in a model organism.

Fugu rubripes and human survival motor neuron genes: structural and functional similarities in comparative genome studies

The compactness of the Fugu rubripes (Fugu) genome has supported its use in comparative genome analysis. Nevertheless, as Fugu is distinct evolution-wise from humans, it is essential to determine the similarity between a Fugu gene and its human counterpart to confirm its potential for comparative genome analysis. We cloned and analyzed the Fugu survival motor neuron gene (fsmn) for similarities with human SMN gene (huSMN). The Fugu genome has a single fsmn that is 13.4 times smaller than huSMN. fsmn and huSMN are highly similar in their genome organization and tissue expression patterns. The functional domains of the Fugu smn and human SMN molecules are also highly conserved. In human MCF-7 cells, expression of fsmn protein resulted in the formation of "gems" in the cytoplasm and nucleus, similar to observations reported for huSMN protein. In these cells, fsmn RNA was also processed correctly and produced alternatively spliced transcripts like huSMN2. These findings indicate close structural and functional similarities between fsmn and huSMN, suggesting that regulation of the two genes may also be similar and supporting the use of fsmn in comparative genome studies for the identification of functional regulatory elements of huSMN.

Genomic organization and expression of CD8alpha and CD8beta genes in fugu Takifugu rubripes

We have cloned cDNAs encoding the alpha and beta chains of CD8 from the tiger pufferfish (fugu), Takifugu rubripes. The cDNA sequences encode a putative leader peptide, extracellular immunoglobulin variable region-like domain, stalk region, transmembrane region, and cytoplasmic tail. A protein tyrosine kinase p56lck binding motif was not found in the putative fugu CD8alpha cytoplasmic tail. O-linked glycosylation sites were found in the stalk of both CD8 chains, suggesting possible stalk formation. Phylogenetic analysis showed that fugu CD8alpha and CD8beta chains cluster with other vertebrate CD8alpha and CD8beta chains, respectively. The fugu CD8 genes comprise six exons separated by five introns. The genes are tandemly aligned 3.6 kb apart and are in the same transcription orientation. Quantitative RT-PCR analysis demonstrated that fugu CD8 is expressed predominantly in lymphoid tissues. In situ hybridization showed that fugu CD8 genes are expressed in thymocytes and lymphocytes within lymphoid organs. Molecular characterization of CD8 in fish provides the basis for development of specific antibodies to identify T-cell subsets, as well as potentially important insights into the evolution of CD8 and the adaptive immunity.

SCL-GFP transgenic zebrafish: in vivo imaging of blood and endothelial development and identification of the initial site of definitive hematopoiesis

The bHLH transcription factor SCL plays a central role in the generation of hematopoietic cells in vertebrates. We modified a PAC containing the whole zebrafish scl locus, inserting GFP into the first coding exon of scl. In germline-transgenic zebrafish generated using this construct, GFP expression completely recapitulates the endogenous expression of scl in blood, endothelium and CNS. We performed in vivo timelapse imaging of blood and endothelial precursor migration at the single-cell level and show that these cells migrate from the posterior lateral plate mesoderm to their site of differentiation in the intermediate cell mass. The anterior lateral plate domain of GFP expression gives rise to primitive macrophages and the blood vessels of the head. In later embryos, GFP expression identifies clusters of hematopoietic cells that develop between the dorsal aorta and posterior cardinal veins after primitive erythrocytes have entered circulation. Two treatments that block definitive hematopoiesis (treatment with dioxin (TCDD), and injection of an antisense morpholino oligonucleotide targeted to runx1) ablate these hematopoietic clusters. This indicates that these clusters represent the first site of definitive hematopoiesis in zebrafish. This site is anatomically homologous to the proposed source of hematopoietic stem cells in amniotes, the aorta-gonad-mesonephros (AGM) region. A second transgenic line, containing the promoter of scl driving GFP, lacks expression in the definitive clusters.

Molecular and biochemical analysis of rainbow trout LCK suggests a conserved mechanism for T-cell signaling in gnathostomes

Two genes were identified in rainbow trout that display high sequence identity to vertebrate Lck. Both of the trout Lck transcripts are associated with lymphoid tissues and were found to be highly expressed in IgM-negative lymphocytes. In vitro analysis of trout lymphocytes indicates that trout Lck mRNA is up-regulated by T-cell mitogens, supporting an evolutionarily conserved function for Lck in the signaling pathways of T-lymphocytes. Here, we describe the generation and characterization of a specific monoclonal antibody raised against the N-terminal domains of recombinant trout Lck that can recognize Lck protein(s) from trout thymocyte lysates that are similar in size ( approximately 57kDa) to mammalian Lck. This antibody also reacted with permeabilized lymphocytes during FACS analysis, indicating its potential usage for cellular analyses of trout lymphocytes, thus representing an important tool for investigations of salmonid T-cell function.

Evolution of the CD4 family: teleost fish possess two divergent forms of CD4 in addition to lymphocyte activation gene-3

The T cell coreceptor CD4 is a transmembrane glycoprotein belonging to the Ig superfamily and is essential for cell-mediated immunity. Two different genes were identified in rainbow trout that resemble mammalian CD4. One (trout CD4) encodes four extracellular Ig domains reminiscent of mammalian CD4, whereas the other (CD4REL) codes for two Ig domains. Structural motifs within the amino acid sequences suggest that the two Ig domains of CD4REL duplicated to generate the four-domain molecule of CD4 and the related gene, lymphocyte activation gene-3. Here we present evidence that both of these molecules in trout are homologous to mammalian CD4 and that teleosts encode an additional CD4 family member, lymphocyte activation gene-3, which is a marker for activated T cells. The syntenic relationships of similar genes in other teleost and non-fish genomes provide evidence for the likely evolution of CD4-related molecules in vertebrates, with CD4REL likely representing the primordial form in fish. Expression of both CD4 genes is highest in the thymus and spleen, and mRNA expression of these genes is limited to surface IgM- lymphocytes. consistent with a role for T cell functionality. Finally, the intracellular regions of both CD4 and CD4REL possess the canonical CXC motif involved in the interaction of CD4 with p56LCK, implying that similar mechanisms for CD4+ T cell activation are present in all vertebrates. Our results therefore raise new questions about T cell development and functionality in lower vertebrates that cannot be answered by current mammalian models and, thus, is of fundamental importance for understanding the evolution of cell-mediated immunity in gnathosomes.

Characterisation of salmon and trout CD8alpha and CD8beta

The genes and corresponding cDNAs of both alpha and beta chains of the Atlantic salmon (Salmo salar) CD8 molecule have been sequenced and characterized. In addition, the cDNAs for alpha and beta chains of brown trout (Salmo trutta) and for the beta chain in rainbow trout (Oncorhynchus mykiss) have been sequenced. The cDNAs code for signal sequences which are preceded by short 5' UTRs. These are followed by typical immunoglobulin superfamily variable sequences all of which contain two conserved cysteines for the intra-chain disulphide bond. The hinge regions display conserved cysteines for dimerisation and several O-glycosylation motifs for each predicted protein. The domain sharing the highest sequence identity with mammals is the single pass transmembrane domain for all sequences. In salmon, each domain is predominantly coded for by a single exon except the cytoplasmic/3' UTR domains, which are coded for by 3 and 2 exons for the alpha and beta genes, respectively. In the alpha gene, the second cytoplasmic exon may be spliced out to form an alternative shorter transcript which if expressed would exhibit a truncated cytoplasmic tail. A splice variant found for the salmon beta gene introduces a stop codon after only 40 amino acids. Overall amino acid identities between salmonid sequences were higher than 90%, whereas they shared only 15-20% identity with species such as, chicken and human. Analysis of the expression patterns of the two salmon genes using quantitative RT-PCR shows a very high expression in the thymus. This is mirrored by the expression of the TCRalpha gene, which is known to be co-expressed with CD8 on mammalian T cells. This is the first report of a sequence for CD8beta in a teleost and together with the CD8alpha sequence, it encodes the ortholog of the CD8 co-receptor molecule on mammalian T cells.

Identification of neuronal enhancers of the proopiomelanocortin gene by transgenic mouse analysis and phylogenetic footprinting

The proopiomelanocortin (POMC) gene is expressed in the pituitary and arcuate neurons of the hypothalamus. POMC arcuate neurons play a central role in the control of energy homeostasis, and rare loss-of-function mutations in POMC cause obesity. Moreover, POMC is the prime candidate gene within a highly significant quantitative trait locus on chromosome 2 associated with obesity traits in several human populations. Here, we identify two phylogenetically conserved neuronal POMC enhancers designated nPE1 (600 bp) and nPE2 (150 bp) located approximately 10 to 12 kb upstream of mammalian POMC transcriptional units. We show that mouse or human genomic regions containing these enhancers are able to direct reporter gene expression to POMC hypothalamic neurons, but not the pituitary of transgenic mice. Conversely, deletion of nPE1 and nPE2 in the context of the entire transcriptional unit of POMC abolishes transgene expression in the hypothalamus without affecting pituitary expression. Our results indicate that the nPEs are necessary and sufficient for hypothalamic POMC expression and that POMC expression in the brain and pituitary is controlled by independent sets of enhancers. Our study advances the understanding of the molecular nature of hypothalamic POMC neurons and will be useful to determine whether polymorphisms in POMC regulatory regions play a role in the predisposition to obesity.

Identification and characterization of a second CD4-like gene in teleost fish

In fish, T cell subdivision is not well studied, although CD8 and CD4 homologues have been reported. This study describes a second teleost CD4-like gene, CD4-like 2 (CD4L-2). Two rainbow trout copies of this gene were found, -2a and -2b, encoding molecules sharing 81% aa identity. The 2a/2b duplication may be related to tetraploid ancestry of salmonid fishes. In the Fugu genome CD4L-2 lies head to tail with an earlier reported, very different CD4-like gene [Suetake, H., Araki, K., Suzuki, Y., 2004. Cloning, expression, and characterization of fugu CD4, the first ectothermic animal CD4. Immunogenetics 56, 368-374], which was designated CD4L-1 in the present article. The flanking genes of the Fugu CD4L-1 and CD4L-2 are reminiscent of the genes surrounding CD4 and LAG-3 in mammals. However, neither synteny nor phylogenetic analysis could decide between CD4 and LAG-3 identity for the fish CD4L genes. CD4L-1 and CD4L-2 share a tyrosine protein kinase p56(lck) binding motif in the cytoplasmic tail with CD4 but not with LAG-3. Trout CD4L-2 expression is highest in the thymus, similar to mammalian and chicken CD4, whereas Fugu CD4L-1 expression was highest in the spleen. However, CD4L-2 encodes only two IG-like domains, whereas CD4L-1, CD4 and LAG-3 encode four. The CD4-like genes 1 and 2 in fish apparently went through an evolution different from that of LAG-3 and CD4 in higher vertebrates.

Comparative genomics using fugu: a tool for the identification of conserved vertebrate cis-regulatory elements

With the imminent completion of the whole genome sequence of humans, increasing attention is being focused on the annotation of cis-regulatory elements in the human genome. Comparative genomics approaches based on evolutionary conservation have proved useful in the detection of conserved cis-regulatory elements. The pufferfish, Fugu rubripes, is an attractive vertebrate model for comparative genomics, by virtue of its compact genome and maximal phylogenetic distance from mammals. Fugu has lost a large proportion of nonessential DNA, and retained single orthologs for many duplicate genes that arose in the fish lineage. Non-coding sequences conserved between fugu and mammals have been shown to be functional cis-regulatory elements. Thus, fugu is a model fish genome of choice for discovering evolutionarily conserved regulatory elements in the human genome. Such evolutionarily conserved elements are likely to be shared by all vertebrates, and related to regulatory interactions fundamental to all vertebrates. The functions of these conserved vertebrate elements can be rapidly assayed in mammalian cell lines or in transgenic systems such as zebrafish/medaka and Xenopus, followed by validation of crucial elements in transgenic rodents.

A novel cell-based binding assay system reconstituting interaction between SARS-CoV S protein and its cellular receptor

Severe acute respiratory syndrome (SARS), a life-threatening disease, is caused by the newly identified virus SARS coronavirus (SARS-CoV). In order to study the spike (S) protein of this highly contagious virus, we established a clonal cell-line, CHO-SG, from the Chinese hamster ovary cells that stably expresses C-terminally EGFP-tagged SARS-CoV S protein (S-EGFP). The ectodomain of the S glycoprotein is localized on the surface of CHO-SG cells with N-acetyl-glucosamine-terminated carbohydrate structure. CHO-SG cells associated tightly with Vero E6 cells, a SARS-CoV receptor (ACE2) expressing cell-line, and the interaction remained stable under highly stringent condition (1M NaCl). This interaction could be blocked by either the serum from a SARS convalescent patient or a goat anti-ACE2 antibody, indicating that the interaction is specific. A binding epitope with lesser degree of glycosylation and native conformation was localized by using rabbit anti-sera raised against five denatured recombinant S protein fragments expressed in Escherichia coli. One of the sera obtained from the fragment encompassing amino acids 48-358 significantly blocked the interaction between CHO-SG and Vero E6 cells. The region is useful for studying neutralizing antibodies in future vaccine development. This paper describes an easy and safe cell-based assay suitable for studying the binding between SARS-CoV S protein and its receptor.

An enigmatic fourth runt domain gene in the fugu genome: ancestral gene loss versus accelerated evolution

Background

The runt domain transcription factors are key regulators of developmental processes in bilaterians, involved both in cell proliferation and differentiation, and their disruption usually leads to disease. Three runt domain genes have been described in each vertebrate genome (the RUNX gene family), but only one in other chordates. Therefore, the common ancestor of vertebrates has been thought to have had a single runt domain gene.

Results

Analysis of the genome draft of the fugu pufferfish (Takifugu rubripes) reveals the existence of a fourth runt domain gene, FrRUNT, in addition to the orthologs of human RUNX1, RUNX2 and RUNX3. The tiny FrRUNT packs six exons and two putative promoters in just 3 kb of genomic sequence. The first exon is located within an intron of FrSUPT3H, the ortholog of human SUPT3H, and the first exon of FrSUPT3H resides within the first intron of FrRUNT. The two gene structures are therefore "interlocked". In the human genome, SUPT3H is instead interlocked with RUNX2. FrRUNT has no detectable ortholog in the genomes of mammals, birds or amphibians. We consider alternative explanations for an apparent contradiction between the phylogenetic data and the comparison of the genomic neighborhoods of human and fugu runt domain genes. We hypothesize that an ancient RUNT locus was lost in the tetrapod lineage, together with FrFSTL6, a member of a novel family of follistatin-like genes.

Conclusions

Our results suggest that the runt domain family may have started expanding in chordates much earlier than previously thought, and exemplify the importance of detailed analysis of whole-genome draft sequence to provide new insights into gene evolution.

Cloning, expression, and characterization of fugu CD4, the first ectothermic animal CD4

We have cloned and sequenced the first ectothermic animal CD4 gene from fugu, Takifugu rubripes, using a public database of the third draft sequence of the fugu genome. The fugu CD4 gene encodes a predicted protein of 463 amino acids containing four extracellular immunoglobulin (Ig)-like domains, a transmembrane region, and a cytoplasmic tail. Fugu CD4 shares low identity of about 15-20% with avian and mammalian CD4 proteins. Unlike avian and mammalian CD4, fugu CD4 lacks the Cys pair of the first Ig-like domain, but has a unique possible disulfide bond in the third domain. These differences suggest that fugu CD4 may have a different structure that could affect binding of major histocompatibility complex class II molecules and subsequent T-cell activation. In the putative fugu cytoplasmic region, the protein tyrosine kinase p56lck binding motif is conserved. The predicted fugu CD4 gene is composed of 12 exons, differing from other CD4 genes, but showing conserved synteny and many conserved sequence motifs in the promoter region. RT-PCR analysis demonstrated that the fugu CD4 gene is expressed predominantly in lymphoid tissues. We also show that fugu CD4 can be expressed on the surface of cells via transfection. Molecular characterization of CD4 in fish provides insights into the evolution of both the CD4 molecule and the immune system.

Identification of 5'-upstream region of pufferfish ribosomal protein L29 gene as a strong constitutive promoter to drive GFP expression in zebrafish

The genomic structure of Tetraodon fluviatilis L29 gene was determined and its promoter activity was analyzed in COS-1 cells and zebrafish embryos. The TfL29 gene comprises four exons and three introns, spanning approximately 1.7kb. The 5(')-upstream 2.2-kb of the first exon contains 10 E-boxes and many putative binding motifs for transcription factors GATA-1, AML-1a, c-Myb, Oct-1, CdxA, and NRF-2. Promoter activity assay showed that the distal 2.2-kb fragment not only had high luciferase activity in COS-1 cells, but also strong and ubiquitous GFP expression in a variety of tissues in zebrafish embryos. On the other hand, there are no TATA or CAAT boxes within a 300-bp region upstream from the transcription initiation site. Although this region has high luciferase activity in COS-1 cells, it is not sufficient to drive GFP expression in zebrafish embryos. In this proximal 300-bp region, there are two E-boxes, two CdxA sites, and one NRF-2 site that is immediately downstream of the transcription start site.

Compact intergenic regions of the pufferfish genome facilitate isolation of gene promoters: characterization ofFugu3′-phosphoadenosine 5′-phosphosulfate synthase 2 (fPapss2) gene promoter function in transgenicXenopus

The highly compact nature of the pufferfish (Fugu rubripes) genome renders it a useful tool not only for annotating coding regions within vertebrate genomes, but also for the identification of sequences important to gene regulation. Indeed, owing to this compaction it will be feasible in many instances to initiate analyses using entire intergenic regions when mapping gene promoters; a strategy that is very rarely feasible with the expanded genomes of other species. Stemming from our interest in studying promoters expressed in chondrocytes, we selected for study the intergenic region upstream of Fugu 3'-phosphoadenosine 5'-phosphosulfate synthase 2, fPapss2, a gene required for the normal development of cartilage extracellular matrix. Functional characterization of the entire fPapss2 5' intergenic region was carried out by monitoring expression of the enhanced green fluorescent protein (EGFP) gene reporter in the developing cartilage of transgenic Xenopus laevis. By evaluating a series of 5' intergenic region deletions we defined a minimal fPapss2 sequence of approximately 300 bp that was essential for EGFP expression in tadpole cartilage. This functional analysis of an entire Fugu intergenic region, combined with the efficiency of Xenopus transgenesis, serves as a model for the rapid characterization of evolutionarily-conserved regulatory regions of other pufferfish genes.

Fugu ESTs: new resources for transcription analysis and genome annotation

The draft Fugu rubripes genome was released in 2002, at which time relatively few cDNAs were available to aid in the annotation of genes. The data presented here describe the sequencing and analysis of 24,398 expressed sequence tags (ESTs) generated from 15 different adult and juvenile Fugu tissues, 74% of which matched protein database entries. Analysis of the EST data compared with the Fugu genome data predicts that approximately 10,116 gene tags have been generated, covering almost one-third of Fugu predicted genes. This represents a remarkable economy of effort. Comparison with the Washington University zebrafish EST assemblies indicates strong conservation within fish species, but significant differences remain. This potentially represents divergence of sequence in the 5' terminal exons and UTRs between these two fish species, although clearly, complete EST data sets are not available for either species. This project provides new Fugu resources, and the analysis adds significant weight to the argument that EST programs remain an essential resource for genome exploitation and annotation. This is particularly timely with the increasing availability of draft genome sequence from different organisms and the mounting emphasis on gene function and regulation.

Characterization of the pufferfish Otx2 cis-regulators reveals evolutionarily conserved genetic mechanisms for vertebrate head specification

The Otx2 gene, containing a highly conserved paired-type homeobox, plays a pivotal role in the development of the rostral head throughout vertebrates. Precise regulation of the temporal and spatial expression of Otx2 is likely to be crucial for proper head specification. However, regulatory mechanisms of Otx2 expression remain largely unknown. In this study, the Otx2 genome of the puffer fish Fugu rubripes, which has been proposed as a model vertebrate owing to its highly compact genome, was cloned. Consistently, Fugu Otx2 possesses introns threefold smaller in size than those of the mouse Otx2 gene. Otx2 mRNA was transcribed after MBT, and expressed in the rostral head region throughout the segmentation and pharyngula periods of wild-type Fugu embryos. To elucidate regulatory mechanisms of Otx2 expression, the expression of Otx2-lacZ reporter genes nearly covering the Fugu Otx2 locus, from -30.5 to +38.5 kb, was analyzed, by generating transgenic mice. Subsequently, seven independent cis-regulators were identified over an expanse of 60 kb; these regulators are involved in the mediation of spatiotemporally distinct subdomains of Otx2 expression. Additionally, these expression domains appear to coincide with local signaling centers and developing sense organs. Interestingly, most domains do not overlap with one another, which implies that cis-regulators for redundant expression may be abolished exclusively in the pufferfish so as to reduce its genome size. Moreover, these cis-regions were also able to direct expression in zebrafish embryos equivalent to that observed in transgenic mice. Further comparative sequence analysis of mouse and pufferfish intergenic regions revealed eight highly conserved elements within these cis-regulators. Therefore, we propose that, in vertebrate evolution, the Otx2 promoter acquires multiple, spatiotemporally specific cis-regulators in order to precisely control highly coordinated processes in head development.

Functional tests of enhancer conservation between distantly related species

Expression patterns of orthologous genes are often conserved, even between distantly related organisms, suggesting that once established, developmental programs can be stably maintained over long periods of evolutionary time. Because many orthologous transcription factors are also functionally conserved, one possible model to account for homologous gene expression patterns, is conservation of specific binding sites within cis-regulatory elements of orthologous genes. If this model is correct, a cis-regulatory element from one organism would be expected to function in a distantly related organism. To test this hypothesis, we fused the green fluorescent protein gene to neuronal and muscular enhancer elements from a variety of Drosophila melanogaster genes, and tested whether these would activate expression in the homologous cell types in Caenorhabditis elegans. Regulatory elements from several genes directed appropriate expression in homologous tissue types, suggesting conservation of regulatory sites. However, enhancers of most Drosophila genes tested were not properly recognized in C. elegans, implying that over this evolutionary distance enough changes occurred in cis-regulatory sequences and/or transcription factors to prevent proper recognition of heterospecific enhancers. Comparisons of enhancer elements of orthologous genes between C. elegans and C. briggsae revealed extensive conservation, as well as specific instances of functional divergence. Our results indicate that functional changes in cis-regulatory sequences accumulate on timescales much shorter than the divergence of arthropods and nematodes, and that mechanisms other than conservation of individual binding sites within enhancer elements are responsible for the conservation of expression patterns of homologous genes between distantly related species.

The regulation of retina specific expression of rhodopsin gene in vertebrates

Rhodopsin genes in most vertebrate species, with the exception of teleost fishes, contain introns. Despite differences in the gene structure, similar regulatory motifs have been identified in fish, amphibian and mammalian rhodopsin promoters, suggesting that rhodopsin gene regulation may be conserved in vertebrates. However, there is no direct evidence to support this notion. To address this, the rhodopsin promoter from the pufferfish, Fugu rubripes, was isolated and tested in transgenic mice and frogs. A 6.5 kb Fugu genomic fragment containing the rhodopsin gene and 4.5 kb 5' flanking region was able to direct expression of the Fugu rhodopsin gene to the retina of transgenic mice. In transgenic tadpoles, photoreceptor rod cell-specific expression of a reporter gene was achieved using only 500 bp Fugu rhodopsin promoter fragment. Mutagenesis of this promoter fragment revealed that a conserved NRE-like motif is crucial for the retina-specific expression. Our investigation suggests that the regulation of retinal specific expression is conserved in the pufferfish, frog and mouse and that the ancestral intron-containing rhodopsin gene has been displaced by an intronless copy in teleosts.

Evolutionary conservation of regulatory elements in vertebrate Hox gene clusters

Comparisons of DNA sequences among evolutionarily distantly related genomes permit identification of conserved functional regions in noncoding DNA. Hox genes are highly conserved in vertebrates, occur in clusters, and are uninterrupted by other genes. We aligned (PipMaker) the nucleotide sequences of the HoxA clusters of tilapia, pufferfish, striped bass, zebrafish, horn shark, human, and mouse, which are separated by approximately 500 million years of evolution. In support of our approach, several identified putative regulatory elements known to regulate the expression of Hox genes were recovered. The majority of the newly identified putative regulatory elements contain short fragments that are almost completely conserved and are identical to known binding sites for regulatory proteins (Transfac database). The regulatory intergenic regions located between the genes that are expressed most anteriorly in the embryo are longer and apparently more evolutionarily conserved than those at the other end of Hox clusters. Different presumed regulatory sequences are retained in either the Aalpha or Abeta duplicated Hox clusters in the fish lineages. This suggests that the conserved elements are involved in different gene regulatory networks and supports the duplication-deletion-complementation model of functional divergence of duplicated genes.

The human L-threonine 3-dehydrogenase gene is an expressed pseudogene

BACKGROUND

L-threonine is an indispensable amino acid. One of the major L-threonine degradation pathways is the conversion of L-threonine via 2-amino-3-ketobutyrate to glycine. L-threonine dehydrogenase (EC 1.1.1.103) is the first enzyme in the pathway and catalyses the reaction: L-threonine + NAD+ = 2-amino-3-ketobutyrate + NADH. The murine and porcine L-threonine dehydrogenase genes (TDH) have been identified previously, but the human gene has not been identified.

RESULTS

The human TDH gene is located at 8p23-22 and has 8 exons spanning 10 kb that would have been expected to encode a 369 residue ORF. However, 2 cDNA TDH transcripts encode truncated proteins of 157 and 230 residues. These truncated proteins are the result of 3 mutations within the gene. There is a SNP, A to G, present in the genomic DNA sequence of some individuals which results in the loss of the acceptor splice site preceding exon 4. The acceptor splice site preceding exon 6 was lost in all 23 individuals genotyped and there is an in-frame stop codon in exon 6 (CGA to TGA) resulting in arginine-214 being replaced by a stop codon. These truncated proteins would be non-functional since they have lost part of the NAD+ binding motif and the COOH terminal domain that is thought to be involved in binding L-threonine. TDH mRNA was present in all tissues examined.

CONCLUSIONS

The human L-threonine 3-dehydrogenase gene is an expressed pseudogene having lost the splice acceptor site preceding exon 6 and codon arginine-214 (CGA) is mutated to a stop codon (TGA).

Fugu and human sequence comparison identifies novel human genes and conserved non-coding sequences

The compact genome of the pufferfish, Fugu rubripes, has been proposed as a 'reference' genome to aid in annotating and analysing the human genome. We have annotated and compared 85 kb of Fugu sequence containing 17 genes with its homologous loci in the human draft genome and identified three 'novel' human genes that were missed or incompletely predicted by the previous gene prediction methods. Two of the novel genes contain zinc finger domains and are designated ZNF366 and ZNF367. They map to human chromosomes 5q13.2 and 9q22.32, respectively. The third novel gene, designated C9orf21, maps to chromosome 9q22.32. This gene is unique to vertebrates, and the protein encoded by it does not contain any known domains. We could not find human homologs for two Fugu genes, a novel chemokine gene and a kinase gene. These genes are either specific to teleosts or lost in the human lineage. The Fugu-human comparison identified several conserved non-coding sequences in the promoter and intronic regions. These sequences, conserved during 450 million years of vertebrate evolution, are likely to be involved in gene regulation. The 85 kb Fugu locus is dispersed over four human loci, occupying about 1.5 Mb. Contiguity is conserved in the human genome between six out of 16 Fugu gene pairs. These contiguous chromosomal segments should share a common evolutionary history dating back to the common ancestor of mammals and teleosts. We propose contiguity as strong evidence to identify orthologous genes in distant organisms. This study confirms the utility of the Fugu as a supplementary tool to uncover and confirm novel genes and putative gene regulatory regions in the human genome.

Sushi gets serious: the draft genome sequence of the pufferfish Fugu rubripes

The publication of the Fugu rubripes draft genome sequence will take this fish from culinary delicacy to potent tool in deciphering the mysteries of human genome function.