From protein sequences to 3D-structures and beyond: the example of the UniProt knowledgebase.

With the dramatic increase in the volume of experimental results in every domain of life sciences, assembling pertinent data and combining information from different fields has become a challenge. Information is dispersed over numerous specialized databases and is presented in many different formats. Rapid access to experiment-based information about well-characterized proteins helps predict the function of uncharacterized proteins identified by large-scale sequencing. In this context, universal knowledgebases play essential roles in providing access to data from complementary types of experiments and serving as hubs with cross-references to many specialized databases. This review outlines how the value of experimental data is optimized by combining high-quality protein sequences with complementary experimental results, including information derived from protein 3D-structures, using as an example the UniProt knowledgebase (UniProtKB) and the tools and links provided on its website ( http://www.uniprot.org/ ). It also evokes precautions that are necessary for successful predictions and extrapolations.

Expression of cysteine proteinase type I and II of Leishmania infantum and their recognition by sera during canine and human visceral leishmaniasis.

In this study, the mature domains of type I (CPB) and type II (CPA) cysteine proteinases (CPs) of Leishmania infantum were expressed and their immunogenic properties defined using sera from active and recovered cases of human visceral leishmaniasis and sera from infected dogs. Immunoblotting and ELISA analysis indicated that a freeze/thaw extract of parasite antigens showed similar and intensive recognition in both active cases of human and dog sera but lower recognition in recovered human individuals. The total IgG of actively infected human sera was higher than in recovered cases when rCPs were used as antigen. In contrast to dog sera, both active and recovered human cases have higher recognition toward rCPB than rCPA. Furthermore, the asymptomatic dogs in contrast to the symptomatic cases exhibited specific lymphocyte proliferation to both crude antigens and rCPs.

Feasibility of RNA studies on illegitimate transcription for molecular characterization of splicing mutations in the ATP7B gene: a case report.

Approximately 520 Wilson disease-causing mutations in the ATP7B gene have been described to date. In this study we report DNA and RNA analyses carried out for molecular characterization of a consensus sequence splicing mutation found in homozygosity in a Swiss Wilson disease patient. RNA analysis of 1946 +6 T→C in both the peripheral lymphoblasts and liver resulted in the production in the propositus of only an alternative transcript lacking exons 6, 7, and 8 resulting most likely in alterations of cell biochemistry and disease. The patient presents an early form of severe hepatic disease characterized by hepatosplenomegaly, reduced hepatic function, anemia and thrombocytopenia indicating that 1946 +6 T→C is a severe mutation. Since identical results were obtained from both peripheral lymphoblasts and liver they also suggest that RNA studies of illegitimate transcripts can be safely used for molecular characterization of ATP7B splicing mutations, thus improving genetic counseling and diagnosis of Wilson disease. Moreover these studies, contribute to reveal the exact molecular mechanisms producing Wilson disease.

The G0/G1 switch gene 2 is a novel PPAR target gene.

PPARs (peroxisome-proliferator-activated receptors) alpha, beta/delta and gamma are a group of transcription factors that are involved in numerous processes, including lipid metabolism and adipogenesis. By comparing liver mRNAs of wild-type and PPARalpha-null mice using microarrays, a novel putative target gene of PPARalpha, G0S2 (G0/G1 switch gene 2), was identified. Hepatic expression of G0S2 was up-regulated by fasting and by the PPARalpha agonist Wy14643 in a PPARalpha-dependent manner. Surprisingly, the G0S2 mRNA level was highest in brown and white adipose tissue and was greatly up-regulated during mouse 3T3-L1 and human SGBS (Simpson-Golabi-Behmel syndrome) adipogenesis. Transactivation, gel shift and chromatin immunoprecipitation assays indicated that G0S2 is a direct PPARgamma and probable PPARalpha target gene with a functional PPRE (PPAR-responsive element) in its promoter. Up-regulation of G0S2 mRNA seemed to be specific for adipogenesis, and was not observed during osteogenesis or myogenesis. In 3T3-L1 fibroblasts, expression of G0S2 was associated with growth arrest, which is required for 3T3-L1 adipogenesis. Together, these data indicate that G0S2 is a novel target gene of PPARs that may be involved in adipocyte differentiation.

A downstream mediator in the growth repression limb of the jasmonate pathway.

Wounding plant tissues initiates large-scale changes in transcription coupled to growth arrest, allowing resource diversion for defense. These processes are mediated in large part by the potent lipid regulator jasmonic acid (JA). Genes selected from a list of wound-inducible transcripts regulated by the jasmonate pathway were overexpressed in Arabidopsis thaliana, and the transgenic plants were then assayed for sensitivity to methyl jasmonate (MeJA). When grown in the presence of MeJA, the roots of plants overexpressing a gene of unknown function were longer than those of wild-type plants. When transcript levels for this gene, which we named JASMONATE-ASSOCIATED1 (JAS1), were reduced by RNA interference, the plants showed increased sensitivity to MeJA and growth was inhibited. These gain- and loss-of-function assays suggest that this gene acts as a repressor of JA-inhibited growth. An alternative transcript from the gene encoding a second protein isoform with a longer C terminus failed to repress jasmonate sensitivity. This identified a conserved C-terminal sequence in JAS1 and related genes, all of which also contain Zim motifs and many of which are jasmonate-regulated. Both forms of JAS1 were found to localize to the nucleus in transient expression assays. Physiological tests of growth responses after wounding were consistent with the fact that JAS1 is a repressor of JA-regulated growth retardation.

Molecular phylogeography of the nose-horned viper (Vipera ammodytes, Linnaeus (1758)): evidence for high genetic diversity and multiple refugia in the Balkan peninsula

The nose-horned viper (Vipera ammodytes) occurs in a large part of the south-eastern Europe and Asia Minor. Phylogenetic relationships were reconstructed for a total of 59 specimens using sequences from three mitochondrial regions (16S and cytochrome b genes, and control region, totalling 2308 bp). A considerable number of clades were observed within this species, showing a large genetic diversity within the Balkan peninsula. Splitting of the basal clades was evaluated to about 4 million years ago. Genetic results are in contradiction with presently accepted taxonomy based on morphological characters: V. a. gregorwallneri and V. a. ruffoi do not display any genetic difference compared with the nominotypic subspecies (V. a. ammodytes), involving that these subspecies can be regarded as synonyms. High genetic divergence in the central part of the Balkan peninsula is not concordant with low morphological differentiation. Finally, the extensive genetic diversity within the Balkan peninsula and the colonisation routes are discussed

A genetic basis for functional hypothalamic amenorrhea.

Background: Functional hypothalamic amenorrhea is a reversible form of gonadotropin-releasing hormone (GnRH) deficiency commonly triggered by stressors such as excessive exercise, nutritional deficits, or psychological distress. Women vary in their susceptibility to inhibition of the reproductive axis by such stressors, but it is unknown whether this variability reflects a genetic predisposition to hypothalamic amenorrhea. We hypothesized that mutations in genes involved in idiopathic hypogonadotropic hypogonadism, a congenital form of GnRH deficiency, are associated with hypothalamic amenorrhea. Methods: We analyzed the coding sequence of genes associated with idiopathic hypogonadotropic hypogonadism in 55 women with hypothalamic amenorrhea and performed in vitro studies of the identified mutations. Results: Six heterozygous mutations were identified in 7 of the 55 patients with hypothalamic amenorrhea: two variants in the fibroblast growth factor receptor 1 gene FGFR1 (G260E and R756H), two in the prokineticin receptor 2 gene PROKR2 (R85H and L173R), one in the GnRH receptor gene GNRHR (R262Q), and one in the Kallmann syndrome 1 sequence gene KAL1 (V371I). No mutations were found in a cohort of 422 controls with normal menstrual cycles. In vitro studies showed that FGFR1 G260E, FGFR1 R756H, and PROKR2 R85H are loss-of-function mutations, as has been previously shown for PROKR2 L173R and GNRHR R262Q. Conclusions: Rare variants in genes associated with idiopathic hypogonadotropic hypogonadism are found in women with hypothalamic amenorrhea, suggesting that these mutations may contribute to the variable susceptibility of women to the functional changes in GnRH secretion that characterize hypothalamic amenorrhea. Our observations provide evidence for the role of rare variants in common multifactorial disease. (Funded by the Eunice Kennedy Shriver National Institute of Child Health and Human Development and others; ClinicalTrials.gov number, NCT00494169.)

S100A1 deficiency results in prolonged ventricular repolarization in response to sympathetic activation.

S100A1 is a Ca(2+)-binding protein and predominantly expressed in the heart. We have generated a mouse line of S100A1 deficiency by gene trap mutagenesis to investigate the impact of S100A1 ablation on heart function. Electrocardiogram recordings revealed that after beta-adrenergic stimulation S100A1-deficient mice had prolonged QT, QTc and ST intervals and intraventricular conduction disturbances reminiscent of 2 : 1 bundle branch block. In order to identify genes affected by the loss of S100A1, we profiled the mutant and wild type cardiac transcriptomes by gene array analysis. The expression of several genes functioning to the electrical activity of the heart were found to be significantly altered. Although the default prediction would be that mRNA and protein levels are highly correlated, comprehensive immunoblot analyses of salient up- or down-regulated candidate genes of any cellular network revealed no significant changes on protein level. Taken together, we found that S100A1 deficiency results in cardiac repolarization delay and alternating ventricular conduction defects in response to sympathetic activation accompanied by a significantly different transcriptional regulation.

Membrane mu poly(A) signal and 3' flanking sequences function as a transcription terminator for immunoglobulin-encoding genes.

Developmentally regulated mechanisms involving alternative RNA splicing and/or polyadenylation, as well as transcription termination, are implicated in controlling the levels of secreted mu (mu s), membrane mu (mu m) and delta immunoglobulin (Ig) heavy chain mRNAs during B cell differentiation (mu gene encodes the mu heavy chain). Using expression vectors constructed with genomic DNA segments composed of the mu m polyadenylation signal region, we analyzed poly(A) site utilization and termination of transcription in stably transfected myeloma cells and in murine fibroblast L cells. We found that the gene segment containing the mu m poly(A) signals, along with 536 bp of downstream flanking sequence, acted as a transcription terminator in both myeloma cells and L cell fibroblasts. Neither a 141-bp DNA fragment (which directed efficient polyadenylation at the mu m site), nor the 536-bp flanking nucleotide sequence alone, were sufficient to obtain a similar regulation. This shows that the mu m poly(A) region plays a central role in controlling developmentally regulated transcription termination by blocking downstream delta gene expression. Because this gene segment exhibited the same RNA processing and termination activities in fibroblasts, it appears that these processes are not tissue-specific.

Mutation patterns of amino acid tandem repeats in the human proteome

Background: Amino acid tandem repeats are found in nearly one-fifth of human proteins. Abnormal expansion of these regions is associated with several human disorders. To gain further insight into the mutational mechanisms that operate in this type of sequence, we have analyzed a large number of mutation variants derived from human expressed sequence tags (ESTs).Results: We identified 137 polymorphic variants in 115 different amino acid tandem repeats. Of these, 77 contained amino acid substitutions and 60 contained gaps (expansions or contractions of the repeat unit). The analysis showed that at least about 21% of the repeats might be polymorphic in humans. We compared the mutations found in different types of amino acid repeats and in adjacent regions. Overall, repeats showed a five-fold increase in the number of gap mutations compared to adjacent regions, reflecting the action of slippage within the repetitive structures. Gap and substitution mutations were very differently distributed between different amino acid repeat types. Among repeats containing gap variants we identified several disease and candidate disease genes.Conclusion: This is the first report at a genome-wide scale of the types of mutations occurring in the amino acid repeat component of the human proteome. We show that the mutational dynamics of different amino acid repeat types are very diverse. We provide a list of loci with highly variable repeat structures, some of which may be potentially involved in disease.

Sequencing human-gibbon breakpoints of synteny reveals mosaic new insertions at rearrangement sites

The gibbon genome exhibits extensive karyotypic diversity with an increased rate of chromosomal rearrangements during evolution. In an effort to understand the mechanistic origin and implications of these rearrangement events, we sequenced 24 synteny breakpoint regions in the white-cheeked gibbon (Nomascus leucogenys, NLE) in the form of high-quality BAC insert sequences (4.2 Mbp). While there is a significant deficit of breakpoints in genes, we identified seven human gene structures involved in signaling pathways (DEPDC4, GNG10), phospholipid metabolism (ENPP5, PLSCR2), beta-oxidation (ECH1), cellular structure and transport (HEATR4), and transcription (ZNF461), that have been disrupted in the NLE gibbon lineage. Notably, only three of these genes show the expected evolutionary signatures of pseudogenization. Sequence analysis of the breakpoints suggested both nonclassical nonhomologous end-joining (NHEJ) and replication-based mechanisms of rearrangement. A substantial number (11/24) of human-NLE gibbon breakpoints showed new insertions of gibbon-specific repeats and mosaic structures formed from disparate sequences including segmental duplications, LINE, SINE, and LTR elements. Analysis of these sites provides a model for a replication-dependent repair mechanism for double-strand breaks (DSBs) at rearrangement sites and insights into the structure and formation of primate segmental duplications at sites of genomic rearrangements during evolution.

Sensitive gene expression profiling of human T cell subsets reveals parallel post-thymic differentiation for CD4+ and CD8+ lineages.

The differentiation of CD4(+) or CD8(+) T cells following priming of naive cells is central in the establishment of the immune response against pathogens or tumors. However, our understanding of this complex process and the significance of the multiple subsets of differentiation remains controversial. Gene expression profiling has opened new directions of investigation in immunobiology. Nonetheless, the need for substantial amount of biological material often limits its application range. In this study, we have developed procedures to perform microarray analysis on amplified cDNA from low numbers of cells, including primary T lymphocytes, and applied this technology to the study of CD4 and CD8 lineage differentiation. Gene expression profiling was performed on samples of 1000 cells from 10 different subpopulations, defining the major stages of post-thymic CD4(+) or CD8(+) T cell differentiation. Surprisingly, our data revealed that while CD4(+) and CD8(+) T cell gene expression programs diverge at early stages of differentiation, they become increasingly similar as cells reach a late differentiation stage. This suggests that functional heterogeneity between Ag experienced CD4(+) and CD8(+) T cells is more likely to be located early during post-thymic differentiation, and that late stages of differentiation may represent a common end in the development of T-lymphocytes.

Numerous potentially functional but non-genic conserved sequences on human chromosome 21.

The use of comparative genomics to infer genome function relies on the understanding of how different components of the genome change over evolutionary time. The aim of such comparative analysis is to identify conserved, functionally transcribed sequences such as protein-coding genes and non-coding RNA genes, and other functional sequences such as regulatory regions, as well as other genomic features. Here, we have compared the entire human chromosome 21 with syntenic regions of the mouse genome, and have identified a large number of conserved blocks of unknown function. Although previous studies have made similar observations, it is unknown whether these conserved sequences are genes or not. Here we present an extensive experimental and computational analysis of human chromosome 21 in an effort to assign function to sequences conserved between human chromosome 21 (ref. 8) and the syntenic mouse regions. Our data support the presence of a large number of potentially functional non-genic sequences, probably regulatory and structural. The integration of the properties of the conserved components of human chromosome 21 to the rapidly accumulating functional data for this chromosome will improve considerably our understanding of the role of sequence conservation in mammalian genomes.

Adaptation of aerobically growing Pseudomonas aeruginosa to copper starvation.

Restricted bioavailability of copper in certain environments can interfere with cellular respiration because copper is an essential cofactor of most terminal oxidases. The global response of the metabolically versatile bacterium and opportunistic pathogen Pseudomonas aeruginosa to copper limitation was assessed under aerobic conditions. Expression of cioAB (encoding an alternative, copper-independent, cyanide-resistant ubiquinol oxidase) was upregulated, whereas numerous iron uptake functions (including the siderophores pyoverdine and pyochelin) were expressed at reduced levels, presumably reflecting a lower demand for iron by respiratory enzymes. Wild-type P. aeruginosa was able to grow aerobically in a defined glucose medium depleted of copper, whereas a cioAB mutant did not grow. Thus, P. aeruginosa relies on the CioAB enzyme to cope with severe copper deprivation. A quadruple cyo cco1 cco2 cox mutant, which was deleted for all known heme-copper terminal oxidases of P. aeruginosa, grew aerobically, albeit more slowly than did the wild type, indicating that the CioAB enzyme is capable of energy conservation. However, the expression of a cioA'-'lacZ fusion was less dependent on the copper status in the quadruple mutant than in the wild type, suggesting that copper availability might affect cioAB expression indirectly, via the function of the heme-copper oxidases.