VarSifter: visualizing and analyzing exome-scale sequence variation data on a desktop computer

VarSifter is a graphical software tool for desktop computers that allows investigators of varying computational skills to easily and quickly sort, filter, and sift through sequence variation data. A variety of filters and a custom query framework allow filtering based on any combination of sample and annotation information. By simplifying visualization and analyses of exome-scale sequence variation data, this program will help bring the power and promise of massively-parallel DNA sequencing to a broader group of researchers. Availability and Implementation: VarSifter is written in Java, and is freely available in source and binary versions, along with a User Guide, at http://research.nhgri.nih.gov/software/VarSifter/.     

MatrixPlot: visualizing sequence constraints.

SUMMARY: MatrixPlot is a program for making high-quality matrix plots, such as mutual information plots of sequence alignments and distance matrices of sequences with known three-dimensional coordinates. The user can add information about the sequences (e.g. a sequence logo profile) along the edges of the plot, as well as zoom in on any region in the plot. AVAILABILITY: MatrixPlot can be obtained on request, and can also be accessed online at http://www. cbs.dtu.dk/services/MatrixPlot. CONTACT: gorodkin@cbs.dtu.dk     

Optimizing copy number variation analysis using genome-wide short sequence oligonucleotide arrays

The detection of copy number variants (CNV) by array-based platforms provides valuable insight into understanding human diversity. However, suboptimal study design and data processing negatively affect CNV assessment. We quantitatively evaluate their impact when short-sequence oligonucleotide arrays are applied (Affymetrix Genome-Wide Human SNP Array 6.0) by evaluating 42 HapMap samples for CNV detection. Several processing and segmentation strategies are implemented, and results are compared to CNV assessment obtained using an oligonucleotide array CGH platform designed to query CNVs at high resolution (Agilent). We quantitatively demonstrate that different reference models (e.g. single versus pooled sample reference) used to detect CNVs are a major source of inter-platform discrepancy (up to 30%) and that CNVs residing within segmental duplication regions (higher reference copy number) are significantly harder to detect (P < 0.0001). After adjusting Affymetrix data to mimic the Agilent experimental design (reference sample effect), we applied several common segmentation approaches and evaluated differential sensitivity and specificity for CNV detection, ranging 39–77% and 86–100% for non-segmental duplication regions, respectively, and 18–55% and 39–77% for segmental duplications. Our results are relevant to any array-based CNV study and provide guidelines to optimize performance based on study-specific objectives.     

VISTA : visualizing global DNA sequence alignments of arbitrary length

Summary: VISTA is a program for visualizing global DNA sequence alignments of arbitrary length. It has a clean output, allowing for easy identification of similarity, and is easily configurable, enabling the visualization of alignments of various lengths at different levels of resolution. It is currently available on the web, thus allowing for easy access by all researchers. Availability: VISTA server is available on the web at http://www-gsd.lbl.gov/vista. The source code is available upon request. Contact: vista@lbl.gov     

Interdependent utilities: how social ranking affects choice behavior

Organization in hierarchical dominance structures is prevalent in animal societies, so a strong preference for higher positions in social ranking is likely to be an important motivation of human social and economic behavior. This preference is also likely to influence the way in which we evaluate our outcome and the outcome of others, and finally the way we choose. In our experiment participants choose among lotteries with different levels of risk, and can observe the choice that others have made. Results show that the relative weight of gains and losses is the opposite in the private and social domain. For private outcomes, experience and anticipation of losses loom larger than gains, whereas in the social domain, gains loom larger than losses, as indexed by subjective emotional evaluations and physiological responses. We propose a theoretical model (interdependent utilities), predicting the implication of this effect for choice behavior. The relatively larger weight assigned to social gains strongly affects choices, inducing complementary behavior: faced with a weaker competitor, participants adopt a more risky and dominant behavior.     

A designed system for assessing how sequence affects alpha to beta conformational transitions in proteins.

The role of amino acid sequence in conformational switching observed in prions and proteins associated with amyloid diseases is not well understood. To study alpha to beta conformational transitions, we designed a series of peptides with structural duality; namely, peptides with sequence features of both an alpha-helical leucine zipper and a beta-hairpin. The parent peptide, Template-alpha, was designed to be a canonical leucine-zipper motif and was confirmed as such using circular dichroism spectroscopy and analytical ultracentrifugation. To introduce beta-structure character into the peptide, glutamine residues at sites away from the leucine-zipper dimer interface were replaced by threonine to give Template-alphaT. Unlike the parent peptide, Template-alphaT underwent a heat-inducible switch to beta-structure, which reversibly formed gels containing amyloid-like fibrils. In contrast to certain other natural proteins where destabilization of the native states facilitate transitions to amyloid, destabilization of the leucine-zipper form of Template-alphaT did not promote a transformation. Cross-linking the termini of the peptides compatible with the alternative beta-hairpin design, however, did promote the change. Furthermore, despite screening various conditions, only the internally cross-linked form of the parent, Template-alpha, peptide formed amyloid-like fibrils. These findings demonstrate that, in addition to general properties of the polypeptide backbone, specific residue placements that favor beta-structure promote amyloid formation.     

Models of deletion for visualizing bacterial variation: an application to tuberculosis spoligotypes

Background  

Molecular typing methods are commonly used to study genetic relationships among bacterial isolates. Many of these methods have become standardized and produce portable data. A popular approach for analyzing such data is to construct graphs, including phylogenies. Inferences from graph representations of data assist in understanding the patterns of transmission of bacterial pathogens, and basing these graph constructs on biological models of evolution of the molecular marker helps make these inferences. Spoligotyping is a widely used method for genotyping isolates of Mycobacterium tuberculosis that exploits polymorphism in the direct repeat region. Our goal was to examine a range of models describing the evolution of spoligotypes in order to develop a visualization method to represent likely relationships among M. tuberculosis isolates.     

3MOTIF: visualizing conserved protein sequence motifs in the protein structure database

SUMMARY: 3MOTIF is a web application that visually maps conserved sequence motifs onto three-dimensional protein structures in the Protein Data Bank (PDB; Berman et al., Nucleic Acids Res., 28, 235-242, 2000). Important properties of motifs such as conservation strength and solvent accessible surface area at each position are visually represented on the structure using a variety of color shading schemes. Users can manipulate the displayed motifs using the freely available Chime plugin. AVAILABILITY: http://motif.stanford.edu/3motif/     

CONAN: copy number variation analysis software for genome-wide association studies

Background  

Genome-wide association studies (GWAS) based on single nucleotide polymorphisms (SNPs) revolutionized our perception of the genetic regulation of complex traits and diseases. Copy number variations (CNVs) promise to shed additional light on the genetic basis of monogenic as well as complex diseases and phenotypes. Indeed, the number of detected associations between CNVs and certain phenotypes are constantly increasing. However, while several software packages support the determination of CNVs from SNP chip data, the downstream statistical inference of CNV-phenotype associations is still subject to complicated and inefficient in-house solutions, thus strongly limiting the performance of GWAS based on CNVs.     

Iron homeostasis in plants： when transcription affects translocation

Despite the fact that iron is one of the most abundant elements of the earth＇s crust, iron deficiencies are serious problems both in human nutrition [ 1 ] and in agriculture [2]. Six to eight percent of the world＇s population is potentially affected by iron deficiency induced anemia, a leading cause of maternal death in African and Asian countries where people rely mostly on plants for their daily intake of iron. Iron can also be a limiting factor in the growth of economically important crop plants because of inadequate soil chemistry, and such deficiencies cannot easily be corrected by amending the soil. Improving the plant＇s ability to absorb iron in adverse conditions and to increase their overall content could offer solutions to these dramatic problems. Therefore understanding the molecular mechanisms regulating iron uptake and homeostasis in plants has potentially important practical applications both in agriculture and human health [3].     

Neisserial surface variation: how and why?

Neisseria gonorrhoeae exhibits striking variability in several of its surface components (pili, Opa proteins and lipooligosaccharide) in vivo and in vitro. Such flagrant variation of this mucosal pathogen's surface components contrasts sharply with changes in single surface components of blood-borne trypanosomes and borreliae. Despite these differences, similar molecular events are sometimes involved.     

Variations in abundance: genome-wide responses to genetic variation and vice versa

How do naturally occurring polymorphisms in DNA sequence relate to variation in gene expression? Recent work to map genetic sources of expression variation has shown a surprising balance between cis and trans effects. Other work suggests some chromosomal clustering of genes by expression pattern. A synthesis of approaches may provide new insight in to adaptive mechanisms in evolution and the population basis of complex traits.