Computing the distribution of the maximum in balls-and-boxes problems with application to clusters of disease cases

Inferring the functional effects of mutation through clusters of mutations in homologous proteins

Inferring functional consequences is a bottleneck in high-throughput cancer mutation discovery and genetic association studies. Most polymorphisms and germline mutations are unlikely to have functionally significant consequences. Most cancer somatic mutations do not contribute to tumorigenesis and are not under selective pressure. Identifying and understanding functionally important mutations can clarify disease biology and lead to new therapeutic and diagnostic opportunities. We investigated the extent to which protein mutations with functional consequences are enriched in clusters at conserved positions across related proteins. We found that disease-causing mutations form clusters more than random mutations or single nucleotide polymorphisms, confirming that mutation hotspots occur at the domain level. In addition to helping to identify functionally significant mutations, analysis of clustered mutations can indicate the mechanism and consequences for protein function. Our analysis focused on somatic cancer mutations suggests functional impact for many, including singleton mutations in FGFR1, FGFR3, GFI1B, PIK3CG, RALB, RAP2B, and STK11. This provides evidence and generates mechanistic hypotheses for the contribution of such mutations to cancer. The same approach can be applied to mutations suspected of involvement in other diseases. An interactive Web application for browsing mutation clusters is available at http://www.mcluster.org.

Subtle evolutionary changes in the distribution of N-glycosylation sequons in the HIV-1 envelope glycoprotein 120

Many viruses are known to undergo rapid evolutionary changes under selective pressures. The HIV-1 envelope glycoprotein 120 (gp120) shows extreme selection for NXS/T sequons, the potential sites of N-glycosylation. Although the average number of sequons in gp120 appears to be relatively stable in the recent past, even slight changes in the distribution of sequons may potentially play crucial roles in protein interaction and viral infection. This study tracked the prevalence and distribution of NXS/T sequons in gp120 over a period of 29 years (from 1981 to 2009). The gp120 showed location specific distribution of sequons with higher density in the outer domain of the molecule. The NXT sequon density decreased in the outer domain (despite the increase in the sequon specific amino acid threonine), but increased in the inner domain. By contrast, the NXS sequon density increased specifically in the outer domain. Related changes were also seen in the distribution probabilities of sequons in two domains. The results indicate that the gp120, chiefly in subtype B, is redistributing NXS/T sequons within the molecule with specific selection for NXS sequons. The subtle evolution of sequons in gp120 may have implications in viral resistance and infection.