Inter-individual variation in gene expression has proven to be in part governed by genetic determinants, which may be trans- or cis-acting. The underlying cause of cis-acting regulatory variation has been identified in only a handful of the hundreds of genes shown to display differential allelic expression. In this report, we describe a systematic effort to map common cis-acting variants in 64 genes, using association methods in HapMap samples. We identified 16 loci (25%), each of which harbors common haplotypes that affect total expression of a gene, and a further 17 loci (27%) with evidence of haplotypes affecting relative allelic expression in heterozygote samples. Our survey suggests that detailed mapping of allele-specific in vivo expression will provide a rich source of regulatory SNPs or haplotypes that should be given high priority in association studies of human phenotypes.
Cis-acting allelic variation in gene regulation is a source of phenotypic variation. Consequently, recent studies have experimentally screened human genes in an attempt to initiate a catalog of genes possessing cis-acting variants. In this study, we use human EST data in dbEST as the source of allelic expression data, and the HapMap database to provide expected allele frequencies in human populations. We demonstrate a greater concordance of allele frequencies estimated from human ESTs in dbEST with those derived from the CEPH HapMap sample representing Caucasians from northern and western Europe, than population samples obtained in Asia and Africa. Deviations between allele frequencies observed in EST databases and the ones obtained from the CEPH HapMap samples may result from common heritable cis-acting variants altering the relative allele distribution in RNA. We provide in silico as well as experimental evidence that this strategy does allow significant enrichment of genes harboring common heritable cis-acting polymorphisms in linkage disequilibrium with expressed alleles.
Diabetes and obesity have long been known to be related. The recently characterized adipocyte hormone resistin (also called FIZZ3/ADSF) has been implicated as a molecular link between impaired glucose tolerance (IGT) and obesity in mice. A search for sequence variants at the human resistin locus identified nine single-nucleotide polymorphisms (SNPs) but no coding variants. An investigation into the association of these SNPs with diabetes and obesity revealed two 5' flanking variants (g.-537 and g.-420), in strong linkage disequilibrium, that are associated with BMI. In nondiabetic individuals from the Quebec City area and the Saguenay-Lac-St-Jean region of Quebec, the g.-537 mutation (allelic frequency = 0.04) was significantly associated with an increase in BMI (P = 0.03 and P = 0.01, respectively). When the data from these two populations were combined and adjusted for age and sex, both the g.-537 (odds ratio [OR] 2.72, 95% CI 1.28-5.81) and the g.-420 variants (1.58, 1.06-2.35) were associated with an increased risk for a BMI > or =30 kg/m(2). In contrast, in case/control and family-based study populations from Scandinavia, we saw no effect on BMI with either of these promoter variants. No association was seen with diabetes in any of the population samples.
Vitamin B 12 (cobalamin) is an essential cofactor of two enzymes, methionine synthase and methylmalonyl-CoA mutase. The conversion of the vitamin to its coenzymes requires a series of biochemical modifications for which several genetic diseases are known, comprising eight complementation groups ( cblA through cblH ). The objective of this study was to clone the gene responsible for the cblA complementation group thought to represent a mitochondrial cobalamin reductase. Examination of bacterial operons containing genes in close proximity to the gene for methylmalonyl-CoA mutase and searching for orthologous sequences in the human genome yielded potential candidates. A candidate gene was evaluated for deleterious mutations in cblA patient cell lines, which revealed a 4-bp deletion in three cell lines, as well as an 8-bp insertion and point mutations causing a stop codon and an amino acid substitution. These data confirm that the identified gene, MMAA , corresponds to the cblA complementation group. It is located on chromosome 4q31.1–2 and encodes a predicted protein of 418 aa. A Northern blot revealed RNA species of 1.4, 2.6, and 5.5 kb predominating in liver and skeletal muscle. The deduced amino acid sequence reveals a domain structure, which belongs to the AAA ATPase superfamily that encompasses a wide variety of proteins including ATP-binding cassette transporter accessory proteins that bind ATP and GTP. We speculate that we have identified a component of a transporter or an accessory protein that is involved in the translocation of vitamin B 12 into mitochondria.
Mutations in the MTR gene, which encodes methionine synthase on human chromosome 1p43, result in the methylcobalamin deficiency G (cblG) disorder, which is characterized by homocystinuria, hyperhomocysteinemia, and hypomethioninemia. To investigate the molecular basis of the disorder, we have characterized the structure of the MTR gene, thereby identifying exon-intron boundaries. This enabled amplification of each of the 33 exons of the gene, from genomic DNA from a panel of 21 patients with cblG. Thirteen novel mutations were identified. These included five deletions (c.12-13delGC, c.381delA, c.2101delT, c.2669-2670delTG, and c.2796-2800delAAGTC) and two nonsense mutations (R585X and E1204X) that would result in synthesis of truncated proteins that lack portions critical for enzyme function. One mutation was identified that resulted in conversion of A to C of the invariant A of the 3' splice site of intron 9. Five missense mutations (A410P, S437Y, S450H, H595P, and I804T) were identified. The latter mutations, as well as the splice-site mutation, were not detected in a panel of 50 anonymous DNA samples, suggesting that these sequence changes are not polymorphisms present in the general population. In addition, a previously described missense mutation, P1173L, was detected in 16 patients in an expanded panel of 24 patients with cblG. Analysis of haplotypes constructed using sequence polymorphisms identified within the MTR gene demonstrated that this mutation, a C-->T transition in a CpG island, has occurred on at least two separate genetic backgrounds.
Autosomal recessive spastic ataxia of Charlevoix-Saguenay (ARSACS or SACS) is a neurodegenerative disease frequent in northeastern Québec. In a previous study, we localized the disease gene to chromosome region 13q11 by identifying excess sharing of a marker allele in patients followed by linkage analysis and haplotyping. To create a detailed physical map of this region, we screened CEPH mega-YACs with 41 chromosome 13 sequence-tagged-sites (STSs) known to map to 13q11-q12. The YAC contig, composed of 27 clones, extends on the genetic map from D13S175 to D13S221, an estimated distance of at least 19.3 cM. A high-resolution BAC and PAC map that includes the ARSACS critical region flanked by D13S1275 and D13S292 was constructed. These YAC and BAC/PAC maps allowed the accurate placement of 29 genes and ESTs previously mapped to the proximal region of chromosome 13q. We confirmed the position of two candidate genes within the critical region and mapped the other 27 genes and ESTs to nearby intervals. Six BAC/PAC clones form a contig between D13S232 and D13S787 for sequencing within the ARSACS critical region.