Phenotype-genotype correlation in Hirschsprung disease is illuminated by comparative analysis of the RET protein sequence.

The ability to discriminate between deleterious and neutral amino acid substitutions in the genes of patients remains a significant challenge in human genetics. The increasing availability of genomic sequence data from multiple vertebrate species allows inclusion of sequence conservation and physicochemical properties of residues to be used for functional prediction. In this study, the RET receptor tyrosine kinase serves as a model disease gene in which a broad spectrum (> or = 116) of disease-associated mutations has been identified among patients with Hirschsprung disease and multiple endocrine neoplasia type 2. We report the alignment of the human RET protein sequence with the orthologous sequences of 12 non-human vertebrates (eight mammalian, one avian, and three teleost species), their comparative analysis, the evolutionary topology of the RET protein, and predicted tolerance for all published missense mutations. We show that, although evolutionary conservation alone provides significant information to predict the effect of a RET mutation, a model that combines comparative sequence data with analysis of physiochemical properties in a quantitative framework provides far greater accuracy. Although the ability to discern the impact of a mutation is imperfect, our analyses permit substantial discrimination between predicted functional classes of RET mutations and disease severity even for a multigenic disease such as Hirschsprung disease.

A common sex-dependent mutation in a RET enhancer underlies Hirschsprung disease risk.

The identification of common variants that contribute to the genesis of human inherited disorders remains a significant challenge. Hirschsprung disease (HSCR) is a multifactorial, non-mendelian disorder in which rare high-penetrance coding sequence mutations in the receptor tyrosine kinase RET contribute to risk in combination with mutations at other genes. We have used family-based association studies to identify a disease interval, and integrated this with comparative and functional genomic analysis to prioritize conserved and functional elements within which mutations can be sought. We now show that a common non-coding RET variant within a conserved enhancer-like sequence in intron 1 is significantly associated with HSCR susceptibility and makes a 20-fold greater contribution to risk than rare alleles do. This mutation reduces in vitro enhancer activity markedly, has low penetrance, has different genetic effects in males and females, and explains several features of the complex inheritance pattern of HSCR. Thus, common low-penetrance variants, identified by association studies, can underlie both common and rare diseases.

Genome-wide association study and mouse model identify interaction between RET and EDNRB pathways in Hirschsprung disease.

Genetic studies of Hirschsprung disease, a common congenital malformation, have identified eight genes with mutations that can be associated with this condition. Mutations at individual loci are, however, neither necessary nor sufficient to cause clinical disease. We conducted a genome-wide association study in 43 Mennonite family trios using 2,083 microsatellites and single-nucleotide polymorphisms and a new multipoint linkage disequilibrium method that searches for association arising from common ancestry. We identified susceptibility loci at 10q11, 13q22 and 16q23; the gene at 13q22 is EDNRB, encoding a G protein-coupled receptor (GPCR) and the gene at 10q11 is RET, encoding a receptor tyrosine kinase (RTK). Statistically significant joint transmission of RET and EDNRB alleles in affected individuals and non-complementation of aganglionosis in mouse intercrosses between Ret null and the Ednrb hypomorphic piebald allele are suggestive of epistasis between EDNRB and RET. Thus, genetic interaction between mutations in RET and EDNRB is an underlying mechanism for this complex disorder.