An enhanced workflow for variant interpretation in UniProtKB/Swiss-Prot improves consistency and reuse in ClinVar.

Personalized genomic medicine depends on integrated analyses that combine genetic and phenotypic data from individual patients with reference knowledge of the functional and clinical significance of sequence variants. Sources of this reference knowledge include the ClinVar repository of human genetic variants, a community resource that accepts submissions from external groups, and UniProtKB/Swiss-Prot, an expert-curated resource of protein sequences and functional annotation. UniProtKB/Swiss-Prot provides knowledge on the functional impact and clinical significance of over 30 000 human protein-coding sequence variants, curated from peer-reviewed literature reports. Here we present a pilot study that lays the groundwork for the integration of curated knowledge of protein sequence variation from UniProtKB/Swiss-Prot with ClinVar. We show that existing interpretations of variant pathogenicity in UniProtKB/Swiss-Prot and ClinVar are highly concordant, with 88% of variants that are common to the two resources having interpretations of clinical significance that agree. Re-curation of a subset of UniProtKB/Swiss-Prot variants according to American College of Medical Genetics and Genomics (ACMG) guidelines using ClinGen tools further increases this level of agreement, mainly due to the reclassification of supposedly pathogenic variants as benign, based on newly available population frequency data. We have now incorporated ACMG guidelines and ClinGen tools into the UniProt Knowledgebase (UniProtKB) curation workflow and routinely submit variant data from UniProtKB/Swiss-Prot to ClinVar. These efforts will increase the usability and utilization of UniProtKB variant data and will facilitate the continuing (re-)evaluation of clinical variant interpretations as data sets and knowledge evolve.

Soil fungal communities of grasslands are environmentally structured at a regional scale in the Alps.

Studying patterns of species distributions along elevation gradients is frequently used to identify the primary factors that determine the distribution, diversity and assembly of species. However, despite their crucial role in ecosystem functioning, our understanding of the distribution of below-ground fungi is still limited, calling for more comprehensive studies of fungal biogeography along environmental gradients at various scales (from regional to global). Here, we investigated the richness of taxa of soil fungi and their phylogenetic diversity across a wide range of grassland types along a 2800 m elevation gradient at a large number of sites (213), stratified across a region of the Western Swiss Alps (700 km(2)). We used 454 pyrosequencing to obtain fungal sequences that were clustered into operational taxonomic units (OTUs). The OTU diversity-area relationship revealed uneven distribution of fungal taxa across the study area (i.e. not all taxa are everywhere) and fine-scale spatial clustering. Fungal richness and phylogenetic diversity were found to be higher in lower temperatures and higher moisture conditions. Climatic and soil characteristics as well as plant community composition were related to OTU alpha, beta and phylogenetic diversity, with distinct fungal lineages suggesting distinct ecological tolerances. Soil fungi, thus, show lineage-specific biogeographic patterns, even at a regional scale, and follow environmental determinism, mediated by interactions with plants.

Protein interaction databases.

Life depends on the interaction of proteins. The availability of the complete human genome sequence has highlighted the need for a tool to analyse protein interactions and several databases have been compiled for this purpose. These databases document, categorize, and analyze interacting proteins and the cellular functions of the interactions.

Characterization of the human ABCG1 gene: liver X receptor activates an internal promoter that produces a novel transcript encoding an alternative form of the protein.

The human ABCG1 gene encodes a member of the ATP-binding cassette (ABC) superfamily of transporter proteins and is highly induced when macrophages are incubated with oxysterols. Using mRNA from oxysterol-treated human THP-1 cells together with 5'-rapid amplification of cDNA ends and polymerase chain reaction, we identified a novel ABCG1 transcript that encodes a putative protein of 786 residues containing a new amino terminus of 203 amino acids. Characterization of the genomic organization and structure of the human ABCG1 gene demonstrates that: (i) the gene consists of 23 exons spanning 98 kilobase pairs (kb) on chromosome 21q22.3, (ii) the 203 amino acids are encoded on three previously unidentified exons, 8-10, and (iii) a promoter, containing a TATA box and two liver X receptor (LXR) alpha response elements (LXREs), is located upstream of exon 8. Northern analysis using exon-specific probes confirms that oxysterol treatment results in >10-fold induction of ABCG1 transcripts that are derived from either exons 8-23 or exons 5, 7, and 11-23. Electromobility shift assays demonstrate that LXRalpha and retinoid X receptor alpha bind to the two LXREs in intron 7. Cells were transiently transfected with reporter luciferase constructs under the control of either (i) 9 kb of genomic DNA corresponding to intron 7 and part of exon 8 and containing either wild-type or mutant LXREs or (ii) two copies of the wild-type or mutant LXRE. In all cases, the wild-type construct was regulated in an LXR- and oxysterol-dependent manner, and this regulation was attenuated when the LXREs were mutated. In conclusion, the human ABCG1 gene contains multiple promoters, spans more than 98 kb and comprises 23 exons that give rise to alternative transcripts encoding proteins with different amino-terminal sequences. Elucidation of the various roles of different ABCG1 isoforms will be important for our understanding of mammalian cholesterol homeostasis.

Complement facilitates early prion pathogenesis.

New-variant Creutzfeldt-Jakob disease and scrapie are typically initiated by extracerebral exposure to the causative agent, and exhibit early prion replication in lymphoid organs. In mouse scrapie, depletion of B-lymphocytes prevents neuropathogenesis after intraperitoneal inoculation, probably due to impaired lymphotoxin-dependent maturation of follicular dendritic cells (FDCs), which are a major extracerebral prion reservoir. FDCs trap immune complexes with Fc-gamma receptors and C3d/C4b-opsonized antigens with CD21/CD35 complement receptors. We examined whether these mechanisms participate in peripheral prion pathogenesis. Depletion of circulating immunoglobulins or of individual Fc-gamma receptors had no effect on scrapie pathogenesis if B-cell maturation was unaffected. However, mice deficient in C3, C1q, Bf/C2, combinations thereof or complement receptors were partially or fully protected against spongiform encephalopathy upon intraperitoneal exposure to limiting amounts of prions. Splenic accumulation of prion infectivity and PrPSc was delayed, indicating that activation of specific complement components is involved in the initial trapping of prions in lymphoreticular organs early after infection.

Three-dimensional cluster analysis identifies interfaces and functional residue clusters in proteins.

Three-dimensional cluster analysis offers a method for the prediction of functional residue clusters in proteins. This method requires a representative structure and a multiple sequence alignment as input data. Individual residues are represented in terms of regional alignments that reflect both their structural environment and their evolutionary variation, as defined by the alignment of homologous sequences. From the overall (global) and the residue-specific (regional) alignments, we calculate the global and regional similarity matrices, containing scores for all pairwise sequence comparisons in the respective alignments. Comparing the matrices yields two scores for each residue. The regional conservation score (C(R)(x)) defines the conservation of each residue x and its neighbors in 3D space relative to the protein as a whole. The similarity deviation score (S(x)) detects residue clusters with sequence similarities that deviate from the similarities suggested by the full-length sequences. We evaluated 3D cluster analysis on a set of 35 families of proteins with available cocrystal structures, showing small ligand interfaces, nucleic acid interfaces and two types of protein-protein interfaces (transient and stable). We present two examples in detail: fructose-1,6-bisphosphate aldolase and the mitogen-activated protein kinase ERK2. We found that the regional conservation score (C(R)(x)) identifies functional residue clusters better than a scoring scheme that does not take 3D information into account. C(R)(x) is particularly useful for the prediction of poorly conserved, transient protein-protein interfaces. Many of the proteins studied contained residue clusters with elevated similarity deviation scores. These residue clusters correlate with specificity-conferring regions: 3D cluster analysis therefore represents an easily applied method for the prediction of functionally relevant spatial clusters of residues in proteins.

Mining literature for protein-protein interactions.

MOTIVATION: A central problem in bioinformatics is how to capture information from the vast current scientific literature in a form suitable for analysis by computer. We address the special case of information on protein-protein interactions, and show that the frequencies of words in Medline abstracts can be used to determine whether or not a given paper discusses protein-protein interactions. For those papers determined to discuss this topic, the relevant information can be captured for the Database of Interacting PROTEINS: Furthermore, suitable gene annotations can also be captured. RESULTS: Our Bayesian approach scores Medline abstracts for probability of discussing the topic of interest according to the frequencies of discriminating words found in the abstract. More than 80 discriminating words (e.g. complex, interaction, two-hybrid) were determined from a training set of 260 Medline abstracts corresponding to previously validated entries in the Database of Interacting Proteins. Using these words and a log likelihood scoring function, approximately 2000 Medline abstracts were identified as describing interactions between yeast proteins. This approach now forms the basis for the rapid expansion of the Database of Interacting Proteins.

DIP: The Database of Interacting Proteins: 2001 update.

The Database of Interacting Proteins (DIP; http://dip.doe-mbi.ucla. edu) is a database that documents experimentally determined protein-protein interactions. Since January 2000 the number of protein-protein interactions in DIP has nearly tripled to 3472 and the number of proteins to 2659. New interactive tools have been developed to aid in the visualization, navigation and study of networks of protein interactions.

Localizing proteins in the cell from their phylogenetic profiles.

We introduce a computational method for identifying subcellular locations of proteins from the phylogenetic distribution of the homologs of organellar proteins. This method is based on the observation that proteins localized to a given organelle by experiments tend to share a characteristic phylogenetic distribution of their homologs-a phylogenetic profile. Therefore any other protein can be localized by its phylogenetic profile. Application of this method to mitochondrial proteins reveals that nucleus-encoded proteins previously known to be destined for mitochondria fall into three groups: prokaryote-derived, eukaryote-derived, and organism-specific (i.e., found only in the organism under study). Prokaryote-derived mitochondrial proteins can be identified effectively by their phylogenetic profiles. In the yeast Saccharomyces cerevisiae, 361 nucleus-encoded mitochondrial proteins can be identified at 50% accuracy with 58% coverage. From these values and the proportion of conserved mitochondrial genes, it can be inferred that approximately 630 genes, or 10% of the nuclear genome, is devoted to mitochondrial function. In the worm Caenorhabditis elegans, we estimate that there are approximately 660 nucleus-encoded mitochondrial genes, or 4% of its genome, with approximately 400 of these genes contributed from the prokaryotic mitochondrial ancestor. The large fraction of organism-specific and eukaryote-derived genes suggests that mitochondria perform specialized roles absent from prokaryotic mitochondrial ancestors. We observe measurably distinct phylogenetic profiles among proteins from different subcellular compartments, allowing the general use of prokaryotic genomes in learning features of eukaryotic proteins.

Protein function in the post-genomic era.

Faced with the avalanche of genomic sequences and data on messenger RNA expression, biological scientists are confronting a frightening prospect: piles of information but only flakes of knowledge. How can the thousands of sequences being determined and deposited, and the thousands of expression profiles being generated by the new array methods, be synthesized into useful knowledge? What form will this knowledge take? These are questions being addressed by scientists in the field known as 'functional genomics'.

Human ERK1 induces filamentous growth and cell wall remodeling pathways in Saccharomyces cerevisiae.

Expression of an activated extracellular signal-regulated kinase 1 (ERK1) construct in yeast cells was used to examine the conservation of function among mitogen-activated protein (MAP) kinases. Sequence alignment of the human MAP kinase ERK1 with all Saccharomyces cerevisiae kinases reveals a particularly strong kinship with Kss1p (invasive growth promoting MAP kinase), Fus3p (pheromone response MAP/ERK kinase), and Mpk1p (cell wall remodeling MAP kinase). A fusion protein of constitutively active human MAP/ERK kinase 1 (MEK) and human ERK1 was introduced under regulated expression into yeast cells. The fusion protein (MEK/ERK) induced a filamentation response element promoter and led to a growth retardation effect concomitant with a morphological change resulting in elongated cells, bipolar budding, and multicell chains. Induction of filamentous growth was also observed for diploid cells following MEK/ERK expression in liquid culture. Neither haploids nor diploids, however, showed marked penetration of agar medium. These effects could be triggered by either moderate MEK/ERK expression at 37 degrees C or by high level MEK/ERK expression at 30 degrees C. The combination of high level MEK/ERK expression and 37 degrees C resulted in cell death. The deleterious effects of MEK/ERK expression and high temperature were significantly mitigated by 1 m sorbitol, which also enhanced the filamentous phenotype. MEK/ERK was able to constitutively activate a cell wall maintenance reporter gene, suggesting misregulation of this pathway. In contrast, MEK/ERK effectively blocked expression from a pheromone-responsive element promoter and inhibited mating. These results are consistent with MEK/ERK promoting filamentous growth and altering the cell wall through its ability to partially mimic Kss1p and stimulate a pathway normally controlled by Mpk1p, while appearing to inhibit the normal functioning of the structurally related yeast MAP kinase Fus3p.

DIP: the database of interacting proteins.

The Database of Interacting Proteins (DIP; http://dip.doe-mbi.ucla.edu) is a database that documents experimentally determined protein-protein interactions. This database is intended to provide the scientific community with a comprehensive and integrated tool for browsing and efficiently extracting information about protein interactions and interaction networks in biological processes. Beyond cataloging details of protein-protein interactions, the DIP is useful for understanding protein function and protein-protein relationships, studying the properties of networks of interacting proteins, benchmarking predictions of protein-protein interactions, and studying the evolution of protein-protein interactions.

Interplays between mouse mammary tumor virus and the cellular and humoral immune response.

Mouse mammary tumor virus has developed strategies to exploit the immune response. It requires vigorous immune stimulation to achieve efficient infection. The infected antigen-presenting cells present a viral superantigen on the cell surface which stimulates strong CD4-mediated T-cell help but CD8 T-cell responses are undetectable. Despite the high frequency of superantigen-reactive T cells, the superantigen-induced immune response is comparable to classical antigen responses in terms of T-cell priming, T-cell-B-cell collaboration as well as follicular and extra-follicular B-cell differentiation. Induction of systemic anergy is observed, similar to classical antigen responses where antigen is administered systemically but does not influence the role of the superantigen-reactive T cells in the maintenance of the chronic germinal center reaction. So far we have been unable to detect a cytotoxic T-cell response to mouse mammary tumor virus peptide antigens or to the superantigen. This might yet represent another step in the viral infection strategy.

Differential reactivity of TCR Vbeta10 alleles to a mouse mammary tumor virus superantigen.

Mouse mammary tumor virus (MMTV) expresses a superantigen (SAg) which plays a critical role in the viral life cycle. We have recently described the new infectious MMTV (SIM) encoding a Vbeta4-specific SAg in mice with a TCR-Vbeta(b) haplotype. We have now compared the SAg activity of this virus in BALB/c mice harboring the TCR-Vbeta(a), TCR-Vbeta(b) or TCR-Vbeta(c) haplotypes which differ by a central deletion in the TCR-Vbeta(a) and TCR-Vbeta(c) locus and by mutations in some of the remaining Vbeta elements. Injection of MMTV (SIM) led to a strong stimulation of Vbeta4+ CD4+ T cells in TCR-Vbeta(b) mice, but only to a weak stimulation of these cells in TCR-Vbeta(a) or TCR-Vbeta(c) mice. A large increase in the percentage of Vbeta10+ cells was observed among CD4+ T cells in mice with the Vbeta(a) or Vbeta(c), but not the Vbeta(b) TCR-Vbeta haplotype. Vbeta10+ cells dominated the response when Vbeta10(a/c) and Vbeta4 subsets were present together. This is the first report of a viral SAg interacting with murine Vbeta10+ cells. Six amino acid differences between Vbeta10(a/c) and Vbeta10(b) could account for the gain of reactivity of Vbeta10(a/c) to the MMTV(SIM) SAg. No mutations were found in the hypervariable region 4 (HV4) of the TCR. Mutations at positions 22 and 28 introduce into Vbeta10(a/c) the same amino acids which are found at these positions in the MMTV(SIM)-reactive Vbeta4. Tridimensional models indicated that these amino acids lie close to HV4 and are likely to be important for the interaction of the SAg with the TCR.

Immune response to mouse mammary tumor virus in mice lacking the alpha/beta interferon or the gamma interferon receptor.

Mouse mammary tumor virus (MMTV) is a retrovirus which induces a strong immune response and a dramatic increase in the number of infected cells through the expression of a superantigen (SAg). Many cytokines are likely to be involved in the interaction between MMTV and the immune system. In particular, alpha/beta interferon (IFN-alpha/beta) and gamma interferon (IFN-gamma) exert many antiviral and immunomodulatory activities and play a critical role in other viral infections. In this study, we have investigated the importance of interferons during MMTV infection by using mice with a disrupted IFN-alpha/beta or IFN-gamma receptor gene. We found that the SAg response to MMTV was not modified in IFN-alpha/betaR(0/0) and IFN-gammaR(0/0) mice. This was true both for the early expansion of B and T cells induced by the SAg and for the deletion of SAg-reactive cells at later stages of the infection. In addition, no increase in the amount of proviral DNA was detected in tissues of IFN-alpha/betaR(0/0) and IFN-gammaR(0/0) mice, suggesting that interferons are not essential antiviral defense mechanisms during MMTV infection. In contrast, IFN-gammaR(0/0) mice had increased amounts of IL-4 mRNA and an altered usage of immunoglobulin isotypes with a reduced frequency of IgG2a- and IgG3-producing cells. This was associated with lower titers of virus-specific antibodies in serum early after infection, although efficient titers were reached later.

IL-4 rapidly produced by V beta 4 V alpha 8 CD4+ T cells instructs Th2 development and susceptibility to Leishmania major in BALB/c mice.

BALB/c mice develop aberrant T helper 2 (Th2) responses and suffer progressive disease after infection with Leishmania major. These outcomes depend on the production of interleukin-4 (IL-4) early after infection. Here we demonstrate that the burst of IL-4 mRNA, peaking in draining lymph nodes of BALB/c mice 16 hr after infection, occurs within CD4+ T cells that express V beta 4 V alpha 8 T cell receptors. In contrast to control and V beta 6-deficient BALB/c mice, V beta 4-deficient BALB/c mice were resistant to infection, demonstrating the role of these cells in Th2 development. The early IL-4 response was absent in these mice, and T helper 1 responses occurred following infection. Recombinant LACK antigen from L. major induced comparable IL-4 production in V beta 4 V alpha 8 CD4+ cells. Thus, the IL-4 required for Th2 development and susceptibility to L. major is produced by a restricted population of V beta 4 V alpha 8 CD4+ T cells after cognate interaction with a single antigen from this complex organism.

New infectious mammary tumor virus superantigen with V beta-specificity identical to staphylococcal enterotoxin B (SEB).

Only few infectious mouse mammary tumor viruses (MMTV) have been characterized which induce a potent superantigen response in vivo. Here we describe the characterization of an MMTV which was isolated from milk of the highly mammary tumor-prone SHN mouse strain. Exposure of newborn mice to milk-borne MMTV (SHN) results in a very slow deletion of V beta 7, 8.1, 8.2 and 8.3 expressing peripheral T cells. Subcutaneous injection of adult mice with this virus induces a rapid and strong stimulation of all four affected V beta-subsets in vivo. Besides the strong T cell effect we observed an early proliferation and activation of the local B cell pool leading to the initial secretion of IgM followed by preferential secretion of IgG2a by day 6. Sequence comparison of the polymorphic C terminus with known open reading frames revealed high homology to the endogenous provirus Mtv-RCS. This is the first report of a virus having a complete overlap in V beta-specificity with a bacterial superantigen stimulating as many as 35% of the whole CD4+ T cell repertoire including V beta 8.2.

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.