Haplotype analysis of the preprotachykinin-1 (TAC1) gene in multiple sclerosis

The potential relevance of chromosome 7q21-22 in susceptibility to multiple sclerosis (MS) has been highlighted in genome-wide linkage screens as well as in association studies of 7q-specific polymorphic microsatellites. Especially, recent, independently performed studies have provided evidence for significant association of the markers D7S554 and D7S3126 with MS in Sardinian, Northern Irish and Spanish-American cohorts. The gene most closely located to these markers is the neuropeptide preprotachykinin-1 (TAC1) gene. Both its position and the array of biological functions exerted by its expression products make it a logical primary choice for further scrutiny as the putative chromosome 7q21-22 MS susceptibility gene. We report identification of eight polymorphisms in this gene by means of a sequencing approach. A Northern Irish case-control was typed for six of these polymorphisms. One of these, an intron 1 single-nucleotide polymorphism (SNP), showed significant association with MS (P=0.009). Two-marker haplotypes composed of allelic combinations of TAC1 promoter-intron 1 SNPs were highly significantly associated with MS and more so with the relapsing-remitting form of this disease. While independent reproduction of these data in other data sets is indicated, our work is suggestive for a role of the TAC1 gene in MS. Genes and Immunity (2005) 6, 265-270. doi:10.1038/sj.gene.6364175 Published online 24 February 2005.     

Using text analysis to identify functionally coherent gene groups

The analysis of large-scale genomic information (such as sequence data or expression patterns) frequently involves grouping genes on the basis of common experimental features. Often, as with gene expression clustering, there are too many groups to easily identify the functionally relevant ones. One valuable source of information about gene function is the published literature. We present a method, neighbor divergence, for assessing whether the genes within a group share a common biological function based on their associated scientific literature. The method uses statistical natural language processing techniques to interpret biological text. It requires only a corpus of documents relevant to the genes being studied (e.g., all genes in an organism) and an index connecting the documents to appropriate genes. Given a group of genes, neighbor divergence assigns a numerical score indicating how "functionally coherent" the gene group is from the perspective of the published literature. We evaluate our method by testing its ability to distinguish 19 known functional gene groups from 1900 randomly assembled groups. Neighbor divergence achieves 79% sensitivity at 100% specificity, comparing favorably to other tested methods. We also apply neighbor divergence to previously published gene expression clusters to assess its ability to recognize gene groups that had been manually identified as representative of a common function.     

Haplotype analysis of the phenylalanine hydroxylase gene in Turkish phenylketonuria families

We have estimated the haplotype distribution of mutant and normal phenylalanine hydroxylase (PAH) alleles for 17 Turkish phenylketonuria (PKU) families: 20 normal and 27 mutated PAH alleles could be identified. Of the latter, the most prevalent were associated with haplotype 6 (29.6%), 1 (18.5%) and 36 (11.1%), while the normal alleles were preferentially associated with haplotype 1 (20%). Of the 19 different haplotypes observed, 5 have not been described previously. The haplotype distribution differed significantly from that of the Northern European population. Two of the eight polymorphic sites were in association with PKU. No deletions of exon sequences were found in the families analysed.     

Haplotype analysis of BRCA1 gene reveals a new gene rearrangement: characterization of a 19.9 KBP deletion

Germ-line mutations in the BRCA1 gene cause hereditary predisposition to breast and ovarian cancer. BRCA1 and BRCA2 mutations account for about 40% of high-risk families. Mutation-screening methods generally focus on genomic DNA and are usually PCR based; they enable the detection of sequence alterations such as point mutations and small deletions and insertions. However, they do not allow the detection of partial or entire exon(s) loss, because the presence of the homologous allele results in a positive PCR signal, giving rise to a false-negative result. Identification of unusual haplotypes in patient samples by an expectation maximization algorithm has recently been suggested as a method for identifying hemizygous regions caused by large intragenic deletions. Using a similar approach, we identified a novel BRCA1 genomic rearrangement in a breast/ovarian cancer family negative at the first mutation screening; we detected a deletion encompassing exons 14-19, probably due to replication slippage between Alu sequences.     

Haplotype analysis of single nucleotide polymorphisms in VEGF gene for vascular dementia.

A case-control study was performed to identify single nucleotide variants of vascular endothelial growth factor (VEGF) gene associated with vascular dementia. Seven SNPs in promoter, 5'-UTR, 3'-UTR, and introns of VEGF gene were identified in 24 Koreans. Three of them, -1154G/A, -7C/T, and 13553C/T, were selected based on allele frequency and linkage disequilibrium (LD), and genotyped in 207 vascular dementia patients and 207 control subjects. Significant association with vascular dementia was not shown (P > 0.05) in the alleles and genotypes of single locus. Subsequent analysis of composing VEGF risk haplotypes associated with vascular dementia was performed with maximum likelihood estimates of their possible haplotypes employing the expectation-maximization (EM) algorithm. Of three-locus haplotypes, only GTC was significantly associated with vascular dementia, conferring a risk of 1.87 (P < 0.05). Of two-locus haplotypes, the risk was observed with the nested forms of the risk haplotype GTC, that is, GT at the loci -1154G/A and -7C/T and TC at the loci -7C/T and 13553C/T (P < 0.05). Our findings suggested some interaction among -1154G/A, -7C/T, and 13553C/T variants in the determination of risk for vascular dementia.     

Haplotype analysis of the estrogen receptor 1 gene in male genital and reproductive abnormalities

BACKGROUND: We have recently suggested that homozygosity for a specific 'AGATA' haplotype within a approximately 50 kb linkage disequilibrium (LD) block of the gene for estrogen receptor alpha (ESR1) may raise the susceptibility to cryptorchidism by enhancing estrogenic effects of environmental endocrine disruptors (EEDs). METHODS: Haplotype analysis of ESR1 was performed in 328 Japanese subjects, i.e. 70 patients with micropenis (MP), 43 patients with hypospadias (HS), 80 patients with spermatogenic failure (SF) and 135 control males. Genotyping was performed by the 5' nuclease assay. RESULTS: The LD block was identified in each of the patient groups and in the control males. The frequency of homozygotes for the specific 'AGATA' haplotype was markedly higher in the HS patients [P = 0.0000033, odds ratio [OR] = 11.26] and slightly higher in the MP patients (P = 0.034, OR = 3.64) than in the control males, and the 'AGATA' haplotype was strongly associated with HS (P = 0.0000022, OR = 11.26) and weakly associated with MP (P = 0.040, OR = 3.64) in a recessive mode. There was no significant difference between the SF patients and the control males. CONCLUSIONS: Our results support the hypothesis that homozygosity for the specific ESR1 'AGATA' haplotype may increase the susceptibility to the development of male genital abnormalities in response to estrogenic EEDs.     

Haplotype and linkage disequilibrium analysis to characterise a region in the calcium channel gene CACNA1A associated with idiopathic generalised epilepsy

Idiopathic generalised epilepsy (IGE) is a common form of epilepsy, including several defined and overlapping syndromes, and likely to be due to the combined actions of mutations in several genes. In a recent study we investigated the calcium channel gene CACNA1A for involvement in IGE, unselected for syndrome, by means of association studies using several polymorphisms within the gene. We reported a highly significant case/control association with a silent single nucleotide polymorphism (SNP) in exon 8 that we confirmed by within-family analyses. In this present study we screened the gene for novel SNPs within 25 kb of exon 8, which have enabled us to define the critical region of CACNA1A in predisposing to IGE. Several intronic SNPs were identified and three, within 1.5 kb of exon 8 and in strong linkage disequilibrium with each other and with the original SNP, were significantly associated with IGE (P=0.00029, P=0.0015 and P=0.010). The associations were not limited to an IGE syndrome or other subgroup. Another SNP, 25 kb away, in intron 6 was also significantly associated with IGE (P=0.0057) but is not in linkage disequilibrium with the SNPs around exon 8. Haplotype predictions revealed even more significant associations (3-marker haplotype: P<10(-6)). Logistic regression showed that all the data can be explained by two of the SNPs, which is consistent with two functionally significant variants being responsible for all five associations, although a single variant cannot be excluded. The functionally significant variant(s) are unlikely to be exonic and suggests an effect on expression or alternative splicing.     

Haplotype analysis of the fragile X syndrome gene FMR1 in the Czech Republic.

We report on the haplotype analysis with polymorphic repeat markers DXS548 and FRAXAC1 next to the FMR1 gene in 37 unrelated fragile X and 36 control chromosomes from Bohemia and Moravia. Our results suggest a significant linkage disequilibrium between fragile X mutations and certain DXS548-FRAXAC1 haplotypes. Allele frequencies obtained differ slightly from those of other European populations with allele 194 being less frequent in our control sample. Rare DXS548 alleles 6.5 (195) and 0 (208) were also present.     

Haplotype analysis of IL-10 gene polymorphism in couples with spontaneous abortions and aborted fetuses

Spontaneous abortion is the loss of pregnancy during an early gestational period. Interleukin-10 is an anti-inflammatory cytokine which plays an important role in successful pregnancy outcome. The aim of the study is to elucidate an association of IL-10 gene promoter polymorphisms (-1082G/A, -819 C/T, -592C/A) in spontaneous abortions from Telangana state of South India. The present population-based retrospective case-control triad study includes a total of 80 case families with spontaneous abortions and 100 control families with medically terminated pregnancies. Peripheral blood from all the couples and fetal tissues of <20 weeks of gestation were collected. Genotype analysis was carried out by a standard amplification refractory mutation system–polymerase chain reaction followed by agarose gel electrophoresis. The strength of the association between IL-10 gene promoter polymorphisms and spontaneous abortions were measured by odd ratios and their respective 95% confidence intervals. Haplotype analysis was carried out for the three polymorphisms to establish an association of specific haplotypes with spontaneous abortions. The increased frequency of AA genotype and A allele of -1082G/A, TT genotype and T allele of -819C/T, and AA genotype and A allele of -592C/A was observed in case fetuses and case mothers compared to their respective controls. Haplotype analysis revealed that A-C-A, G-C-A haplotypes in fetuses and haplotypes A-C-C, G-T-C, A-T-A, and G-C-A in mothers were associated with increased risk of spontaneous abortions. IL-10 gene promoter polymorphisms may act as a major genetic regulator in the etiology of spontaneous abortions with maternal genome imprinting effects.     

Haplotype analysis to determine the position of a mutation among closely linked DNA markers

Positional cloning involves first finding linkage between aninherited phenotype (such as a disease) and a DNA marker, followedby the use of a variety of physical and genetic mapping techniquesto move from linkage to mutation. If there is a founder effectwithin a population, crossovers are often rare between the mutationcausing the phenotype and closely situated markers and increasingdisequilibrium may be observed as the site of the mutation isapproached. Standard coefficients of disequilibrium may, however,be insensitive to the relative position of close markers andthe mutation, because they depend upon allele frequencies inthe normal population compared to those of the founder chromosome.Using cystic fibrosis in European populations as a model system,alternative methods for determining the position of a mutationare discussed. These include hapiotype parsimony and three-wayinterval likelihood analysis. Both methods predict the locationof the major CF mutation accurately from a real set of morethan 600 European CF chromosomes.     

Haplotype tagging efficiency in worldwide populations in CTLA4 gene

The cytotoxic T lymphocyte antigen 4 (CTLA4) acts as a potent negative regulator of T-cell response, and has been suggested as a pivotal candidate gene for autoimmune disorders such as Graves' disease, type 1 diabetes and autoimmune hypothyroidism, among others. Several single-nucleotide polymorphisms (SNPs) have been proposed as the susceptibility variants, or to be in strong linkage disequilibrium (LD) with the variant. Nevertheless, contradictory results have been found, which may be due to lack of knowledge of the genetic structure of CTLA4 and its geographic variation. We have typed 17 SNPs throughout the CTLA4 gene region in order to analyze the haplotype diversity and LD structure in a worldwide population set (1262 individuals from 44 populations) to understand the variation pattern of the region. Allele and haplotype frequency differentiation between populations is consistent with genomewide averages and points to a lack of strong population-specific selection pressures. LD is high and its pattern is not significantly different within or between continents. However, haplotype composition is significantly different between geographical groups. A continent-specific set of haplotype tagging SNPs has been designed to be used for future association studies. These are portable among populations, although their efficiency might vary depending on the population haplotype spectrum.     

Haplotype analysis for CF-linked DNA polymorphisms in Switzerland

A total of 295 patients, parents and unaffected sibs from 106 CF-families in central and northeastern Switzerland were investigated with probes 7C22(D7S16), metH, metD, pKM19, pXV-2c and pJ3.11(D7S8) for eight DNA polymorphisms (RFLP's). Linkage disequilibrium to the CF locus and haplotype frequencies were compared to those in other populations. They are comparable to other Caucasian populations and, for pKM 19 and pXV-2c, very close to the findings in Italy. The prevalence of certain haplotypes among the CF and the normal allele-bearing chromosomes indicate that the majority of the CF cases are probably the result of one ancient mutation in a common ancestor, but that there may be allelic heterogeneity accounting for an important proportion of patients, that may differ between countries or regions. Informative family constellations for the different polymorphisms in Switzerland and strategies for carrier detection and prenatal diagnosis are discussed. Haplotype analyses for each country and its ethnic subgroups are recommended.     

Multivariate QTL linkage analysis suggests a QTL for platelet count on chromosome 19q

Platelet count is a highly heritable trait with genetic factors responsible for around 80% of the phenotypic variance. We measured platelet count longitudinally in 327 monozygotic and 418 dizygotic twin pairs at 12, 14 and 16 years of age. We also performed a genome-wide linkage scan of these twins and their families in an attempt to localize QTLs that influenced variation in platelet concentrations. Suggestive linkage was observed on chromosome 19q13.13-19q13.31 at 12 (LOD = 2.12, P = 0.0009), 14 (LOD = 2.23, P = 0.0007) and 16 (LOD = 1.01, P = 0.016) years of age and multivariate analysis of counts at all three ages increased the LOD to 2.59 (P = 0.0003). A possible candidate in this region is the gene for glycoprotein VI, a receptor involved in platelet aggregation. Smaller linkage peaks were also seen at 2p, 5p, 5q, 10p and 15q. There was little evidence for linkage to the chromosomal regions containing the genes for thrombopoietin (3q27) and the thrombopoietin receptor (1q34), suggesting that polymorphisms in these genes do not contribute substantially to variation in platelet count between healthy individuals.     

GenBlastA: enabling BLAST to identify homologous gene sequences

BLAST is an extensively used local similarity search tool for identifying homologous sequences. When a gene sequence (either protein sequence or nucleotide sequence) is used as a query to search for homologous sequences in a genome, the search results, represented as a list of high-scoring pairs (HSPs), are fragments of candidate genes rather than full-length candidate genes. Relevant HSPs ("signals"), which represent candidate genes in the target genome sequences, are buried within a report that contains also hundreds to thousands of random HSPs ("noises"). Consequently, BLAST results are often overwhelming and confusing even to experienced users. For effective use of BLAST, a program is needed for extracting relevant HSPs that represent candidate homologous genes from the entire HSP report. To achieve this goal, we have designed a graph-based algorithm, genBlastA, which automatically filters HSPs into well-defined groups, each representing a candidate gene in the target genome. The novelty of genBlastA is an edge length metric that reflects a set of biologically motivated requirements so that each shortest path corresponds to an HSP group representing a homologous gene. We have demonstrated that this novel algorithm is both efficient and accurate for identifying homologous sequences, and that it outperforms existing approaches with similar functionalities.     

Candidate gene analysis of SPARCL1 gene in patients with multiple sclerosis

Recently, proteomic analysis in cerebrospinal fluid (CSF) from patients with MS identified four proteins which are present in MS but not in normal human CSF, including SPARCL1, an extracellular matrix-associated protein member of the SPARC family. One hundred eighty-six patients with MS and 185 age-matched controls were genotyped for A/G single nucleotide polymorphism (SNP) in exon 1 (rs1049539), C/G SNP in exon 4 (rs1049544), resulting in a substitution of an aspartate with an histidine, and A/G substitution in the exon 5 (rs1130643), leading to the substitution of alanine with threonine. No significant differences in either allelic or genotypic frequency of the three SNPs were found (P>0.05), even in stratifying MS patients according to the course of the disease. Stratifying according to gender, a trend towards a decreased frequency of the C/C genotype of the rs1049544 was observed in male patients as compared with male controls (30.2% versus 44.0%; P=0.217). Despite proteomic studies in CSF from MS patients suggested an important role for SPARCL1 in the development of the disease, SPARCL1 gene does not appear to act as susceptibility factor for MS in the population investigated here. However, the frequency of the C/C genotype of rs1049544 was decreased in male patients, possibly conferring a lower risk of developing MS in male population. Further studies are needed to clarify this issue.