MGMT methylation analysis of glioblastoma on the Infinium methylation BeadChip identifies two distinct CpG regions associated with gene silencing and outcome, yielding a prediction model for comparisons across datasets, tumor grades, and CIMP-status.

The methylation status of the O(6)-methylguanine-DNA methyltransferase (MGMT) gene is an important predictive biomarker for benefit from alkylating agent therapy in glioblastoma. Recent studies in anaplastic glioma suggest a prognostic value for MGMT methylation. Investigation of pathogenetic and epigenetic features of this intriguingly distinct behavior requires accurate MGMT classification to assess high throughput molecular databases. Promoter methylation-mediated gene silencing is strongly dependent on the location of the methylated CpGs, complicating classification. Using the HumanMethylation450 (HM-450K) BeadChip interrogating 176 CpGs annotated for the MGMT gene, with 14 located in the promoter, two distinct regions in the CpG island of the promoter were identified with high importance for gene silencing and outcome prediction. A logistic regression model (MGMT-STP27) comprising probes cg1243587 and cg12981137 provided good classification properties and prognostic value (kappa = 0.85; log-rank p < 0.001) using a training-set of 63 glioblastomas from homogenously treated patients, for whom MGMT methylation was previously shown to be predictive for outcome based on classification by methylation-specific PCR. MGMT-STP27 was successfully validated in an independent cohort of chemo-radiotherapy-treated glioblastoma patients (n = 50; kappa = 0.88; outcome, log-rank p < 0.001). Lower prevalence of MGMT methylation among CpG island methylator phenotype (CIMP) positive tumors was found in glioblastomas from The Cancer Genome Atlas than in low grade and anaplastic glioma cohorts, while in CIMP-negative gliomas MGMT was classified as methylated in approximately 50 % regardless of tumor grade. The proposed MGMT-STP27 prediction model allows mining of datasets derived on the HM-450K or HM-27K BeadChip to explore effects of distinct epigenetic context of MGMT methylation suspected to modulate treatment resistance in different tumor types.

Clear detection of ADIPOQ locus as the major gene for plasma adiponectin: results of genome-wide association analyses including 4659 European individuals.

OBJECTIVE: Plasma adiponectin is strongly associated with various components of metabolic syndrome, type 2 diabetes and cardiovascular outcomes. Concentrations are highly heritable and differ between men and women. We therefore aimed to investigate the genetics of plasma adiponectin in men and women. METHODS: We combined genome-wide association scans of three population-based studies including 4659 persons. For the replication stage in 13795 subjects, we selected the 20 top signals of the combined analysis, as well as the 10 top signals with p-values less than 1.0 x 10(-4) for each the men- and the women-specific analyses. We further selected 73 SNPs that were consistently associated with metabolic syndrome parameters in previous genome-wide association studies to check for their association with plasma adiponectin. RESULTS: The ADIPOQ locus showed genome-wide significant p-values in the combined (p=4.3 x 10(-24)) as well as in both women- and men-specific analyses (p=8.7 x 10(-17) and p=2.5 x 10(-11), respectively). None of the other 39 top signal SNPs showed evidence for association in the replication analysis. None of 73 SNPs from metabolic syndrome loci exhibited association with plasma adiponectin (p>0.01). CONCLUSIONS: We demonstrated the ADIPOQ gene as the only major gene for plasma adiponectin, which explains 6.7% of the phenotypic variance. We further found that neither this gene nor any of the metabolic syndrome loci explained the sex differences observed for plasma adiponectin. Larger studies are needed to identify more moderate genetic determinants of plasma adiponectin.

Stochastic models and numerical algorithms for a class of regulatory gene networks.

Regulatory gene networks contain generic modules, like those involving feedback loops, which are essential for the regulation of many biological functions (Guido et al. in Nature 439:856-860, 2006). We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady-state distribution of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops.

Clinical trial substantiates the predictive value of O-6-methylguanine-DNA methyltransferase promoter methylation in glioblastoma patients treated with temozolomide.

PURPOSE: In the setting of a prospective clinical trial, we determined the predictive value of the methylation status of the O-6-methylguanine-DNA methyltransferase (MGMT) promoter for outcome in glioblastoma patients treated with the alkylating agent temozolomide. Expression of this excision repair enzyme has been associated with resistance to alkylating chemotherapy. EXPERIMENTAL DESIGN: The methylation status of MGMT in the tumor biopsies was evaluated in 38 patients undergoing resection for newly diagnosed glioblastoma and enrolled in a Phase II trial testing concomitant and adjuvant temozolomide and radiation. The epigenetic silencing of the MGMT gene was determined using methylation-specific PCR. RESULTS: Inactivation of the MGMT gene by promoter methylation was associated with longer survival (P = 0.0051; Log-rank test). At 18 months, survival was 62% (16 of 26) for patients testing positive for a methylated MGMT promoter but reached only 8% (1 of 12) in absence of methylation (P = 0.002; Fisher's exact test). In the presence of other clinically relevant factors, methylation of the MGMT promoter remains the only significant predictor (P = 0.017; Cox regression). CONCLUSIONS: This prospective clinical trial identifies MGMT-methylation status as an independent predictor for glioblastoma patients treated with a methylating agent. The association of the epigenetic inactivation of the DNA repair gene MGMT with better outcome in this homogenous cohort may have important implications for the design of future trials and supports efforts to deplete MGMT by O-6-benzylguanine, a noncytotoxic substrate of this enzyme.

In vivo gene transfer into the ocular ciliary muscle mediated by ultrasound and microbubbles.

This study aimed to assess application of ultrasound (US) combined with microbubbles (MB) to transfect the ciliary muscle of rat eyes. Reporter DNA plasmids encoding for Gaussia luciferase, ^6;-galactosidase or the green fluorescent protein (GFP), alone or mixed with 50% Artison MB, were injected into the ciliary muscle, with or without US exposure (US set at 1 MHz, 2 W/cm(2), 50% duty cycle for 2 min). Luciferase activity was measured in ocular fluids at 7 and 30 days after sonoporation. At 1 week, the US+MB treatment showed a significant increase in luminescence compared with control eyes, injected with plasmid only, with or without MB (×2.6), and, reporter proteins were localized in the ciliary muscle by histochemical analysis. At 1 month, a significant decrease in luciferase activity was observed in all groups. A rise in lens and ciliary muscle temperature was measured during the procedure but did not result in any observable or microscopic damages at 1 and 8 days. The feasibility to transfer gene into the ciliary muscle by US and MB suggests that sonoporation may allow intraocular production of proteins for the treatment of inflammatory, angiogenic and/or degenerative retinal diseases.

Members of the PHO1 gene family show limited functional redundancy in phosphate transfer to the shoot, and are regulated by phosphate deficiency via distinct pathways.

The PHO1 family comprises 11 members in Arabidopsis thaliana. In order to decipher the role of these genes in inorganic phosphate (Pi) transport and homeostasis, complementation of the pho1 mutant, deficient in loading Pi to the root xylem, was determined by the expression of the PHO1 homologous genes under the control of the PHO1 promoter. Only PHO1 and the homologue PHO1;H1 could complement pho1. The PHO1;H1 promoter was active in the vascular cylinder of roots and shoots. Expression of PHO1;H1 was very low in Pi-sufficient plants, but was strongly induced under Pi-deficient conditions. T-DNA knock-out mutants of PHO1;H1 neither showed growth defects nor alteration in Pi transport dynamics, or Pi content, compared with wild type. However, the double mutant pho1/pho1;h1 showed a strong reduction in growth and in the capacity to transfer Pi from the root to the shoot compared with pho1. Grafting experiments revealed that phenotypes associated with the pho1 and pho1/pho1;h1 mutants were linked to the lack of gene expression in the root. The increased expression of PHO1;H1 under Pi deficiency was largely controlled by the transcription factor PHR1 and was suppressed by the phosphate analogue phosphite, whereas the increase of PHO1 expression was independent of PHR1 and was not influenced by phosphite. Together, these data reveal that although transfer of Pi to the root xylem vessel is primarily mediated by PHO1, the homologue PHO1;H1 also contributes to Pi loading to the xylem, and that the two corresponding genes are regulated by Pi deficiency by distinct signal transduction pathways.

Comparative genomics suggests that the human pathogenic fungus Pneumocystis jirovecii acquired obligate biotrophy through gene loss.

Pneumocystis jirovecii is a fungal parasite that colonizes specifically humans and turns into an opportunistic pathogen in immunodeficient individuals. The fungus is able to reproduce extracellularly in host lungs without eliciting massive cellular death. The molecular mechanisms that govern this process are poorly understood, in part because of the lack of an in vitro culture system for Pneumocystis spp. In this study, we explored the origin and evolution of the putative biotrophy of P. jirovecii through comparative genomics and reconstruction of ancestral gene repertoires. We used the maximum parsimony method and genomes of related fungi of the Taphrinomycotina subphylum. Our results suggest that the last common ancestor of Pneumocystis spp. lost 2,324 genes in relation to the acquisition of obligate biotrophy. These losses may result from neutral drift and affect the biosyntheses of amino acids and thiamine, the assimilation of inorganic nitrogen and sulfur, and the catabolism of purines. In addition, P. jirovecii shows a reduced panel of lytic proteases and has lost the RNA interference machinery, which might contribute to its genome plasticity. Together with other characteristics, that is, a sex life cycle within the host, the absence of massive destruction of host cells, difficult culturing, and the lack of virulence factors, these gene losses constitute a unique combination of characteristics which are hallmarks of both obligate biotrophs and animal parasites. These findings suggest that Pneumocystis spp. should be considered as the first described obligate biotrophs of animals, whose evolution has been marked by gene losses.

Contrasting diffusion of Quaternary gene pools across Europe: the case of the arctic-alpine Gentiana nivalis L. (Gentianaceae).

The fate of European arctic-alpine species during Pleistocene climatic oscillations still remains debated. Did these cold-adapted species invade much of the continental steppe or did they remain restricted to warmer slopes of inner mountain massifs? To examine this question, we investigated the phylogeography of Gentiana nivalis, a typical European arctic-alpine plant species. Genome fingerprinting analyses revealed that four genetic pools are actually unevenly distributed across the continent. One cluster covers almost all mountain massifs as well as northern areas, and thus coincides with a scenario of past distribution covering a large part of the European glacial steppe. In contrast, the three other lineages are strongly restricted spatially to western, central, and eastern Alps, respectively, thus arguing towards a scenario of in situ glacial survival. The coexistence of lineages with such contrasting demographic histories in Europe challenges our classical view of refugia and corroborates several hypotheses of biogeographers from the twentieth century.

Differential regulation of Na-K-ATPase isoform gene expression by T3 during rat brain development.

A fetal rat telencephalon organotypic cell culture system was found to reproduce the developmental pattern of Na-K-adenosinetriphosphatase (ATPase) gene expression observed in vivo [Am. J. Physiol. 258 (Cell Physiol. 27): C1062-C1069, 1990]. We have used this culture system to study the effects of triiodothyronine (T3; 0.003-30 nM) on mRNA abundance and basal transcription rates of Na-K-ATPase isoforms. Steady-state mRNA levels were low at culture day 6 (corresponding to the day of birth) but distinct for each isoform alpha 3 much greater than beta 1 = beta 2 greater than alpha 2 greater than alpha 1. At culture day 6, T3 did not modify mRNA abundance of any isoform. At culture day 12 (corresponding to day 7 postnatal), T3 increased the mRNA level of alpha 2 (4- to 7-fold), beta 2 (4- to 5-fold), alpha 1 (3- to 6-fold), and beta 1 (1.5-fold), whereas alpha 3 mRNA levels remained unchanged. Interestingly, the basal transcription rate for each isoform differed strikingly (alpha 2 greater than alpha 1 much greater than beta 1 = beta 2 greater than alpha 3) but remained stable throughout 12 days of culture and was not regulated by T3. Thus we observed an inverse relationship between rate of transcription and rate of mRNA accumulation for each alpha-isoform, suggesting that alpha 1- and alpha 2-mRNA are turning over rapidly whereas alpha 3-mRNA is turning over slowly. Our data indicate that one of the mechanisms by which T3 selectively controls Na-K-ATPase gene expression during brain development in vitro occurs at the posttranscriptional level.

Feasibility of RNA studies on illegitimate transcription for molecular characterization of splicing mutations in the ATP7B gene: a case report.

Approximately 520 Wilson disease-causing mutations in the ATP7B gene have been described to date. In this study we report DNA and RNA analyses carried out for molecular characterization of a consensus sequence splicing mutation found in homozygosity in a Swiss Wilson disease patient. RNA analysis of 1946 +6 T→C in both the peripheral lymphoblasts and liver resulted in the production in the propositus of only an alternative transcript lacking exons 6, 7, and 8 resulting most likely in alterations of cell biochemistry and disease. The patient presents an early form of severe hepatic disease characterized by hepatosplenomegaly, reduced hepatic function, anemia and thrombocytopenia indicating that 1946 +6 T→C is a severe mutation. Since identical results were obtained from both peripheral lymphoblasts and liver they also suggest that RNA studies of illegitimate transcripts can be safely used for molecular characterization of ATP7B splicing mutations, thus improving genetic counseling and diagnosis of Wilson disease. Moreover these studies, contribute to reveal the exact molecular mechanisms producing Wilson disease.

In vitro emergence of rifampicin resistance in Propionibacterium acnes and molecular characterization of mutations in the rpoB gene.

OBJECTIVES: Activity of rifampicin against Propionibacterium acnes biofilms was recently demonstrated, but rifampicin resistance has not yet been described in this organism. We investigated the in vitro emergence of rifampicin resistance in P. acnes and characterized its molecular background. METHODS: P. acnes ATCC 11827 was used (MIC 0.007 mg/L). The mutation rate was determined by inoculation of 10(9) cfu of P. acnes on rifampicin-containing agar plates incubated anaerobically for 7 days. Progressive emergence of resistance was studied by serial exposure to increasing concentrations of rifampicin in 72 h cycles using a low (10(6) cfu/mL) and high (10(8) cfu/mL) inoculum. The stability of resistance was determined after three subcultures of rifampicin-resistant isolates on rifampicin-free agar. For resistant mutants, the whole rpoB gene was amplified, sequenced and compared with a P. acnes reference sequence (NC006085). RESULTS: P. acnes growth was observed on rifampicin-containing plates with mutation rates of 2 ± 1 cfu  ×  10(-9) (4096× MIC) and 12 ± 5 cfu  ×  10(-9) (4 × MIC). High-level rifampicin resistance emerged progressively after 4 (high inoculum) and 13 (low inoculum) cycles. In rifampicin-resistant isolates, the MIC remained >32 mg/L after three subcultures. Mutations were detected in clusters I (amino acids 418-444) and II (amino acids 471-486) of the rpoB gene after sequence alignment with a Staphylococcus aureus reference sequence (CAA45512). The five following substitutions were found: His-437 → Tyr, Ser-442 → Leu, Leu-444 → Ser, Ile-483 → Val and Ser-485 → Leu. CONCLUSION: The rifampicin MIC increased from highly susceptible to highly resistant values. The resistance remained stable and was associated with mutations in the rpoB gene. To our knowledge, this is the first report of the emergence of rifampicin resistance in P. acnes.

Evolution of the P-type II ATPase gene family in the fungi and presence of structural genomic changes among isolates of Glomus intraradices.

BACKGROUND: The P-type II ATPase gene family encodes proteins with an important role in adaptation of the cell to variation in external K+, Ca2+ and Na2+ concentrations. The presence of P-type II gene subfamilies that are specific for certain kingdoms has been reported but was sometimes contradicted by discovery of previously unknown homologous sequences in newly sequenced genomes. Members of this gene family have been sampled in all of the fungal phyla except the arbuscular mycorrhizal fungi (AMF; phylum Glomeromycota), which are known to play a key-role in terrestrial ecosystems and to be genetically highly variable within populations. Here we used highly degenerate primers on AMF genomic DNA to increase the sampling of fungal P-Type II ATPases and to test previous predictions about their evolution. In parallel, homologous sequences of the P-type II ATPases have been used to determine the nature and amount of polymorphism that is present at these loci among isolates of Glomus intraradices harvested from the same field. RESULTS: In this study, four P-type II ATPase sub-families have been isolated from three AMF species. We show that, contrary to previous predictions, P-type IIC ATPases are present in all basal fungal taxa. Additionally, P-Type IIE ATPases should no longer be considered as exclusive to the Ascomycota and the Basidiomycota, since we also demonstrate their presence in the Zygomycota. Finally, a comparison of homologous sequences encoding P-type IID ATPases showed unexpectedly that indel mutations among coding regions, as well as specific gene duplications occur among AMF individuals within the same field. CONCLUSION: On the basis of these results we suggest that the diversification of P-Type IIC and E ATPases followed the diversification of the extant fungal phyla with independent events of gene gains and losses. Consistent with recent findings on the human genome, but at a much smaller geographic scale, we provided evidence that structural genomic changes, such as exonic indel mutations and gene duplications are less rare than previously thought and that these also occur within fungal populations.

Gene by environment interaction: effects of a single-gene and social-environment on reproductive phenotypes of fire ant queens

1. The gene Pgm-3 (or a closely linked gene) influences the phenotype and reproductive success of queens in multiple-queen (polygynous) colonies but not single-queen (monogynous) colonies of the Fire Ant Solenopsis invicta. 2. We investigated the mechanisms of differential phenotypic expression of Pgm-3 in these alternate social forms. Mature winged queens with the homozygous genotype Pgm-3(a/a) averaged 26% heavier than queens with the genotypes Pgm-3(a/b) and Pgm 3(b/b) in the polygynous form. Heterozygotes were slightly heavier (2%) than Pgm-3(b/b) queens in this form, demonstrating that the allele Pgm-3(a) is not completely recessive in its effects on weight. 3. There was no significant difference in weight among queens of the three Pgm-3 genotypes in the monogynous form, with the mean weight of monogynous queens slightly greater than that of polygynous Pgm-3(a/a) queens. Differences in weight between queens of the two social forms and among queens of the three genotypes in the polygynous form are not evident at the pupal stage and thus appear to develop during sexual maturation of the adults. This suggests that some component of the social environment of polygynous colonies inhibits weight gains during queen maturation and that Pgm-(3a/a) queens are relatively less sensitive to this factor. 4. To test whether the high cumulative queen pheromone level characteristic of polygynous colonies is the factor responsible for the differential queen maturation, we compared phenotypes of winged queens reared in split colonies in which pheromone levels were manipulated by adjusting queen number. Queens produced in colony fragments made monogynous were heavier than those produced in polygynous fragments, a finding consistent with the hypothesis that pheromone level affects the reproductive development of queens. However, genotype-specific differences in weights of queens were similar between the two treatments, suggesting that pheromone level was not the key factor of the social environment responsible for the gene-environment interaction. 5. To test whether limited food availability to winged queens associated with the high brood/worker ratios in polygynous colonies is the factor responsible for this interaction, similar split-colony experiments were performed. Elevated brood/worker ratios decreased the weight of winged queens but there was no evidence that this treatment intensified differential weight gains among queens with different Pgm-3 genotypes. Manipulation of the amount of food provided to colonies had no effect on queen weight. 6. The combined data indicate that cumulative pheromone level and brood/worker ratio are two of the factors responsible for the differences in reproductive phenotypes between monogynous and polygynous winged queens but that these factors are not directly responsible for inducing the phenotypic effects of Pgm-3 in polygynous colonies.

The vitamin D receptor gene bAt (CCA) haplotype impairs the response to pegylated-interferon/ribavirin-based therapy in chronic hepatitis C patients.

BACKGROUND: Chronic hepatitis C infection is a major cause of end-stage liver disease. Therapy outcome is influenced by 25-OH vitamin D deficiency. To further address this observation, our study investigates the impact of the vitamin D receptor (NR1I1) haplotype and combined effects of plasma vitamin D levels in a well-described cohort of hepatitis C patients. METHODS: A total of 155 chronic hepatitis C patients were recruited from the Swiss Hepatitis C Cohort Study for NR1I1 genotyping and plasma 25-OH vitamin D level measurement. NR1I1 genotype data and combined effects of plasma 25-OH vitamin D level were analysed regarding therapy response (sustained virological response). RESULTS: A strong association was observed between therapy non-response and the NR1I1 CCA (bAt) haplotype consisting of rs1544410 (BsmI) C, rs7975232 (ApaI) C and rs731236 (TaqI) A alleles. Of the HCV patients carrying the CCA haplotype, 50.3% were non-responders (odds ratio [OR] 1.69, 95% CI 1.07, 2.67; P=0.028). A similar association was observed for the combinational CCCCAA genotype (OR 2.94, 95% CI 1.36, 6.37; P=0.007). The combinational CCCCAA genotype was confirmed as an independent risk factor for non-response in multivariate analysis (OR 2.50, 95% CI 1.07, 5.87; P=0.034). Analysing combined effects, a significant impact of low 25-OH vitamin D levels on sustained virological response were only seen in patients with the unfavourable NR1I1 CCA (bAt) haplotype (OR for non-SVR 3.55; 95% CI 1.005, 12.57; P=0.049). CONCLUSIONS: NR1I1 vitamin D receptor polymorphisms influence response to pegylated-interferon/ribavirin-based therapy in chronic hepatitis C and exert an additive genetic predisposition to previously described low 25-OH vitamin D serum levels.

Validation of real-time methylation-specific PCR to determine O6-methylguanine-DNA methyltransferase gene promoter methylation in glioma.

Epigenetic silencing of the DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) by promoter methylation predicts successful alkylating agent therapy, such as with temozolomide, in glioblastoma patients. Stratified therapy assignment of patients in prospective clinical trials according to tumor MGMT status requires a standardized diagnostic test, suitable for high-throughput analysis of small amounts of formalin-fixed, paraffin-embedded tumor tissue. A direct, real-time methylation-specific PCR (MSP) assay was developed to determine methylation status of the MGMT gene promoter. Assay specificity was obtained by selective amplification of methylated DNA sequences of sodium bisulfite-modified DNA. The copy number of the methylated MGMT promoter, normalized to the beta-actin gene, provides a quantitative test result. We analyzed 134 clinical glioma samples, comparing the new test with the previously validated nested gel-based MSP assay, which yields a binary readout. A cut-off value for the MGMT methylation status was suggested by fitting a bimodal normal mixture model to the real-time results, supporting the hypothesis that there are two distinct populations within the test samples. Comparison of the tests showed high concordance of the results (82/91 [90%]; Cohen's kappa = 0.80; 95% confidence interval, 0.82-0.95). The direct, real-time MSP assay was highly reproducible (Pearson correlation 0.996) and showed valid test results for 93% (125/134) of samples compared with 75% (94/125) for the nested, gel-based MSP assay. This high-throughput test provides an important pharmacogenomic tool for individualized management of alkylating agent chemotherapy.