A human model for multigenic inheritance: Phenotypic expression in Hirschsprung disease requires both the RET gene and a new 9q31 locus

Reduced penetrance in genetic disorders may be either dependent or independent of the genetic background of gene carriers. Hirschsprung disease (HSCR) demonstrates a complex pattern of inheritance with ≈50% of familial cases being heterozygous for mutations in the receptor tyrosine kinase RET. Even when identified, the penetrance of RET mutations is only 50–70%, gender-dependent, and varies with the extent of aganglionosis. We searched for additional susceptibility genes which, in conjunction with RET, lead to phenotypic expression by studying 12 multiplex HSCR families. Haplotype analysis and extensive mutation screening demonstrated three types of families: six families harboring severe RET mutations (group I); and the six remaining families, five of which are RET-linked families with no sequence alterations and one RET-unlinked family (group II). Although the presence of RET mutations in group I families is sufficient to explain HSCR inheritance, a genome scan reveals a new susceptibility locus on 9q31 exclusively in group II families. As such, the gene at 9q31 is a modifier of HSCR penetrance. These observations imply that identification of new susceptibility factors in a complex disease may depend on classification of families by mutational type at known susceptibility genes.

Estimating the probability for a protein to have a new fold: A statistical computational model

Structural genomics aims to solve a large number of protein structures that represent the protein space. Currently an exhaustive solution for all structures seems prohibitively expensive, so the challenge is to define a relatively small set of proteins with new, currently unknown folds. This paper presents a method that assigns each protein with a probability of having an unsolved fold. The method makes extensive use of protomap, a sequence-based classification, and scop, a structure-based classification. According to protomap, the protein space encodes the relationship among proteins as a graph whose vertices correspond to 13,354 clusters of proteins. A representative fold for a cluster with at least one solved protein is determined after superposition of all scop (release 1.37) folds onto protomap clusters. Distances within the protomap graph are computed from each representative fold to the neighboring folds. The distribution of these distances is used to create a statistical model for distances among those folds that are already known and those that have yet to be discovered. The distribution of distances for solved/unsolved proteins is significantly different. This difference makes it possible to use Bayes' rule to derive a statistical estimate that any protein has a yet undetermined fold. Proteins that score the highest probability to represent a new fold constitute the target list for structural determination. Our predicted probabilities for unsolved proteins correlate very well with the proportion of new folds among recently solved structures (new scop 1.39 records) that are disjoint from our original training set.