The exploration of copy-number variation (CNV), notably of somatic cells, is an understudied aspect of genome biology. Any differences in the genetic makeup between twins derived from the same zygote represent an irrefutable example of somatic mosaicism. We studied 19 pairs of monozygotic twins with either concordant or discordant phenotype by using two platforms for genome-wide CNV analyses and showed that CNVs exist within pairs in both groups. These findings have an impact on our views of genotypic and phenotypic diversity in monozygotic twins and suggest that CNV analysis in phenotypically discordant monozygotic twins may provide a powerful tool for identifying disease-predisposition loci. Our results also imply that caution should be exercised when interpreting disease causality of de novo CNVs found in patients based on analysis of a single tissue in routine disease-related DNA diagnostics.

Two major types of genetic variation are known: single nucleotide polymorphisms (SNPs), and a more recently discovered structural variation, involving changes in copy number (CNVs) of kilobase- to megabase-sized chromosomal segments. It is unknown whether CNVs arise in somatic cells, but it is, however, generally assumed that normal cells are genetically identical. We tested 34 tissue samples from three subjects and, having analyzed for each tissue < or =10(-6) of all cells expected in an adult human, we observed at least six CNVs, affecting a single organ or one or more tissues of the same subject. The CNVs ranged from 82 to 176 kb, often encompassing known genes, potentially affecting gene function. Our results indicate that humans are commonly affected by somatic mosaicism for stochastic CNVs, which occur in a substantial fraction of cells. The majority of described CNVs were previously shown to be polymorphic between unrelated subjects, suggesting that some CNVs previously reported as germline might represent somatic events, since in most studies of this kind, only one tissue is typically examined and analysis of parents for the studied subjects is not routinely performed. A considerable number of human phenotypes are a consequence of a somatic process. Thus, our conclusions will be important for the delineation of genetic factors behind these phenotypes. Consequently, biobanks should consider sampling multiple tissues to better address mosaicism in the studies of somatic disorders.

We have simulated an odor ligand's dynamic behavior in the binding region of an olfactory receptor (OR). Our short timescale computational studies (up to 200 ps) have helped identify unprecedented postdocking ligand behavior of ligands. From in vacuo molecular dynamics simulations of interactions between models of rat OR I7 and 10 aldehyde ligands, we have identified a dissociative pathway along which the ligand exits and enters the OR-binding pocket--a transit event. The ligand's transit through the receptor's binding region may mark the beginning of a signal transduction cascade leading to odor recognition. We have graphically traced the rotameric changes in key OR amino acid side chains during the transit. Our results have helped substantiate or refute previously held notions of amino acid contribution to ligand stability in the binding pocket. Our observations of ligand activity when compared to those of experimental (electroolfactogram response) OR-activation studies provide a view to predicting the stability of ligands in the binding pocket as a precursor to OR activation by the ligand.

The Olfactory Receptor Database (ORDB; http://senselab.med.yale.edu/senselab/ordb) is a central repository of olfactory receptor (OR) and olfactory receptor-like gene and protein sequences. To deal with the very large OR gene family, we have constructed an algorithm that automatically downloads sequences from web sources such as GenBank and SWISS-PROT into the database. The algorithm uses hypertext markup language (HTML) parsing techniques that extract information relevant to ORDB. The information is then correlated with the metadata in the ORDB knowledge base to encode the unstructured text extracted into the structured format compliant with the database architecture, entity attribute value with classes and relationship (EAV/CR), which supports the SenseLab project as a whole. Three population methods: batch, automatic and semi-automatic population are discussed. The data is imported into the database using extensible markup language (XML).

DNA polymerase beta (Pol beta) is one of the key enzymes in the base excision repair pathway. The amino-terminal 8 kDa domain of Pol beta has an activity for excising a 5'-deoxyribose phosphate (dRP) group from preincised apurine/apyrimidine (AP) sites. Recent biochemical studies have identified the catalytic center of the 8 kDa domain and provided new insight into the mechanism of DNA repair by DNA polymerase beta. By incorporating both structural and biochemical data, we present here a reaction mechanism for the 5'-dRP excision activity of the 8 kDa domain. This mechanism focuses on a catalytic groove near the helix-hairpin-helix (HhH) motif of the 8 kDa domain. Our model shows that the dRP group of the AP site can be stabilized in the catalytic groove through extensive interactions with the residues of the groove and be positioned close to the active center, Lys72, which catalyzes a beta-elimination reaction by forming a Schiff base with the C1' of the dRP group.