ICGE: an R package for detecting relevant clusters and atypical units in gene expression

Background: Gene expression technologies have opened up new ways to diagnose and treat cancer and other diseases. Clustering algorithms are a useful approach with which to analyze genome expression data. They attempt to partition the genes into groups exhibiting similar patterns of variation in expression level. An important problem associated with gene classification is to discern whether the clustering process can find a relevant partition as well as the identification of new genes classes. There are two key aspects to classification: the estimation of the number of clusters, and the decision as to whether a new unit (gene, tumor sample ... ) belongs to one of these previously identified clusters or to a new group. Results: ICGE is a user-friendly R package which provides many functions related to this problem: identify the number of clusters using mixed variables, usually found by applied biomedical researchers; detect whether the data have a cluster structure; identify whether a new unit belongs to one of the pre-identified clusters or to a novel group, and classify new units into the corresponding cluster. The functions in the ICGE package are accompanied by help files and easy examples to facilitate its use. Conclusions: We demonstrate the utility of ICGE by analyzing simulated and real data sets. The results show that ICGE could be very useful to a broad research community.

Identification of a biomarker panel for colorectal cancer diagnosis

Background: Malignancies arising in the large bowel cause the second largest number of deaths from cancer in the Western World. Despite progresses made during the last decades, colorectal cancer remains one of the most frequent and deadly neoplasias in the western countries. Methods: A genomic study of human colorectal cancer has been carried out on a total of 31 tumoral samples, corresponding to different stages of the disease, and 33 non-tumoral samples. The study was carried out by hybridisation of the tumour samples against a reference pool of non-tumoral samples using Agilent Human 1A 60- mer oligo microarrays. The results obtained were validated by qRT-PCR. In the subsequent bioinformatics analysis, gene networks by means of Bayesian classifiers, variable selection and bootstrap resampling were built. The consensus among all the induced models produced a hierarchy of dependences and, thus, of variables. Results: After an exhaustive process of pre-processing to ensure data quality–lost values imputation, probes quality, data smoothing and intraclass variability filtering–the final dataset comprised a total of 8, 104 probes. Next, a supervised classification approach and data analysis was carried out to obtain the most relevant genes. Two of them are directly involved in cancer progression and in particular in colorectal cancer. Finally, a supervised classifier was induced to classify new unseen samples. Conclusions: We have developed a tentative model for the diagnosis of colorectal cancer based on a biomarker panel. Our results indicate that the gene profile described herein can discriminate between non-cancerous and cancerous samples with 94.45% accuracy using different supervised classifiers (AUC values in the range of 0.997 and 0.955).

Second-to-Fourth Digit Ratio Has a Non-Monotonic Impact on Altruism

Gene-culture co-evolution emphasizes the joint role of culture and genes for the emergence of altruistic and cooperative behaviors and behavioral genetics provides estimates of their relative importance. However, these approaches cannot assess which biological traits determine altruism or how. We analyze the association between altruism in adults and the exposure to prenatal sex hormones, using the second-to-fourth digit ratio. We find an inverted U-shaped relation for left and right hands, which is very consistent for men and less systematic for women. Subjects with both high and low digit ratios give less than individuals with intermediate digit ratios. We repeat the exercise with the same subjects seven months later and find a similar association, even though subjects' behavior differs the second time they play the game. We then construct proxies of the median digit ratio in the population (using more than 1000 different subjects), show that subjects' altruism decreases with the distance of their ratio to these proxies. These results provide direct evidence that prenatal events contribute to the variation of altruistic behavior and that the exposure to fetal hormones is one of the relevant biological factors. In addition, the findings suggest that there might be an optimal level of exposure to these hormones from social perspective.

Fine mapping and functional analysis of the multiple sclerosis risk gene CD6

CD6 has recently been identified and validated as risk gene for multiple sclerosis (MS), based on the association of a single nucleotide polymorphism (SNP), rs17824933, located in intron 1. CD6 is a cell surface scavenger receptor involved in T-cell activation and proliferation, as well as in thymocyte differentiation. In this study, we performed a haptag SNP screen of the CD6 gene locus using a total of thirteen tagging SNPs, of which three were non-synonymous SNPs, and replicated the recently reported GWAS SNP rs650258 in a Spanish-Basque collection of 814 controls and 823 cases. Validation of the six most strongly associated SNPs was performed in an independent collection of 2265 MS patients and 2600 healthy controls. We identified association of haplotypes composed of two non-synonymous SNPs [rs11230563 (R225W) and rs2074225 (A257V)] in the 2nd SRCR domain with susceptibility to MS (Pmax(T) permutation=161024). The effect of these haplotypes on CD6 surface expression and cytokine secretion was also tested. The analysis showed significantly different CD6 expression patterns in the distinct cell subsets, i.e. – CD4+ naı¨ve cells, P = 0.0001; CD8+ naı¨ve cells, P,0.0001; CD4+ and CD8+ central memory cells, P = 0.01 and 0.05, respectively; and natural killer T (NKT) cells, P = 0.02; with the protective haplotype (RA) showing higher expression of CD6. However, no significant changes were observed in natural killer (NK) cells, effector memory and terminally differentiated effector memory T cells. Our findings reveal that this new MS-associated CD6 risk haplotype significantly modifies expression of CD6 on CD4+ and CD8+ T cells.

Preprocess and data analysis techniques for affymetrix DNA microarrays using bioconductor: a case study in Alzheimer disease

DNA microarray, or DNA chip, is a technology that allows us to obtain the expression level of many genes in a single experiment. The fact that numerical expression values can be easily obtained gives us the possibility to use multiple statistical techniques of data analysis. In this project microarray data is obtained from Gene Expression Omnibus, the repository of National Center for Biotechnology Information (NCBI). Then, the noise is removed and data is normalized, also we use hypothesis tests to find the most relevant genes that may be involved in a disease and use machine learning methods like KNN, Random Forest or Kmeans. For performing the analysis we use Bioconductor, packages in R for the analysis of biological data, and we conduct a case study in Alzheimer disease. The complete code can be found in https://github.com/alberto-poncelas/ bioc-alzheimer

Design,optimization and in vivo evaluation of non-viral vectors for gene therapy. Applications in ocular disorders

227 págs.

Behavioral, neurochemical and morphological changes induced by the overexpression of munc18-1a in brain of mice: relevance to schizophrenia

Overexpression of the mammalian homolog of the unc-18 gene (munc18-1) has been described in the brain of subjects with schizophrenia. Munc18-1 protein is involved in membrane fusion processes, exocytosis and neurotransmitter release. A transgenic mouse strain that overexpresses the protein isoform munc18-1a in the brain was characterized. This animal displays several schizophrenia-related behaviors, supersensitivity to hallucinogenic drugs and deficits in prepulse inhibition that reverse after antipsychotic treatment. Relevant brain areas (that is, cortex and striatum) exhibit reduced expression of dopamine D-1 receptors and dopamine transporters together with enhanced amphetamine-induced in vivo dopamine release. Magnetic resonance imaging demonstrates decreased gray matter volume in the transgenic animal. In conclusion, the mouse overexpressing brain munc18-1a represents a new valid animal model that resembles functional and structural abnormalities in patients with schizophrenia.

Role of G protein coupled inwardly rectifier K+ channels (GIRK) on the neurobiology depression

353 págs.

Control of gene expression by modulated self-assembly

Numerous transcription factors self-assemble into different order oligomeric species in a way that is actively regulated by the cell. Until now, no general functional role has been identified for this widespread process. Here, we capture the effects of modulated self-assembly in gene expression with a novel quantitative framework. We show that this mechanism provides precision and flexibility, two seemingly antagonistic properties, to the sensing of diverse cellular signals by systems that share common elements present in transcription factors like p53, NF-kappa B, STATs, Oct and RXR. Applied to the nuclear hormone receptor RXR, this framework accurately reproduces a broad range of classical, previously unexplained, sets of gene expression data and corroborates the existence of a precise functional regime with flexible properties that can be controlled both at a genome-wide scale and at the individual promoter level.

Prognostic Role of Host Cyclooxygenase and Cytokine Genotypes in a Caucasian Cohort of Patients with Gastric Adenocarcinoma

11 p.

The promoter of cell growth- and RNA protection-associated SND1 gene is activated by endoplasmic reticulum stress in human hepatoma cells

Background: Staphyloccocal nuclease domain-containing protein 1 (SND1) is involved in the regulation of gene expression and RNA protection. While numerous studies have established that SND1 protein expression is modulated by cellular stresses associated with tumor growth, hypoxia, inflammation, heat- shock and oxidative conditions, little is known about the factors responsible for SND1 expression. Here, we have approached this question by analyzing the transcriptional response of human SND1 gene to pharmacological endoplasmic reticulum (ER) stress in liver cancer cells. Results: We provide first evidence that SND1 promoter activity is increased in human liver cancer cells upon exposure to thapsigargin or tunicamycin or by ectopic expression of ATF6, a crucial transcription factor in the unfolded protein response triggered by ER stress. Deletion analysis of the 5'-flanking region of SND1 promoter identified maximal activation in fragment (-934, +221), which contains most of the predicted ER stress response elements in proximal promoter. Quantitative real- time PCR revealed a near 3 fold increase in SND1 mRNA expression by either of the stress- inducers; whereas SND1 protein was maximally upregulated (3.4-fold) in cells exposed to tunicamycin, a protein glycosylation inhibitor. Conclusion: Promoter activity of the cell growth- and RNA-protection associated SND1 gene is up-regulated by ER stress in human hepatoma cells.

The Interplay between Natural Selection and Susceptibility to Melanoma on Allele 374F of SLC45A2 Gene in a South European Population

We aimed to study the selective pressures interacting on SLC45A2 to investigate the interplay between selection and susceptibility to disease. Thus, we enrolled 500 volunteers from a geographically limited population (Basques from the North of Spain) and by resequencing the whole coding region and intron 5 of the 34 most and the 34 least pigmented individuals according to the reflectance distribution, we observed that the polymorphism Leu374Phe (L374F, rs16891982) was statistically associated with skin color variability within this sample. In particular, allele 374F was significantly more frequent among the individuals with lighter skin. Further genotyping an independent set of 558 individuals of a geographically wider population with known ancestry in the Spanish population also revealed that the frequency of L374F was significantly correlated with the incident UV radiation intensity. Selection tests suggest that allele 374F is being positively selected in South Europeans, thus indicating that depigmentation is an adaptive process. Interestingly, by genotyping 119 melanoma samples, we show that this variant is also associated with an increased susceptibility to melanoma in our populations. The ultimate driving force for this adaptation is unknown, but it is compatible with the vitamin D hypothesis. This shows that molecular evolution analysis can be used as a useful technology to predict phenotypic and biomedical consequences in humans.

Physico-chemical interactions between compartment-forming lipids and other prebiotically relevant biomolecules

Lipids are essential constituents of contemporary living cells, serving as structural molecules that are necessary to form membranous compartments. Amphiphilic lipid-like molecules may also have contributed to prebiotic chemical evolution by promoting the synthesis, aggregation and cooperative encapsulation of other biomolecules. The resulting compartments would allow systems of molecules to be maintained that represent microscopic experiments in a natural version of combinatorial chemistry. Here we address these possibilities and describe recent results related to interactions between amphiphiles and other biomolecules during early evolution toward the first living cells.

Lymphangioleiomyomatosis Biomarkers Linked to Lung Metastatic Potential and Cell Stemness

Lymphangioleiomyomatosis (LAM) is a rare lung-metastasizing neoplasm caused by the proliferation of smooth muscle-like cells that commonly carry loss-of-function mutations in either the tuberous sclerosis complex 1 or 2 (TSC1 or TSC2) genes. While allosteric inhibition of the mechanistic target of rapamycin (mTOR) has shown substantial clinical benefit, complementary therapies are required to improve response and/or to treat specific patients. However, there is a lack of LAM biomarkers that could potentially be used to monitor the disease and to develop other targeted therapies. We hypothesized that the mediators of cancer metastasis to lung, particularly in breast cancer, also play a relevant role in LAM. Analyses across independent breast cancer datasets revealed associations between low TSC1/2 expression, altered mTOR complex 1 (mTORC1) pathway signaling, and metastasis to lung. Subsequently, immunohistochemical analyses of 23 LAM lesions revealed positivity in all cases for the lung metastasis mediators fascin 1 (FSCN1) and inhibitor of DNA binding 1 (ID1). Moreover, assessment of breast cancer stem or luminal progenitor cell biomarkers showed positivity in most LAM tissue for the aldehyde dehydrogenase 1 (ALDH1), integrin-beta 3 (ITGB3/CD61), and/or the sex-determining region Y-box 9 (SOX9) proteins. The immunohistochemical analyses also provided evidence of heterogeneity between and within LAM cases. The analysis of Tsc2-deficient cells revealed relative over-expression of FSCN1 and ID1; however, Tsc2-deficient cells did not show higher sensitivity to ID1-based cancer inhibitors. Collectively, the results of this study reveal novel LAM biomarkers linked to breast cancer metastasis to lung and to cell stemness, which in turn might guide the assessment of additional or complementary therapeutic opportunities for LAM.

Exploring Genetic Factors Involved in Huntington Disease Age of Onset: E2F2 as a New Potential Modifier Gene

Age of onset (AO) of Huntington disease (HD) is mainly determined by the length of the CAG repeat expansion (CAGexp) in exon 1 of the HTT gene. Additional genetic variation has been suggested to contribute to AO, although the mechanism by which it could affect AO is presently unknown. The aim of this study is to explore the contribution of candidate genetic factors to HD AO in order to gain insight into the pathogenic mechanisms underlying this disorder. For that purpose, two AO definitions were used: the earliest age with unequivocal signs of HD (earliest AO or eAO), and the first motor symptoms age (motor AO or mAO). Multiple linear regression analyses were performed between genetic variation within 20 candidate genes and eAO or mAO, using DNA and clinical information of 253 HD patients from REGISTRY project. Gene expression analyses were carried out by RT-qPCR with an independent sample of 35 HD patients from Basque Country Hospitals. We found suggestive association signals between HD eAO and/or mAO and genetic variation within the E2F2, ATF7IP, GRIN2A, GRIN2B, LINC01559, HIP1 and GRIK2 genes. Among them, the most significant was the association between eAO and rs2742976, mapping to the promoter region of E2F2 transcription factor. Furthermore, rs2742976 T allele patient carriers exhibited significantly lower lymphocyte E2F2 gene expression, suggesting a possible implication of E2F2-dependent transcriptional activity in HD pathogenesis. Thus, E2F2 emerges as a new potential HD AO modifier factor.