Floating Search Methodology for Combining Classification Models for Site Recognition in DNA Sequences.

Recognition of the functional sites of genes, such as translation initiation sites, donor and acceptor splice sites and stop codons, is a relevant part of many current problems in bioinformatics. The best approaches use sophisticated classifiers, such as support vector machines. However, with the rapid accumulation of sequence data, methods for combining many sources of evidence are necessary as it is unlikely that a single classifier can solve this problem with the best possible performance. A major issue is that the number of possible models to combine is large and the use of all of these models is impractical. In this paper we present a methodology for combining many sources of information to recognize any functional site using "floating search", a powerful heuristics applicable when the cost of evaluating each solution is high. We present experiments on four functional sites in the human genome, which is used as the target genome, and use another 20 species as sources of evidence. The proposed methodology shows significant improvement over state-of-the-art methods. The results show an advantage of the proposed method and also challenge the standard assumption of using only genomes not very close and not very far from the human to improve the recognition of functional sites.

Improved Pathogenic Variant Localization via a Hierarchical Model of Sub-regional Intolerance.

Different parts of a gene can be of differential importance to development and health. This regional heterogeneity is also apparent in the distribution of disease-associated mutations, which often cluster in particular regions of disease-associated genes. The ability to precisely estimate functionally important sub-regions of genes will be key in correctly deciphering relationships between genetic variation and disease. Previous methods have had some success using standing human variation to characterize this variability in importance by measuring sub-regional intolerance, i.e., the depletion in functional variation from expectation within a given region of a gene. However, the ability to precisely estimate local intolerance was restricted by the fact that only information within a given sub-region is used, leading to instability in local estimates, especially for small regions. We show that borrowing information across regions using a Bayesian hierarchical model stabilizes estimates, leading to lower variability and improved predictive utility. Specifically, our approach more effectively identifies regions enriched for ClinVar pathogenic variants. We also identify significant correlations between sub-region intolerance and the distribution of pathogenic variation in disease-associated genes, with AUCs for classifying de novo missense variants in Online Mendelian Inheritance in Man (OMIM) genes of up to 0.86 using exonic sub-regions and 0.91 using sub-regions defined by protein domains. This result immediately suggests that considering the intolerance of regions in which variants are found may improve diagnostic interpretation. We also illustrate the utility of integrating regional intolerance into gene-level disease association tests with a study of known disease-associated genes for epileptic encephalopathy.

Clinical utility of next-generation sequencing in the diagnosis of hereditary haemolytic anaemias.

Hereditary haemolytic anaemias are genetically and phenotypically heterogeneous disorders characterized by increased red cell destruction, with consequences ranging from innocuous to severe life-threatening anaemia. Diagnostic laboratories endeavour to assist clinicians reach the exact patient diagnosis by using tests principally based on morphological and biochemical techniques. However, these routine studies may be inconclusive, particularly in newborn infants and when transfusions have recently been administered. Large numbers and size of the potentially involved genes also impose a practical challenge for molecular diagnosis using routine sequencing approaches. To overcome these diagnostic shortcomings, we have utilized next-generation sequencing to provide a high-throughput, highly sensitive assay. We developed a panel interrogating 28 genes encoding cytoskeletal proteins and enzymes with sequencing coverage of the coding regions, splice site junctions, deep intronic and regulatory regions. We then evaluated 19 samples, including infants with unexplained extreme hyperbilirubinaemia and patients with transfusion-dependent haemolytic anaemia. Where possible, inheritance patterns of pathogenic mutations were determined by sequencing of immediate relatives. We conclude that this next-generation sequencing panel could be a cost-effective approach to molecular diagnosis of hereditary haemolytic anaemia, especially when the family history is uninformative or when routine laboratory testing fails to identify the causative haemolytic process.

Modularization of genetic elements promotes synthetic metabolic engineering.

In the context of emerging synthetic biology, metabolic engineering is moving to the next stage powered by new technologies. Systematical modularization of genetic elements makes it more convenient to engineer biological systems for chemical production or other desired purposes. In the past few years, progresses were made in engineering metabolic pathway using synthetic biology tools. Here, we spotlighted the topic of implementation of modularized genetic elements in metabolic engineering. First, we overviewed the principle developed for modularizing genetic elements and then discussed how the genetic modules advanced metabolic engineering studies. Next, we picked up some milestones of engineered metabolic pathway achieved in the past few years. Last, we discussed the rapid raised synthetic biology field of "building a genome" and the potential in metabolic engineering.

The first echinoderm gamma-interferon-inducible lysosomal thiol reductase (GILT) identified from sea cucumber (Stichopus monotuberculatus).

Gamma-interferon-inducible lysosomal thiol reductase (GILT) has been described as a key enzyme that facilitating the processing and presentation of major histocompatibility complex class II-restricted antigen in mammals. In this study, the first echinoderm GILT named StmGILT was identified from sea cucumber (Stichopus monotuberculatus). The StmGILT cDNA is 1529 bp in length, containing a 5'-untranslated region (UTR) of 87 bp, a 3'-UTR of 674 bp and an open reading frame (ORF) of 768 bp that encoding a protein of 255 amino acids with a deduced molecular weight of 27.82 kDa and a predicted isoelectric point of 4.73. The putative StmGILT protein possesses all the main characteristics of known GILT proteins, including a signature sequence, a reductase active site CXXC, twelve conserved cysteines, and two potential N-linked glycosylation sites. For the gene structure, StmGILT contains four exons separated by three introns. In the promoter region of StmGILT gene, an NF-κB binding site and an IFN-γ activation site were found. The thiol reductase activity of recombinant StmGILT protein was also demonstrated in this study. In addition, the highest level of mRNA expression was noticed in coelomocytes of S. monotuberculatus. In in vitro experiments performed in coelomocytes, the expression of StmGILT mRNA was significantly up-regulated by lipopolysaccharides (LPS), inactivated bacteria or polyriboinosinic polyribocytidylic acid [poly (I:C)] challenge, suggested that the sea cucumber GILT might play critical roles in the innate immune defending against bacterial and viral infections.

Two NF-κB inhibitor-alpha (IκBα) genes from rock bream (Oplegnathus fasciatus): molecular characterization, genomic organization and mRNA expression analysis after immune stimulation.

IkBa is a member of IkB family, which sequesters NF-kB in an inactivate form in the cytoplasm and blocks the translocation of NF-kB to nucleus. The IkBa paralogs of rock bream (OfIkBa-A and OfIkBa-B) encoded IkBa proteins with typical features including, highly conserved IkB degradation motif, six ankyrin repeats and a PEST sequence. However, their amino acid identity and similarity were only 55.6 and 69.7%, respectively suggesting that these two genes could be the two different isoforms of IkBa. The number and size of the exons of OfIkBa-A and OfIkBa-B were conserved well with all the compared vertebrate species, although they have significantly different genomic sizes. Phylogenetic analysis revealed that OfIkBa-A and OfIkBa-B proteins cluster with IkBa family members; however, they were grouped with different subclades in IkBa family. Tissue specific expression of OfIkBa mRNA was constitutively detected in all the tested tissues, and they showed the higher transcription level in heart, liver, gill and peripheral blood cells, respectively. The injection of flagellin stimulated the mRNA expression of OfIkBa paralogs in head kidney and intestine. Moreover, the OfIkBa mRNA expression in gill and liver was significantly upregulated by LPS, poly I:C and Edwardsiella tarda challenges. The transcription of OfIkBa was up-regulated in early-phase of injection and then rapidly restored. These results suggest that the OfIkBa paralogs might be involved in rapid immune responsive reactions in rock bream against bacterial and viral pathogens.

A universal small molecule, inorganic phosphate, restricts the substrate specificity of Dicer-2 in small RNA biogenesis.

The enzyme Dicer is central to the production of small silencing RNAs such as microRNAs (miRNAs) and small interfering RNAs (siRNAs). Like other insects, Drosophila melanogaster uses different Dicers to make siRNAs and miRNAs: Dicer-1 produces miRNAs from pre-miRNAs, whereas Dicer-2 generates siRNAs from long double-stranded RNA (dsRNA). How do the 2 Dicers achieve their substrate specificity? Here, we review recent findings that inorganic phosphate restricts the substrate specificity of Dicer-2 to long dsRNA. Inorganic phosphate inhibits Dicer-2 from binding and cleaving pre-miRNAs, without affecting the processing of long dsRNA. Crystal structures of a fragment of human Dicer in complex with an RNA duplex identify a phosphate-binding pocket that recognizes both the 5'-monophosphate of a substrate RNA and inorganic phosphate. We propose that inorganic phosphate occupies the phosphate-binding pocket in the fly Dicer-2, blocking binding of pre-miRNA and restricting pre-miRNA processing to Dicer-1. Thus, a small molecule can alter the substrate specificity of a nucleic acid-processing enzyme.

Genomic organization and functional diversification of two warm-temperature-acclimation-associated 65-kDa protein genes in rockbream (Oplegnathus fasciatus; Perciformes).

Two paralogue genes of warm-temperature-acclimation-associated 65-kDa protein were characterized and their mRNA expression patterns during various experimental stimulations were examined in the rockbream (Oplegnathus fasciatus; Perciformes). Rockbream Wap65 isoforms (rbWap65-1 and rbWap65-2) share basically common structural features with other teleostean orthologues and human hemopexin (HPX) at both amino acid (conserved cysteine and histidine residues) and genomic levels (ten-exon structure), although the rbWap65-2 reveals more homologous characteristics to human HPX than does rbWap65-1 isoform. Southern blot analysis indicates that each rbWap65 isoform exists as a single copy gene in the rockbream genome. Both rbWap65 genes were predicted to possess various transcription factor (TF) binding motifs related with stress and innate immunity in their 5ʹ-upstream regions, in which inflammation-related motifs were more highlighted in the rbWap65-2 than in rbWap65-1. Based on the RT-PCR assay, the liver-predominant expression pattern was more apparent in rbWap65-1 than rbWap65-2 isoform. During thermal elevation, clear upregulation was found only for the rbWap65-1. In contrast, immune stimulations (bacterial challenges, viral infection and iron overload) activated more preferentially the rbWap65-2 isoform in overall, although the inducibility was affected by the kinds of stimulators and tissue types. Taken together, our data suggest that the two paralogue rbWap65 isoforms have experienced subfunctionalization and/or neofunctionalization during their evolutionary history, in which the rbWap65-2 has retained closer, functional orthology to the human HPX while the rbWap65-1 have been diversified to be more related with thermal acclimation physiology.

Structural and functional characterization of the Senegalese sole (Solea senegalensis) Mx promoter.

Mx proteins are one of the most studied interferon-stimulated genes (ISGs). The antiviral activity against different fish viruses has been demonstrated for diverse fish Mx proteins, including the Senegalese sole (Solea senegalensis) Mx protein (SsMx). The aim of the current study is to characterize the structure and functional activity of the SsMx promoter. Several polyclonal cell populations expressing the luciferase reporter gene under the control of the SsMx promoter have been used to determine the ability of this promoter to drive the expression of the luciferase gene after poly I:C stimulation. In addition, the implication of each interferon-stimulated response element (ISRE) in the activation of the promoter has also been analysed. The genomic structure of the Senegalese sole and Japanese flounder Mx promoters (containing three ISREs) differs from the rest of the fish Mx promoters described to date. The ISRE1, the one closest to the start codon, is the main ISRE involved in the SsMx promoter activity, whereas ISRE2 and ISRE3 show a minor additive effect on this activity. Another feature differing SsMx promoter from the rest of the fish Mx promoters is the presence of a 24-bp GC island close to the ATG codon, including one Sp1 binding site, which may constitute the transcriptional start site. Furthermore, the SsMx promoter contains a gamma interferon activation site (GAS) element.

An evolutionarily ancient NO synthase (NOS) in shrimp.

Nitric oxide (NO) is a well known essential molecule that is involved in multiple functions such as neuron transduction, cardiac disease, immune responses, etc.; nitric oxide synthase (NOS) is a critical enzyme that catalyzes the synthesis of it. A very few crustacean NOS molecules were biochemically characterized so far. In the present study, we cloned and characterized a NOS cDNA from haemocytes of tiger shrimp (Penaeus monodon) (PmNOS). The full-length of PmNOS cDNA contained 3997 bp, including a 5'UTR of 249 bp, ORF of 3582 bp and a 3'UTR of 166 bp. The putative peptide was 1193 amino acid residues in length, with an estimated molecular weight of 134.7 kDa and pI 6.7. Structurally, PmNOS contained oxygenase and reductase domains at N-terminal and C-terminal, respectively, and connected with a calmodulin binding motif. The deduced amino acid sequence of PmNOS shared 98% identical to the Chinese shrimp (Fenneropenaeus chinensis) NOS. Phylogenetically, PmNOS clustered with invertebrate NOS, but not clustered with iNOS, eNOS or nNOS found in vertebrates. PmNOS mRNA was expressed in many tissues or organs including thoracic and ventral nerves, midgut, gill, eyestalk, haemocytes, subcuticular epithelium and heart, but not found in hepatopancreas, muscle and lymphoid organ. But there was no significant difference in PmNOS mRNA expression after stimulation with LPS either by different concentration or time course or against CpG-ODN 2006. The enzyme activities of rPmNOS or crude homogenates from different tissues were detected, and were shown its highest activity in thoracic and ventral nerves, moderate in midgut and haemocytes but the lowest activity were seen in muscle. The addition of NOS antibody against NADPH binding domain leads to less activity which suggested that NADPH was an essential cofactor for PmNOS catalytic activity. The calcium dependency of PmNOS was ascertained using calmodulin inhibitor, Trifluroperazine. To confirm the population of haemocyte which produce NOS, the florescence test was assayed, and it implicated that the production of NO was catalyzed by subset of granulocytic NOS. Since the MW range, inducible/noninducible transcript, calcium-dependent activity and tissue distribution, we suggest that PmNOS may recognize as an ancient NOS evolutionarily.

Evolution of the tbx6/16 subfamily genes in vertebrates: insights from zebrafish.

In any comparative studies striving to understand the similarities and differences of the living organisms at the molecular genetic level, the crucial first step is to establish the homology (orthology and paralogy) of genes between different organisms. Determination of the homology of genes becomes complicated when the genes have undergone a rapid divergence in sequence or when the involved genes are members of a gene family that has experienced a differential gain or loss of its constituents in different taxonomic groups. Organisms with duplicated genomes such as teleost fishes might have been especially prone to these problems because the functional redundancies provided by the duplicate copies of genes would have allowed a rapid divergence or loss of genes during evolution. In this study, we will demonstrate that much of the ambiguities in the determination of the homology between fish and tetrapod genes resulting from the problems like these can be eliminated by complementing the sequence-based phylogenies with nonsequence information, such as the exon-intron structure of a gene or the composition of a gene's genomic neighbors. We will use the Tbx6/16 subfamily genes of zebrafish (tbx6, tbx16, tbx24, and mga genes), which have been well known for the ambiguity of their evolutionary relationships to the Tbx6/16 subfamily genes of tetrapods, as an illustrative example. We will show that, despite the similarity of sequence and expression to the tetrapod Tbx6 genes, zebrafish tbx6 gene is actually a novel T-box gene more closely related to the tetrapod Tbx16 genes, whereas the zebrafish tbx24 gene, hitherto considered to be a novel gene due to the high level of sequence divergence, is actually an ortholog of tetrapod Tbx6 genes. We will also show that, after their initial appearance by the multiplication of a common ancestral gene at the beginning of vertebrate evolution, the Tbx6/16 subfamily of vertebrate T-box genes might have experienced differential losses of member genes in different vertebrate groups and gradual pooling of member gene's functions in surviving members, which might have prevented the revelation of the true identity of member genes by way of the comparison of sequence and function.

Gene expression levels are correlated with synonymous codon usage, amino acid composition, and gene architecture in the red flour beetle, Tribolium castaneum.

Gene expression levels correlate with multiple aspects of gene sequence and gene structure in phylogenetically diverse taxa, suggesting an important role of gene expression levels in the evolution of protein-coding genes. Here we present results of a genome-wide study of the influence of gene expression on synonymous codon usage, amino acid composition, and gene structure in the red flour beetle, Tribolium castaneum. Consistent with the action of translational selection, we find that synonymous codon usage bias increases with gene expression. However, the correspondence between tRNA gene copy number and optimal codons is weak. At the amino acid level, translational selection is suggested by the positive correlation between tRNA gene numbers and amino acid usage, which is stronger for highly expressed genes. In addition, there is a clear trend for increased use of metabolically cheaper, less complex amino acids as gene expression increases. tRNA gene numbers also correlate negatively with amino acid size/complexity (S/C) score indicating the coupling between translational selection and selection to minimize the use of large/complex amino acids. Interestingly, the analysis of 10 additional genomes suggests that the correlation between tRNA gene numbers and amino acid S/C score is widespread and might be explained by selection against negative consequences of protein misfolding. At the level of gene structure, three major trends are detected: 1) complete coding region length increases across low and intermediate expression levels but decreases in highly expressed genes; 2) the average intron size shows the opposite trend, first decreasing with expression, followed by a slight increase in highly expressed genes; and 3) intron density remains nearly constant across all expression levels. These changes in gene architecture are only in partial agreement with selection favoring reduced cost of biosynthesis.

Characterization and developmental expression of genes encoding the early carotenoid biosynthetic enzymes in Citrus paradisi Macf.

In the present study, the full-length cDNA sequences of PSY, PDS, and ZDS, encoding the early carotenoid biosynthetic enzymes in the carotenoid pathway of grapefruit (Citrus paradisi), were isolated and characterized for the first time. CpPSY contained a 1311-bp open reading frame (ORF) encoding a polypeptide of 436 amino acids, CpPDS contained a 1659-bp ORF encoding a polypeptide of 552 amino acids, and CpZDS contained a 1713-bp ORF encoding a polypeptide of 570 amino acids. Phylogenetic analysis indicated that CpPSY shares homology with PSYs from Citrus, tomato, pepper, Arabidopsis, and the monocot PSY1 group, while CpPDS and CpZDS are most closely related to orthologs from Citrus and tomato. Expression analysis revealed fluctuations in CpPSY, CpPDS, and CpZDS transcript abundance and a non-coordinated regulation between the former and the two latter genes during fruit development in albedo and juice vesicles of white ('Duncan') and red ('Flame') grapefruits. A 3× higher upregulation of CpPSY expression in juice vesicles of red-fleshed 'Flame' as compared to white-fruited 'Duncan' was observed in the middle stages of fruit development, which correlates with the well documented accumulation pattern of lycopene in red grapefruit. Together with previous data, our results suggest that the primary mechanism controlling lycopene accumulation in red grapefruit involves the transcriptional upregulation of CpPSY, which controls the flux into the carotenoid pathway, and the downregulated expression of CpLCYB2, which controls the step of cyclization of lycopene in chromoplasts during fruit ripening. A correlation between CpPSY expression and fruit color evolution in red grapefruit is demonstrated.

Zebrafish type XVII collagen: gene structures, expression profiles, and morpholino "knock-down" phenotypes.

The human COL17A1 gene encodes type XVII collagen (also known as the 180-kDa bullous pemphigoid antigen), an integral component of hemidesmosomes, attachment complexes providing integrity to the dermal-epidermal junction. Zebrafish, a useful model system to study skin development, displays fully developed hemidesmosomes at approximately 5 days post-fertilization (dpf). We have identified two COL17A1 orthologues in the zebrafish genome, col17a1a and col17a1b, which are expressed in the skin and the neural system, respectively. The proteins coded by these genes have structural module organizations homologous to the human type XVII collagen. "Knock-down" of the expression of col17a1a with a specific morpholino targeting the 5' UTR of the gene resulted in a blistering phenotype and in perturbations in the basement membrane zone. "Knock-down" of col17a1b expression resulted in ablation or in marked reduction of neuromasts in the lateral line. Thus, zebrafish has two COL17A1 orthologues which may have evolved tissue-specific functions during vertebrate development. Collectively, zebrafish provides a model system to study the molecular aspects of skin development and offers insights into the corresponding human diseases.

Aspergillus fumigatus survival in alkaline and extreme zinc-limiting environments relies on the induction of a zinc homeostasis system encoded by the zrfC and aspf2 genes.

Aspergillus fumigatus has three zinc transporter-encoding genes whose expression is regulated by both pH and the environmental concentration of zinc. We have previously reported that the zrfA and zrfB genes of A. fumigatus are transcribed at higher levels and are required for fungal growth under acidic zinc-limiting conditions whereas they are dispensable for growth in neutral or alkaline zinc-limiting media. Here we report that the transporter of the zinc uptake system that functions in A. fumigatus growing in neutral or alkaline environments is encoded by zrfC. The transcription of zrfC occurs divergently with respect to the adjacent aspf2 gene, which encodes an immunodominant antigen secreted by A. fumigatus. The two genes-zrfC and aspf2-are required to different extents for fungal growth in alkaline and extreme zinc-limiting media. Indeed, these environmental conditions induce the simultaneous transcription of both genes mediated by the transcriptional regulators ZafA and PacC. ZafA upregulates the expression of zrfC and aspf2 under zinc-limiting conditions regardless of the ambient pH, whereas PacC represses the expression of these genes under acidic growth conditions. Interestingly, the mode of action of PacC for zrfC-aspf2 transcription contrasts with the more widely accepted model for PacC function, according to which under alkaline growth conditions PacC would activate the transcription of alkaline-expressed genes but would repress the transcription of acid-expressed genes. In sum, this report provides a good framework for investigating several important aspects of the biology of species of Aspergillus, including the repression of alkaline genes by PacC at acidic pH and the interrelationship that must exist between tissue pH, metal availability in the host tissue, and fungal virulence.

Association of HPA axis genes with suicidal behaviour in schizophrenia.

Family, adoption and twin studies show that genetics influences suicidal behaviour, but do not indicate specific susceptibility variants. Stress response is thought to be mediated by the corticotrophin-releasing hormone (CRH), which is known to be a regulator of the hypothalamic-pituitary-adrenal pathway (HPA). Alterations in HPA system have been related to impulsivity, aggression and suicidal behaviour, common feature in schizophrenia. CRH is the hypothalamic factor that stimulates the pituitary gland. To search for markers conferring genetic susceptibility to suicide, we typed six HPA axis genes (CRH, CRHR1, CRHR2, CRHBP, MC2R, NC3R1) in a cohort of 231 subjects with schizophrenia in which 81 attempted suicide. The genotype analyses yielded significant association between CRH binding protein (CRHBP) and suicide attempt (P = 0.035). The genotype analysis for quantitative measures of suicidal behaviour showed no association. The interaction analysis showed a significant interaction between CRH receptor type 1 (CRHR1) and CRH binding protein (CRHBP) in influencing suicide attempt and the severity of suicidal behaviour. Current results show that genetic variation in HPA axis genes could be associated with suicidal behaviour in schizophrenia. This is to our knowledge the first study on suicidal behaviour investigating the interaction among the HPA axis genes.

Calcitonin receptor-stimulating peptide: Its evolutionary and functional relationship with calcitonin/calcitonin gene-related peptide based on gene structure.

This review focuses on the evolutionary and functional relationship of calcitonin receptor-stimulating peptide (CRSP) with calcitonin (CT)/calcitonin gene-related peptide (CGRP) in mammals. CRSP shows high sequence identity with CGRP, but distinct biological properties. CRSP genes (CRSPs) have been identified in mammals such as pigs and dogs of the Laurasiatheria, but not in primates and rodents of the Euarchontoglires or in non-placental mammals. CRSPs have genomic organizations highly similar to those of CT/CGRP genes (CT/CGRPs), which are located along with CGRPs in a locus between CYP2R1 and INSC, while the other members of the CGRP superfamily, adrenomedullin and amylin, show genomic organizations and locations distinct from CT, CGRP, and CRSP. Thus, we categorized these three peptides into the CT/CGRP/CRSP family. Non-placental mammals having one and placental mammals having multiple CT/CGRP/CRSP family genes suggests that multiplicity of CT/CGRP started at an early stage of mammalian evolution. In the placental mammals, Laurasiatheria generally possesses multiple CRSPs and only one CT/CGRP, while Euarchontoglires possesses CT/CGRP and CGRPbeta but no CRSP, indicating an increase in the diversity and multiplicity of this family of genes in mammalian evolution. Phylogenetic analysis suggests that some CRSPs have been generated very recently in mammalian evolution. Taken together, the increase in the number and complexity of the CT/CGRP/CRSP family genes may have due to evolutionary pressure to facilitate adaptation during mammalian evolution. In this regard, it is important to elucidate the physiological roles of CT, CGRP and CRSP from the viewpoint of the CT/CGRP/CRSP family even in Euarchontoglires.

Expression of the gene for ancestral glycoprotein hormone beta subunit in the nerve cord of amphioxus.

Amphioxus belongs to the subphylum cephalochordata, a clade of chordates phylogenetically placed at the most basal position. Despite many studies on the endocrine system of amphioxus, there were no confident lines of evidence on the presence of pituitary hormones, whereas recent amphioxus genome analysis reported that amphioxus has no pituitary hormone except for thyrostimulin, which is a glycoprotein hormone in the pituitary, brain, and other organs of vertebrates. In the present study, we cloned cDNA for one glycoprotein hormone beta subunit (GPB) from amphioxus, AmpGPB5, and phylogenetically indicated that AmpGPB5 is the ancestral molecule of glycoprotein hormone beta subunits of vertebrates including pituitary glycoprotein hormones. Synteny analyses showed conservation of chromosomal location of genes near GPB genes from amphioxus through human. The AmpGPB5 gene was expressed in a restricted region of the dorsal part of the nerve cord, glandular atrial cells of gills, and pre-vitellogenic oocytes in amphioxus. However, expression was not detected in the Hatschek's pit which is considered to be a primitive pituitary gland. On the basis of present results, we hypothesize that a portion of vertebrate pituitary hormones might be derived from an ancestral glycoprotein hormone of amphioxus that functions as a neuroendocrine hormone.

An analysis of the human chemokine CXC receptor 4 gene.

In this article we analyze some of the structural characteristics of the coding section and the intron of the human chemokine CXC receptor 4 (a 7-transmembrane receptor) pre-mRNA. In the coding sequence the frequencies of the individual nucleotides do not depart significantly from 0.25, while in the intron the frequencies of the As and Gs are significantly lower and higher, respectively, than expected from a random distribution. Analysis of the pattern of association of nucleotides into triplets or couples shows that some triplets or couples occur with frequencies significantly higher or lower than expected when assuming a random association of nucleotides. In particular, in the intron combinations of the same nucleotide are over-represented. 7-or-more nucleotide repeats occur in both the coding section and the intron with frequencies which exceed the confidence limits for a random distribution. For the coding sequence this is possibly explained by the alternans of relatively similar hydrophobic-coding sections and relatively similar intervening intracellular and extracellular hydrophilic-coding sections. 7-or-more nucleotide repeats in reverse order and in reverse/complemented order occur in the intron, but not in the coding section, with frequencies which significantly exceed a random distribution. The numerous intronic repeats in reverse/complemented order may be of relevance for the secondary structure of the intron and might be one important element of the integrated splicing code.

Two macrophage colony-stimulating factor genes exist in fish that differ in gene organization and are differentially expressed.

Mammalian macrophage CSF (MCSF; CSF-1) is the primary regulator of the mononuclear phagocyte lineage. We, for the first time, report the complete sequencing of five MCSF cDNAs from three fish species, rainbow trout, zebrafish, and goldfish. Despite the difference in the lengths of the MCSF transcripts, all of the fish MCSF molecules encode a signal peptide, a CSF-1 domain, a transmembrane domain, and an intracellular region. Each fish MCSF gene has a unique exon/intron structure. The primordial MCSF gene may have had a nine exon/eight intron structure. In this model, insertion of an intron in exon 6 in primitive fish created the fish type I MCSF, while the loss of this exon or part of the original exon 6 created the fish type II MCSF. Investigation of alternative splicing variants in trout suggests that no mammalian equivalent splice variants exist. The two trout MCSF genes are differentially expressed in vivo and contributed differently to the high-level expression of MCSF in spleen and head kidney. In contrast to the up-regulation of MCSF by PMA in mammals, in trout MCSF1 expression is down-regulated by PMA treatment. As in mammals, recombinant trout MCSF1 can promote the growth of head kidney leukocytes, and it up-regulates the expression of CXCR3 in head kidney macrophages, with the latter suggesting a role of MCSF in the trafficking of macrophages to sites of inflammation or injury where the CXCR3 ligands are expressed. Thus MCSF has an important role in the immune system of fish as in mammals.