MyHits: improvements to an interactive resource for analyzing protein sequences.

The MyHits web site (http://myhits.isb-sib.ch) is an integrated service dedicated to the analysis of protein sequences. Since its first description in 2004, both the user interface and the back end of the server were improved. A number of tools (e.g. MAFFT, Jacop, Dotlet, Jalview, ESTScan) were added or updated to improve the usability of the service. The MySQL schema and its associated API were revamped and the database engine (HitKeeper) was separated from the web interface. This paper summarizes the current status of the server, with an emphasis on the new services.

Improving the sensitivity of the sequence profile method.

The sequence profile method (Gribskov M, McLachlan AD, Eisenberg D, 1987, Proc Natl Acad Sci USA 84:4355-4358) is a powerful tool to detect distant relationships between amino acid sequences. A profile is a table of position-specific scores and gap penalties, providing a generalized description of a protein motif, which can be used for sequence alignments and database searches instead of an individual sequence. A sequence profile is derived from a multiple sequence alignment. We have found 2 ways to improve the sensitivity of sequence profiles: (1) Sequence weights: Usage of individual weights for each sequence avoids bias toward closely related sequences. These weights are automatically assigned based on the distance of the sequences using a published procedure (Sibbald PR, Argos P, 1990, J Mol Biol 216:813-818). (2) Amino acid substitution table: In addition to the alignment, the construction of a profile also needs an amino acid substitution table. We have found that in some cases a new table, the BLOSUM45 table (Henikoff S, Henikoff JG, 1992, Proc Natl Acad Sci USA 89:10915-10919), is more sensitive than the original Dayhoff table or the modified Dayhoff table used in the current implementation. Profiles derived by the improved method are more sensitive and selective in a number of cases where previous methods have failed to completely separate true members from false positives.

Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila.

Ionotropic glutamate receptors (iGluRs) mediate neuronal communication at synapses throughout vertebrate and invertebrate nervous systems. We have characterized a family of iGluR-related genes in Drosophila, which we name ionotropic receptors (IRs). These receptors do not belong to the well-described kainate, AMPA, or NMDA classes of iGluRs, and they have divergent ligand-binding domains that lack their characteristic glutamate-interacting residues. IRs are expressed in a combinatorial fashion in sensory neurons that respond to many distinct odors but do not express either insect odorant receptors (ORs) or gustatory receptors (GRs). IR proteins accumulate in sensory dendrites and not at synapses. Misexpression of IRs in different olfactory neurons is sufficient to confer ectopic odor responsiveness. Together, these results lead us to propose that the IRs comprise a novel family of chemosensory receptors. Conservation of IR/iGluR-related proteins in bacteria, plants, and animals suggests that this receptor family represents an evolutionarily ancient mechanism for sensing both internal and external chemical cues.

The complete amino acid sequence for brain beta spectrin (beta fodrin): relationship to globin sequences.

The amino acid sequence of mouse brain beta spectrin (beta fodrin), deduced from the nucleotide sequence of complementary DNA clones, reveals that this non-erythroid beta spectrin comprises 2363 residues, with a molecular weight of 274,449 Da. Brain beta spectrin contains three structural domains and we suggest the position of several functional domains including f-actin, synapsin I, ankyrin and spectrin self association sites. Analysis of deduced amino acid sequences indicated striking homology and similar structural characteristics of brain beta spectrin repeats beta 11 and beta 12 to globins. In vitro analysis has demonstrated that heme is capable of specific attachment to brain spectrin, suggesting possible new functions in electron transfer, oxygen binding, nitric oxide binding or heme scavenging.

Global Genomic Diversity of Human Papillomavirus 6 Based on 724 Isolates and 190 Complete Genome Sequences.

Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution. IMPORTANCE: This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

DNA binding properties of peroxisome proliferator-activated receptor subtypes on various natural peroxisome proliferator response elements. Importance of the 5'-flanking region.

The three subtypes of the peroxisome proliferator-activated receptors (PPARalpha, beta/delta, and gamma) form heterodimers with the 9-cis-retinoic acid receptor (RXR) and bind to a common consensus response element, which consists of a direct repeat of two hexanucleotides spaced by one nucleotide (DR1). As a first step toward understanding the molecular mechanisms determining PPAR subtype specificity, we evaluated by electrophoretic mobility shift assays the binding properties of the three PPAR subtypes, in association with either RXRalpha or RXRgamma, on 16 natural PPAR response elements (PPREs). The main results are as follows. (i) PPARgamma in combination with either RXRalpha or RXRgamma binds more strongly than PPARalpha or PPARbeta to all natural PPREs tested. (ii) The binding of PPAR to strong elements is reinforced if the heterodimerization partner is RXRgamma. In contrast, weak elements favor RXRalpha as heterodimerization partner. (iii) The ordering of the 16 natural PPREs from strong to weak elements does not depend on the core DR1 sequence, which has a relatively uniform degree of conservation, but correlates with the number of identities of the 5'-flanking nucleotides with respect to a consensus element. This 5'-flanking sequence is essential for PPARalpha binding and thus contributes to subtype specificity. As a demonstration of this, the PPARgamma-specific element ARE6 PPRE is able to bind PPARalpha only if its 5'-flanking region is exchanged with that of the more promiscuous HMG PPRE.

The extended GLUT-family of sugar/polyol transport facilitators: nomenclature, sequence characteristics, and potential function of its novel members (review).

During the last 2 years, several novel genes that encode glucose transporter-like proteins have been identified and characterized. Because of their sequence similarity with GLUT1, these genes appear to belong to the family of solute carriers 2A (SLC2A, protein symbol GLUT). Sequence comparisons of all 13 family members allow the definition of characteristic sugar/polyol transporter signatures: (1) the presence of 12 membrane-spanning helices, (2) seven conserved glycine residues in the helices, (3) several basic and acidic residues at the intracellular surface of the proteins, (4) two conserved tryptophan residues, and (5) two conserved tyrosine residues. On the basis of sequence similarities and characteristic elements, the extended GLUT family can be divided into three subfamilies, namely class I (the previously known glucose transporters GLUT1-4), class II (the previously known fructose transporter GLUT5, the GLUT7, GLUT9 and GLUT11), and class III (GLUT6, 8, 10, 12, and the myo-inositol transporter HMIT1). Functional characteristics have been reported for some of the novel GLUTs. Like GLUT1-4, they exhibit a tissue/cell-specific expression (GLUT6, leukocytes, brain; GLUT8, testis, blastocysts, brain, muscle, adipocytes; GLUT9, liver, kidney; GLUT10, liver, pancreas; GLUT11, heart, skeletal muscle). GLUT6 and GLUT8 appear to be regulated by sub-cellular redistribution, because they are targeted to intra-cellular compartments by dileucine motifs in a dynamin dependent manner. Sugar transport has been reported for GLUT6, 8, and 11; HMIT1 has been shown to be a H+/myo-inositol co-transporter. Thus, the members of the extended GLUT family exhibit a surprisingly diverse substrate specificity, and the definition of sequence elements determining this substrate specificity will require a full functional characterization of all members.

HIV-1 co/super-infection in intravenous drug users.

BACKGROUND: The frequency of HIV-1 co/super-infection is unknown despite their implications for public health and vaccine development. This issue was addressed during an epidemic of both CRF11 and B subtype among intravenous drug users (IVDUs). METHODS: Bulk sequencing of reverse transcriptase, protease and C2V3 regions and subtype-specific nested polymerase chain reaction (PCR) in plasma and proviral DNA were performed using baseline and follow-up samples collected in recently infected IVDUs between 1998-2002 and in IVDUs with chronic infection living in the same area and presenting an unexpected rise of viremia (> 1 log10). RESULTS: In 58 recently infected patients, three B/CRF-11 co-infections, 25 B, 28 CRF-11 and two other subtypes were detected at baseline. In the three co-infected patients, both CRF-11 and B were detected in plasma and proviral DNA and persisted during follow-up. B- and CFR-11-specific PCR performed on follow-up samples of 40 of 58 recently infected patients (median follow-up, 14.5 months) revealed a transient B super-infection in a patient initially infected by CRF-11. Five of 156 chronic IVDUs (total follow-up: 346 years) had an unexpected rise of viremia. In two of them, aviremic without treatment for years after an initial B infection, a symptomatic CRF-11 super-infection occurred and was associated with high viral load and a fall of CD4 cell count. CONCLUSIONS: In recently infected IVDUs, co-infection B/CRF-11 is relatively frequent (5%). In chronically infected IVDUs super-infection may be transient and may occur in patients controlling efficiently HIV infection by the initial strain.

Analyses of six homologous proteins of Protochlamydia amoebophila UWE25 encoded by large GC-rich genes (lgr): a model of evolution and concatenation of leucine-rich repeats.

BACKGROUND: Along the chromosome of the obligate intracellular bacteria Protochlamydia amoebophila UWE25, we recently described a genomic island Pam100G. It contains a tra unit likely involved in conjugative DNA transfer and lgrE, a 5.6-kb gene similar to five others of P. amoebophila: lgrA to lgrD, lgrF. We describe here the structure, regulation and evolution of these proteins termed LGRs since encoded by "Large G+C-Rich" genes. RESULTS: No homologs to the whole protein sequence of LGRs were found in other organisms. Phylogenetic analyses suggest that serial duplications producing the six LGRs occurred relatively recently and nucleotide usage analyses show that lgrB, lgrE and lgrF were relocated on the chromosome. The C-terminal part of LGRs is homologous to Leucine-Rich Repeats domains (LRRs). Defined by a cumulative alignment score, the 5 to 18 concatenated octacosapeptidic (28-meric) LRRs of LGRs present all a predicted alpha-helix conformation. Their closest homologs are the 28-residue RI-like LRRs of mammalian NODs and the 24-meres of some Ralstonia and Legionella proteins. Interestingly, lgrE, which is present on Pam100G like the tra operon, exhibits Pfam domains related to DNA metabolism. CONCLUSION: Comparison of the LRRs, enable us to propose a parsimonious evolutionary scenario of these domains driven by adjacent concatenations of LRRs. Our model established on bacterial LRRs can be challenged in eucaryotic proteins carrying less conserved LRRs, such as NOD proteins and Toll-like receptors.

Cerebellar cleft: confirmation of the neuroimaging pattern.

We recently described the neuroimaging and clinical findings in 6 children with cerebellar clefts and proposed that they result from disruptive changes following prenatal cerebellar hemorrhage. We now report an additional series of 9 patients analyzing the clinical and neuroimaging findings. The clefts were located in the left cerebellar hemisphere in 5 cases, in the right in 3, and bilaterally in one child who had bilateral cerebellar hemorrhages as a preterm infant at 30 weeks gestation. In one patient born at 24 weeks of gestation a unilateral cerebellar hemorrhage has been found at the age of 4 months. Other findings included disordered alignment of the folia and fissures, an irregular gray/white matter junction, and abnormal arborization of the white matter in all cases. Supratentorial abnormalities were found in 4 cases. All but 2 patients were born at term. We confirm the distinct neuroimaging pattern of cerebellar clefts. Considering the documented fetal cerebellar hemorrhage in our first series, we postulate that cerebellar clefts usually represent residual disruptive changes after a prenatal cerebellar hemorrhage. Exceptionally, as now documented in 2 patients, cerebellar clefts can be found after neonatal cerebellar hemorrhages in preterm infants. The short-term outcome in these children was variable.

The murine interleukin-2 receptor gamma chain gene: organization, chromosomal localization and expression in the adult thymus.

Defects in the interleukin-2 receptor gamma (IL-2R gamma) chain in the man result in an X-linked severe combined immunodeficiency, SCIDX1, characterized by an absence of T-cell differentiation. This phenotype may result from pertubations in IL-2, IL-4-, IL-7- or IL-15-mediated signaling, as the IL-2R gamma chain forms an integral component of these receptor systems. We have isolated and characterized cDNA and genomic clones for the murine IL-2R gamma. The gene (Il2rg) is well conserved between mouse and man with respect to overall structure and size, and contains regions of high conservation in the promoter region as well. Il2rg maps to mouse X chromosome region 40, in a region of synteny with human Xq12-13.1. We have also explored the expression of the IL-2R gamma during thymocyte development. IL-2R gamma transcripts are detected in the earliest thymocyte precursor cells and persist throughout intrathymic development into the mature peripheral compartment. Genomic clones for the murine IL-2R gamma will allow for further studies on the regulation and function of this gene in vivo.

Copper acquisition by the SenC protein regulates aerobic respiration in Pseudomonas aeruginosa PAO1.

Aerobic respiration of Pseudomonas aeruginosa involves four terminal oxidases belonging to the heme-copper family (that is, three cytochrome c oxidases and one quinol oxidase) plus one copper-independent, cyanide-insensitive quinol oxidase (CIO). The PA0114 gene encoding an SCO1/SenC-type protein, which is known to be important for copper delivery to cytochrome c in yeast, Rhodobacter spp. and Agrobacterium tumefaciens, was found to be important for copper acquisition and aerobic respiration in P. aeruginosa. A PA0114 (senC) mutant grew poorly in low-copper media and had low cytochrome cbb(3)-type oxidase activity, but expressed CIO at increased levels, by comparison with the wild-type PAO1. Addition of copper reversed these phenotypes, suggesting that periplasmic copper capture by the SenC protein helps P. aeruginosa to adapt to copper deprivation.

Postglacial recolonization at a snail's pace (Trochulus villosus): confronting competing refugia hypotheses using model selection.

The localization of Last Glacial Maximum (LGM) refugia is crucial information to understand a species' history and predict its reaction to future climate changes. However, many phylogeographical studies often lack sampling designs intensive enough to precisely localize these refugia. The hairy land snail Trochulus villosus has a small range centred on Switzerland, which could be intensively covered by sampling 455 individuals from 52 populations. Based on mitochondrial DNA sequences (COI and 16S), we identified two divergent lineages with distinct geographical distributions. Bayesian skyline plots suggested that both lineages expanded at the end of the LGM. To find where the origin populations were located, we applied the principles of ancestral character reconstruction and identified a candidate refugium for each mtDNA lineage: the French Jura and Central Switzerland, both ice-free during the LGM. Additional refugia, however, could not be excluded, as suggested by the microsatellite analysis of a population subset. Modelling the LGM niche of T. villosus, we showed that suitable climatic conditions were expected in the inferred refugia, but potentially also in the nunataks of the alpine ice shield. In a model selection approach, we compared several alternative recolonization scenarios by estimating the Akaike information criterion for their respective maximum-likelihood migration rates. The 'two refugia' scenario received by far the best support given the distribution of genetic diversity in T. villosus populations. Provided that fine-scale sampling designs and various analytical approaches are combined, it is possible to refine our necessary understanding of species responses to environmental changes.

OpenFluDB, a database for human and animal influenza virus.

Although research on influenza lasted for more than 100 years, it is still one of the most prominent diseases causing half a million human deaths every year. With the recent observation of new highly pathogenic H5N1 and H7N7 strains, and the appearance of the influenza pandemic caused by the H1N1 swine-like lineage, a collaborative effort to share observations on the evolution of this virus in both animals and humans has been established. The OpenFlu database (OpenFluDB) is a part of this collaborative effort. It contains genomic and protein sequences, as well as epidemiological data from more than 27,000 isolates. The isolate annotations include virus type, host, geographical location and experimentally tested antiviral resistance. Putative enhanced pathogenicity as well as human adaptation propensity are computed from protein sequences. Each virus isolate can be associated with the laboratories that collected, sequenced and submitted it. Several analysis tools including multiple sequence alignment, phylogenetic analysis and sequence similarity maps enable rapid and efficient mining. The contents of OpenFluDB are supplied by direct user submission, as well as by a daily automatic procedure importing data from public repositories. Additionally, a simple mechanism facilitates the export of OpenFluDB records to GenBank. This resource has been successfully used to rapidly and widely distribute the sequences collected during the recent human swine flu outbreak and also as an exchange platform during the vaccine selection procedure. Database URL: http://openflu.vital-it.ch.