Computational problems in perfect phylogeny haplotyping: typing without calling the allele.

A haplotype is an m-long binary vector. The XOR-genotype of two haplotypes is the m-vector of their coordinate-wise XOR. We study the following problem: Given a set of XOR-genotypes, reconstruct their haplotypes so that the set of resulting haplotypes can be mapped onto a perfect phylogeny (PP) tree. The question is motivated by studying population evolution in human genetics, and is a variant of the perfect phylogeny haplotyping problem that has received intensive attention recently. Unlike the latter problem, in which the input is "full" genotypes, here we assume less informative input, and so may be more economical to obtain experimentally. Building on ideas of Gusfield, we show how to solve the problem in polynomial time, by a reduction to the graph realization problem. The actual haplotypes are not uniquely determined by that tree they map onto, and the tree itself may or may not be unique. We show that tree uniqueness implies uniquely determined haplotypes, up to inherent degrees of freedom, and give a sufficient condition for the uniqueness. To actually determine the haplotypes given the tree, additional information is necessary. We show that two or three full genotypes suffice to reconstruct all the haplotypes, and present a linear algorithm for identifying those genotypes.

The novel ss469415590 variant predicts virological response to therapy in patients with chronic hepatitis C virus type 1 infection.

BACKGROUND: A novel dinucleotide variant TT/∆G (ss469415590) has been associated with hepatitis C virus clearance. AIM: To assess the role of the ss469415590 variant, compared with the known IL28B polymorphisms (rs8099917, rs12979860 and rs12980275) for predicting virological response to therapy in chronic hepatitis C, and its association with the CXCL10 chemokine serum levels - a surrogate marker of interferon-stimulated genes activation. METHODS: Multivariate analysis of factors predicting rapid and sustained virological response in 280 consecutive, treatment-naïve, nondiabetic, Caucasian patients with chronic hepatitis C treated with peginterferon alpha and ribavirin. RESULTS: In hepatitis C virus genotype 1, the OR (95% CI) for rapid and sustained virological response for the wild-type ss469415590 TT was 9.88 (1.99-48.99) and 7.25 (1.91-27.51), respectively, similar to those found for rs12979860 CC [9.55 (1.93-47.37) and 6.30 (1.71-23.13)] and for rs12980275 AA [9.62 (1.94-47.77] and 7.83 (2.02-30.34)], but higher than for rs8099917 TT [4.8 (1.73-13.33) and 4.75 (2.05-10.98)]. In hepatitis C virus genotype 1, mean (SD) CXCL10 levels in patients with the TT/TT, TT/∆G and ∆G/∆G variants were, respectively, 355.1 (240.6), 434.4 (247.4) and 569.9 (333.3) (P = 0.04). In patients with genotypes 2 and 3 no significant association was found for TT/∆G with viral response. The predictive value of ss469415590 was stronger in patients with advanced fibrosis. CONCLUSIONS: The novel IL28B variants at marker ss469415590 predict response to IFN alpha in chronic hepatitis C patients, especially in those with advanced fibrosis. Their determination may be superior to that of known IL28B variants for patient management using IFN-based regimens.

MRI monitoring of focal cerebral ischemia in peroxisome proliferator-activated receptor (PPAR)-deficient mice.

Peroxisome proliferator-activated receptors (PPARs) are a potential target for neuroprotection in focal ischemic stroke. These nuclear receptors have major effects in lipid metabolism, but they are also involved in inflammatory processes. Three PPAR isotypes have been identified: alpha, beta (or delta) and gamma. The development of PPAR transgenic mice offers a promising tool for prospective therapeutic studies. This study used MRI to assess the role of PPARalpha and PPARbeta in the development of stroke. Permanent middle cerebral artery occlusion induced focal ischemia in wild-type, PPARalpha-null mice and PPARbeta-null mice. T(2)-weighted MRI was performed with a 7 T MRI scan on day 0, 1, 3, 7 and 14 to monitor lesion growth in the various genotypes. General Linear Model statistical analysis found a significant difference in lesion volume between wild-type and PPAR-null mice for both alpha and beta isotypes. These data validate high-resolution MRI for monitoring cerebral ischemic lesions, and confirm the neuroprotective role of PPARalpha and PPARbeta in the brain.

Elderly age is not a negative predictive factor for virological response to therapy with pegylated interferon-α and ribavirin in chronic hepatitis C virus patients.

BACKGROUND & AIMS: Age is frequently discussed as negative host factor to achieve a sustained virological response (SVR) to antiviral therapy of chronic hepatitis C. However, elderly patients often show advanced fibrosis/cirrhosis as known negative predictive factor. The aim of this study was to assess age as an independent predictive factor during antiviral therapy. METHODS: Overall, 516 hepatitis C patients were treated with pegylated interferon-α and ribavirin, thereof 66 patients ≥60 years. We analysed the impact of host factors (age, gender, fibrosis, haemoglobin, previous hepatitis C treatment) and viral factors (genotype, viral load) on SVR per therapy course by performing a generalized estimating equations (GEE) regression modelling, a matched pair analysis and a classification tree analysis. RESULTS: Overall, SVR per therapy course was 42.9 and 26.1%, respectively, in young and elderly patients with hepatitis C virus (HCV) genotypes 1/4/6. The corresponding figures for HCV genotypes 2/3 were 74.4 and 84%. In the GEE model, age had no significant influence on achieving SVR. In matched pair analysis, SVR was not different in young and elderly patients (54.2 and 55.9% respectively; P = 0.795 in binominal test). In classification tree analysis, age was not a relevant splitting variable. CONCLUSIONS: Age is not a significant predictive factor for achieving SVR, when relevant confounders are taken into account. As life expectancy in Western Europe at age 60 is more than 20 years, it is reasonable to treat chronic hepatitis C in selected elderly patients with relevant fibrosis or cirrhosis but without major concomitant diseases, as SVR improves survival and reduces carcinogenesis.

Dyslipidemia in HIV-infected individuals: from pharmacogenetics to pharmacogenomics.

HIV-infected individuals may have accelerated atherogenesis and an increased risk for premature coronary artery disease. Dyslipidemia represents a key pro-atherogenic mechanism. In HIV-infected patients, dyslipidemia is typically attributed to the adverse effects of antiretroviral therapy. Nine recent genome-wide association studies have afforded a comprehensive, unbiased inventory of common SNPs at 36 genetic loci that are reproducibly associated with dyslipidemia in the general population. Genome-wide association study-validated SNPs have now been demonstrated to contribute to dyslipidemia in the setting of HIV infection and antiretroviral therapy. In a Swiss HIV-infected study population, a similar proportion of serum lipid variability was explained by antiretroviral therapy and by genetic background. In the individual patient, both antiretroviral therapy and the cumulative effect of SNPs contribute to the risk of high low-density lipoprotein cholesterol, low high-density lipoprotein cholesterol and hypertriglyceridemia. Genetic variants presumably contribute to additional major metabolic complications in HIV-infected individuals, including diabetes mellitus and coronary artery disease. In an effort to explain an increasing proportion of the heritability of complex metabolic traits, ongoing large-scale gene resequencing studies are focusing on the effects of rare SNPs and structural genetic variants.

Adaptive diversity in heterogeneous environments for populations regulated by a mixture of soft and hard selection

The stable co-existence of two haploid genotypes or two species is studied in a spatially heterogeneous environment submitted to a mixture of soft selection (within-patch regulation) and hard selection (outside-patch regulation) and where two kinds of resource are available. This is analysed both at an ecological time-scale (short term) and at an evolutionary time-scale (long term). At an ecological scale, we show that co-existence is very unlikely if the two competitors are symmetrical specialists exploiting different resources. In this case, the most favourable conditions are met when the two resources are equally available, a situation that should favour generalists at an evolutionary scale. Alternatively, low within-patch density dependence (soft selection) enhances the co-existence between two slightly different specialists of the most available resource. This results from the opposing forces that are acting in hard and soft regulation modes. In the case of unbalanced accessibility to the two resources, hard selection favours the most specialized genotype, whereas soft selection strongly favours the less specialized one. Our results suggest that competition for different resources may be difficult to demonstrate in the wild even when it is a key factor in the maintenance of adaptive diversity. At an evolutionary scale, a monomorphic invasive evolutionarily stable strategy (ESS) always exists. When a linear trade-off exists between survival in one habitat versus that in another, this ESS lies between an absolute adjustment of survival to niche size (for mainly soft-regulated populations) and absolute survival (specialization) in a single niche (for mainly hard-regulated populations). This suggests that environments in agreement with the assumptions of such models should lead to an absence of adaptive variation in the long term.

Role of hepatitis C virus genotype 3 in liver fibrosis progression: a systematic review and meta-analysis.

The progression of liver fibrosis in chronic hepatitis C has long been considered to be independent from viral genotypes. However, recent studies suggest an association between Hepatitis C virus (HCV) genotype 3 and accelerated liver disease progression. We completed a systematic review and meta-analysis of studies evaluating the association between HCV genotypes and fibrosis progression. PubMed, Embase and ISI Web of Knowledge databases were searched for cohort, cross-sectional and case-control studies on treatment-naïve HCV-infected adults in which liver fibrosis progression rate (FPR) was assessed by the ratio of fibrosis stage in one single biopsy to the duration of infection (single-biopsy studies) or from the change in fibrosis stage between two biopsies (paired biopsies studies). A random effect model was used to derive FPR among different HCV genotypes. Eight single-biopsy studies (3182 patients, mean/median duration of infection ranging from 9 to 21 years) and eight paired biopsies studies (mean interval between biopsies 2-12 years) met the selection criteria. The odds ratio for the association of genotype 3 with accelerated fibrosis progression was 1.52 (95% CI 1.12-2.07, P = 0.007) in single-biopsy studies and 1.37 (95% CI 0.87-2.17, P = 0.17) in paired biopsy studies. In conclusion, viral genotype 3 was associated with faster fibrosis progression in single-biopsy studies. This observation may have important consequences on the clinical management of genotype 3-infected patients. The association was not significant in paired biopsies studies, although the latter may be limited by important indication bias, short observation time and small sample size.

Treatment of Recurrent Hepatitis C after Liver Transplantation with PEG-Interferon and Ribavirin.

Background and aim: Recurrent hepati tis C is a major cause of morbidity and mortality after li ver transpl ant ati on (LT), and optimal treatm ent algorithms have yet to be defined. Here, we present our experience of 22 patients with recurrent hepatitis C treated in our institution .Patients and methods: Twenty-two patients with hi stology-proven recurrent hepati tis Cafter LT were treated since 2003. Treatment was ini ti ated with pegylated interferon-a2a 135 IIg per week and ribavirin 400 mg per day in the majority of patients, and subsequent doses were adapted individllally based on on-treatment virologieal responses and c1inical and/or biochemical si de effeets.Results: On an intention-to-treat basis, ustained virological re ponse(SVR) was achieved in 12/21 (54.5%) patie nts (5/12 [41 .6%], 2/3 [67%], 4/5 [80%] and 1/2 [50%] of patients infected with genotypes 1,2,3 and 4, respectively). Two patients experieneed relap e and 6 did not respond to treatm ent (NR). Treatment duration ranged from 24 to 90 weeks. It was stopped prematurely due to adverse events in 6/22 (27.2%) patients (with SVR achieved in 2 patients, NR in 2 patients, and death of 2 patients: one patient awaiting retransplantation and a second patient with HCV-HJV co-infection and fibrosing cholestat ic hepatiti s, nine months after transplantation). Of note, SVR was achi eved in a patient \Vi th combined liver and kidney transplantation. Importantly, SVR \Vas ach ieved in some patients despite the lack ofan early virological response or HCV RNA negativity at week 24. Darbepoetin a and fil ~,'rasti m were used in 36% and 18%, respectively.Conclusion: Individually adapted treatment of recurrent hepatitis C canachieve SVR in a substantial proponion ofLT patients. Conventional stopping rules do not apply in this setting so that prolonged therapy may be useful in selected patients.

Association of genetic variation with systolic and diastolic blood pressure among African Americans: the Candidate Gene Association Resource study.

The prevalence of hypertension in African Americans (AAs) is higher than in other US groups; yet, few have performed genome-wide association studies (GWASs) in AA. Among people of European descent, GWASs have identified genetic variants at 13 loci that are associated with blood pressure. It is unknown if these variants confer susceptibility in people of African ancestry. Here, we examined genome-wide and candidate gene associations with systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the Candidate Gene Association Resource (CARe) consortium consisting of 8591 AAs. Genotypes included genome-wide single-nucleotide polymorphism (SNP) data utilizing the Affymetrix 6.0 array with imputation to 2.5 million HapMap SNPs and candidate gene SNP data utilizing a 50K cardiovascular gene-centric array (ITMAT-Broad-CARe [IBC] array). For Affymetrix data, the strongest signal for DBP was rs10474346 (P= 3.6 × 10(-8)) located near GPR98 and ARRDC3. For SBP, the strongest signal was rs2258119 in C21orf91 (P= 4.7 × 10(-8)). The top IBC association for SBP was rs2012318 (P= 6.4 × 10(-6)) near SLC25A42 and for DBP was rs2523586 (P= 1.3 × 10(-6)) near HLA-B. None of the top variants replicated in additional AA (n = 11 882) or European-American (n = 69 899) cohorts. We replicated previously reported European-American blood pressure SNPs in our AA samples (SH2B3, P= 0.009; TBX3-TBX5, P= 0.03; and CSK-ULK3, P= 0.0004). These genetic loci represent the best evidence of genetic influences on SBP and DBP in AAs to date. More broadly, this work supports that notion that blood pressure among AAs is a trait with genetic underpinnings but also with significant complexity.

Fatty acid amide hydrolase deficiency enhances intraplaque neutrophil recruitment in atherosclerotic mice.

OBJECTIVE: Endocannabinoid levels are elevated in human and mouse atherosclerosis, but their causal role is not well understood. Therefore, we studied the involvement of fatty acid amide hydrolase (FAAH) deficiency, the major enzyme responsible for endocannabinoid anandamide degradation, in atherosclerotic plaque vulnerability. METHODS AND RESULTS: We assessed atherosclerosis in apolipoprotein E-deficient (ApoE(-/-)) and ApoE(-/-)FAAH(-/-) mice. Before and after 5, 10, and 15 weeks on high-cholesterol diet, we analyzed weight, serum cholesterol, and endocannabinoid levels, and atherosclerotic lesions in thoracoabdominal aortas and aortic sinuses. Serum levels of FAAH substrates anandamide, palmitoylethanolamide (PEA), and oleoylethanolamide (OEA) were 1.4- to 2-fold higher in case of FAAH deficiency. ApoE(-/-)FAAH(-/-) mice had smaller plaques with significantly lower content of smooth muscle cells, increased matrix metalloproteinase-9 expression, and neutrophil content. Circulating and bone marrow neutrophil counts were comparable between both genotypes, whereas CXC ligand1 levels were locally elevated in aortas of FAAH-deficient mice. We observed enhanced recruitment of neutrophils, but not monocytes, to large arteries of ApoE(-/-) mice treated with FAAH inhibitor URB597. Spleens of ApoE(-/-)FAAH(-/-) mice had reduced CD4+FoxP3+regulatory T-cell content, and in vitro stimulation of splenocytes revealed significantly elevated interferon-γ and tumor necrosis factor-α production in case of FAAH deficiency. CONCLUSIONS: Increased anandamide and related FAAH substrate levels are associated with the development of smaller atherosclerotic plaques with high neutrophil content, accompanied by an increased proinflammatory immune response.

Therapeutic monitoring of antidepressants in the era of pharmacogenetics studies.

As for other drugs, there is a large interindividual variability of the plasma concentrations of antidepressants for a given dose. Within the last 2 decades, a very large number of pharmacogenetic studies have made it possible to understand the importance of genetic factors on the disposition of drugs in the organism, many of them at the levels of drug metabolism. Polymorphism of CYP2D6 and of other drug-metabolizing enzymes may thus lead to very large differences in drug exposure between patients and possibly also to toxicity or ineffective drug concentrations in some subjects. In consequence, dose recommendations of antidepressants based on genotypes, justified by the principle of administering bioequivalent individualized drug doses, are now proposed. However, blood (and thus possibly brain) concentrations also depend on other factors than the genetic makeup of the patients. Therapeutic drug monitoring of antidepressants allows us to take into account the influence of factors such as comedications, diet, smoking habit, impaired organ function, and compliance. Therapeutic drug monitoring and genotyping are thus complementary, and their combined use contributes to improve pharmacotherapy with antidepressants and other drugs.

Validation of a Low-Cost Human Papillomavirus Genotyping Assay Based on PGMY PCR and Reverse Blotting Hybridization with Reusable Membranes.

The genotyping of human papillomaviruses (HPV) is essential for the surveillance of HPV vaccines. We describe and validate a low-cost PGMY-based PCR assay (PGMY-CHUV) for the genotyping of 31 HPV by reverse blotting hybridization (RBH). Genotype-specific detection limits were 50 to 500 genome equivalents per reaction. RBH was 100% specific and 98.61% sensitive using DNA sequencing as the gold standard (n = 1,024 samples). PGMY-CHUV was compared to the validated and commercially available linear array (Roche) on 200 samples. Both assays identified the same positive (n = 182) and negative samples (n = 18). Seventy-six percent of the positives were fully concordant after restricting the comparison to the 28 genotypes shared by both assays. At the genotypic level, agreement was 83% (285/344 genotype-sample combinations; κ of 0.987 for single infections and 0.853 for multiple infections). Fifty-seven of the 59 discordant cases were associated with multiple infections and with the weakest genotypes within each sample (P < 0.0001). PGMY-CHUV was significantly more sensitive for HPV56 (P = 0.0026) and could unambiguously identify HPV52 in mixed infections. PGMY-CHUV was reproducible on repeat testing (n = 275 samples; 392 genotype-sample combinations; κ of 0.933) involving different reagents lots and different technicians. Discordant results (n = 47) were significantly associated with the weakest genotypes in samples with multiple infections (P < 0.0001). Successful participation in proficiency testing also supported the robustness of this assay. The PGMY-CHUV reagent costs were estimated at $2.40 per sample using the least expensive yet proficient genotyping algorithm that also included quality control. This assay may be used in low-resource laboratories that have sufficient manpower and PCR expertise.

The ecological significance of arbuscular mycorrhizal fungal effects on clonal reproduction in plants

The population ecology of clonal plants depends on the number and distribution of ramets formed during growth. Variation in clonal reproduction has previously been explained by variation in effects of abiotic resource heterogeneity and by plant genotypic variation. Different co-occurring species of the mutualistic arbuscular mycorrhizal fungi (AMF) have been shown to differentially alter growth traits of Prunella vulgaris which we hypothesize would lead to changes in clonal reproduction. Two experiments were carried out to test whether different co-occurring mycorrhizal fungi significantly influence clonal reproduction of P. vulgaris whether this effect also occurs when P. vulgaris is growing in an artificial plant community and how the effects compare with plant genotype effects on clonal growth of P. vulgaris. In the first experiment the number of ramets of P. vulgaris growing in a plant community of simulated calcareous grassland was significantly affected by inoculation with different mycorrhizal fungi. The number of ramets produced by P. vulgaris differed by a factor of up to 1.8 with different mycorrhizal fungi. The fungal effects on the number of new ramets were independent of their effects on the biomass of P. vulgaris. In a second experiment 17 different genotypes of P. vulgaris were inoculated with different mycorrhizal fungi. There were significant main effects of genotypes and mycorrhizal fungi on clonal reproduction of P. vulgaris. The effect of different mycorrhizal fungi contributed more than the effect of plant genotype to variation in size and ramet production. However mean stolon length and spacer length which determine the spatial arrangement of ramets were only significantly affected by plant genotype. There were no mycorrhizal fungal X plant genotype interactions on clonal growth of P. vulgaris indicating that there is no obvious evidence that selection pressures would favor further coevolution between P. vulgaris and mycorrhizal fungal species. In natural communities plants can be colonized by several different AMF at the same time. The effect of the mixed AMF treatment on the growth and clonal reproduction of P. vulgaris could not be predicted from the responses of the plants to the single AMF To what extent however the patterns of colonization by different AMF differ among plants in a natural community is unknown. Since the effects of AMF on growth and clonal reproduction occur on a population of P. vulgaris in a microcosm plant community and because the effects are also as great as those caused by plant genotypic variation we conclude that the effects are strong enough to potentially affect population size and variation of clonal plants in communities.

Staphylococcus aureus host range and human-bovine host shift.

Staphylococcus aureus is a major agent of bovine mastitis. The concomitant emergence of pig-associated methicillin-resistant S. aureus (MRSA) in human carriage and infection requires a reexamination of the host range and specificity of human- and cow-associated S. aureus strains, something which has not been systematically studied previously. The genetic relatedness of 500 S. aureus isolates from bovine mastitis cases, 57 isolates from nasal carriage of farmers, and 133 isolates from nonfarmers was determined by amplified fragment length polymorphism (AFLP) analysis and spa typing. Multilocus sequence typing (MLST) was conducted on a subset of isolates to match AFLP clusters with MLST clonal complexes (CCs). This data set allowed us to study host range and host specificity and to estimate the extent of bovine-to-human transmission. The genotype compositions of S. aureus isolates from farmers and nonfarmers were very similar, while the mastitis isolates were quite distinct. Overall, transmission was low, but specific genotypes did show increased cow-to-human transmission. Unexpectedly, more than one-third of mastitis isolates belonged to CC8, a lineage which has not been considered to be bovine mastitis associated, but it is well known from human carriage and infection (i.e., USA300). Despite the fact that we did detect some transmission of other genotypes from cows to farmers, no transmission of CC8 isolates to farmers was detected, except for one tentative case. This was despite the close genetic relatedness of mastitis CC8 strains to nonfarmer carriage strains. These results suggest that the emergence of the new bovine-adapted genotype was due to a recent host shift from humans to cows concurrent with a loss of the ability to colonize humans. More broadly, our results indicate that host specificity is a lineage-specific trait that can rapidly evolve.

Identification and validation of copy number variants using SNP genotyping arrays from a large clinical cohort.

BACKGROUND: Genotypes obtained with commercial SNP arrays have been extensively used in many large case-control or population-based cohorts for SNP-based genome-wide association studies for a multitude of traits. Yet, these genotypes capture only a small fraction of the variance of the studied traits. Genomic structural variants (GSV) such as Copy Number Variation (CNV) may account for part of the missing heritability, but their comprehensive detection requires either next-generation arrays or sequencing. Sophisticated algorithms that infer CNVs by combining the intensities from SNP-probes for the two alleles can already be used to extract a partial view of such GSV from existing data sets. RESULTS: Here we present several advances to facilitate the latter approach. First, we introduce a novel CNV detection method based on a Gaussian Mixture Model. Second, we propose a new algorithm, PCA merge, for combining copy-number profiles from many individuals into consensus regions. We applied both our new methods as well as existing ones to data from 5612 individuals from the CoLaus study who were genotyped on Affymetrix 500K arrays. We developed a number of procedures in order to evaluate the performance of the different methods. This includes comparison with previously published CNVs as well as using a replication sample of 239 individuals, genotyped with Illumina 550K arrays. We also established a new evaluation procedure that employs the fact that related individuals are expected to share their CNVs more frequently than randomly selected individuals. The ability to detect both rare and common CNVs provides a valuable resource that will facilitate association studies exploring potential phenotypic associations with CNVs. CONCLUSION: Our new methodologies for CNV detection and their evaluation will help in extracting additional information from the large amount of SNP-genotyping data on various cohorts and use this to explore structural variants and their impact on complex traits.