The Structure-Function Linkage Database

The Structure–Function Linkage Database (SFLD, http://sfld.rbvi.ucsf.edu/) is a manually curated classification resource describing structure–function relationships for functionally diverse enzyme superfamilies. Members of such superfamilies are diverse in their overall reactions yet share a common ancestor and some conserved active site features associated with conserved functional attributes such as a partial reaction. Thus, despite their different functions, members of these superfamilies ‘look alike’, making them easy to misannotate. To address this complexity and enable rational transfer of functional features to unknowns only for those members for which we have sufficient functional information, we subdivide superfamily members into subgroups using sequence information, and lastly into families, sets of enzymes known to catalyze the same reaction using the same mechanistic strategy. Browsing and searching options in the SFLD provide access to all of these levels. The SFLD offers manually curated as well as automatically classified superfamily sets, both accompanied by search and download options for all hierarchical levels. Additional information includes multiple sequence alignments, tab-separated files of functional and other attributes, and sequence similarity networks. The latter provide a new and intuitively powerful way to visualize functional trends mapped to the context of sequence similarity

The International Gene Trap Consortium Website: a portal to all publicly available gene trap cell lines in mouse

Gene trapping is a method of generating murine embryonic stem (ES) cell lines containing insertional mutations in known and novel genes. A number of international groups have used this approach to create sizeable public cell line repositories available to the scientific community for the generation of mutant mouse strains. The major gene trapping groups worldwide have recently joined together to centralize access to all publicly available gene trap lines by developing a user-oriented Website for the International Gene Trap Consortium (IGTC). This collaboration provides an impressive public informatics resource comprising ∼45 000 well-characterized ES cell lines which currently represent ∼40% of known mouse genes, all freely available for the creation of knockout mice on a non-collaborative basis. To standardize annotation and provide high confidence data for gene trap lines, a rigorous identification and annotation pipeline has been developed combining genomic localization and transcript alignment of gene trap sequence tags to identify trapped loci. This information is stored in a new bioinformatics database accessible through the IGTC Website interface. The IGTC Website () allows users to browse and search the database for trapped genes, BLAST sequences against gene trap sequence tags, and view trapped genes within biological pathways. In addition, IGTC data have been integrated into major genome browsers and bioinformatics sites to provide users with outside portals for viewing this data. The development of the IGTC Website marks a major advance by providing the research community with the data and tools necessary to effectively use public gene trap resources for the large-scale characterization of mammalian gene function

Genetic Variation in the Proximal Promoter of ABC and SLC Superfamilies: Liver and Kidney Specific Expression and Promoter Activity Predict Variation

Membrane transporters play crucial roles in the cellular uptake and efflux of an array of small molecules including nutrients, environmental toxins, and many clinically used drugs. We hypothesized that common genetic variation in the proximal promoter regions of transporter genes contribute to observed variation in drug response. A total of 579 polymorphisms were identified in the proximal promoters (−250 to +50 bp) and flanking 5′ sequence of 107 transporters in the ATP Binding Cassette (ABC) and Solute Carrier (SLC) superfamilies in 272 DNA samples from ethnically diverse populations. Many transporter promoters contained multiple common polymorphisms. Using a sliding window analysis, we observed that, on average, nucleotide diversity (π) was lowest at approximately 300 bp upstream of the transcription start site, suggesting that this region may harbor important functional elements. The proximal promoters of transporters that were highly expressed in the liver had greater nucleotide diversity than those that were highly expressed in the kidney consistent with greater negative selective pressure on the promoters of kidney transporters. Twenty-one promoters were evaluated for activity using reporter assays. Greater nucleotide diversity was observed in promoters with strong activity compared to promoters with weak activity, suggesting that weak promoters are under more negative selective pressure than promoters with high activity. Collectively, these results suggest that the proximal promoter region of membrane transporters is rich in variation and that variants in these regions may play a role in interindividual variation in drug disposition and response

Comparison of methods for genomic localization of gene trap sequences

<p>Abstract</p> <p>Background</p> <p>Gene knockouts in a model organism such as mouse provide a valuable resource for the study of basic biology and human disease. Determining which gene has been inactivated by an untargeted gene trapping event poses a challenging annotation problem because gene trap sequence tags, which represent sequence near the vector insertion site of a trapped gene, are typically short and often contain unresolved residues. To understand better the localization of these sequences on the mouse genome, we compared stand-alone versions of the alignment programs BLAT, SSAHA, and MegaBLAST. A set of 3,369 sequence tags was aligned to build 34 of the mouse genome using default parameters for each algorithm. Known genome coordinates for the cognate set of full-length genes (1,659 sequences) were used to evaluate localization results.</p> <p>Results</p> <p>In general, all three programs performed well in terms of localizing sequences to a general region of the genome, with only relatively subtle errors identified for a small proportion of the sequence tags. However, large differences in performance were noted with regard to correctly identifying exon boundaries. BLAT correctly identified the vast majority of exon boundaries, while SSAHA and MegaBLAST missed the majority of exon boundaries. SSAHA consistently reported the fewest false positives and is the fastest algorithm. MegaBLAST was comparable to BLAT in speed, but was the most susceptible to localizing sequence tags incorrectly to pseudogenes.</p> <p>Conclusion</p> <p>The differences in performance for sequence tags and full-length reference sequences were surprisingly small. Characteristic variations in localization results for each program were noted that affect the localization of sequence at exon boundaries, in particular.</p

Comparison of methods for genomic localization of gene trap sequences.

BackgroundGene knockouts in a model organism such as mouse provide a valuable resource for the study of basic biology and human disease. Determining which gene has been inactivated by an untargeted gene trapping event poses a challenging annotation problem because gene trap sequence tags, which represent sequence near the vector insertion site of a trapped gene, are typically short and often contain unresolved residues. To understand better the localization of these sequences on the mouse genome, we compared stand-alone versions of the alignment programs BLAT, SSAHA, and MegaBLAST. A set of 3,369 sequence tags was aligned to build 34 of the mouse genome using default parameters for each algorithm. Known genome coordinates for the cognate set of full-length genes (1,659 sequences) were used to evaluate localization results.ResultsIn general, all three programs performed well in terms of localizing sequences to a general region of the genome, with only relatively subtle errors identified for a small proportion of the sequence tags. However, large differences in performance were noted with regard to correctly identifying exon boundaries. BLAT correctly identified the vast majority of exon boundaries, while SSAHA and MegaBLAST missed the majority of exon boundaries. SSAHA consistently reported the fewest false positives and is the fastest algorithm. MegaBLAST was comparable to BLAT in speed, but was the most susceptible to localizing sequence tags incorrectly to pseudogenes.ConclusionThe differences in performance for sequence tags and full-length reference sequences were surprisingly small. Characteristic variations in localization results for each program were noted that affect the localization of sequence at exon boundaries, in particular

Enhancing Data Sharing in Collaborative Research Projects with DASH T.E. Ferrin, C.C. Huang, D.M. Greenblatt, D. Stryke, K.M. Giacomini, and J.H. Morris Pacific Symposium on Biocomputing 10:260-271(2005) ENHANCING DATA SHARING IN COLLABORATIVE RESEARCH PR

We describe a software framework, called DASH, that enables the facile access, maintenance, curation and sharing of computational biology data among collaborating research scientists. The DASH event-based framework enables members of team-based research projects to describe the multistep computational processing pipelines frequently required to generate data for sharing, monitors multiple distributed data stores for changes, and will then automatically invoke the appropriate processing pipeline(s). These pipelines can be used to communicate the results of data analyses to collaborators using mechanisms such as Web Services. We describe the overall design of the DASH system and the application of a simple DASH prototype to a collaborative pharmacogenomics research project involving several dozen researchers located at several different sites—the UCSF Pharmacogenetics of Membrane Transporters project. 1

Longitudinal comparison of the developing gut virome in infants and their mothers

The virome of the human gut and its development in early life are poorly understood. Here we performed viral metagenomic sequencing on stool samples from a multiethnic, socioeconomically diverse cohort of 53 infants collected longitudinally over their first 3 years of life and their mothers to investigate and compare their viromes. The asymptomatic infant virome consisted of bacteriophages, dietary/environmental viruses, and human pathogenic viruses, in contrast to the material virome, in which sequence reads from human pathogenic viruses were absent or present at extremely low levels. Picornaviruses and phages in the family Microviridae (microviruses) dominated the infant virome, while microviruses and tomato mosaic virus dominated the maternal virome. As the infants aged, the human pathogenic and dietary/environmental virus components remained distinct from the materal virome, while the phage component evolved to become more similar. However, the composition of the evolving infant virome was not determined by the mother and was still maturing to the adult virome at three years of age. Importance The development of the human gut virome in early childhood is poorly understood. Here we use viral metagenomic sequencing in a cohort of 53 infants to the characterize their gut viromes and compare them to their mothers’.. This study finds that the infant virome consists of phages and human pathogenic viruses in asymptomatic individuals and is still maturing into the adult virome at three years of age

Longitudinal comparison of the developing gut virome in infants and their mothers.

The human gut virome and its early life development are poorly understood. Prior studies have captured single-point assessments with the evolution of the infant virome remaining largely unexplored. We performed viral metagenomic sequencing on stool samples collected longitudinally from a cohort of 53 infants from age 2&nbsp;weeks to 3 years (80.7 billion reads), and from their mothers (9.8 billion reads) to examine and compare viromes. The asymptomatic infant virome consisted of bacteriophages, nonhuman dietary/environmental viruses, and human-host viruses, predominantly picornaviruses. In contrast, human-host viruses were largely absent from the maternal virome. Previously undescribed, sequence-divergent vertebrate viruses were detected in the maternal but not infant virome. As infants aged, the phage component evolved to resemble the maternal virome, but by age 3, the human-host component remained dissimilar from the maternal virome. Thus, early life virome development is determined predominantly by dietary, infectious, and environmental factors rather than direct maternal acquisition