Exome sequencing in 32 patients with anophthalmia/microphthalmia and developmental eye defects.

Anophthalmia/microphthalmia (A/M) is a genetically heterogeneous birth defect for which the etiology is unknown in more than 50% of patients. We used exome sequencing with the ACE Exome(TM) (Personalis, Inc; 18 cases) and UCSF Genomics Core (21 cases) to sequence 28 patients with A/M and four patients with varied developmental eye defects. In the 28 patients with A/M, we identified de novo mutations in three patients (OTX2, p.(Gln91His), RARB, p.Arg387Cys and GDF6, p.Ala249Glu) and inherited mutations in STRA6 in two patients. In patients with developmental eye defects, a female with cataracts and cardiomyopathy had a de novo COL4A1 mutation, p.(Gly773Arg), expanding the phenotype associated with COL4A1 to include cardiomyopathy. A male with a chorioretinal defect, microcephaly, seizures and sensorineural deafness had two PNPT1 mutations, p.(Ala507Ser) and c.401-1G>A, and we describe eye defects associated with this gene for the first time. Exome sequencing was efficient for identifying mutations in pathogenic genes for which there is no clinical testing available and for identifying cases that expand phenotypic spectra, such as the PNPT1 and COL4A1-associated disorders described here.

The direct cellular target of topically applied pimecrolimus may not be infiltrating lymphocytes.

BACKGROUND: Topically applied calcineurin inhibitors have been shown to be effective in the treatment of atopic dermatitis. When systemically administered, these agents cause immunosuppression via inhibition of calcineurin in lymphocytes. As topical agents, the mechanism of action is poorly defined. OBJECTIVES: To test the hypothesis that skin-infiltrating lymphocytes are directly targeted when calcineurin inhibitors are applied to the skin. METHODS: Ten patients with atopic dermatitis were treated with 1% pimecrolimus cream twice daily to target lesions. Skin biopsies were performed before and 48 h after beginning therapy. We assessed the cellular localization of NFAT1 and NFAT2 as a surrogate measure of intracellular calcineurin activity (e.g. increasing cytoplasmic localization with increasing calcineurin inhibition). RESULTS: All patients showed a clinical response, at both 48 h and 2 weeks. As previously described, NFAT2 localized to the follicular keratinocytes, and its activation was partially inhibited by topical pimecrolimus. NFAT1 was found to be expressed by follicular and interfollicular keratinocytes, and its mostly nuclear localization was not affected by topical pimecrolimus therapy. Both NFAT1 and NFAT2 were found in the infiltrating lymphocytes. However, using both manual counting as well as an automated method to assess nuclear intensity of NFAT staining, we found that the proportion of infiltrating leucocytes with nuclear ('activated') NFAT did not change following therapy with pimecrolimus. CONCLUSIONS: Our results suggest that topical pimecrolimus does not act primarily by inhibiting the calcineurin/NFAT axis in lymphocytes but may instead act by other mechanisms, possibly by decreasing NFAT2 activity in follicular keratinocytes.

Automated diagnosis of data-model conflicts using metadata.

The authors describe a methodology for helping computational biologists diagnose discrepancies they encounter between experimental data and the predictions of scientific models. The authors call these discrepancies data-model conflicts. They have built a prototype system to help scientists resolve these conflicts in a more systematic, evidence-based manner. In computational biology, data-model conflicts are the result of complex computations in which data and models are transformed and evaluated. Increasingly, the data, models, and tools employed in these computations come from diverse and distributed resources, contributing to a widening gap between the scientist and the original context in which these resources were produced. This contextual rift can contribute to the misuse of scientific data or tools and amplifies the problem of diagnosing data-model conflicts. The authors' hypothesis is that systematic collection of metadata about a computational process can help bridge the contextual rift and provide information for supporting automated diagnosis of these conflicts. The methodology involves three major steps. First, the authors decompose the data-model evaluation process into abstract functional components. Next, they use this process decomposition to enumerate the possible causes of the data-model conflict and direct the acquisition of diagnostically relevant metadata. Finally, they use evidence statically and dynamically generated from the metadata collected to identify the most likely causes of the given conflict. They describe how these methods are implemented in a knowledge-based system called GRENDEL and show how GRENDEL can be used to help diagnose conflicts between experimental data and computationally built structural models of the 30S ribosomal subunit.

RNA secondary structure as a reusable interface to biological information resources.

The dissemination of biological information has become critically dependent on the Internet and World Wide Web (WWW), which enable distributed access to information in a platform independent manner. The mode of interaction between biologists and on-line information resources, however, has been mostly limited to simple interface technologies such has hypertext links, tables and forms. The introduction of platform-independent runtime environments facilitates the development of more sophisticated WWW-based user interfaces. Until recently, most such interfaces have been tightly coupled to the underlying computation engines, and not separated as reusable components. We believe that many subdisciplines of biology have intuitive and familiar graphical representations of knowledge that can serve as multipurpose user interface elements. We call such graphical idioms "domain graphics". In order to illustrate the power of such graphics, we have built a reusable interface based on the standard two dimensional (2D) layout of RNA secondary structure. The interface can be used to represent any pre-computed layout of RNA, and takes as a parameters the sets of actions to be performed as a user interacts with the interface. It can provide to any associated application program information about the base, helix, or subsequence selected by the user. We show the versatility of this interface by using it as a special purpose interface to BLAST, Medline and the RNA MFOLD search/compute engines. These demonstrations are available at: http://www-smi.stanford.edu/projects/helix/pubs/ gene-combis-96/

Standardized representations of the literature: combining diverse sources of ribosomal data.

We are building a knowledge base (KB) of published structural data on the 30s ribosomal subunit in prokaryotes. Our KB is distinguished by a standardized representation of biological experiments and their results, in a reusable format. It can be accessed by computer programs that exploit the rich interconnections within the data. The KB is designed to support the construction of 3D models of the 30S subunit, as well as the analysis and extension of relevant functional and phylogenetic information. Most published information about the structure of the ubiquitous ribosome focuses on E. coli as a model system. At the same time, thousands of RNA sequences for the ribosome have been gathered and cataloged. The volume and complexity of these data can complicate attempts to separate structural data peculiar to E. coli from data of universal relevance. We have written an application that dynamically queries the KB and the Ribosome Database Project, a repository of ribosomal RNA sequences from other organisms, in order to assess the relevance of structural data to particular organisms. The application uses the RDP alignment to determine whether a set of data refer primarily to conserved, mismatched, or gapped positions. For a set of 16 representative articles evaluated over 211 sequences, 73% of observations have unambiguous translations from E. coli to the other organisms, 21% have somewhat ambiguous translations, and 6% have no translations. There is a wide variation in these numbers over different articles and organisms, confirming that some articles report structural information specific to E. coli while others report information that is quite general.

RIBOWEB: linking structural computations to a knowledge base of published experimental data.

The world wide web (WWW) has become critical for storing and disseminating biological data. It offers an additional opportunity, however, to support distributed computation and sharing of results. Currently, computational analysis tools are often separated from the data in a manner that makes iterative hypothesis testing cumbersome. We hypothesize that the cycle of scientific reasoning (using data to build models, and evaluating models in light of data) can be facilitated with resources that link computations with semantic models of the data. Riboweb is an on-line knowledge-based resource that supports the creation of three-dimensional models of the 30S ribosomal subunit. It has three components: (I) a knowledge base containing representations of the essential physical components and published structural data, (II) computational modules that use the knowledge base to build or analyze structural models, and (III) a web-based user interface that supports multiple users, sessions and computations. We have built a prototype of Riboweb, and have used it to refine a rough model of the central domain of the 30S subunit from E. coli. procedure. Our results suggest that sophisticated and integrated computational capabilities can be delivered to biologists using this simple three-component architecture.

Computational methods for defining the allowed conformational space of 16S rRNA based on chemical footprinting data.

Structural models for 16S ribosomal RNA have been proposed based on combinations of crosslinking, chemical protection, shape, and phylogenetic evidence. These models have been based for the most part on independent data sets and different sets of modeling assumptions. In order to evaluate such models meaningfully, methods are required to explicitly model the spatial certainty with which individual structural components are positioned by specific data sets. In this report, we use a constraint satisfaction algorithm to explicitly assess the location of the secondary structural elements of the 16S RNA, as well as the certainty with which these elements can be positioned. The algorithm initially assumes that these helical elements can occupy any position and orientation and then systematically eliminates those positions and orientations that do not satisfy formally parameterized interpretations of structural constraints. Using a conservative interpretation of the hydroxyl radical footprinting data, the positions of the ribosomal proteins as defined by neutron diffraction studies, and the secondary structure of 16S rRNA, the location of the RNA secondary structural elements can be defined with an average precision of 25 A (ranging from 12.8 to 56.3 A). The uncertainty in individual helix positions is both heterogeneous and dependent upon the number of constraints imposed on the helix. The topology of the resulting model is consistent with previous models based on independent approaches. The result of our computation is a conservative upper bound on the possible positions of the RNA secondary structural elements allowed by this data set, and provides a suitable starting point for refinement with other sources of data or different sets of modeling assumptions.

Constraining volume by matching the moments of a distance distribution.

The problem of computing a molecular structure from a set of distances arises in the interpretation of NMR data as well as other experimental methods that yield distance information. Techniques for computing structures must find conformations consistent with the distance data. There are often other constraints on the structure that must be satisfied as well. One of the most problematic constraints is the constraint on the total volume occupied by the atoms. In this paper, we use the first two moments (mean and variance) of an estimated distance distribution to constrain the volume of a computed structure. We show that a probabilistic algorithm for matching the first two moments of the estimated distance distribution significantly improves the quality of the solution, especially when the distance information alone is not sufficient to define the structure precisely. We also show that our method is not sensitive to small errors in the estimates of mean and variance of the distance distribution. Finally, we demonstrate the use of this constraint in computing a low-resolution structure of the 30S prokaryotic ribosomal subunit. Quantitative analysis of our results allows us to assess the information content contained in constraints on volume, and to show that in some cases addition of a volume constraint adds information roughly equivalent to doubling the number of input distances. Our results also demonstrate the flexibility of probabilistic representations of structural constraints, and the importance of including volume information to constrain structural computations-especially in the case of sparse data.

Three-dimensional transformation of capsids associated with genome packaging in a bacterial virus.

Three-dimensional structures of the empty procapsid and the mature capsid of the Salmonella bacteriophage P22 have been determined to a resolution of 28 A using electron cryomicroscopy and computer image processing. The coat subunits in both the structures are arranged as pentamers and hexamers on a T = 7 icosahedral lattice. The two structures display significant differences in shape, size and intersubunit interactions. The empty procapsid is spherical in contrast to the distinctly larger and polyhedral mature capsid. The empty procapsid structure exhibits holes at all the quasi sixfold positions that are absent in the mature capsid. These holes may be the exit ports for scaffolding subunits. Detailed comparisons of the two structures indicate that extensive structural changes take place during maturation in all seven quasi-equivalent subunits. These changes cause flattening of the icosahedral facets, capsid expansion and closing of the holes. This process results in a stable and impenetrable capsid that protects the bacterial genome.