Low hanging fruit: a subset of human cSNPs is both highly non-uniform and predictable

Relative effects of mutability and selection on single nucleotide polymorphisms in transcribed regions of the human genome

MOTIVATION

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation in humans. However, the factors that affect SNP density are poorly understood. The goal of this study was to estimate the relative effects of mutability and selection on SNP density in transcribed regions of human genes. It is important for prediction of the regions that harbor functional polymorphisms.

RESULTS

We used frequency-validated SNPs resulting from single-nucleotide substitutions. SNPs were subdivided into five functional categories: (i) 5' untranslated region (UTR) SNPs, (ii) 3' UTR SNPs, (iii) synonymous SNPs, (iv) SNPs producing conservative missense mutations, and (v) SNPs producing radical missense mutations. Each of these categories was further subdivided into nine mutational categories on the basis of the single-nucleotide substitution type. Thus, 45 functional/mutational categories were analyzed. The relative mutation rate in each mutational category was estimated on the basis of published data. The proportion of segregating sites (PSSs) for each functional/mutational category was estimated by dividing the observed number of SNPs by the number of potential sites in the genome for a given functional/mutational category. By analyzing each functional group separately, we found significant positive correlations between PSSs and relative mutation rates (Spearman's correlation coefficient, at least r = 0.96, df = 9, P < 0.001). We adjusted the PSSs for the mutation rate and found that the functional category had a significant effect on SNP density (F = 5.9, df = 4, P = 0.001), suggesting that selection affects SNP density in transcribed regions of the genome. We used analyses of variance and covariance to estimate the relative effects of selection (functional category) and mutability (relative mutation rate) on the PSSs and found that approximately 87% of variation in PSS was due to variation in the mutation rate and approximately 13% was due to selection, suggesting that the probability that a site located in a transcribed region of a gene is polymorphic mostly depends on the mutability of the site.

Arginine (CGC) codon targeting in the human prostacyclin receptor gene (PTGIR) and G-protein coupled receptors (GPCR)

The human prostacyclin receptor (hIP) has recently been recognized as an important seven transmembrane G-protein coupled receptor that plays critical roles in atheroprevention and cardioprotection. To date, four non-synonymous genetic variants have been identified, two of which occur at the same Arg amino acid position (R212H, R212C). This observation instigated further genetic screening for prostacyclin receptor variants on 1455 human genomic samples. A total of 31 distinct genetic variants were detected, with 6 (19%) involving Arg residues. Distinct differences in location and frequencies of genetic variants were noted between Caucasian, Asian, Hispanic and African Americans, with the most changes noted in the Asian cohort. From the sequencing results, three Arg-targeted changes at the same 212 position within the third cytoplasmic loop of the human prostacyclin (hIP) receptor were detected: 1) R212C (CGC-->TGC), 2) R212H (CGC-->CAC), and 3) R212R (CGC-->CGT). Three additional Arg codon variants (all exhibiting the same CGC to TGC change) were also detected, R77C, R215C, and R279C. Analysis (GPCR and SNP databases) of 200 other GPCRs, with recorded non-synonymous mutations, confirmed a high frequency of Arg-targeted missense mutations, particularly within the important cytoplasmic domain. Preferential nucleotide changes (at Arg codons), were observed involving cytosine (C) to thymine (T) (pyrimidine to pyrimidine), as well as guanine (G) to adenine (A) (purine to purine) (p<0.001, Pearson's goodness-of-fit test). Such targeting of Arg residues, leading to significant changes in coding amino acid size and/or charge, may have potentially-important structural and evolutionary implications on the hIP and GPCRs in general. In the case of the human prostacyclin receptor, such alterations may reduce the cardio-, vasculo-, and cytoprotective effects of prostacyclin.

Precise and parallel characterization of coding polymorphisms, alternative splicing, and modifications in human proteins by mass spectrometry

The human proteome is a highly complex extension of the genome wherein a single gene often produces distinct protein forms due to alternative splicing, RNA editing, polymorphisms, and posttranslational modifications. Such biological variation compounded by the high sequence identity within gene families currently overwhelms the complete and routine characterization of mammalian proteins by MS. A new data base of human proteins (and their possible variants) was created and searched using tandem mass spectrometric data from intact proteins. This first application of top down MS/MS to wild-type human proteins demonstrates both gene-specific identification and the unambiguous characterization of multifaceted mass shifts (Deltam values). Such Deltam values found from the precise identification of 45 protein forms from HeLa cells reveal 34 coding single nucleotide polymorphisms, two protein forms from alternative splicing, and 12 diverse modifications (not including simple N-terminal processing), including a previously unknown phosphorylation at 10% occupancy. Automated protein identification was achieved with a median expectation value of 10(-13) and often occurred simultaneously with dissection of diverse sources of protein variability as they occur in combination. Top down MS therefore has a bright future for enabling precise annotation of gene products expressed from the human genome by non-mass spectrometrists.