Identification of 96 single nucleotide polymorphisms in eight genes involved in iron metabolism: efficiency of bioinformatic extraction compared with a systematic sequencing approach

Single nucleotide polymorphisms (SNPs) can significantly contribute to the characterization of the genes predisposing to iron overloads or deficiencies. We report an SNP survey of coding and non-coding regions of eight genes involved in iron metabolism, by two successive methods. First, we made use of the public domain sequence data, by using assembled expressed sequence tags, non-redundant sequences, and SNP database screening. We extracted 77 potential SNPs of which only 31 could be further validated by sequencing DNA from 44 unrelated multi-ethnic individuals. Our results indicate that a bioinformatic approach may be effective only in those cases where candidate SNPs are extracted from two different data sources or in cases of experimentally confirmed SNPs. Second, additional systematic sequencing of DNA from 24 unrelated Breton subjects increased the number of SNPs over a total length of 86 kb to 96. The average distance between the SNPs and minor allele frequencies were higher than reported by others authors; this discrepancy may reflect the nature of the genes studied and the ethnic homogeneity of our test population.     

Identification of 187 single nucleotide polymorphisms (SNPs) among 41 candidate genes for ischemic heart disease in the Japanese population

To investigate whether common variants in the human genetic background are associated with pathogenesis of ischemic heart diseases, we systematically surveyed 41 possible candidate genes for single-nucleotide polymorphisms (SNPs) by directly sequencing 96 independent alleles at each locus, derived from 48 unrelated Japanese patients with myocardial infarction, including 25.8 kb 5' flanking regions, 56.8 kb exonic and 35.4 kb intronic sequences, and 1.8 kb 3' flanking regions. In this genomic DNA of nearly 120 kb, we identified 187 SNPs: 55 in 5' flanking regions, seven in 5' untranslated regions (UTRs), 52 in coding elements, 64 in introns, eight in 3' UTRs, and one in a 3' flanking region. Among the 52 coding SNPs, 26 were non-synonymous changes. Allelic frequencies of some of the polymorphisms were significantly different from those reported in European populations. For example, the Q506R substitution in the coagulation factor V gene, the so-called "Leiden mutation", has a reported frequency of 2.3% in Europeans, but we detected the Leiden mutation in none of the Japanese genomes that we investigated. The allelic frequencies of the -33A>G SNP in the thrombomodulin gene were also very different; this allele occurred at a 12% frequency in the Japanese patients that we examined, although it had been detected in none of 82 Caucasians reported previously. These data support the hypothesis that some SNPs are specific to particular ethnic groups.     

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.     

The analysis of nucleotide substitutions, gaps, and recombination events between RHD and RHCE genes through complete sequencing

We determined the entire nucleotide sequences of all introns within the RHD and RHCE genes by amplifying genomic DNA using long PCR methods. The RHD and RHCE genes were 57,295 and 57,831 bp in length, respectively. Aligning both genes revealed 138 gaps (insertions and deletions) below 100 bp, 1116 substitutions in all introns and all exons (coding region), and 5 gaps of over 100 bp. Homologies (%) between the RH genes were 93.8% over all introns and coding exons and 91.7% over all exons and introns. Various short tandem repeats (STRs) and many interspersed nuclear elements were identified in both genes. The proportions of Alu sequences in the RHD and RHCE genes were 25.9 and 25.7%, respectively and these Alu sequences were concentrated in several regions. We confirmed multiple recombinations in introns 1 and 2. Such multiple recombination, which probably arose due to the concentrations of Alu sequences and the high level of the homology (%), is one of most important factors in the formation and evolution of RH gene. The variability of the Rh system may be generated because of these features of RH genes. Apparent mutational hotspots and regions with low of K values (the numbers of substitutions per nucleotide site) caused by recombinations as well as true mutational hotspots may be found in human genome. Accordingly, in searching for and identifying single nucleotide polymorphisms (SNPs) especially in noncoding regions, apparent mutational hotspots and areas of low K values by recombination should be noted since the unequal distribution of SNPs will reduce the power of SNPs as genetic maker. Combining the complete sequences' data of both RH genes with serological findings will provide beneficial information with which to elucidate the mechanism of recombination, mutation, polymorphism, and evolution of other genes containing the RH gene as well as to analyze Rh variants and develop new methods of Rh genotyping.     

Variable SNP density in aspartyl-protease genes of the malaria parasite Plasmodium falciparum

An analysis of the diversity of the aspartyl proteases of Plasmodium falciparum, known as plasmepsins (PMs), was completed in view of their possible role as drug targets. DNA sequence polymorphisms were identified in nine pm genes including their non-coding (introns and 5' flanking) sequences. All genes contained at least one single nucleotide polymorphism (SNP). Extensive microsatellite diversity was observed predominantly in non-coding sequences. All but one non-synonymous polymorphism (a conservative substitution) were mapped to the surface of the predicted protein, contradicting a possible role in enzymatic activity. The distribution of SNPs was found to be non-random among pm genes, with pm6 and pm10 having significantly higher SNP densities, suggesting they were under selection. For pm6 the majority of the SNPs were in introns and some of these may contribute to splice site variation. SNPs were found at a high density in both the coding and non-coding sequences of pm10. Recombination was important in generating additional diversity at this locus. Although direct selection for pm10 mutations could not be ruled out, the presence of balancing selection and a high density of SNPs in non-coding sequence led us to propose that another gene under selection may be influencing the diversity in the region. By sequencing short DNA tags in a 200 kb region flanking pm10 we show that a cluster of antigen genes, known to be under diversifying selection, may contribute to the observed diversity. We discuss the importance of diversity and local selection effects when choosing drug targets for intervention strategies.     

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

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

Complex high-resolution linkage disequilibrium and haplotype patterns of single-nucleotide polymorphisms in 2.5 Mb of sequence on human chromosome 21.

One approach to identify potentially important segments of the human genome is to search for DNA regions with nonrandom patterns of human sequence variation. Previous studies have investigated these patterns primarily in and around candidate gene regions. Here, we determined patterns of DNA sequence variation in 2.5 Mb of finished sequence from five regions on human chromosome 21. By sequencing 13 individual chromosomes, we identified 1460 single-nucleotide polymorphisms (SNPs) and obtained unambiguous haplotypes for all chromosomes. For all five chromosomal regions, we observed segments with high linkage disequilibrium (LD), extending from 1.7 to>81 kb (average 21.7 kb), disrupted by segments of similar or larger size with no significant LD between SNPs. At least 25% of the contig sequences consisted of segments with high LD between SNPs. Each of these segments was characterized by a restricted number of observed haplotypes,with the major haplotype found in over 60% of all chromosomes. In contrast, the interspersed segments with low LD showed significantly more haplotype patterns. The position and extent of the segments of high LD with restricted haplotype variability did not coincide with the location of coding sequences. Our results indicate that LD and haplotype patterns need to be investigated with closely spaced SNPs throughout the human genome, independent of the location of coding sequences, to reliably identify regions with significant LD useful for disease association studies.     

人群变异的分子基础:从单核苷酸多态性到单氨基酸多态性

单核苷酸多态性可以划分为位于基因编码区的SNP和非编码区的SNP两大种类;而在基因编码区的SNP还可以进一步划分为两个亚类:不改变氨基酸序列的同义SNP和改变氨基酸序列的非同义SNP.显然,非同义SNP将导致氨基酸序列的改变,即形成单氨基酸多态性.基于蛋白质组学方法,对亚洲人群血浆样本中的SAP进行了系统研究,发现某一特定SAP在纯合人群和杂合人群中可能与生理或病理性状有着不同的关联.更为重要的是,近期有研究发现,在生物体中广泛存在着RNA序列与DNA序列不一致的现象.导致这种差异的主要原因是在转录水平上存在着规模化的RNA编辑(被称为RNA编辑组,RNA editome).该发现表明,个体拥有的SAP中可能有一部分与基因组SNP无关,而是源于RNA编辑组.进一步推论,可能在翻译水平上存在着不依赖DNA和RNA序列的全新的SAP.     

Conserved worldwide linkage disequilibrium in the human factor XI gene

We have identified, in four diverse human populations, five common single-nucleotide polymorphisms (SNPs) in the coding region of the gene for the blood coagulation protease factor XI. Each SNP has an allele frequency >5% in at least one population. Three of the SNPs (C472T, A844G, and T1234C), spread out over approximately 10 kb of genomic DNA, are in marked linkage disequilibrium (LD) with one another (P < 10(-4)). Interestingly, haplotypes associated with the linked SNPs are conserved across all populations studied, despite significantly different allele frequencies between populations. The presence of such common, widely dispersed haplotypes could complicate the interpretation of LD studies and emphasizes the need for a better understanding of general patterns of LD to facilitate identification of genes for common disorders.     

Effects of GWAS-Associated Genetic Variants on lncRNAs within IBD and T1D Candidate Loci

Long non-coding RNAs are a new class of non-coding RNAs that are at the crosshairs in many human diseases such as cancers, cardiovascular disorders, inflammatory and autoimmune disease like Inflammatory Bowel Disease (IBD) and Type 1 Diabetes (T1D). Nearly 90% of the phenotype-associated single-nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWAS) lie outside of the protein coding regions, and map to the non-coding intervals. However, the relationship between phenotype-associated loci and the non-coding regions including the long non-coding RNAs (lncRNAs) is poorly understood. Here, we systemically identified all annotated IBD and T1D loci-associated lncRNAs, and mapped nominally significant GWAS/ImmunoChip SNPs for IBD and T1D within these lncRNAs. Additionally, we identified tissue-specific cis-eQTLs, and strong linkage disequilibrium (LD) signals associated with these SNPs. We explored sequence and structure based attributes of these lncRNAs, and also predicted the structural effects of mapped SNPs within them. We also identified lncRNAs in IBD and T1D that are under recent positive selection. Our analysis identified putative lncRNA secondary structure-disruptive SNPs within and in close proximity (+/−5 kb flanking regions) of IBD and T1D loci-associated candidate genes, suggesting that these RNA conformation-altering polymorphisms might be associated with diseased-phenotype. Disruption of lncRNA secondary structure due to presence of GWAS SNPs provides valuable information that could be potentially useful for future structure-function studies on lncRNAs.     

Characteristics of X- and Y-chromosome specific regions of the amelogenin gene and a PCR-based method for sex identification in red deer (<Emphasis Type="Italic">Cervus elaphus</Emphasis>)

The present study attempts to analyse sequences of the X- and Y-chromosome specific regions of the amelogenin (AMEL) gene in red deer. To this end, primers specific for each form of the gene (AMELX and AMELY) were designed based on bovine genomic sequences and the homologous regions of the genes were sequenced. The obtained sequence of AMELX gene showed high similarity with the corresponding region in cattle (91%) and humans (77%), but this similarity was slightly lower among AMELY genes and showed 87 and 73% of identical nucleotides, respectively. In addition, three single nucleotide polymorphisms (SNPs) were found in the AMELX gene of the female red deer investigated. Comparative analysis of the homologous fragments of the red deer AMELX and AMELY genes confirmed the deletion of an AMELY gene fragment in relation to AMELX. Homology of both sequences was 82% of identical nucleotides in the coding region and 74% in 3′ non-coding sequence. The sequences studied showed considerable similarity to homologous fragments of the human and bovine gene, but the structural differences observed lead us to design PCR-based method for sex identification in red deer, based on the presented sequences.     

A genome-wide survey of copy number variation regions in various chicken breeds by array comparative genomic hybridization method

The discovery of copy number variation (CNV) in the genome has provided new insight into genomic polymorphism. Studies with chickens have identified a number of large CNV segments using a 385k comparative genomic hybridization (CGH) chip (mean length >140 kb). We present a detailed CNV map for local Chinese chicken breeds and commercial chicken lines using an Agilent 400k array CGH platform with custom-designed probes. We identified a total of 130 copy number variation regions (CNVRs; mean length = 25.70 kb). Of these, 104 (80.0%) were novel segments reported for the first time in chickens. Among the 104 novel CNVRs, 56 (53.8%) of the segments were non-coding sequences, 65 (62.5%) showed the gain of DNA and 40 (38.5%) showed the loss of DNA (one locus showed both loss and gain). Overlapping with the formal selective sweep data and the quantitative trait loci data, we identified four loci that might be considered to be high-confidence selective segments that arose during the domestication of chickens. Compared with the CNVRs reported previously, genes for the positive regulation of phospholipase A2 activity were discovered to be significantly over-represented in the novel CNVRs reported here by gene ontology analysis. Availability of our results should facilitate further research in the study of the genetic variability in chicken breeds.     

Single-nucleotide polymorphism detection by denaturing high-performance liquid chromatography and direct sequencing in genes in the MHC class III region encoding novel cell surface molecules

The class III region of the human major histocompatibility complex (MHC) contains approximately 59 genes, many of which encode polypeptides with a variety of different functions. Eight of these genes are of particular interest because they encode novel surface molecules that could be involved in immune and/or inflammatory responses and are excellent candidates as disease susceptibility loci. These molecules are members of two different superfamilies, the immunoglobulin superfamily (1C7, G6B, and G6F genes) and the leucocyte antigen-6 superfamily (G6C, G6D, G6E, G5C, and G5B genes). Some level of variation was found when overlapping genomic DNAs from different haplotypes were compared. The present work describes a systematic search for single-nucleotide polymorphisms (SNPs) in these genes using direct sequencing and denaturing high-performance liquid chromatography (DHPLC) in 24 unrelated healthy individuals. We validated the DHPLC methodology by first studying the 1C7 gene. This gene was directly sequenced in all 24 samples, and DHPLC was found to resolve all the polymorphic sites present in the heterozygote samples tested. We screened the rest of the genes by DHPLC only, and only those chromatograms that revealed a polymorphic profile were sequenced. We detected one SNP every 489 bp in the 18 kb of DNA studied, corresponding to theta = 4.61x10-4. The diversity in noncoding regions is 1 SNP/560 bp, but a higher frequency was detected in coding regions with 1 SNP/423 bp corresponding to theta =5.33x10-4. Of the coding SNPs, 63.6% caused amino acid substitutions. The power of this study is emphasized by the fact that of the 37 SNPs/indels detected, only 6 can be found in the SNP database at the NCBI.     

First-generation SNP/InDel markers tagging loci for pathogen resistance in the potato genome

A panel of 17 tetraploid and 11 diploid potato genotypes was screened by comparative sequence analysis of polymerase chain reaction (PCR) products for single nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (InDels), in regions of the potato genome where genes for qualitative and/or quantitative resistance to different pathogens have been localized. Most SNP and InDel markers were derived from bacterial artificial chromosome (BAC) insertions that contain sequences similar to the family of plant genes for pathogen resistance having nucleotide-binding-site and leucine-rich-repeat domains (NBS-LRR-type genes). Forty-four such NBS-LRR-type genes containing BAC-insertions were mapped to 14 loci, which tag most known resistance quantitative trait loci (QTL) in potato. Resistance QTL not linked to known resistance-gene-like (RGL) sequences were tagged with other markers. In total, 78 genomic DNA fragments with an overall length of 31 kb were comparatively sequenced in the panel of 28 genotypes. 1498 SNPs and 127 InDels were identified, which corresponded, on average, to one SNP every 21 base pairs and one InDel every 243 base pairs. The nucleotide diversity of the tetraploid genotypes (pi = 0.72 x 10(-3)) was lower when compared with diploid genotypes (pi = 2.31 x 10(-3)). RGL sequences showed higher nucleotide diversity when compared with other sequences, suggesting evolution by divergent selection. Information on sequences, sequence similarities, SNPs and InDels is provided in a database that can be queried via the Internet.     

Single Nucleotide polymorphisms and their relationship to codon usage bias in the Pacific oyster Crassostrea gigas

DNA sequence polymorphism and codon usage bias were investigated in a set of 41 nuclear loci in the Pacific oyster Crassostrea gigas. Our results revealed a very high level of DNA polymorphism in oysters, in the order of magnitude of the highest levels reported in animals to date. A total of 290 single nucleotide polymorphisms (SNPs) were detected, 76 of which being localised in exons and 214 in non-coding regions. Average density of SNPs was estimated to be one SNP every 60 bp in coding regions and one every 40 bp in non-coding regions. Non-synonymous substitutions contributed substantially to the polymorphism observed in coding regions. The non-synonymous to silent diversity ratio was 0.16 on average, which is fairly higher to the ratio reported in other invertebrate species recognised to display large population sizes. Therefore, purifying selection does not appear to be as strong as it could have been expected for a species with a large effective population size. The level of non-synonymous diversity varied greatly from one gene to another, in accordance with varying selective constraints. We examined codon usage bias and its relationship with DNA polymorphism. The table of optimal codons was deduced from the analysis of an EST dataset, using EST counts as a rough assessment of gene expression. As recently observed in some other taxa, we found a strong and significant negative relationship between codon bias and non-synonymous diversity suggesting correlated selective constraints on synonymous and non-synonymous substitutions. Codon bias as measured by the frequency of optimal codons for expression might therefore provide a useful indicator of the level of constraint upon proteins in the oyster genome.     

Accepted mutations in a gene family: Evolutionary diversification of duplicated DNA

Summary We report and compare the DNA sequences of 14 silkmoth (Antheraea polyphemus) chorion genes, derived from either cDNA or chromosomal DNA clones. Seven of these genes are members of the A multigene family, and seven are members of the B family. Where available, the previously reported (Jones and Kafatos 1980) intronic and extragenic flanking DNA sequences are also considered. Closely related sequences are compared, revealing the types of spontaneous mutations that were fixed during paralogous evolution. Segmental mutations (i.e. mutations other than substitutions) are nearly always interpretable as small duplications or deletions. related to small direct repeats. Segmental mutations are strongly constrained in the coding regions, although they do occur. Nucleotide substitutions also appear to be under selective constraints: relatively few substitutions leading to amino acid replacements are accepted, silent substitutions leading to some codons (especially purine-terminated ones) are disfavored, and different compositional biases are maintained in different parts of the sequences. Other sequence differences can be interpreted as indicative of neutral drift, including most differences in non-coding regions and most T/C transitions in third-base positions. In the non-coding regions, which are thought to be only loosely constrained by selection, transitions are observed more frequently than might be expected: they account for 52% of all substitutions, and they appear to be favored two to threefold over transversions when allowance is made for the skewed base composition of these regions.     

Linkage disequilibrium and allele-frequency distributions for 114 single-nucleotide polymorphisms in five populations

Single-nucleotide polymorphisms (SNPs) may be extremely important for deciphering the impact of genetic variation on complex human diseases. The ultimate value of SNPs for linkage and association mapping studies depends in part on the distribution of SNP allele frequencies and intermarker linkage disequilibrium (LD) across populations. Limited information is available about these distributions on a genomewide scale, particularly for LD. Using 114 SNPs from 33 genes, we compared these distributions in five American populations (727 individuals) of African, European, Chinese, Hispanic, and Japanese descent. The allele frequencies were highly correlated across populations but differed by >20% for at least one pair of populations in 35% of SNPs. The correlation in LD was high for some pairs of populations but not for others (e.g., Chinese American or Japanese American vs. any other population). Regardless of population, average minor-allele frequencies were significantly higher for SNPs in noncoding regions (20%-25%) than for SNPs in coding regions (12%-16%). Interestingly, we found that intermarker LD may be strongest with pairs of SNPs in which both markers are nonconservative substitutions, compared to pairs of SNPs where at least one marker is a conservative substitution. These results suggest that population differences and marker location within the gene may be important factors in the selection of SNPs for use in the study of complex disease with linkage or association mapping methods.     

Genetic diversity of tyrosine hydroxylase (TH) and dopamine β-hydroxylase (DBH) genes in cattle breeds

DNA from four cattle breeds was used to re-sequence all of the exons and 56% of the introns of the bovine tyrosine hydroxylase (TH) gene and 97% and 13% of the bovine dopamine β-hydroxylase (DBH) coding and non-coding sequences, respectively. Two novel single nucleotide polymorphisms (SNPs) and a microsatellite motif were found in the TH sequences. The DBH sequences contained 62 nucleotide changes, including eight non-synonymous SNPs (nsSNPs) that are of particular interest because they may alter protein function and therefore affect the phenotype. These DBH nsSNPs resulted in amino acid substitutions that were predicted to destabilize the protein structure. Six SNPs (one from TH and five from DBH non-synonymous SNPs) were genotyped in 140 animals; all of them were polymorphic and had a minor allele frequency of > 9%. There were significant differences in the intra- and inter-population haplotype distributions. The haplotype differences between Brahman cattle and the three B. t. taurus breeds (Charolais, Holstein and Lidia) were interesting from a behavioural point of view because of the differences in temperament between these breeds.     

Characteristics of the tomato nuclear genome as determined by sequencing undermethylated EcoRI digested fragments

A collection of 9,990 single-pass nuclear genomic sequences, corresponding to 5 Mb of tomato DNA, were obtained using methylation filtration (MF) strategy and reduced to 7,053 unique undermethylated genomic islands (UGIs) distributed as follows: (1) 59% non-coding sequences, (2) 28% coding sequences, (3) 12% transposons—96% of which are class I retroelements, and (4) 1% organellar sequences integrated into the nuclear genome over the past approximately 100 million years. A more detailed analysis of coding UGIs indicates that the unmethylated portion of tomato genes extends as far as 676 bp upstream and 766 bp downstream of coding regions with an average of 174 and 171 bp, respectively. Based on the analysis of the UGI copy distribution, the undermethylated portion of the tomato genome is determined to account for the majority of the unmethylated genes in the genome and is estimated to constitute 61±15 Mb of DNA (~5% of the entire genome)—which is significantly less than the 220 Mb estimated for gene-rich euchromatic arms of the tomato genome. This result indicates that, while most genes reside in the euchromatin, a significant portion of euchromatin is methylated in the intergenic spacer regions. Implications of the results for sequencing the genome of tomato and other solanaceous species are discussed.