Cubic exact solutions for the estimation of pairwise haplotype frequencies: implications for linkage disequilibrium analyses and a web tool 'CubeX'

<p>Abstract</p> <p>Background</p> <p>The frequency of a haplotype comprising one allele at each of two loci can be expressed as a cubic equation (the 'Hill equation'), the solution of which gives that frequency. Most haplotype and linkage disequilibrium analysis programs use iteration-based algorithms which substitute an estimate of haplotype frequency into the equation, producing a new estimate which is repeatedly fed back into the equation until the values converge to a maximum likelihood estimate (expectation-maximisation).</p> <p>Results</p> <p>We present a program, "CubeX", which calculates the biologically possible exact solution(s) and provides estimated haplotype frequencies, D', r²and <it>χ</it>²values for each. CubeX provides a "complete" analysis of haplotype frequencies and linkage disequilibrium for a pair of biallelic markers under situations where sampling variation and genotyping errors distort sample Hardy-Weinberg equilibrium, potentially causing more than one biologically possible solution. We also present an analysis of simulations and real data using the algebraically exact solution, which indicates that under perfect sample Hardy-Weinberg equilibrium there is only one biologically possible solution, but that under other conditions there may be more.</p> <p>Conclusion</p> <p>Our analyses demonstrate that lower allele frequencies, lower sample numbers, population stratification and a possible |D'| value of 1 are particularly susceptible to distortion of sample Hardy-Weinberg equilibrium, which has significant implications for calculation of linkage disequilibrium in small sample sizes (eg HapMap) and rarer alleles (eg paucimorphisms, q < 0.05) that may have particular disease relevance and require improved approaches for meaningful evaluation.</p

A genetic algorithm based method for stringent haplotyping of family data

<p>Abstract</p> <p>Background</p> <p>The linkage phase, or haplotype, is an extra level of information that in addition to genotype and pedigree can be useful for reconstructing the inheritance pattern of the alleles in a pedigree, and computing for example Identity By Descent probabilities. If a haplotype is provided, the precision of estimated IBD probabilities increases, as long as the haplotype is estimated without errors. It is therefore important to only use haplotypes that are strongly supported by the available data for IBD estimation, to avoid introducing new errors due to erroneous linkage phases.</p> <p>Results</p> <p>We propose a genetic algorithm based method for haplotype estimation in family data that includes a stringency parameter. This allows the user to decide the error tolerance level when inferring parental origin of the alleles. This is a novel feature compared to existing methods for haplotype estimation. We show that using a high stringency produces haplotype data with few errors, whereas a low stringency provides haplotype estimates in most situations, but with an increased number of errors.</p> <p>Conclusion</p> <p>By including a stringency criterion in our haplotyping method, the user is able to maintain the error rate at a suitable level for the particular study; one can select anything from haplotyped data with very small proportion of errors and a higher proportion of non-inferred haplotypes, to data with phase estimates for every marker, when haplotype errors are tolerable. Giving this choice makes the method more flexible and useful in a wide range of applications as it is able to fulfil different requirements regarding the tolerance for haplotype errors, or uncertain marker-phases.</p

Intriguing Balancing Selection on the Intron 5 Region of LMBR1 in Human Population

Background: The intron 5 of gene LMBR1 is the cis-acting regulatory module for the sonic hedgehog (SHH) gene. Mutation in this non-coding region is associated with preaxial polydactyly, and may play crucial roles in the evolution of limb and skeletal system. Methodology/Principal Findings: We sequenced a region of the LMBR1 gene intron 5 in East Asian human population, and found a significant deviation of Tajima’s D statistics from neutrality taking human population growth into account. Data from HapMap also demonstrated extended linkage disequilibrium in the region in East Asian and European population, and significantly low degree of genetic differentiation among human populations. Conclusion/Significance: We proposed that the intron 5 of LMBR1 was presumably subject to balancing selection during the evolution of modern human

Prognostic impact of urokinase-type plasminogen activator receptor (uPAR) in cytosols and pellet extracts derived from primary breast tumours

Using a previously developed enzyme-linked immunosorbent assay (ELISA), the levels of the receptor for urokinase-type plasminogen activator (uPAR) were determined in cytosols and corresponding membrane pellets derived from 878 primary breast tumours. The levels of uPAR in the pellet extracts were more than 3-fold higher than those measured in the cytosols (P< 0.001). Moreover, the uPAR levels in the two types of extracts were weakly, though significantly, correlated with each other (rS= 0.20, P< 0.001). In Cox univariate analysis, high cytosolic levels of uPAR were significantly associated with reduced overall survival (OS) and relapse-free survival (RFS). The levels of uPAR in pellet extracts appeared not to be related with patient survival. In multivariate analysis, elevated levels of uPAR measured in cytosols and pellet extracts were found to be independent predictors of poor OS, not RFS. The prediction of poor prognosis on the basis of high uPAR levels emphasizes its important role in plasmin-mediated degradation of extracellular matrix proteins during cancer invasion and metastasis. © 2001 Cancer Research Campaign http://www.bjcancer.co

Signaling pathways downstream of P2 receptors in human neutrophils

Extracellular nucleotides stimulate human neutrophils by activating the purinergic P2Y2 receptor. However, it is not completely understood which types of G proteins are activated downstream of this P2 receptor subtype. We investigated the G-protein coupling to P2Y2 receptors and several subsequent signaling events. Treatment of neutrophils with pertussis toxin (PTX), a Gi protein inhibitor, caused only ∼75% loss of nucleotide-induced Ca2+ mobilization indicating that nucleotides cause Ca2+ mobilization both through Gi-dependent and Gi-independent pathways. However, the PLC inhibitor U73122 almost completely inhibited Ca2+ mobilization in both nucleotide- and fMLP-stimulated neutrophils, strongly supporting the view that both the PTX-sensitive and the PTX-insensitive mechanism of Ca2+ increase require activation of PLC. We investigated the dependence of ERK phosphorylation on the Gi pathway. Treatment of neutrophils with PTX caused almost complete inhibition of ERK phosphorylation in nucleotide or fMLP activated neutrophils. U73122 caused inhibition of nucleotide- or fMLP-stimulated ERK phosphorylation, suggesting that although pertussis toxin-insensitive pathways cause measurable Ca2+ mobilization, they are not sufficient for causing ERK phosphorylation. Since PLC activation leads to intracellular Ca2+ increase and PKC activation, we investigated if these intracellular events are necessary for ERK phosphorylation. Exposure of cells to the Ca2+ chelator BAPTA had no effect on nucleotide- or fMLP-induced ERK phosphorylation. However, the PKC inhibitor GF109203X was able to almost completely inhibit nucleotide- or fMLP-induced ERK phosphorylation. We conclude that the P2Y2 receptor can cause Ca2+ mobilization through a PTX-insensitive but PLC-dependent pathway and ERK phosphorylation is highly dependent on activation of the Gi proteins

The Use of Haplotypes in the Identification of Interaction between SNPs

Although haplotypes can provide great insight into the complex relationships between functional polymorphisms at a locus, their use in modern association studies has been limited. This is due to our inability to directly observe haplotypes in studies of unrelated individuals, but also to the extra complexity involved in their analysis and the difficulty in identifying which is the truly informative haplotype. Using a series of simulations, we tested a number of different models of a haplotype carrying two functional single nucleotide polymorphisms (SNPs) to assess the ability of haplotypic analysis to identify functional interactions between SNPs at the same locus. We found that, when phase is known, analysis of the haplotype is more powerful than analysis of the individual SNPs. The difference between the two approaches becomes less either as an increasing number of non-informative SNPs are included, or when the haplotypic phase is unknown, while in both cases the SNP association becomes progressively better at identifying the association. Our results suggest that when novel genotyping and bioinformatics methods are available to reconstruct haplotypic phase, this will permit the emergence of a new wave of haplotypic analysis able to consider interactions between SNPs with increased statistical power.</p

Predictors of gastrointestinal lesions on endoscopy in iron deficiency anemia without gastrointestinal symptoms

<p>Abstract</p> <p>Background</p> <p>Iron deficiency anaemia (IDA) due to occult gastrointestinal (GI) blood loss usually remains unnoticed until patient become symptomatic. There is sparse data in IDA patients without gastrointestinal symptoms. This study was designed to find out the frequency and predictors of endoscopic lesions in IDA without gastrointestinal symptoms. Cross-sectional study performed on a convenience sample of consecutive subjects.</p> <p>Methods</p> <p>Ninety five consecutive patients with laboratory based diagnosis of IDA having no gastrointestinal symptoms were interviewed and their clinical and biochemical variables were recorded. All the study patients underwent esophago-gastroduodenoscopy (EGD) and colonoscopy. Endoscopic findings were documented as presence/absence of bleeding related lesion and presence/absence of cause of IDA. Multiple logistic regressions were performed to identify variables significantly related to outcome variables.</p> <p>Results</p> <p>Possible cause of anaemia was found in 71% and bleeding related lesions were found in 53% of patients. Upper gastrointestinal tract lesions were found in 41% of patients with bleeding related lesions. On multivariable logistic regression; advancing age, low mean corpuscular volume (MCV ≤ 60 fl), and positive fecal occult blood test were predictive factors for bleeding related GI lesions and cause of IDA</p> <p>Conclusion</p> <p>Clinical and Biochemical markers can predict gastrointestinal lesions on endoscopy in IDA patients without gastrointestinal symptoms. High proportion of upper gastrointestinal involvement warrants EGD as initial endoscopic procedure however, this needs validation by further studies.</p

Examination of NRCAM, LRRN3, KIAA0716, and LAMB1 as autism candidate genes

BACKGROUND: A substantial body of research supports a genetic involvement in autism. Furthermore, results from various genomic screens implicate a region on chromosome 7q31 as harboring an autism susceptibility variant. We previously narrowed this 34 cM region to a 3 cM critical region (located between D7S496 and D7S2418) using the Collaborative Linkage Study of Autism (CLSA) chromosome 7 linked families. This interval encompasses about 4.5 Mb of genomic DNA and encodes over fifty known and predicted genes. Four candidate genes (NRCAM, LRRN3, KIAA0716, and LAMB1) in this region were chosen for examination based on their proximity to the marker most consistently cosegregating with autism in these families (D7S1817), their tissue expression patterns, and likely biological relevance to autism. METHODS: Thirty-six intronic and exonic single nucleotide polymorphisms (SNPs) and one microsatellite marker within and around these four candidate genes were genotyped in 30 chromosome 7q31 linked families. Multiple SNPs were used to provide as complete coverage as possible since linkage disequilibrium can vary dramatically across even very short distances within a gene. Analyses of these data used the Pedigree Disequilibrium Test for single markers and a multilocus likelihood ratio test. RESULTS: As expected, linkage disequilibrium occurred within each of these genes but we did not observe significant LD across genes. None of the polymorphisms in NRCAM, LRRN3, or KIAA0716 gave p < 0.05 suggesting that none of these genes is associated with autism susceptibility in this subset of chromosome 7-linked families. However, with LAMB1, the allelic association analysis revealed suggestive evidence for a positive association, including one individual SNP (p = 0.02) and three separate two-SNP haplotypes across the gene (p = 0.007, 0.012, and 0.012). CONCLUSIONS: NRCAM, LRRN3, KIAA0716 are unlikely to be involved in autism. There is some evidence that variation in or near the LAMB1 gene may be involved in autism

GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers

We describe methods with enhanced power and specificity to identify genes targeted by somatic copy-number alterations (SCNAs) that drive cancer growth. By separating SCNA profiles into underlying arm-level and focal alterations, we improve the estimation of background rates for each category. We additionally describe a probabilistic method for defining the boundaries of selected-for SCNA regions with user-defined confidence. Here we detail this revised computational approach, GISTIC2.0, and validate its performance in real and simulated datasets