Metric for Measuring the Effectiveness of Clustering of DNA Microarray Expression

BACKGROUND: The recent advancement of microarray technology with lower noise and better affordability makes it possible to determine expression of several thousand genes simultaneously. The differentially expressed genes are filtered first and then clustered based on the expression profiles of the genes. A large number of clustering algorithms and distance measuring matrices are proposed in the literature. The popular ones among them include hierarchal clustering and k-means clustering. These algorithms have often used the Euclidian distance or Pearson correlation distance. The biologists or the practitioners are often confused as to which algorithm to use since there is no clear winner among algorithms or among distance measuring metrics. Several validation indices have been proposed in the literature and these are based directly or indirectly on distances; hence a method that uses any of these indices does not relate to any biological features such as biological processes or molecular functions. RESULTS: In this paper we have proposed a metric to measure the effectiveness of clustering algorithms of genes by computing inter-cluster cohesiveness and as well as the intra-cluster separation with respect to biological features such as biological processes or molecular functions. We have applied this metric to the clusters on the data set that we have created as part of a larger study to determine the cancer suppressive mechanism of a class of chemicals called retinoids. We have considered hierarchal and k-means clustering with Euclidian and Pearson correlation distances. Our results show that genes of similar expression profiles are more likely to be closely related to biological processes than they are to molecular functions. The findings have been supported by many works in the area of gene clustering. CONCLUSION: The best clustering algorithm of genes must achieve cohesiveness within a cluster with respect to some biological features, and as well as maximum separation between clusters in terms of the distribution of genes of a behavioral group across clusters. We claim that our proposed metric is novel in this respect and that it provides a measure of both inter and intra cluster cohesiveness. Best of all, computation of the proposed metric is easy and it provides a single quantitative value, which makes comparison of different algorithms easier. The maximum cluster cohesiveness and the maximum intra-cluster separation are indicated by the metric when its value is 0. We have demonstrated the metric by applying it to a data set with gene behavioral groupings such as biological process and molecular functions. The metric can be easily extended to other features of a gene such as DNA binding sites and protein-protein interactions of the gene product, special features of the intron-exon structure, promoter characteristics, etc. The metric can also be used in other domains that use two different parametric spaces; one for clustering and the other one for measuring the effectiveness

Expression of MAGE-C1/CT7 and MAGE-C2/CT10 Predicts Lymph Node Metastasis in Melanoma Patients

MAGE-C1/CT7 and MAGE-C2/CT10 are members of the large MAGE family of cancer-testis (CT) antigens. CT antigens are promising targets for immunotherapy in cancer because their expression is restricted to cancer and germ line cells and a proportion of cancer patients presents with immune responses against CT antigens, which clearly demonstrates their immunogenicity. This study investigates the expression of MAGE-C1/CT7 and MAGE-C2/CT10 in primary and metastatic melanoma. Immunohistochemical staining of tissue microarrays that consisted of 59 primary malignant melanomas of the skin, 163 lymph node and distant melanoma metastases and 68 melanoma cell lines was performed. We found MAGE-C1/CT7 expression in 15 out of 50 (24%) primary melanomas and 15 out of 50 (24%) cell lines, whereas MAGE-C2/CT10 was detected in 17 out of 51 (33%) primary melanomas and 14 out of 68 (17%) cell lines. MAGE-C1/CT7 and MAGE-C2/CT10 were both detected in 40% of melanoma metastases. Patients with MAGE-C1/CT7 or MAGE-C2/CT10 positive primary melanoma had significantly more lymph node metastases (p = 0.005 and p<0.001, resp.). Prediction of lymph node metastasis by MAGE-C1/CT7 and MAGE-C2/CT10 was independent of tumor cell proliferation rate (Ki67 labeling index) in a multivariate analysis (p = 0.01). Our results suggest that the expression of MAGE-C1/CT7 and MAGE-C2/CT10 in primary melanoma is a potent predictor of sentinel lymph node metastasis

Laser excitation of the 1s-hyperfine transition in muonic hydrogen

The CREMA collaboration is pursuing a measurement of the ground-state hyperfine splitting (HFS) in muonic hydrogen ($\mu$p) with 1 ppm accuracy by means of pulsed laser spectroscopy to determine the two-photon-exchange contribution with $2\times10^{-4}$ relative accuracy. In the proposed experiment, the $\mu$p atom undergoes a laser excitation from the singlet hyperfine state to the triplet hyperfine state, {then} is quenched back to the singlet state by an inelastic collision with a H$_2$ molecule. The resulting increase of kinetic energy after the collisional deexcitation is used as a signature of a successful laser transition between hyperfine states. In this paper, we calculate the combined probability that a $\mu$p atom initially in the singlet hyperfine state undergoes a laser excitation to the triplet state followed by a collisional-induced deexcitation back to the singlet state. This combined probability has been computed using the optical Bloch equations including the inelastic and elastic collisions. Omitting the decoherence effects caused by {the laser bandwidth and }collisions would overestimate the transition probability by more than a factor of two in the experimental conditions. Moreover, we also account for Doppler effects and provide the matrix element, the saturation fluence, the elastic and inelastic collision rates for the singlet and triplet states, and the resonance linewidth. This calculation thus quantifies one of the key unknowns of the HFS experiment, leading to a precise definition of the requirements for the laser system and to an optimization of the hydrogen gas target where $\mu$p is formed and the laser spectroscopy will occur.Comment: 21 pages, 4 figure

The next generation of laser spectroscopy experiments using light muonic atoms

Precision spectroscopy of light muonic atoms provides unique information about the atomic and nuclear structure of these systems and thus represents a way to access fundamental interactions, properties and constants. One application comprises the determination of absolute nuclear charge radii with unprecedented accuracy from measurements of the 2S$\,$-$\,$2P Lamb shift. Here, we review recent results of nuclear charge radii extracted from muonic hydrogen and helium spectroscopy and present experiment proposals to access light muonic atoms with $Z \geq 3$. In addition, our approaches towards a precise measurement of the Zemach radii in muonic hydrogen ($\mu$p) and helium ($\mu$$^{3}$He$^{+}$) are discussed. These results will provide new tests of bound-state quantum-electrodynamics in hydrogen-like systems and can be used as benchmarks for nuclear structure theories.Comment: 17 pages, 8 figure

Rational Design of Protein Stability: Effect of (2S,4R)-4-Fluoroproline on the Stability and Folding Pathway of Ubiquitin

BACKGROUND: Many strategies have been employed to increase the conformational stability of proteins. The use of 4-substituted proline analogs capable to induce pre-organization in target proteins is an attractive tool to deliver an additional conformational stability without perturbing the overall protein structure. Both, peptides and proteins containing 4-fluorinated proline derivatives can be stabilized by forcing the pyrrolidine ring in its favored puckering conformation. The fluorinated pyrrolidine rings of proline can preferably stabilize either a C(γ)-exo or a C(γ)-endo ring pucker in dependence of proline chirality (4R/4S) in a complex protein structure. To examine whether this rational strategy can be generally used for protein stabilization, we have chosen human ubiquitin as a model protein which contains three proline residues displaying C(γ)-exo puckering. METHODOLOGY/PRINCIPAL FINDINGS: While (2S,4R)-4-fluoroproline ((4R)-FPro) containing ubiquitinin can be expressed in related auxotrophic Escherichia coli strain, all attempts to incorporate (2S,4S)-4-fluoroproline ((4S)-FPro) failed. Our results indicate that (4R)-FPro is favoring the C(γ)-exo conformation present in the wild type structure and stabilizes the protein structure due to a pre-organization effect. This was confirmed by thermal and guanidinium chloride-induced denaturation profile analyses, where we observed an increase in stability of -4.71 kJ·mol(-1) in the case of (4R)-FPro containing ubiquitin ((4R)-FPro-ub) compared to wild type ubiquitin (wt-ub). Expectedly, activity assays revealed that (4R)-FPro-ub retained the full biological activity compared to wt-ub. CONCLUSIONS/SIGNIFICANCE: The results fully confirm the general applicability of incorporating fluoroproline derivatives for improving protein stability. In general, a rational design strategy that enforces the natural occurring proline puckering conformation can be used to stabilize the desired target protein

Diffusion of muonic hydrogen in hydrogen gas and the measurement of the 1s hyperfine splitting of muonic hydrogen

725039; and the Swiss National Science Foundation through the projects SNF 200021_165854 and SNF 200020_197052. Publisher Copyright: © 2022 Cosmos Publishing Cooperative. All rights reserved.The CREMA collaboration is pursuing a measurement of the ground-state hyperfine splitting (HFS) in muonic hydrogen (µp) with 1 ppm accuracy by means of pulsed laser spectroscopy. In the proposed experiment, the µp atom is excited by a laser pulse from the singlet to the triplet hyperfine sub-levels, and is quenched back to the singlet state by an inelastic collision with a H2 molecule. The resulting increase of kinetic energy after this cycle modifies the µp atom diffusion in the hydrogen gas and the arrival time of the µp atoms at the target walls. This laser-induced modification of the arrival times is used to expose the atomic transition. In this paper we present the simulation of the µp diffusion in the H2 gas which is at the core of the experimental scheme. These simulations have been implemented with the Geant4 framework by introducing various low-energy processes including the motion of the H2 molecules, i.e. the effects related with the hydrogen target temperature. The simulations have been used to optimize the hydrogen target parameters (pressure, temperatures and thickness) and to estimate signal and background rates. These rates allow to estimate the maximum time needed to find the resonance and the statistical accuracy of the spectroscopy experiment.publishersversionpublishe

Gene Expression Profiles in Parkinson Disease Prefrontal Cortex Implicate FOXO1 and Genes under Its Transcriptional Regulation

Parkinson disease (PD) is a complex neurodegenerative disorder with largely unknown genetic mechanisms. While the degeneration of dopaminergic neurons in PD mainly takes place in the substantia nigra pars compacta (SN) region, other brain areas, including the prefrontal cortex, develop Lewy bodies, the neuropathological hallmark of PD. We generated and analyzed expression data from the prefrontal cortex Brodmann Area 9 (BA9) of 27 PD and 26 control samples using the 44K One-Color Agilent 60-mer Whole Human Genome Microarray. All samples were male, without significant Alzheimer disease pathology and with extensive pathological annotation available. 507 of the 39,122 analyzed expression probes were different between PD and control samples at false discovery rate (FDR) of 5%. One of the genes with significantly increased expression in PD was the forkhead box O1 (FOXO1) transcription factor. Notably, genes carrying the FoxO1 binding site were significantly enriched in the FDR–significant group of genes (177 genes covered by 189 probes), suggesting a role for FoxO1 upstream of the observed expression changes. Single-nucleotide polymorphisms (SNPs) selected from a recent meta-analysis of PD genome-wide association studies (GWAS) were successfully genotyped in 50 out of the 53 microarray brains, allowing a targeted expression–SNP (eSNP) analysis for 52 SNPs associated with PD affection at genome-wide significance and the 189 probes from FoxO1 regulated genes. A significant association was observed between a SNP in the cyclin G associated kinase (GAK) gene and a probe in the spermine oxidase (SMOX) gene. Further examination of the FOXO1 region in a meta-analysis of six available GWAS showed two SNPs significantly associated with age at onset of PD. These results implicate FOXO1 as a PD–relevant gene and warrant further functional analyses of its transcriptional regulatory mechanisms

Co-Localization of the Oncogenic Transcription Factor MYCN and the DNA Methyl Binding Protein MeCP2 at Genomic Sites in Neuroblastoma

MYCN is a transcription factor that is expressed during the development of the neural crest and its dysregulation plays a major role in the pathogenesis of pediatric cancers such as neuroblastoma, medulloblastoma and rhabdomyosarcoma. MeCP2 is a CpG methyl binding protein which has been associated with a number of cancers and developmental disorders, particularly Rett syndrome.Using an integrative global genomics approach involving chromatin immunoprecipitation applied to microarrays, we have determined that MYCN and MeCP2 co-localize to gene promoter regions, as well as inter/intragenic sites, within the neuroblastoma genome (MYCN amplified Kelly cells) at high frequency (70.2% of MYCN sites were also positive for MeCP2). Intriguingly, the frequency of co-localization was significantly less at promoter regions exhibiting substantial hypermethylation (8.7%), as determined by methylated DNA immunoprecipitation (MeDIP) applied to the same microarrays. Co-immunoprecipitation of MYCN using an anti-MeCP2 antibody indicated that a MYCN/MeCP2 interaction occurs at protein level. mRNA expression profiling revealed that the median expression of genes with promoters bound by MYCN was significantly higher than for genes bound by MeCP2, and that genes bound by both proteins had intermediate expression. Pathway analysis was carried out for genes bound by MYCN, MeCP2 or MYCN/MeCP2, revealing higher order functions.Our results indicate that MYCN and MeCP2 protein interact and co-localize to similar genomic sites at very high frequency, and that the patterns of binding of these proteins can be associated with significant differences in transcriptional activity. Although it is not yet known if this interaction contributes to neuroblastoma disease pathogenesis, it is intriguing that the interaction occurs at the promoter regions of several genes important for the development of neuroblastoma, including ALK, AURKA and BDNF