Identification of 142 single nucleotide polymorphisms in 41 candidate genes for rheumatoid arthritis in the Japanese population

Single-nucleotide polymorphisms (SNPs) can make an important contribution to our understanding of genetic backgrounds that may influence medical conditions and ethnic diversity. We undertook a systematic survey of genomic DNA for SNPs located not only in coding sequences but also in non-coding regions (e.g., introns and 5' flanking regions) of selected genes. Using DNA samples from 48 Japanese patients with rheumatoid arthritis (RA) as templates, we surveyed 41 genes that represent candidates for RA, screening a total of 104 kb of DNA (30 kb of coding sequences and 74 kb of non-coding DNA). Within this 104 kb of genomic sequences we identified 163 polymorphisms (1 per 638 bases on average), of which 142 were single-nucleotide substitutions and the remainder, insertions or deletions. Of the coding SNPs, 52% were non-synonymous substitutions, and non-conservative amino acid changes were observed in a quarter of those. Sixty-nine polymorphisms showed high frequencies for minor alleles (more than 15%) and 20 revealed low frequencies (<5%). Our results indicated a greater average distance between SNPs than others have reported, but this disparity may reflect the type of genes surveyed and/or the relative ethnic homogeneity of our test population.     

Sequence diversity in 36 candidate genes for cardiovascular disorders.

Two strategies involving whole-genome association studies have been proposed for the identification of genes involved in complex diseases. The first one seeks to characterize all common variants of human genes and to test their association with disease. The second one seeks to develop dense maps of single-nucleotide polymorphisms (SNPs) and to detect susceptibility genes through linkage disequilibrium. We performed a molecular screening of the coding and/or flanking regions of 36 candidate genes for cardiovascular diseases. All polymorphisms identified by this screening were further genotyped in 750 subjects of European descent. In the whole set of genes, the lengths explored spanned 53.8 kb in the 5' regions, 68.4 kb in exonic regions, and 13 kb in the 3' regions. The strength of linkage disequilibrium within candidate regions suggests that genomewide maps of SNPs might be efficient ways to identify new disease-susceptibility genes, provided that the maps are sufficiently dense. However, the relatively large number of polymorphisms within coding and regulatory regions of candidate genes raises the possibility that several of them might be functional and that the pattern of genotype-phenotype association might be more complex than initially envisaged, as actually has been observed in some well-characterized genes. These results argue in favor of both genomewide association studies and detailed studies of the overall sequence variation of candidate genes, as complementary approaches.     

Evolutionary constraint in flanking regions of avian genes

An important comprehension from comparative genomic analysis is that sequence conservation beyond neutral expectations is frequently found outside protein-coding regions, indicating important functional roles of noncoding DNA. Understanding the causes of constraint on noncoding sequence evolution forms an important area of research, not least in light of the importance for understanding the evolution of gene expression. We aligned all orthologous genes of chicken and zebra finch together with 5 kb of their upstream and downstream noncoding sequences, to study the evolution of gene flanking sequences in the avian genome. Using ancestral repeats as a neutral reference, we detected significant evolutionary constraint in the 3' flanking region, highest directly after termination (60%) and then gradually decreasing to about 20% 5 kb downstream. Constraint was higher in annotated 3' untranslated regions (UTRs) than in non-UTRs at the same distance from the stop codon and higher in sequences annotated as microRNA (miRNA)-binding sites than in non-miRNA-binding sites within 3' UTRs. Constraint was also higher when estimated for a smaller data set of genes from more closely related songbird species, indicating turnover of functional elements during avian evolution. On the 5' flanking side constraint was readily seen within the first 125 bp immediately upstream of the start codon (34%) and was about 10% for remaining sequence within 5 kb upstream. Analysis of chicken polymorphism data gave further support for the highest constraint directly before and after the translated region. Finally, we found that genes evolving under the highest constraint measured by d(N)/d(S) also had the highest level of constraint in the 3' flanking region. This study broadens the insights into gene flanking sequence evolution by adding new findings from a vertebrate lineage other than mammals.     

Fibroblast growth factor receptor 2 (FGFR2): genomic sequence and variations

Fibroblast growth factor receptors (FGFRs) play an important role in development and tumorigenesis. Mutations in FGFR2 cause more than five craniosynostosis syndromes. The FGFR2 genomic structure is the largest of the FGFR family. We have refined and extended the genomic organization of the FGFR2 gene by sequencing more than 119 kb of PACs, cosmids, and PCR products and assembling a region of approximately 175 kb. Although the gene structure has been reported to include only 20 exons, we have verified the presence of at least 22 exons, some of which are alternatively spliced. The sizes of six exons differed from those reported previously. Comparison of our sequence and those in the NCBI database detected more than 300 potential single nucleotide polymorphisms (SNPs). However, sequencing regions containing 52 of these potential SNPs verified only 14 in PCR products generated from 16 CEPH alleles. In contrast, direct sequencing of the CEPH DNAs revealed 21 other polymorphisms. Only one SNP was found in the 2,926 bp of coding sequence. Twenty-seven SNPs, two insertion polymorphisms and five microsatellite polymorphisms are contained in approximately 16.6 kb of non-coding sequence. These data yield an average of one polymorphism for approximately 488 bp of non-coding sequence examined. This collection of SNP, insertion, and repeat polymorphisms will aid future association studies between the FGFR2 gene and human disease and will enhance mutation detection.     

Substantial linkage disequilibrium across the insulin-degrading enzyme locus but no association with late-onset Alzheimer's disease

Insulin-degrading enzyme (IDE; insulysin; EC 3.4.24.56) is a 110-kDa neutral metallopeptidase that can degrade a number of peptides including beta-amyloid. The gene encoding IDE is located on chromosome 10 close to a region of linkage for late-onset Alzheimer's disease (LOAD) and thus is a functional and positional candidate for this disorder. We analysed all of the coding exons, untranslated regions and 1000 bp of 5'-flanking sequence of IDE by using denaturing high-performance liquid chromatography and sequencing. We detected eight single nucleotide polymorphisms (SNPs), three in the 5' flanking sequence and five in the coding sequence, of which three were found at lower than 5% frequency. None of them changed the amino acid sequence. We genotyped the five SNPs with allele frequencies of more than 5% in 133 Caucasian LOAD cases and 135 controls collected in the UK and 95 cases and 117 controls collected at the Mayo Clinic, Rochester, USA. Two of the SNPs were analysed in a further independent case-control sample (Washington University, St. Louis: 86 cases, 94 controls). No significant association was found with any individual SNP in any of the samples or with any haplotypes. Analysis of the marker D10S583, which maps 36 kb upstream of IDE, also failed to show association in 134 cases and 111 matched controls from the UK ( P=0.63). Strong linkage disequilibrium was detected between the five SNPs that spanned the whole of the 120-kb genomic region of IDE and one major and a number of minor haplotypes were detected in the populations studied. We conclude that IDE does not make a substantial contribution to the aetiology of LOAD and therefore cannot account for the linkage between LOAD and 10q.     

DNA variation at the rp49 gene region of Drosophila simulans: evolutionary inferences from an unusual haplotype structure.

An approximately 1.3-kb region including the rp49 gene plus its 5' and 3' flanking regions was sequenced in 24 lines of Drosophila simulans (10 from Spain and 14 from Mozambique). Fifty-four nucleotide and 8 length polymorphisms were detected. All nucleotide polymorphisms were silent: 52 in noncoding regions and 2 at synonymous sites in the coding region. Estimated silent nucleotide diversity was similar in both populations (pi = 0.016, for the total sample). Nucleotide variation revealed an unusual haplotype structure showing a subset of 11 sequences with a single polymorphism. This haplotype was present at intermediate frequencies in both the European and the African samples. The presence of such a major haplotype in a highly recombining region is incompatible with the neutral equilibrium model. This haplotype structure in both a derived and a putatively ancestral population can be most parsimoniously explained by positive selection. As the rate of recombination in the rp49 region is high, the target of selection should be close to or within the region studied.     

HLA-E, HLA-F, and HLA-G polymorphism: genomic sequence defines haplotype structure and variation spanning the nonclassical class I genes

Despite several studies that defined the polymorphism of the nonclassical human leukocyte antigen-E (HLA-E), HLA-F, and HLA-G genes, most polymorphisms thus far examined in correlative studies were derived from the coding sequences of these genes. In addition, some discrepancies and ambiguities in the available data have persisted in current databases. To expand the data available and to resolve some of the discrepant data, we have defined protocols that allow for the amplification of 6 to 7 kb of contiguous genomic sequence for each gene, including all of the coding and intron sequences, approximately 2 kb of 5' flanking promoter sequence, and 1 kb of 3' flanking sequence. Using long-range polymerase chain reaction (PCR) protocols, generating either one or two PCR products depending on the locus, amplified genomic DNA was directly sequenced to completion using a set of about 30 primers over each locus to yield contiguous sequence data from both strands. Using this approach, we sequenced 33 genomic DNAs, from Asian, African American, and Caucasian samples. The results of this analysis confirmed several previously reported coding sequence variants, identified several new allelic variants, and also defined extensive variation in intron and flanking sequences. It was possible to construct haplotype maps and to identify tagging single nucleotide polymorphisms that can be used to detect the composite variation spanning all three genes.     

Homozygosity for multiple contiguous single-nucleotide polymorphisms as an indicator of large heterozygous deletions: identification of a novel heterozygous 8-kb intragenic deletion (IVS7-19 to IVS15-17) in a patient with glycogen storage disease type II

Current methods for detection of mutations by polymerase chain reaction (PCR) and sequence analysis frequently are not able to detect heterozygous large deletions. We report the successful use of a novel approach to identify such deletions, based on detection of apparent homozygosity of contiguous single-nucleotide polymorphisms (SNPs). The sequence analysis of genomic DNA PCR products containing all coding exons and flanking introns identified only a single heterozygous mutation (IVS18+2t-->a) in a patient with classic infantile-onset autosomal recessive glycogen storage disease type II (GSDII). Apparent homozygosity for multiple contiguous SNPs detected by this sequencing suggested presence of a large deletion as the second mutation; primers flanking the region of homozygous SNPs permitted identification and characterization by PCR of a large genomic deletion (8.26 kb) extending from IVS7 to IVS15. The data clearly demonstrate the utility of SNPs as markers for large deletions in autosomal recessive diseases when only a single mutation is found, thus complementing currently standard DNA PCR sequence methods for identifying the molecular basis of disease.     

Two human gamma-crystallin genes are linked and riddled with Alu-repeats

A human genomic cosmid clone, pHcos gamma-1, has been isolated containing two closely linked gamma-crystallin genes, oriented in the same direction. The sequence of these genes and their 5' and 3' flanking regions has been determined. The coding regions of both genes are interrupted by two introns. The first introns (94 and 100 bp, respectively) are located in the 5' region of the genes. The second introns (2.82 and 0.95 kb, respectively) divide the genes into two halves, each encoding a structural domain of the gamma-crystallin protein. The coding regions of the two genes show 80% homology. Due to a mutation in the splice acceptor site of the second intron of the first gene, the coding region of its third exon is 3 bp longer than that of the second gene. In the flanking regions several conserved sequence elements were found, including those elements that are known to be necessary for the correct expression of eukaryotic genes. The flanking and intronic regions of the genes contain 'simple sequence' DNA and Alu repeats. The Alu repeats are usually clustered, contain truncated elements, and are often located near simple sequence DNA.     

PON基因簇序列变异筛查研究

摘要：系统筛查PON1、PON2及PON3基因编码、剪接及侧翼序列,以期发现所有潜在功能多态基因座,为进一步探讨PON基因家族与心血管疾病的关系做准备。随机选择48例冠心病患者作为筛查对象, 以PCR产物直接测序检测DNA序列变异。扩增片断涵盖整个外显子, 其两侧部分内含子区域及5’和3’侧翼序列。（1）13.9kb测序范围内共发现31个多态性基因座,均为单核甘酸多态（SNP）,其中17个SNP为首次报道。（2）国人中SNP构成和等位基因频率与高加索人群存在显著差异。（3）一个基因内部两个或多个多态性基因座间存在完全或近乎完全连锁不平衡相当常见。中国汉族人群中PON基因簇多个潜在功能多态基因座的识别及这些基因座间的强连锁不平衡状态,为在国人中探讨PON基因簇与心血管疾病关系提供了重要的基础数据。     

Genomic heterogeneity in the density of noncoding single-nucleotide and microsatellite polymorphisms in Plasmodium falciparum

The density and distribution of single-nucleotide polymorphisms (SNPs) across the genome has important implications for linkage disequilibrium mapping and association studies, and the level of simple-sequence microsatellite polymorphisms has important implications for the use of oligonucleotide hybridization methods to genotype SNPs. To assess the density of these types of polymorphisms in P. falciparum, we sampled introns and noncoding DNA upstream and downstream of coding regions among a variety of geographically diverse parasites. Across 36,229 base pairs of noncoding sequence representing 41 genetic loci, a total of 307 polymorphisms including 248 polymorphic microsatellites and 39 SNPs were identified. We found a significant excess of microsatellite polymorphisms having a repeat unit length of one or two, compared to those with longer repeat lengths, as well as a nonrandom distribution of SNP polymorphisms. Almost half of the SNPs localized to only three of the 41 genetic loci sampled. Furthermore, we find significant differences in the frequency of polymorphisms across the two chromosomes (2 and 3) examined most extensively, with an excess of SNPs and a surplus of polymorphic microsatellites on chromosome 3 as compared to chromosome 2 (P=0.0001). Furthermore, at some individual genetic loci we also find a nonrandom distribution of polymorphisms between coding and flanking noncoding sequences, where completely monomorphic regions may flank highly polymorphic genes. These data, combined with our previous findings of nonrandom distribution of SNPs across chromosome 2, suggest that the Plasmodium falciparum genome may be a mosaic with regard to genetic diversity, containing chromosomal regions that are highly polymorphic interspersed with regions that are much less polymorphic.     

Comparing the frequencies of restriction fragment length polymorphisms for dystrophin gene in Chinese with those from Japanese and Caucasian populations

The restriction fragment length polymorphisms distribution and frequency of dystrophin gene in Chinese were studied by using 14 subclones of the entire 14kb cDNA for the dystrophin as hybridization probes. Allelic fragments were detected in hybridization patterns of PvuⅡ/la, Taq Ⅰ/2b-3, Taq Ⅰ/5b-7, and Xba Ⅰ/10. Among them, the allelic fragments (26kb and 3.8kb) in PvuⅡ/2b-3 pattern and the allelic fragments (10.0kb and 8.4kb) in Taq Ⅰ/5b-7 patterns had never been reported previously. Compared with the data from Caucasians and Japanese, it indicated that there was a significant difference (P<0.01) of the allelic fragment frequency in Taq Ⅰ/2b-3 and Xba Ⅰ/10 patterns between Chinese and Caucasians. The frequencies of allelic fragments A2 (5.6kb) in Taq Ⅰ/8 and A2 (10.Tkb) in EcoR Ⅴ/9 were high in Caucasians, yet had not been detected in Chinese, the differences were also highly significant. But in Chinese and Caucasians, the B1B2 allelic frequencies in Taq Ⅰ/5b-7 are the same. As to the frequency of the allelic fragments A1A2 and B1B2 in Pvu Ⅱ/la, there was no significant difference between Chinese and Japanese.     

The Bos taurus casein genes. Isolation and characterization of the kappa-casein gene

A region spanning 25 kb of genomic DNA containing the kappa-casein gene, has been isolated from two genomic libraries in EMBL3 and EMBL4 phage vectors. Five phage clones containing kappa-casein gene have been found. Gene organisation has been determined using restriction mapping and a partial sequencing the 5' and 3' flanking regions. The kappa-casein gene includes 5 exons, the first of them coding for 64 nucleotides from the 5' untranslated mRNA zone. The gene is 12.5 kb long, which is almost 16 times longer than the corresponding mRNA. The first intron spans 2.5 kb, the second is the largest one and spans 5.5 kb. The 5' flanking region sequence has been analysed; it contains a TATA box from -30 to -25 bp, somewhat different from the canonic sequence, and a CAAT box at -80 bp.     

Cap-independent translational enhancement of turnip crinkle virus genomic and subgenomic RNAs

The presence of translational control elements and cap structures has not been carefully investigated for members of the Carmovirus genus, a group of small icosahedral plant viruses with positive-sense RNA genomes. In this study, we examined both the 5' and 3' untranslated regions (UTRs) of the turnip crinkle carmovirus (TCV) genomic RNA (4 kb) as well as the 5' UTR of the coat protein subgenomic RNA (1.45 kb) for their roles in translational regulation. All three UTRs enhanced translation of the firefly luciferase reporter gene to different extents. Optimal translational efficiency was achieved when mRNAs contained both 5' and 3' UTRs. The synergistic effect due to the 5'-3' cooperation was at least fourfold greater than the sum of the contributions of the individual UTRs. The observed translational enhancement of TCV mRNAs occurred in a cap-independent manner, a result consistent with the demonstration, using a cap-specific antibody, that the 5' end of the TCV genomic RNA was uncapped. Finally, the translational enhancement activity within the 5' UTR of 1.45-kb subgenomic RNA was shown to be important for the translation of coat protein in protoplasts and for virulent infection in Arabidopsis plants.