Sex-specific regulation of neuronal functions in Caenorhabditis elegans: the sex-determining protein TRA-1 represses goa-1/Gα(i/o).

Females and males differ substantially in various neuronal functions in divergent, sexually dimorphic animal species, including humans. Despite its developmental, physiological and medical significance, understanding the molecular mechanisms by which sex-specific differences in the anatomy and operation of the nervous system are established remains a fundamental problem in biology. Here, we show that in Caenorhabditis elegans (nematodes), the global sex-determining factor TRA-1 regulates food leaving (mate searching), male mating and adaptation to odorants in a sex-specific manner by repressing the expression of goa-1 gene, which encodes the Gα(i/o) subunit of heterotrimeric G (guanine-nucleotide binding) proteins triggering physiological responses elicited by diverse neurotransmitters and sensory stimuli. Mutations in tra-1 and goa-1 decouple behavioural patterns from the number of X chromosomes. TRA-1 binds to a conserved binding site located in the goa-1 coding region, and downregulates goa-1 expression in hermaphrodites, particularly during embryogenesis when neuronal development largely occurs. These data suggest that the sex-determination machinery is an important modulator of heterotrimeric G protein-mediated signalling and thereby various neuronal functions in this organism and perhaps in other animal phyla.

Characterization of all possible single-nucleotide change caused amino acid substitutions in the kinase domain of Bruton tyrosine kinase.

Knowledge about features distinguishing deleterious and neutral variations is crucial for interpretation of novel variants. Bruton tyrosine kinase (BTK) contains the highest number of unique disease-causing variations among the human protein kinases, still it is just 10% of all the possible single-nucleotide substitution-caused amino acid variations (SNAVs). In the BTK kinase domain (BTK-KD) can appear altogether 1,495 SNAVs. We investigated them all with bioinformatic and protein structure analysis methods. Most disease-causing variations affect conserved and buried residues disturbing protein stability. Minority of exposed residues is conserved, but strongly tied to pathogenicity. Sixty-seven percent of variations are predicted to be harmful. In 39% of the residues, all the variants are likely harmful, whereas in 10% of sites, all the substitutions are tolerated. Results indicate the importance of the entire kinase domain, involvement in numerous interactions, and intricate functional regulation by conformational change. These results can be extended to other protein kinases and organisms.

Abnormal cerebellar development and ataxia in CARP VIII morphant zebrafish.

Congenital ataxia and mental retardation are mainly caused by variations in the genes that affect brain development. Recent reports have shown that mutations in the CA8 gene are associated with mental retardation and ataxia in humans and ataxia in mice. The gene product, carbonic anhydrase-related protein VIII (CARP VIII), is predominantly present in cerebellar Purkinje cells, where it interacts with the inositol 1,4,5-trisphosphate receptor type 1, a calcium channel. In this study, we investigated the effects of the loss of function of CARP VIII during embryonic development in zebrafish using antisense morpholino oligonucleotides against the CA8 gene. Knockdown of CA8 in zebrafish larvae resulted in a curved body axis, pericardial edema and abnormal movement patterns. Histologic examination revealed gross morphologic defects in the cerebellar region and in the muscle. Electron microscopy studies showed increased neuronal cell death in developing larvae injected with CA8 antisense morpholinos. These data suggest a pivotal role for CARP VIII during embryonic development. Furthermore, suppression of CA8 expression leads to defects in motor and coordination functions, mimicking the ataxic human phenotype. This work reveals an evolutionarily conserved function of CARP VIII in brain development and introduces a novel zebrafish model in which to investigate the mechanisms of CARP VIII-related ataxia and mental retardation in humans.

Carbonic anhydrase related proteins: molecular biology and evolution.

The catalytically inactive isoforms of α-carbonic anhydrases are known as carbonic anhydrase related proteins (CARPs). The CARPs occur independently or as domains of other proteins in animals (both vertebrates and invertebrates) and viruses. The catalytic inactivity of CARPs is due to the lack of histidine residues required for the coordination of the zinc atom. The phylogenetic analysis shows that these proteins are highly conserved across the species. The three CARPs in vertebrates are known as CARP VIII, X and XI. CARPs orthologous to CARP VIII are found in deuterostome invertebrates, whereas protostomes only possess orthologs of CARP X. The CA-like domains of receptor-type protein tyrosine phosphatases (PTPR) are found only in PTPRG and PTPRZ. Most of these CARPs are predominantly expressed in central nervous system. Among the three vertebrate CA isoforms, CARP VIII is functionally associated with motor coordination in human, mouse and zebrafish and certain types of cancers in humans. Vertebrate expression studies show that CARP X is exclusively expressed in the brain. CARP XI is only found in tetrapods and is highly expressed in the central nervous system (CNS) of humans and mice and is also associated with several cancers. CARP VIII, PTPRZ and PTPRG have been shown to coordinate the function of other proteins by protein-protein interaction, and viral CARPs participate in attachment to host cells, but the precise biological function of CARPs X and XI is still unknown. The findings so far suggest many novel functions for the CARP subfamily, most likely related to binding to other proteins.

Co-acting gene networks predict TRAIL responsiveness of tumour cells with high accuracy.

BACKGROUND: Identification of differentially expressed genes from transcriptomic studies is one of the most common mechanisms to identify tumor biomarkers. This approach however is not well suited to identify interaction between genes whose protein products potentially influence each other, which limits its power to identify molecular wiring of tumour cells dictating response to a drug. Due to the fact that signal transduction pathways are not linear and highly interlinked, the biological response they drive may be better described by the relative amount of their components and their functional relationships than by their individual, absolute expression. RESULTS: Gene expression microarray data for 109 tumor cell lines with known sensitivity to the death ligand cytokine tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) was used to identify genes with potential functional relationships determining responsiveness to TRAIL-induced apoptosis. The machine learning technique Random Forest in the statistical environment "R" with backward elimination was used to identify the key predictors of TRAIL sensitivity and differentially expressed genes were identified using the software GeneSpring. Gene co-regulation and statistical interaction was assessed with q-order partial correlation analysis and non-rejection rate. Biological (functional) interactions amongst the co-acting genes were studied with Ingenuity network analysis. Prediction accuracy was assessed by calculating the area under the receiver operator curve using an independent dataset. We show that the gene panel identified could predict TRAIL-sensitivity with a very high degree of sensitivity and specificity (AUC=0·84). The genes in the panel are co-regulated and at least 40% of them functionally interact in signal transduction pathways that regulate cell death and cell survival, cellular differentiation and morphogenesis. Importantly, only 12% of the TRAIL-predictor genes were differentially expressed highlighting the importance of functional interactions in predicting the biological response. CONCLUSIONS: The advantage of co-acting gene clusters is that this analysis does not depend on differential expression and is able to incorporate direct- and indirect gene interactions as well as tissue- and cell-specific characteristics. This approach (1) identified a descriptor of TRAIL sensitivity which performs significantly better as a predictor of TRAIL sensitivity than any previously reported gene signatures, (2) identified potential novel regulators of TRAIL-responsiveness and (3) provided a systematic view highlighting fundamental differences between the molecular wiring of sensitive and resistant cell types.

Are proposed early genetic codes capable of encoding viable proteins?

Proteins are elaborate biopolymers balancing between contradicting intrinsic propensities to fold, aggregate, or remain disordered. Assessing their primary structural preferences observable without evolutionary optimization has been reinforced by the recent identification of de novo proteins that have emerged from previously non-coding sequences. In this paper we investigate structural preferences of hypothetical proteins translated from random DNA segments using the standard genetic code and three of its proposed evolutionarily predecessor models encoding 10, 6, and 4 amino acids, respectively. Our only main assumption is that the disorder, aggregation, and transmembrane helix predictions used are able to reflect the differences in the trends of the protein sets investigated. We found that the 10-residue code encodes proteins that resemble modern proteins in their predicted structural properties. All of the investigated early genetic codes give rise to proteins with enhanced disorder and diminished aggregation propensities. Our results suggest that an ancestral genetic code similar to the proposed 10-residue one is capable of encoding functionally diverse proteins but these might have existed under conditions different from today's common physiological ones. The existence of a protein functional repertoire for the investigated earlier stages which is quite distinct as it is today can be deduced from the presented results.

Analysis of evolution of carbonic anhydrases IV and XV reveals a rich history of gene duplications and a new group of isozymes.

Carbonic anhydrase (CA) isozymes CA IV and CA XV are anchored on the extracellular cell surface via glycosylphosphatidylinositol (GPI) linkage. Analysis of evolution of these isozymes in vertebrates reveals an additional group of GPI-linked CAs, CA XVII, which has been lost in mammals. Our work resolves nomenclature issues in GPI-linked fish CAs. Review of expression data brings forth previously unreported tissue and cancer types in which human CA IV is expressed. Analysis of collective glycosylation patterns of GPI-linked CAs suggests functionally important regions on the protein surface.

Conserved and quickly evolving immunome genes have different evolutionary paths.

Genetic, transcript, and protein level variations have important functional and evolutionary consequences. We performed systematic data collection and analysis of copy-number variations, single-nucleotide polymorphisms, disease-causing variations, messenger RNA splicing variants, and protein posttranslational modifications for the genes and proteins essential for human immune system. Information about polymorphic and evolutionarily fixed genetic variations was used to group immunome genes to the most conserved and the most quickly changing ones under directed selection during the recent immunome evolution. Gene Ontology terms related to adaptive immunity are associated with gene groups subject to recent directing selection. In addition, several other characteristics of the immunome genes and proteins in these two categories have statistically significant differences. The presented findings question the usability of directed mouse genes as models for human diseases and conditions and shed light on the fine tuning of human immunity and its diverse functions.

Systematic classification of primary immunodeficiencies based on clinical, pathological, and laboratory parameters.

The classification of diseases has several important applications ranging from diagnosis and choice of treatment to demographics. To date, classifications have been successfully created manually, often within international consortia. Some groups of diseases, such as primary immunodeficiencies (PIDs), are especially hard to nosologically cluster due, on one hand, to the presence of a wide variety of disorders and, in contrast, because of overlapping characteristics. More than 200 PIDs affecting components of the innate and adaptive immune systems have been described. Clinical, pathological, and laboratory characteristics were collected and used to group PIDs. A consensus of at least five independent methods provided a novel classification of 11 groups, which revealed previously unknown features and relationships of PIDs. Comparison of the classification to independent features, including the severity and therapy of the diseases, functional classification of proteins, and network vulnerability, indicated a strong statistical support. The method can be applied to any group of diseases.

Identification of candidate disease genes by integrating Gene Ontologies and protein-interaction networks: case study of primary immunodeficiencies.

Disease gene identification is still a challenge despite modern high-throughput methods. Many diseases are very rare or lethal and thus cannot be investigated with traditional methods. Several in silico methods have been developed but they have some limitations. We introduce a new method that combines information about protein-interaction network properties and Gene Ontology terms. Genes with high-calculated network scores and statistically significant gene ontology terms based on known diseases are prioritized as candidate genes. The method was applied to identify novel primary immunodeficiency-related genes, 26 of which were found. The investigation uses the protein-interaction network for all essential immunome human genes available in the Immunome Knowledge Base and an analysis of their enriched gene ontology annotations. The identified disease gene candidates are mainly involved in cellular signaling including receptors, protein kinases and adaptor and binding proteins as well as enzymes. The method can be generalized for any disease group with sufficient information.

Aspartic Aminopeptidase Is a Novel Biomarker of Aggressive Chronic Lymphocytic Leukemia.

Treatment of chronic lymphocytic leukemia has advanced substantially as our understanding of the kinase signal transduction pathways driven by the B cell receptor (BcR) has developed. Particularly, understanding the role of Bruton tyrosine kinase and phosphatidyl inositol 3 kinase delta in driving prosurvival signal transduction in chronic lymphocytic leukemia (CLL) cells and their targeting with pharmacological inhibitors (ibrutinib and idelalisib, respectively) has improved patient outcomes significantly. The kinase signaling pathway induced by the BcR is highly complex and has multiple interconnecting branches mediated by tyrosine and serine/threonine kinases activated downstream of the BcR. There is a high level of redundancy in the biological responses, with several BcR-signaling kinases driving nuclear factor kappa B activation or inducing antiapoptotic Bcl-2 genes. Accordingly, common gene targets of BcR-signaling kinases may serve as biomarkers indicating enhanced BCR-signaling and aggressive disease progression. This study used a gene expression correlation analysis of malignant B cell lines and primary CLL cells to identify genes whose expression correlated with BCR-signaling kinases overexpressed and/or overactivated in CLL, namely: AKT1, AKT2, BTK, MAPK1, MAPK3, PI3KCD and ZAP70. The analysis identified a 32-gene signature with a strong prognostic potential and DNPEP, the gene coding for aspartic aminopeptidase, as a predictor of aggressive CLL. DNPEP gene expression correlated with MAPK3, PI3KCD, and ZAP70 expression and, in the primary CLL test dataset, showed a strong prognostic potential. The inhibition of DNPEP with a pharmacological inhibitor enhanced the cytotoxic potential of idelalisib and ibrutinib, indicating a biological functionality of DNPEP in CLL. DNPEP, as an aminopeptidase, contributes to the maintenance of the free amino acid pool in CLL cells found to be an essential process for the survival of many cancer cell types, and thus, these results warrant further research into the exploitation of aminopeptidase inhibitors in the treatment of drug-resistant CLL.

Immunome knowledge base (IKB): an integrated service for immunome research.

BACKGROUND: Functioning of the immune system requires the coordinated expression and action of many genes and proteins. With the emergence of high-throughput technologies, a great amount of molecular data is available for the genes and proteins of the immune system. However, these data are scattered into several databases and literature and therefore integration is needed. DESCRIPTION: The Immunome Knowledge Base (IKB) is a dedicated resource for immunological information. We identified and collected genes that are essential for the immunome. Nucleotide and protein sequences, as well as information about the related pseudogenes are available for 893 human essential immunome genes. To allow the study of the evolution of the immune system, data for the orthologs of human genes was collected. In addition to the human immunome, ortholog groups of 1811 metazoan immunity genes are available with information about the evidence of their immunity function. IKB combines three previous databases and several additional data items in an integrated system. CONCLUSION: IKB provides in one single service access to several databases and resources and contains plenty of new data about immune system. The most recent addition is variation data on genomic, transcriptomic and proteomic levels for all the immunome genes and proteins. In the future, more data will be added on the function of these genes. The service has a free and public web interface.

PseudoGeneQuest - service for identification of different pseudogene types in the human genome.

BACKGROUND: Pseudogenes, nonfunctional copies of genes, evolve fast due the lack of evolutionary pressures and thus appear in several different forms. PseudoGeneQuest is an online tool to search the human genome for a given query sequence and to identify different types of pseudogenes as well as novel genes and gene fragments. DESCRIPTION: The service can detect pseudogenes, that have arisen either by retrotransposition or segmental genome duplication, many of which are not listed in the public pseudogene databases. The service has a user-friendly web interface and uses a powerful computer cluster in order to perform parallel searches and provide relatively fast runtimes despite exhaustive database searches and analyses. CONCLUSION: PseudoGeneQuest is a versatile tool for detecting novel pseudogene candidates from the human genome. The service searches human genome sequences for five types of pseudogenes and provides an output that allows easy further analysis of observations. In addition to the result file the system provides visualization of the results linked to Ensembl Genome Browser. PseudoGeneQuest service is freely available.

Distribution and evolution of short tandem repeats in closely related bacterial genomes.

Simultaneous identification and comparison of perfect and imperfect microsatellites within a genome is a valuable tool both to overcome the lack of a consensus definition of SSRs and to assess repeat history. Detailed analysis of the overall distribution of perfect and imperfect microsatellites in closely related bacterial taxa is expected to give new insight into the evolution of prokaryotic genomes. We have performed a genome-wide analysis of microsatellite distribution in four Escherichia coli and seven Chlamydial strains. Chlamydial strains generally have a higher density of SSRs and show greater intra-group differences of SSR distribution patterns than E. coli genomes. In most investigated genomes the distribution of the total lengths of matching perfect and imperfect trinucleotide repeats are highly similar, with the notable exception of C. muridarum. Closely related strains show more similar repeat distribution patterns than strains separated by a longer divergence time. The discrepancy between the preferred classes of perfect and imperfect repeats in C. muridarum implies accelerated evolution of SSRs in this particular strain. Our results suggest that microsatellites, although considerably less abundant than in eukaryotic genomes, may nevertheless play an important role in the evolution of prokaryotic genomes and several gene families.

Repression of Mcl-1 expression by the CDC7/CDK9 inhibitor PHA-767491 overcomes bone marrow stroma-mediated drug resistance in AML.

Acute myeloid leukaemia (AML) is an aggressive cancer with 50-75% of patients relapsing even after successful chemotherapy. The role of the bone marrow microenvironment (BMM) in protecting AML cells from chemotherapeutics and causing consequent relapse is increasingly recognised. However the role that the anti-apoptotic Bcl-2 proteins play as effectors of BMM-mediated drug resistance are less understood. Here we show that bone marrow mesenchymal stromal cells (BMSC) provide resistance to AML cells against BH3-mimetics, cytarabine and daunorubicin, but this is not mediated by Bcl-2 and/or Bcl-XL as previously thought. Instead, BMSCs induced Mcl-1 expression over Bcl-2 and/or Bcl-XL in AML cells and inhibition of Mcl-1 with a small-molecule inhibitor, A1210477, or repressing its expression with the CDC7/CDK9 dual-inhibitor, PHA-767491 restored sensitivity to BH3-mimetics. Furthermore, combined inhibition of Bcl-2/Bcl-XL and Mcl-1 could revert BMSC-mediated resistance against cytarabine + daunorubicin. Importantly, the CD34+/CD38- leukemic stem cell-encompassing population was equally sensitive to the combination of PHA-767491 and ABT-737. These results indicate that Bcl-2/Bcl-XL and Mcl-1 act in a redundant fashion as effectors of BMM-mediated AML drug resistance and highlight the potential of Mcl-1-repression to revert BMM-mediated drug resistance in the leukemic stem cell population, thus, prevent disease relapse and ultimately improve patient survival.

Sex-specific regulation of aging in Caenorhabditis elegans.

A fascinating aspect of sexual dimorphism in various animal species is that the two sexes differ substantially in lifespan. In humans, for example, women's life expectancy exceeds that of men by 3-7 years. Whether this trait can be attributed to dissimilar lifestyles or genetic (regulatory) factors remains to be elucidated. Herein, we demonstrate that in the nematode Caenorhabditis elegans, the significantly longer lifespan of hermaphrodites-which are essentially females capable of sperm production-over males is established by TRA-1, the terminal effector of the sex-determination pathway. This transcription factor directly controls the expression of daf-16/FOXO, which functions as a major target of insulin/IGF-1 signaling (IIS) and key modulator of aging across diverse animal phyla. TRA-1 extends hermaphrodite lifespan through promoting daf-16 activity. Furthermore, TRA-1 also influences reproductive growth in a DAF-16-dependent manner. Thus, the sex-determination machinery is an important regulator of IIS in this organism. These findings provide a mechanistic insight into how longevity and development are specified unequally in the two genders. As TRA-1 is orthologous to mammalian GLI (glioma-associated) proteins, a similar sex-specific mechanism may also operate in humans to determine lifespan.

Divergent microsatellite evolution in the human and chimpanzee lineages.

Comparison of the complete human genome sequence to one of its closest relatives, the chimpanzee genome, provides a unique opportunity for exploring recent evolutionary events affecting the microsatellites in these species. A simple assumption on microsatellite distribution is that the total length of perfect repeats is constant compared to that of imperfect ones regardless of the repeat sequence. In this paper, we show that this is valid for most of the chimpanzee genome but not for a number of human chromosomes. Our results suggest accelerated evolution of microsatellites in the human genome relative to the chimpanzee lineage.

DoOP: Databases of Orthologous Promoters, collections of clusters of orthologous upstream sequences from chordates and plants.

DoOP (http://doop.abc.hu/) is a database of eukaryotic promoter sequences (upstream regions) aiming to facilitate the recognition of regulatory sites conserved between species. The annotated first exons of human and Arabidopsis thaliana genes were used as queries in BLAST searches to collect the most closely related orthologous first exon sequences from Chordata and Viridiplantae species. Up to 3000 bp DNA segments upstream from these first exons constitute the clusters in the chordate and plant sections of the Database of Orthologous Promoters. Release 1.0 of DoOP contains 21,061 chordate clusters from 284 different species and 7548 plant clusters from 269 different species. The database can be used to find and retrieve promoter sequences of a given gene from various species and it is also suitable to see the most trivial conserved sequence blocks in the orthologous upstream regions. Users can search DoOP with either sequence or text (annotation) to find promoter clusters of various genes. In addition to the sequence data, the positions of the conserved sequence blocks derived from multiple alignments, the positions of repetitive elements and the positions of transcription start sites known from the Eukaryotic Promoter Database (EPD) can be viewed graphically.

Identification and characterization of a novel zebrafish (Danio rerio) pentraxin-carbonic anhydrase.

BACKGROUND: Carbonic anhydrases (CAs) are ubiquitous, essential enzymes which catalyze the conversion of carbon dioxide and water to bicarbonate and H+ ions. Vertebrate genomes generally contain gene loci for 15-21 different CA isoforms, three of which are enzymatically inactive. CA VI is the only secretory protein of the enzymatically active isoforms. We discovered that non-mammalian CA VI contains a C-terminal pentraxin (PTX) domain, a novel combination for both CAs and PTXs. METHODS: We isolated and sequenced zebrafish (Danio rerio) CA VI cDNA, complete with the sequence coding for the PTX domain, and produced the recombinant CA VI-PTX protein. Enzymatic activity and kinetic parameters were measured with a stopped-flow instrument. Mass spectrometry, analytical gel filtration and dynamic light scattering were used for biophysical characterization. Sequence analyses and Bayesian phylogenetics were used in generating hypotheses of protein structure and CA VI gene evolution. A CA VI-PTX antiserum was produced, and the expression of CA VI protein was studied by immunohistochemistry. A knock-down zebrafish model was constructed, and larvae were observed up to five days post-fertilization (dpf). The expression of ca6 mRNA was quantitated by qRT-PCR in different developmental times in morphant and wild-type larvae and in different adult fish tissues. Finally, the swimming behavior of the morphant fish was compared to that of wild-type fish. RESULTS: The recombinant enzyme has a very high carbonate dehydratase activity. Sequencing confirms a 530-residue protein identical to one of the predicted proteins in the Ensembl database (ensembl.org). The protein is pentameric in solution, as studied by gel filtration and light scattering, presumably joined by the PTX domains. Mass spectrometry confirms the predicted signal peptide cleavage and disulfides, and N-glycosylation in two of the four observed glycosylation motifs. Molecular modeling of the pentamer is consistent with the modifications observed in mass spectrometry. Phylogenetics and sequence analyses provide a consistent hypothesis of the evolutionary history of domains associated with CA VI in mammals and non-mammals. Briefly, the evidence suggests that ancestral CA VI was a transmembrane protein, the exon coding for the cytoplasmic domain was replaced by one coding for PTX domain, and finally, in the therian lineage, the PTX-coding exon was lost. We knocked down CA VI expression in zebrafish embryos with antisense morpholino oligonucleotides, resulting in phenotype features of decreased buoyancy and swim bladder deflation in 4 dpf larvae. DISCUSSION: These findings provide novel insights into the evolution, structure, and function of this unique CA form.

KinMutBase: a registry of disease-causing mutations in protein kinase domains.

A large number of disease-causing mutations have been identified from several protein kinases. KinMutBase is a comprehensive knowledge base for human disease-related mutations in protein kinase domains (http://bioinf.uta.fi/KinMutBase/). The latest version contains 582 different mutations for 1,790 cases in 1,322 families. KinMutBase entries are described on the DNA, mRNA, and protein level. Numbers for affected patients and families are also provided. KinMutBase has extensive amount of links and cross-references to literature, other databases, and information sources. There are numerous interactive pages about sequences, structures, mutation statistics, and diseases. Detailed statistical study was done on frequencies of different types of mutations both on the DNA and protein level in serine/threonine kinase (PSK) and tyrosine kinase (PTK). Three-dimensional structures indicate clustering of disease-related mutations mainly to conserved subdomains, and substrate and coligand binding amino acids, although mutations appear throughout the sequences. CpG containing codons, especially for arginine, constitute the majority of mutational hotspots. There are certain clear differences in mutation patterns and types between PSKs and PTKs.