A combined long-range phasing and long haplotype imputation method to impute phase for SNP genotypes

Background

Knowing the phase of marker genotype data can be useful in genome-wide association studies, because it makes it possible to use analysis frameworks that account for identity by descent or parent of origin of alleles and it can lead to a large increase in data quantities via genotype or sequence imputation. Long-range phasing and haplotype library imputation constitute a fast and accurate method to impute phase for SNP data.

Methods

A long-range phasing and haplotype library imputation algorithm was developed. It combines information from surrogate parents and long haplotypes to resolve phase in a manner that is not dependent on the family structure of a dataset or on the presence of pedigree information.

Results

The algorithm performed well in both simulated and real livestock and human datasets in terms of both phasing accuracy and computation efficiency. The percentage of alleles that could be phased in both simulated and real datasets of varying size generally exceeded 98% while the percentage of alleles incorrectly phased in simulated data was generally less than 0.5%. The accuracy of phasing was affected by dataset size, with lower accuracy for dataset sizes less than 1000, but was not affected by effective population size, family data structure, presence or absence of pedigree information, and SNP density. The method was computationally fast. In comparison to a commonly used statistical method (fastPHASE), the current method made about 8% less phasing mistakes and ran about 26 times faster for a small dataset. For larger datasets, the differences in computational time are expected to be even greater. A computer program implementing these methods has been made available.

Conclusions

The algorithm and software developed in this study make feasible the routine phasing of high-density SNP chips in large datasets.     

Nine new polymorphic microsatellite loci for the amplification of archived otolith DNA of Atlantic cod,Gadus morhua L.

Nine out of 22 microsatellite primers tested were successfully amplified on three samples of cod Gadus morhua L. (two contemporary and one archived otolith samples). All loci were polymorphic (5–23 alleles/locus). The average observed heterozygosity across loci and samples was 0.625, ranging from 0.294 to 0.895 at each locus. All loci were under Hardy–Weinberg equilibrium, except PGmo56 that showed significant excess of heterozygotes in all studied samples. The isolated loci were suitable for degraded DNA and therefore useful for conducting a long‐term temporal study with DNA obtained from archived otoliths of cod.     

Molecular mechanisms of homologous and heterologous desensitization mediated by vasopressin in smooth muscle cells

Arginine vasopressin (AVP) has been shown previously to enhance phosphatidylinositol (PI) turnover and mobilize calcium in the rat aortic smooth muscle cell-line (A10; ATCC CRL 1476) via the V1 receptor (Aiyar, N., Nambi, P., Stassen, F. L., and Crooke, S. T. (1986) Life Sci. 39, 37-45). Exposure of A10 cells to AVP for periods ranging from 5 min to 2 h resulted in 30-40% loss in AVP-binding sites and an inhibition of the production of inositol di- and trisphosphates and the mobilization of calcium when the cells were rechallenged by addition of AVP. We now report that during the same time course AVP induces a dose- and time-dependent decrease in labeled PI, phosphatidylinositol-4-phosphate and phosphatidylinositol-4,5-bisphosphate reaching a minimum after 30 min of incubation. After 2 h of exposure to AVP, the levels of labeled PI, phosphatidylinositol-4-phosphate and phosphatidylinositol-4,5-bisphosphate increased to a new basal level approximately 30% less than the untreated cultures. The decrease in inositol lipid labeling mediated by AVP was inhibited when the V1 antagonist SK&F 100273 was included in the incubations with AVP. No decrease was observed when the V2 agonist 1-deamino, [8-D-arginine]vasopressin was used for pretreatment of the cells. Furthermore, when PI kinase activity was measured in cell extracts from untreated and AVP-treated (2 h) cells a significant decrease (p less than 0.05) was observed in the absence, but not in the presence, of added PI in the AVP-treated cells as compared with the control cells. Thrombin also stimulates PI metabolism and calcium mobilization in these cells and brought about both a prolonged decrease in inositol lipids and inhibition of PI kinase activity. AVP pretreatment affected the release of intracellular Ca2+ induced by AVP, thrombin, and ATP, differently. The time of AVP pretreatment required to induce half-maximal inhibition of intracellular Ca2+ release in response to AVP, thrombin, and ATP was approximately 8, 24, and 30 min, respectively. Consequently, we suggest that the reduction in response to AVP with short term preincubation is due to homologous desensitization as reflected by 30-40% decrease in V1 receptors. Subsequently, a decrease in inositol lipid pools and PI kinase activity results in heterologous desensitization in response to AVP, thrombin, and ATP.     

Through regulation of TCR expression levels, an Idd7 region gene(s) interactively contributes to the impaired thymic deletion of autoreactive diabetogenic CD8+ T cells in nonobese diabetic mice

When expressed in NOD, but not C57BL/6 (B6) genetic background mice, the common class I variants encoded by the H2g7 MHC haplotype aberrantly lose the ability to mediate the thymic deletion of autoreactive CD8+ T cells contributing to type 1 diabetes (T1D). This indicated some subset of the T1D susceptibility (Idd) genes located outside the MHC of NOD mice interactively impair the negative selection of diabetogenic CD8+ T cells. In this study, using both linkage and congenic strain analyses, we demonstrate contributions from a polymorphic gene(s) in the previously described Idd7 locus on the proximal portion of Chromosome 7 predominantly, but not exclusively, determines the extent to which H2g7 class I molecules can mediate the thymic deletion of diabetogenic CD8+ T cells as illustrated using the AI4 TCR transgenic system. The polymorphic Idd7 region gene(s) appears to control events that respectively result in high vs low expression of the AI4 clonotypic TCR alpha-chain on developing thymocytes in B6.H2g7 and NOD background mice. This expression difference likely lowers levels of the clonotypic AI4 TCR in NOD, but not B6.H2g7 thymocytes, below the threshold presumably necessary to induce a signaling response sufficient to trigger negative selection upon Ag engagement. These findings provide further insight to how susceptibility genes, both within and outside the MHC, may interact to elicit autoreactive T cell responses mediating T1D development in both NOD mice and human patients.     

5′-O-Aliphatic and amino acid ester prodrugs of (−)-β-d-(2R,4R)-dioxolane-thymine (DOT): Synthesis,anti-HIV activity,cytotoxicity and stability studies

A series of (?)-β-d-(2R,4R)-dioxolane-thymine-5′-O-aliphatic acid esters as well as amino acid esters were synthesized as prodrugs of (?)-β-d-(2R,4R)-dioxolane-thymine (DOT). The compounds were evaluated for anti-HIV activity against HIV-1_LAI in human peripheral blood mononuclear (PBM) cells as well as for their cytotoxicity in PBM, CEM and Vero cells. Improved anti-HIV potency in vitro was observed for the compound 2–4 (5′-O-aliphatic acid esters) without increase in cytotoxicity in comparison to the parent drug. Chemical and enzymatic hydrolysis of the prodrugs was also studied, in which the prodrugs exhibited good chemical stability with the half-lives from 3 h to 54 h at pH 2.0 and 7.4 phosphate buffer. However, the prodrugs were relatively labile to porcine esterase with the half-lives from 12.3 to 48.0 min.     

Adverse environmental conditions influence age-related innate immune responsiveness

Hutchinson-Gilford progeria syndrome (HGPS) is a rare premature aging disorder that belongs to a group of conditions called laminopathies which affect nuclear lamins. Mutations in two genes, LMNA and ZMPSTE24, have been found in patients with HGPS. The p.G608G LMNA mutation is the most commonly reported mutation. The aim of this work was to compile a comprehensive literature review of the clinical features and genetic mutations and mechanisms of this syndrome as a contribution to health care workers. This review shows the necessity of a more detailed clinical identification of Hutchinson-Gilford progeria syndrome and the need for more studies on the pharmacologic and pharmacogenomic approach to this syndrome.     

2D difference gel electrophoresis of prepubertal and pubertal rat mammary gland proteomes

Rat mammary gland proteomes at day 21 (prepubertal) and day 50 (late puberty) were compared by 2D difference gel electrophoresis. Two-hundred fifty-one spots were significantly different ( p < 0.05) in abundance. Peptide mass fingerprint analysis of a subset of these proteins identified two significantly over-represented classes including structural and blood proteins (increased), and metabolism-relevant proteins (reduced) in day 50 relative to day 21 glands. This is a first report of mammary gland proteome differences at these important breast cancer-relevant time-points.     

Evaluating statistical methods using plasmode data sets in the age of massive public databases: an illustration using false discovery rates

Plasmode is a term coined several years ago to describe data sets that are derived from real data but for which some truth is known. Omic techniques, most especially microarray and genomewide association studies, have catalyzed a new zeitgeist of data sharing that is making data and data sets publicly available on an unprecedented scale. Coupling such data resources with a science of plasmode use would allow statistical methodologists to vet proposed techniques empirically (as opposed to only theoretically) and with data that are by definition realistic and representative. We illustrate the technique of empirical statistics by consideration of a common task when analyzing high dimensional data: the simultaneous testing of hundreds or thousands of hypotheses to determine which, if any, show statistical significance warranting follow-on research. The now-common practice of multiple testing in high dimensional experiment (HDE) settings has generated new methods for detecting statistically significant results. Although such methods have heretofore been subject to comparative performance analysis using simulated data, simulating data that realistically reflect data from an actual HDE remains a challenge. We describe a simulation procedure using actual data from an HDE where some truth regarding parameters of interest is known. We use the procedure to compare estimates for the proportion of true null hypotheses, the false discovery rate (FDR), and a local version of FDR obtained from 15 different statistical methods.